|
Здравствуйте. Необходимо методом Рунге-Кутта 4-го порядка решить систему уравнений вида:
dy/dt=-2y+100x
dz/dt=-10z+25gz+1000y+18ax
da/dt=-35xa-46a+41gz
db/dt=-15b+73a
b+a+g=1
x+y+z=1
При этом когда t=o: x=1, y=0, z=0, a=0, b=0, g=0.
В Maple возникает ошибка:
> sys:={diff(y(t),t)=-2*y(t)+100*x(t),diff(z(t),t)=-10*z(t)+25*g(t)*z(t)+1000*y(t)+18*a(t)*x(t),diff(a(t),t)=-35*x(t)*a(t)-46*a(t)+41*g(t)*z(t),diff(b(t),t)=-15*b(t)+73*a(t),b(t)+a(t)+g(t)=1,x(t)+y(t)+z(t)=1};
> fens:={x(t),y(t),z(t),a(t),b(t),g(t)};
> F:=dsolve({sys,x(0)=1,y(0)=0,z(0)=0,a(0)=0,b(0)=0,g(0)=0},fens,numeric);
Error, (in dsolve/numeric/process_input) system must be entered as a set/list of expressions/equations
> mathod=rkf45;
Не подскажите, с чем это связано?
|
