As delivered, the Smath integrator accuracy is limited to "n subspace 10000".
The attached document exposes two examples that the native integrator fails.
It suggests to code for "n at least 1 000 000".
On your fast PC, will takes ~ 35 min to calculate the document.
Accordingly to Wolfram [Mathematica], they will fail no integration.
In Engineering works, this is not an option.
Conclusively: would be nice Smath implements the equivalent Mathcad/Mathsoft.
various integration options ... I don't have those codes [Mathsoft proprietary].

Jean

Integrate Adaptive Simpson.sm (30 KiB) downloaded 65 time(s).

Romberg.sm (13 KiB) downloaded 63 time(s).

Hello Martin,

The Romberg integration transform improves significantly the accuracy
of the result when applicable. Mathcad 11 and earlier performs two tests
to validate application: a trapezoidal on some length and a Richardson
interpolation. If the two tests coincide on some pre-decided basis,
then Romberg is accepted. Otherwise, it performs a "Romberg open end"
[I don't know what it means], to decide an adaptive method [for which
I lost the code also].
From what you are saying, Maxima works on timing => looks weird, isn't it.
The attached exemplifies applicability. All that is because Simpson is
universally a good integrator, but NOT for all integrands.

Integrate Romberg Demo.sm (52 KiB) downloaded 60 time(s).

Timeout-problem resolved. http://smath.info/bts/Issues/IssueDetail.aspx?id=189


I am open to proposals on what integrator to use. The Maxima plugin has the option to search for a symbolic solution if symbolic optimization is requested. If numeric optimization is requested and no unresolved constants remain (units are ok) in the expression, then we could use whatever numeric procedure, be it from Maxima or whatever library.

Wrote

I am open to proposals on what integrator to use. The Maxima plugin has the option to search for a symbolic solution if symbolic optimization is requested. If numeric optimization is requested and no unresolved constants remain (units are ok) in the expression, then we could use whatever numeric procedure, be it from Maxima or whatever library.

The best library is the Mathcad/Mathsoft, don't know about Mathcad/PTC.
From all that I posted, easy to manage 'singular & infinite limit(s)'.
For those cases that need be checked for convergence [up to 1 000 000],
it is just a matter of increasing the "Integral accuracy" in options.

WolframAlpha is the piece of gold as it gives the indefinite integral
on top of the definite integral result [cumulative within limits].

Cheers, Jean