ó k™xTc@ srdZddlmZddlZd„Zed„ZedkrnddlZedej dkƒd GHndS( s Approximating Log(BesselK[1/4,x]) using a Pade form. Coefficients were fitted with gnuplot on numbers obtained from mathematica. f(x) = (a0+a1*x+a2*x*x+a3*x*x*x+a4*x**4)/(1+b1*x+b2*x*x+b3*x*x*x+b4*x**4+b5*x**5) The total range has been split in 3 regions: 0s