Hi,

It seems unlikely that speed of execution will play a part in the
outcome of this program. 

I've not used Ocaml but one possibility is that string_of_float commands
modify their arguments? Perhaps when the numbers in piA or piB are very
small they get truncated to zero to match what is actually printed to
the screen. This would avoid underflow errors at least.

Also might be worth using exp1m function instead of exp to compute the
exponential.

Best,
Bojan

-- 
Bojan Nikolic          ||          http://www.bnikolic.co.uk



_______________________________________________
Help-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/help-gsl

Reply via email to