Как записать в Matlab систему ДУ 2-го порядка , Не могу понятным для Matlab образом записать систему ДУ 2-го порядка с нач. усл.

Всем привет!

Имею систему ДУ второго порядка с н.у. Посчитал её в MatCad'е - все замечательно. Теперь тоже самое нужно проделать в Matlab. Скачал кучу самоучителей, в т.ч. 1000-страничный учебник Ануфриева. Нигде не нашел подробного описания записи системы ДУ.
В интернетах нашел что-то более внятное. Записал - не считает, выдает ошибки записи. Ниже приведена моя система ДУ и программа моего кривого М-файла.
Прикреплённый файл_______.doc (26,5 Кбайт, скачиваний: 174)
Введенные обозначения:
x=x(1)
x'=x(2)
y=y(1)
y'=y(2)
Начальные условия:
x(0)=0
x'(0)=0
y(0)=0
y'(0)=0
Промежуток интегрирования от 0 до 1.

function f=rod(t,x,y)
lambda=2*10^7;l=100;d=6;theta=1;F=30;E=2*10^6;
k=l^2/(d*pi^2);
f=[x(2);y(2);(4*E*F*lambda*theta-8*lambda^2*k*t^2+...
    2*E*F*h*lambda+4*E*F*k*lambda*y(1)+...
    4*E*F*k*lambda*x(1)+8*E*F*k*lambda*t*y(2)+...
    8*E*F*k*lambda*t*x(2))/(E^2*F^2*h-...
    4*E*F*k*lambda*t^2);...
    (8*k*lambda^2*t^2+4*E*F*lambda*theta-...
    2*E*F*h*lambda+4*E*F*k*lambda*y(1)+...
    4*E*F*k*lambda*x(1)+8*E*F*k*lambda*t*y(2)+...
    8*E*F*k*lambda*t*x(2))/(E^2*F^2*h-...
    4*E*F*k*lambda*t^2)];
function dynpoint=s(t,x,y)
[t,x,y]=ode45(@rod,[0 1],[0 0 0 0]);
plot(t,x(:,1),'b.:')
hold on
plot(t,y(:,1),'r.-')
hold off

Буду признателен за помощь и указания на ошибки.

Ищи в хелпе
ode45

Так подправил.
Считает. Но результат странный. Параметра h не было задал 1. Копаться в страшной формул возьмусь только за отдельную плату)
Прикреплённый файлdypoint.m (474 байт, скачиваний: 208)
Прикреплённый файлrod.m (877 байт, скачиваний: 194)

LuckyboyLuckyboy
здравствуйте..можете помочь разобраться в разспознавании фигурок на матлабе? Как лучше делать? Через тулбокс сначала?