Advanced Search

Question: How to solve complex-valued IVP?

Question:How to solve complex-valued IVP?

adel-00 135 Maple
February 22 2016
1

Hi experts,

attached the following code,, but i dont know what is the wrong with it, looking forward for helpful advise.

dsys :={diff(u(t),t)=-(N1+M*cos(2*I*omega*t))*u(t)-1+v(t)*exp(-2*I*omega*t)+w(t)*exp(2*I*omega*t), diff(v(t),t)=-(N1+I*Delta-2*M*exp(2*I*omega*t))*v(t)-(N1+u(t))*exp(2*I*omega*t)-2*M, diff(w(t),t)=-(N1-I*Delta-2*M*exp(-2*I*omega*t))*w(t)-(N1+u(t))*exp(-2*I*omega*t)-2*M}:
res:=dsolve(dsys union {u(0)=-1,v(0)=0,w(0)=0},numeric,output=listprocedure,maxfun=0):
plots[odeplot](res,[[t,(Re(w(t)))]],0..10,axes=boxed,titlefont=[SYMBOL,14],font=[1,1,18],color=black,linestyle=1,tickmarks=[3, 4],font=[1,1,14],thickness=2,titlefont=[SYMBOL,12]);
Warning, computation interrupted
Please Wait...