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

k = 2*10^9; d = 10; M = 2.19*10^9; P =
. ^& t( N" `! S+ p+ ^0 O
2.19*10^9; t = 1; h = 0.025; q = 1000;: |' ~& q7 r! s
w = 0.01; r = 0.15 X8 r3 v( Y4 M# v! |
A = (3*M*P)/(M*P + 3*k*M + 3*k*P);! L' {3 x' J" B; C5 y$ U
- @9 c7 K6 ^$ K3 y& Z7 F
B = (3*M)/(P + 3*k);
. X$ U+ r' P: H4 i+ I9 d$ w
" R9 a: D* @! ~! M; n
Q = (3*k*M + 3*k*P + 7*M*P)/(M*P + 3*k*M + 3*k*P);/ `& W$ m4 i% [9 _5 N+ c. J/ {
s = (7*M + 6*k)/(P + 3*k);9 r% R; i! c }7 E" N- l0 ]
a = E^((M*P*t + 3*k*M*t + 3*k*P*t)/(3*k*d + P*d))4 U) f: f L r5 _& D
G = (z - h)/Sqrt[(r^2 + z^2 - 2*h*z + h^2)] - (z - h)/6 ^7 b: l- K) U; A- \
Sqrt[(r^2 + h^2 + 2*h*z + z^2)];" `/ | C& C, o" T7 A
H = (z^3 - h^3 - 3*z^2*h + 3*h^2*z)/Power[z^2 - 2*h*z + h^2 + r^2, (4 z; S, x* z; _4 [; Y- B
3)^-1];
2 \5 P0 Z, [) n" D8 R
U = (z^3 + 2*h*z^2 + z*h^2)/Power[z^2 + 2*h*z + h^2 + r^2, (3)^-1];. O1 C& a1 w( p7 r4 n( l4 f
" c, L! i+ P- {6 f& `0 m. l
J = (18*h^2*z^3 + 18*h^3*z^2 + 6*h^4*z + 6*h*z^4)/Power[2 Y6 k) o( C. B# P. Y
z^2 + 2*h*z + h^2 + r^2, (5)^-1];
/ P7 t* v4 B i- O7 d/ r" h5 z% m
M = 1 - h^3/Power[h^2 + w^2, (3)^-1];
K- q) O8 U) B3 z3 B
L[r_, z_] = ( G2 r0 [3 ]/ Z, i, C
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 - 3 v" H" q7 h0 E7 X
4*q*M*a*U*Q + 4*q*M*s*a*U + 4*q*M*s;
' ?! X0 n. J7 M6 F
ContourPlot[c[r, z], {r, 0, 0.0001}, {z, 0.025, 0.035}]

复制代码