function d=delta(t,j,k)
% this is the triangular function used in the wavelet generation of
% Brownian motion
if (j>0)
    d=del(t*2^j-k);
else
    d=t.*(t<1).*(t>0);
end
