FUNCTION inverf,z
z2=!pi*z^2
s=(((((34807/364953600*z2+4369.0/11612160)*z2+127.0/30640)*z2 + $
  7.0/960)*z2+1.0/24)*z2+1/2.0)*z*sqrt(!pi)

return,s
end