2014 dxdy logo

Научный форум dxdy

Математика, Физика, Computer Science, Machine Learning, LaTeX, Механика и Техника, Химия,
Биология и Медицина, Экономика и Финансовая Математика, Гуманитарные науки




 
 Решение нелинейных диф уравнений
Сообщение19.04.2013, 13:55 
Помогите, пожалуйста, решить следующую систему дифф. уравнений в MatLam. Необходимо получить график в трехмерном пространстве (x,y,z). У меня показывает ошибки в командном окне. Может я что-то неправильно набрала M-файле?

M-файл

function dy=pr8(t,y,H)
dy=zeros(3,2,1);
dy(1)=(1/(0.275*(5.88*((7.5/H)^3)-1.78*((7.5/H)^2)+0.541*(7.5/H))*(0.037569-y(1)*(0.290*(3.442*((7.5/H)^3)-2.024*((7.5/H)^2)+0.4795*(7.5/H)+1)+0.02186))+y(2)*(0.029*(3.964*((7.5/H)^3)-2.399*((7.5/H)^2)+0.5157*(7.5/H)+1)+0.02696189)-0.494*y(1)*y(2)*sign(y(1),y(2));
dy(2)=(1/(0.016*(10*(7.5/H)^3)-6.8*(7.5/H)^2)+2.35*(7.5/H)+1)*(0.018671793+y(1)*(0.068*(3.756*((7.5/H)^3)-2.254*((7.5/H)^2)+0.5449*(7.5/H)+1)-0.01079567))+y(2)*(-0.031*(4.04*((7.5/H)^3)-2.43*((7.5/H)^2)+0.555*(7.5/H)+1)-0.03161575595)+0.139*y(1)*y(2)+sign(y(1),Y(2));
dy(3)=(0.01/(0.176*(-0.559*((7.5/H)^3)+3.62*((7.5/H)^2)-0.406*(7.5/H)+1)*(y(3)^2))-((0.0247*y(3))/(0.176*(-0.559*((7.5/H)^3)+3.62*((7.5/H)^2)-0.406*(7.5/H)+1))-((0.275*(5.88*((7.5/H)^3)-1.78*((7.5/H)^2)+0.541*(7.5/H)+1))/(0.176*(-0.559*((7.5/H)^3)+3.62*((7.5/H)^2)-0.406*(7.5/H)+1)))*y(1)*y(2)*y(3)+dy(1)*y(1)*y(3);


Командное окно

>> [T,Y,H]=ode45(@pr8,[0 500],[0 200 100], [0 100]);
??? Error using ==> feval
Error: File: pr8.m Line: 3 Column: 264
Unbalanced or unexpected parenthesis or bracket.

Error in ==> odearguments at 111
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.

Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...

-- 19.04.2013, 15:00 --

В MatLab, a не в matLam. Я извиняюсь

 
 
 [ 1 сообщение ] 


Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group