0
$\begingroup$

I am trying to plot the roots of the following equation as mentioned in the code below But it is not giving me output. If anyone can resolve this problem for me is most welcome.

del = -1.5;
del0 = 1.5;
ome = 2;
A1 = 0;
B1 = 1;
g0 = 2;
G1 = 5;
k1 = 0.1;
kex = 0.1;
X1 = g0/2*(1 - del0/ome)*A1;
Y1 = P0/Sqrt[2]*g0/ome*B1;
C1 = (G1/2)^2 + ome^2;
kL = (k1 + kex)/2 - del0*(k1 - kex)/(2*ome);
sol = Solve[{I*del*a1 + 
   I*g0*(1 - del0/ome)*A1*(X1*Dagger[a1]*a1 + Y1*Re[a1])/C1*ome + 
   I*P0/Sqrt[2]*(1 - del0/(2*ome)) + 
   I*P0*g0*B1/(Sqrt[2]*ome)*(X1*Dagger[a1]*a1 + Y1*Re[a1])/C1*
    ome - kL*a1/2 - (k1 - kex)*g0*
    B1/ome*(X1*Dagger[a1]*a1 + Y1*Re[a1])/C1*a1 == 0}, 
a1]; // FullSimplify
 Plot[Evaluate@(a1 /. sol), {P0, 0, 2}, 
 PlotStyle -> (Directive[Thick, ColorData[97][#]] & /@ {1, 3, 2}), 
  PlotLegends -> Placed[Automatic, {0.5, 0.625}]]
$\endgroup$
  • $\begingroup$ Try replacing Plot[Evaluate@(a1/.sol)...] with Plot[ReIm[Evaluate@(a1/.sol)]...] and you should see a plot of the real and imaginary components. $\endgroup$ – Bill Feb 10 at 8:34
  • $\begingroup$ @Bill In my code how can i write a1 dagger in mathmatica? If you can edit and run my code please. $\endgroup$ – vini Feb 10 at 18:36
  • $\begingroup$ Does Print[a1*Dagger[a1]/.sol] give you what you are asking for? If not then please try to explain more details and I will see what I can do. $\endgroup$ – Bill Feb 12 at 19:19

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.