Blog List View

请教一下为什么ContourPlot出不了图,代码如下,图如下 

2016-03-29 22:09:32
861 次阅读


QQ截图20160329214818.png

k = 2*10^9; d = 10; M = 2.19*10^9; P = 
 2.19*10^9; t = 1; h = 0.025; q = 1000;
w = 0.01; r = 0.1
A = (3*M*P)/(M*P + 3*k*M + 3*k*P);

B = (3*M)/(P + 3*k);

Q = (3*k*M + 3*k*P + 7*M*P)/(M*P + 3*k*M + 3*k*P);
s = (7*M + 6*k)/(P + 3*k);
a = E^((M*P*t + 3*k*M*t + 3*k*P*t)/(3*k*d + P*d))
G = (z - h)/Sqrt[(r^2 + z^2 - 2*h*z + h^2)] - (z - h)/
   Sqrt[(r^2 + h^2 + 2*h*z + z^2)];
H = (z^3 - h^3 - 3*z^2*h + 3*h^2*z)/Power[z^2 - 2*h*z + h^2 + r^2, (
  3)^-1];
U = (z^3 + 2*h*z^2 + z*h^2)/Power[z^2 + 2*h*z + h^2 + r^2, (3)^-1];

J = (18*h^2*z^3 + 18*h^3*z^2 + 6*h^4*z + 6*h*z^4)/Power[
  z^2 + 2*h*z + h^2 + r^2, (5)^-1];
M = 1 - h^3/Power[h^2 + w^2, (3)^-1];
L[r_, z_] = 
  4*q*M*A*G - 4*q*M*a*G*A + 4*q*M*a*G*B + 4*q*M*H + 4*q*M*U*Q - 
   4*q*M*a*U*Q + 4*q*M*s*a*U + 4*q*M*s;
ContourPlot[L[r, z], {r, 0, 0.0001}, {z, 0.025, 0.035}]


发表评论