%   function for noneq
    function val=noneqfun(z)
    global t
    val=exp(-(z+t)).*besselj(0,i.*(sqrt(4.*z.*t)));