Skip to main content

k0ef_small

pxfm::bessel::k0ef

Function k0ef_small 

Source 
fn k0ef_small(x: f32) -> f32
Expand description

K0(x) + log(x) * I0(x) = P(x^2) hence K0(x) = P(x^2) - log(x)*I0(x)

Polynomial generated by Wolfram Mathematica:

<<FunctionApproximations`
ClearAll["Global`*"]
f[x_]:=BesselK[0,x]+Log[x]BesselI[0,x]
g[z_]:=f[Sqrt[z]]
{err, approx}=MiniMaxApproximation[g[z],{z,{0.000000001,1},6,0},WorkingPrecision->60]
poly=Numerator[approx][[1]];
coeffs=CoefficientList[poly,z];
TableForm[Table[Row[{"'",NumberForm[coeffs[[i+1]],{50,50}, ExponentFunction->(Null&)],"',"}],{i,0,Length[coeffs]-1}]]