2
$\begingroup$

I have a problem with ParametricNDSolve:

a = 1; d = 0.13; s = 0.4; f = .1;

eq = c'[t] + I a (Abs[c[t]])^2 c[t] + d (Abs[c[t]])^2 c[t] -s (Abs[c[t]])^2 c[t] - f Exp[-I w t] == 0;

r =  ParametricNDSolve[{eq, c[0] == 0}, c[t], {t, 0, 25}, {w}];
Plot[Evaluate[{Re@#, Im@#} &[c[t] /. r, {w, 1, 10, .1}]], {t, 0, 25}]

i have no idea what is wrong here; i checked the documentation center and some posts on similar topic - and i'm still stuck. Any help appreciated! Thanks in advance!

Improve this question
$\endgroup$

2 Answers 2

Reset to default
1
$\begingroup$ 
a = 1; d = 0.13; s = 0.4; f = .1;
eq = c'[t] + I a (Abs[c[t]])^2 c[t] + d (Abs[c[t]])^2 c[t] - 
    s (Abs[c[t]])^2 c[t] - f Exp[-I w t] == 0;
r = ParametricNDSolve[{eq, c[0] == 0}, c, {t, 0, 25}, {w}];
Plot[Table[Through@{Re, Im}@c[w][t] /. r, {w, 1, 10, .1}], {t, 0, 25}, Evaluated -> True]

enter image description here

Update: various views

pre0 = Plot[Table[Re@c[w][t] /. r, {w, 1, 10, .1}], {t, 0, 25}, 
   Evaluated -> True, ImageSize -> 500,  
   PlotLabel -> Style["Plot@Table[Re@c[w][t]/.r,{w,1,10,.1}],{t,0,25}]", 16]];
pim0 = Plot[Table[Im@c[w][t] /. r, {w, 1, 10, .1}], {t, 0, 25}, 
   Evaluated -> True, ImageSize -> 500, 
   PlotLabel -> Style["Plot@Table[Im@c[w][t]/.r,{w,1,10,.1}],{t,0,25}]", 16]];
pre1 = Plot[Table[Re@c[w][t] /. r, {t, 0, 25, 5}], {w, 1, 10}, 
   Evaluated -> True, ImageSize -> 500, 
   PlotLabel -> Style["Plot[Re@c[w][t]/.r,{t,0,25,5}],{w,1,10}]", 16]];
pim1 = Plot[Table[Im@c[w][t] /. r, {t, 0, 25, 5}], {w, 1, 10}, 
   Evaluated -> True, ImageSize -> 500, 
   PlotLabel -> Style["Plot[Table[Im@c[w][t]/.r,{t,0,25,5}],{w,1,10}]", 16]];

Row[{pre0, pim0}, Spacer[5]]

enter image description here

Row[{pre1, pim1}, Spacer[5]]

enter image description here

Plot3D[Re@c[w][t] /. r, {t, 0, 10}, {w, 1, 5}, Evaluated -> True,
  PlotPoints -> 100, BoxRatios -> 1, ImageSize -> 500,
  ColorFunction -> "TemperatureMap", Mesh -> None, Lighting -> "Neutral"]

enter image description here

Improve this answer
$\endgroup$
2
  • $\begingroup$ thank you guys! and btw. is there a way to plot the c dependence of w without transforming the equation explicitly? $\endgroup$
    – Bar Lac
    Nov 6, 2014 at 20:35
  • $\begingroup$ @BarLac, do the plots in my updated post answer your question re "the c dependence of w without transforming the equation explicitly"? $\endgroup$
    – kglr
    Nov 11, 2014 at 19:20
0
$\begingroup$ 
r = ParametricNDSolve[{eq, c[0] == 0}, c, {t, 0, 25}, w];

Plot[Flatten[{Re@#, Im@#} &[Table[c[w][t] /. r, {w, 1, 4}]]], {t, 0, 
  25}, Evaluated -> True]
Improve this answer
$\endgroup$
1
  • $\begingroup$ ok, thanks; also i'd like to make a fourier transform for this equation, but somehow FourierTransform[eq, t, [Omega]] is stuck on evaluation (putting a or w instead of omega doesn't work) $\endgroup$
    – Bar Lac
    Nov 11, 2014 at 15:57

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.