请教ContourPlot出不了图这是为什么，代码如下，图如下 - 数学建模社区-数学中国

k = 2*10^9; d = 10; M = 2.19*10^9; P = + f: V( ^& k3 f' V, _" l- {
2.19*10^9; t = 1; h = 0.025; q = 1000;- C2 g) D1 }1 t; h. I
w = 0.01; r = 0.1
& ~4 w0 o1 K4 x" B$ q8 V
A = (3*M*P)/(M*P + 3*k*M + 3*k*P);
4 e. |# Q5 p$ d2 _/ V1 c' P1 n1 P
/ y1 ?2 U6 t( R! a( G" F
B = (3*M)/(P + 3*k);
. `, W2 Q, Q' |' c
L- f( d. a" l
Q = (3*k*M + 3*k*P + 7*M*P)/(M*P + 3*k*M + 3*k*P);% s' @! ~6 N: I, i2 _2 G6 c- y* y
s = (7*M + 6*k)/(P + 3*k);
; |' p6 I+ ], B& [2 p3 N
a = E^((M*P*t + 3*k*M*t + 3*k*P*t)/(3*k*d + P*d))
# K, r: R& U: e% _6 ~ H* u
G = (z - h)/Sqrt[(r^2 + z^2 - 2*h*z + h^2)] - (z - h)/0 S0 P/ K( L$ B- {% u) {- k- A
Sqrt[(r^2 + h^2 + 2*h*z + z^2)];0 H* u" N) ~! Y
H = (z^3 - h^3 - 3*z^2*h + 3*h^2*z)/Power[z^2 - 2*h*z + h^2 + r^2, () ]8 U+ k/ s- J1 q
3)^-1];7 \& o4 O* l) \( W# [6 [5 u
U = (z^3 + 2*h*z^2 + z*h^2)/Power[z^2 + 2*h*z + h^2 + r^2, (3)^-1];6 @' u& b' Y$ {2 _
! b; b8 C7 N/ S
J = (18*h^2*z^3 + 18*h^3*z^2 + 6*h^4*z + 6*h*z^4)/Power[
; o& w/ |4 r) ~
z^2 + 2*h*z + h^2 + r^2, (5)^-1];
4 s! g4 u M# T9 H3 r
M = 1 - h^3/Power[h^2 + w^2, (3)^-1];% j$ A4 _6 B0 w# _
L[r_, z_] =
# u4 X/ {" [) w7 ?# J
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 - 7 i% G) A1 }( W
4*q*M*a*U*Q + 4*q*M*s*a*U + 4*q*M*s;# w9 Y8 f9 ^% E$ g& B7 D" g" F% A
ContourPlot[c[r, z], {r, 0, 0.0001}, {z, 0.025, 0.035}]

复制代码