Интерпретация вывода dsolve для использования с ODE45 [MATLAB]

я пытался решить дифференциальное уравнение второго порядка, используя этот 1, однако Я не мог понять это правильно и не нашел ничего полезного в Интернете, но я считаю, что добился прогресса.

Я использовал dsolve;

syms x(t) v(t) fi(t)

[x(t), v(t)] = dsolve(diff(x) == v, diff(v) == fi/m, x(0) == [-L, -L], v(0) == [5, 10] )

что дает мне;

   x(t) =



int(fi(x)/5, x, 0, t, 'IgnoreSpecialCases', true, 'IgnoreAnalyticConstraints', true) - 5





v(t) =

C2 + t*(int(fi(x)/5, x, 0, t, 'IgnoreSpecialCases', true, 'IgnoreAnalyticConstraints', true) - 5) + int(-(x*fi(x))/5, x, 0, t, 'IgnoreSpecialCases', true, 'IgnoreAnalyticConstraints', true)

Теперь мне нужна помощь в обработке результата, и мне интересно, могу ли я использовать этот результат, чтобы получить что-то из ode45? Также я хочу построить решение как эталонную траекторию для моделирования 500 частиц, движущихся через силовое поле.

используя ODE45:

function dxdt = solution(t,y0)
frprintf('Second stop')
.....
dxdt = [x, v]
end

вызов из основного файла:

t:dt:t_f
y0 = [x0,v0]
fprintf('first stop')
[x, v] = ode45(@solution, y0, t)

Я настроил его так, что если код будет работать гладко, он будет печатать «первая остановка, вторая остановка, третья остановка» и «четвертая остановка», он печатает только первую остановку, и вот где я получаю ошибку.