| Математический форум 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/ |
|