> (%i1)
> quad_qags(lambda([x],quad_qags(lambda([y],exp(x-y^2)),y,1,2+x)[1]),x,0,1);
> (%o1) [.2384683615524653, 2.647530656401834e-15, 21, 0]
>
>
> I have an interest in determining knot energies, and double integrals
> arise naturally in these calculations. My tests show that for double
> integrals quad_qags usually provides a better estimate than romberg.
>
> Regards,
> James M Rotenberry
----------------------------------
Hi again,
The expression "g" can be used if double quoted (two single quotes):
(%i17) quad_qags(lambda([x], first(quad_qags(lambda([y],
''g),y,0,x/2))),x,1,3);
(%o17) [0.819302, 9.09608385E-15, 21, 0]
Best Wishes,
Ted