Численное решение системы ОДУ в Maple. Помогите!

[reterty]

Сразу оговорюсь: я не спец в теории ДУ (тем более, использования численных методов).
Недавно знакомый-инженер попросил помочь решить ему численно систему уравнений (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

Что не так с этой системой?

[Pphantom]

Очевидно, это не система первого порядка, и о ее переделке в таковую вы должны были позаботиться сами.