function s = psi_doubleprime(t,beta)
  
  s = -.5 * (t + beta^2).^(-1.5);