function s = psi_doubleprime(t,beta)
  
  s = -.5 ./ sqrt(t + beta^2).^3;