Использование функции ode45 в Matlab

[Kpoxa]

Помогите пожалуйста!!! У меня имеется 3 уравнения, их необходимо решить с помощью функции ode45 в Matlab...
Вот и сама система уравнений:
m1*(d^2 x1/dt^2)= -k1x1-k2(x1-x2)
m2*(d^2 x2/dt^2)= -k2(x2-x1)-k3(x2-x3)
m3*(d^2 x3/dt^2)= -k4x3-k3(x3-x2)
Необходимо найти 3 значения: x1, x2 и x3.
Кто знает, как пользоваться этой функцией, подскажите!!!

[photon]

Читаем help, там естъ разобранный пример для:

$y_1'=y_2$
$y_2'=\mu(1-y_1^2)y_2-y_1$

Цитата:

The code below represents the van der Pol system in the function, vdp1. The vdp1 function assumes that $\mu=1$. The variables $y_1$ and $y_2$ are the entries $y(1)$ and $y(2)$ of a two-element vector.

function dydt = vdp1(t,y)
dydt = [y(2); (1-y(1)^2)*y(2)-y(1)];

For the van der Pol system, you can use ode45 on time interval [0 20] with initial values $y(1) = 2$ and $y(2) = 0$.

[t,y] = ode45(@vdp1,[0 20],[2; 0]);

а дальше делаете просто по аналогии для своей системы