Сразу оговорюсь: я не спец в теории ДУ (тем более, использования численных методов). Недавно знакомый-инженер попросил помочь решить ему численно систему уравнений (3-х) динамики. Решил использовать Maple. Набрал такие коды: restart:A:=42403: > B:=269.95: > C:=79051: > J1:=0.55: > J2:=1.02: > a:=196200: > b:=14342220: > c:=1589220:
> dsys := {A-B*diff(phi1(t),t)=J1*diff(phi1(t),t,t)+C*(phi1(t)-phi2(t)),J2*diff(phi2(t),t,t)-C*(phi1(t)-phi2(t))=-(a*Q(t)+b*((Q(t))^3)+c*diff(Q(t),t)*((Q(t))))/(157.08*Q(t)/0.0875), diff(phi2(t),t)= 157.08*Q(t)/0.0875,Q(0)=0, D(phi1)(0) = 157.08, Q(5)=0.0875, D(phi1)(5)=154.98}: > dsn := dsolve(dsys,numeric):
Программа после их исполнения выдает следующее сообщение: Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system
Что не так с этой системой?
|