# Why won't mathematica plot the 3D Contour of this function?

I want to find the 3D contour of the equation1:

C1 = 10^(-10);
C2 = 0.1*C1;
R = 50;
Tb = 0.1;
Geb = 5*10^-15;
Z0 = 50;
L[Te_] := 10^-9 + 10^-9*(Te - 0.1);
Zlcr[Te_, w_] := (1/R + 1/(I*L[Te]*w) + I*C1*w)^-1;
Zload[Te_, w_] := -I*w*C2 + Zlcr[Te, w];
y[Te_, w_] := (Abs[\[CapitalGamma][Te, w]])^2;
eqn1[w_, Te_, Pprobe_] := (1 -   y[Te, w]) Pprobe  ==  (Te - Tb) Geb
ContourPlot3D[
eqn1[w, Te, Pprobe], {w, 0, 5*10^9}, {Te, 0, 1}, {Pprobe, 0, 10^-14}]

• I don't understand this part: == (Te - Tb) Geb! Why do use == inside definition of a function? Jul 15, 2014 at 10:05
• (Te-Tb)Geb is simply equal to (Te - 0.1)* (5*10^(-15)) Jul 15, 2014 at 10:07
• Yes that is correct but I think this definition (1 - y[Te, w]) Pprobe == (Te - Tb) Geb is wrong. Are (1 - y[Te, w]) Pprobe and (Te - Tb) Geb equal? Jul 15, 2014 at 10:09
• Use Evaluate[eqn1[w, Te, Pprobe]] inside ContourPlot3D. Jul 15, 2014 at 12:22

exp = (1 - y[Te, w]) Pprobe == (Te - Tb) Geb;