Hello experts in Maple,
I'm going to solve IVP numerically by this command bellow:
var:={u(t),v(t),w(t)}:
dsys1 :=diff(u(t),t)=-(Gamma+AM*cos(q*t+phi))*u(t)+(Delta-AM*sin(q*t+phi))*v(t)+Omega*w(t), diff(v(t),t)=-(Delta+AM*sin(q*t+phi))*u(t)-(Gamma-AM*cos(q*t+phi))*v(t), diff(w(t),t)=-0.5-2*Gamma*w(t)-Omega*u(t):
dsol1 :=dsolve({dsys1,ini1},numeric,method=lsode,var, abserr=1e-9, relerr=1e-8,output=listprocedure):
dsolu:=subs(dsol1,u(t)):dsolv:=subs(dsol1,v(t)):dsolw:=subs(dsol1,w(t)):...