10-30-2017, 01:47 PM
hi,
this is only a first version and has'n control for singular values. I hope someone could help to do better, thank you.
Fresnel_S(x) and Fresnel_C(x) are the functions for S(x) and C(x)
Fresnel_Sn(x) and Fresnel_Cn(x) are the correspondent functions, normalised with a π/2 factor by some authors, like Abramowitz and Stegun (see here).
Enjoy!
Salvo Micciché
this is only a first version and has'n control for singular values. I hope someone could help to do better, thank you.
Fresnel_S(x) and Fresnel_C(x) are the functions for S(x) and C(x)
Fresnel_Sn(x) and Fresnel_Cn(x) are the correspondent functions, normalised with a π/2 factor by some authors, like Abramowitz and Stegun (see here).
Enjoy!
Salvo Micciché
Code:
EXPORT Fresnel_S(x)
// Fresnel integral S(x)
BEGIN
local t;
RETURN int(SIN(t^2),t,0,x);
END;
EXPORT Fresnel_Sn(x)
// Fresnel integral S(x) normalized π/2
BEGIN
local t;
RETURN int(SIN((π/2)*t^2),t,0,x);
END;
EXPORT Fresnel_C(x)
// Fresnel integral C(x)
BEGIN
local t;
RETURN int(COS(t^2),t,0,x);
END;
EXPORT Fresnel_Cn(x)
// Fresnel integral C(x) normalized π/2
BEGIN
local t;
RETURN int(COS((π/2)*t^2),t,0,x);
END;