# I want to plot root of the following complex equation

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}]]

• Try replacing Plot[Evaluate@(a1/.sol)...] with Plot[ReIm[Evaluate@(a1/.sol)]...] and you should see a plot of the real and imaginary components. – Bill Feb 10 at 8:34
• @Bill In my code how can i write a1 dagger in mathmatica? If you can edit and run my code please. – vini Feb 10 at 18:36
• 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. – Bill Feb 12 at 19:19