I1(x)≈x2(1+12[x2]2+[x2]4P([x2]2));x<7.75