(y,err)=hyp3f0(a,b,c,x) returns (y,err) with the hypergeometric function 3F0 in y and an error estimate in err.
scipy.special.hyp1f2
scipy.special.pbdv