share/contrib/fresnel/fresnelFG.mac

   1 fresnelFG(z):=block([fc,fs,cz2,sz2,zp2,FLOAT2BF:true],
   2    z:bfloat(z),
   3    zp2:bfloat(%pi)*z*z/2.0b0,
   4    cz2:cos(zp2),
   5    sz2:sin(zp2),
   6    fc:fresnelC(z),
   7    fs:fresnelS(z),
   8    [z,
   9     ((2 * fc * sz2 - sz2 - 2 * cz2 * fs + cz2)/(2 * sz2^2 + 2 * cz2^2)),
  10     (( - 2 * fs * sz2 + sz2 - 2 * cz2 * fc + cz2)/(2 * sz2^2 + 2 * cz2^2))]
  11 )$