Математический форум Math Help Planet
http://mathhelpplanet.com/

Дифференциальные уравнения второго порядка
http://mathhelpplanet.com/viewtopic.php?f=38&t=32908
Страница 1 из 1

Автор:  Eugene85 [ 29 апр 2014, 20:00 ]
Заголовок сообщения:  Дифференциальные уравнения второго порядка

Есть уравнение x''+(x^2-1)x'+x=0. Как представить это уравнение в виде системы дифференциальных уравнений первого порядка? Распишите пожалуйста как это делать.

Автор:  Eugene85 [ 29 апр 2014, 20:38 ]
Заголовок сообщения:  Re: Дифференциальные уравнения второго порядка

Код:
>> global Q1 Q2 C1x C1y C2x C2y
>> Q1=-50; Q2=-0.; C1x=5; C1y=0; C2x=0; C2y=10;
>> x0=0; y0=0; vx0=0; vy0=4.3; T1=4000;
>> tol=1e-6;
>> [t,h]=ode45(@pointq12,[0,T1],[x0,y0,vx0,vy0], ...
odeset ('RelTol',tol));

Функция:
function f=pointq12(t,x)
global Q1 Q2 C1x C1y C2x C2y
f=zeros(3,1);
f(1)=x(3);
f(2)=x(4);
f(3)=[Q1*(x(1)-C1x)/(sqrt((x(1)-C1x)^2+(x(2)-C1y)^2))^3+...
Q2*(x(1)-C2x)/(sqrt((x(1)-C2x)^2+(x(2)-C2y)^2))^3....
+Q1*(x(2)-C1y)/(sqrt((x(1)-C1x)^2+(x(2)-C1y)^2))^3....
Q2*(x(2)-C2y)/(sqrt((x(1)-C2x)^2+(x(2)-C2y)^2))^3];

В результате выполнения всего этого выдает ошибку:

??? In an assignment A(I) = B, the number of elements in B and
I must be the same.

Error in ==> pointq12 at 6
f(3)=[Q1*(x(1)-C1x)/(sqrt((x(1)-C1x)^2+(x(2)-C1y)^2))^3+...

Error in ==> odearguments at 109
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, ...

Подскажите, что не так?

Страница 1 из 1 Часовой пояс: UTC + 3 часа [ Летнее время ]
Powered by phpBB © 2000, 2002, 2005, 2007 phpBB Group
http://www.phpbb.com/