> We can work around some of the problems and produce an answer close to
> what Mathematica produces by using bfloats.
>
> fpprec:32;
> x=1;
> y=1b-15;
> elliptic_kc(((x-y)/(x+y))^2);
>
> 17.9265353780153059...
>
> Ray
I suspected that to be the problem when I tried Scilab and saw the truncated results. I tried "fpprec" but with "1e-15", not "1b-15", hence the difference that lead me to that conclusion.
The numbers are needed for the current work and they're the worst case scenario; I didn't choose them for a willing head-ache :)
Thank you,
Vlad.