Question: how can I solve a system of differential equations in Maple?

Hello, I have tried to solve a system of ODE about the gyroscope for a week. I have entered the following commands : >eq1:=(l^2+r^2/2)*(sin(theta(t))^2*diff(diff(psi(t),t),t)+2*diff(psi(t),t)*diff(theta(t),t)*cos(theta(t)))-r^2*(cos(theta(t))*diff(psi(t),t)+omega(t))*diff(theta(t),t)=0;>eq2:=(l^2+r^2/2)*(diff(diff(theta(t),t),t)-(diff(psi(t),t))^2*sin(theta(t))*cos(theta(t)))+r^2*(diff(psi(t),t)*cos(theta(t))+omega(t))*diff(psi(t),t)*sin(theta(t))=g*l*sin(theta(t)); >eq3:=diff(diff(psi(t),t)*cos(theta(t))+omega,t)=0; >l:=0.05;r:=0.063/2;g:=0.81;omega:=evalf(2*Pi*10000/3600); >ic1:=theta(0)=evalf(9*Pi/180),psi(0)=0,(D(theta(t)))(0)=0,D(psi(t))(0)=0; >Sol:=dsolve({eq1,eq2,eq3} union {ic1},{theta(t),psi(t)},numeric); And maple answers me "Error (in DEtools/convertsys) unable to convert to an explicit first-order system". What could I do? My purpose is to use ODEplot after that.