Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
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 (31kb) downloaded 28 time(s).
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,980 Was thanked: 1120 time(s) in 718 post(s)
|
|
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(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 (53kb) downloaded 21 time(s).
|
|
|
|
Rank: Advanced Member Groups: Registered
Joined: 15/04/2012(UTC) Posts: 1,980 Was thanked: 1120 time(s) in 718 post(s)
|
Timeout-problem resolved. http://smath.info/bts/Is.../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. |
|
|
|
|
Rank: Guest
Groups: Registered
Joined: 04/07/2015(UTC) Posts: 6,868 Was thanked: 980 time(s) in 808 post(s)
|
Originally Posted by: mkraska 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
|
|
|
|
Forum Jump
You cannot post new topics in this forum.
You cannot reply to topics in this forum.
You cannot delete your posts in this forum.
You cannot edit your posts in this forum.
You cannot create polls in this forum.
You cannot vote in polls in this forum.