@clement.erdman лучше всего воспользоваться функцией ode45(), она помогает решить дифуры и системы дифференциальных уравнений в подавляющем большинстве случаев.
Привожу пример решения системы линейных дифференциальных уравнений и полученные решения
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 |
% система линейных дифференциальный уравнений: % y1'=0.4*t; % y2'=3-t; tspan = [0 5]; % диапазон получения решения y0 = [0 1]; % начальные значения [t,y] = ode45(@(t,y) [0.4*t; 3-t], tspan, y0); plot(t,y,'-o'); grid on; title('Решение системы дифференциальных уравнений'); xlabel('t'); ylabel('y'); legend('y1', 'y2', 'location', 'best'); saveas(gca, 'ODE1.jpg'); |