Hallo Scilab experts, If I use ode() to simulate the two-body problem sun-earth (state vector of 12 barycentric coordinates), there is a large discrepancy between the numerical and the analytical solution: In a simulation time of 550 (analytically calculated) orbital periods, the Earth completes 551 orbits. The circular orbit shape is fine. I used ode( ) with the default solver and the default tolerances. By narrowing the tolerances, the result becomes not more accurate. Can another solver help here that improves the numerical precision?
Kind regards Jens This email and any attachments are intended solely for the use of the individual or entity to whom it is addressed and may be confidential and/or privileged. If you are not one of the named recipients or have received this email in error, (i) you should not read, disclose, or copy it, (ii) please notify sender of your receipt by reply email and delete this email and all attachments, (iii) Dassault Systèmes does not accept or assume any liability or responsibility for any use of or reliance on this email. Please be informed that your personal data are processed according to our data privacy policy as described on our website. Should you have any questions related to personal data protection, please contact 3DS Data Protection Officer https://www.3ds.com/privacy-policy/contact/
_______________________________________________ users mailing list - users@lists.scilab.org Click here to unsubscribe: <mailto:users-unsubscr...@lists.scilab.org> https://lists.scilab.org/mailman/listinfo/users
