Система ДУ - подскажите как повысить точность?

[Vaselisa]

Добрый день, решаю систему ДУ в программе, через встроенную функцию ode15s:

T=[1:1:86399];
options=odeset('BDF','on', 'AbsTol',1e-55, 'RelTol', 1e-13);
 [t,Y] =  ode15s('fun',T,y0, options);

строю график - функция растет, хотя не должна.
Подскажите пожалуйста, как повысить точность расчета? (Если убираю из уравнения движения ЦМ несферичность - радиус-вектор не растет.)

[Pphantom]

Догадаться по названию функции, что "программа" - это MATLAB/Octave/etc , еще более-менее реально. Но вот осознать, какую систему ДУ Вы решаете, при чем тут "несферичность" и т.п., уже совершенно невозможно.

[Vaselisa]

пытаюсь решить систему уравнений движения центра масс в абсолютной геоцентрической системе координат

[Pphantom]

Vaselisa в сообщении #1113022 писал(а):

пытаюсь решить систему уравнений движения центра масс в абсолютной геоцентрической системе координат

Центра масс чего?