function s = psi_prime(t,beta)
     s = 1 ./ sqrt(t + beta^2);