I was attempting to perform some numerical quadrature, and I thought it might be nice to use the NAG library. According to HyperDoc, this command:
d01ajf(0.0 ,1.0 ,0.0 ,1.0e-4 ,800 ,200 ,-1 ,(X/(X^7+1)::_expression_ Float) :: ASP1(F))
should integrate x/(^7+1) between 0 and 1. But in fact it doesn't - it just hangs Axiom. What do I need to do to use these routines, and is there an easy front end command to them?
Thanks,
Alasdair
--
Blog:
http://amca01.wordpress.com