function y = tlhtdens( a,b,t)
% tilted line a+bt hitting time density
y=(abs(a)/sqrt(2*pi))./t.^(3/2).*exp(-0.5*(a+b*t).^2./t);
end

