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

k = 2*10^9; d = 10; M = 2.19*10^9; P =
$ c$ Z( b9 q' Q
2.19*10^9; t = 1; h = 0.025; q = 1000;
# W. ~' o' o0 f& T2 q9 e# b
w = 0.01; r = 0.1$ `) D) I) @4 N! L& r! [: Y V7 C' {4 ~
A = (3*M*P)/(M*P + 3*k*M + 3*k*P);
G/ X8 a& L O3 _- A
3 \* n: x0 j9 c1 e# ]# b
B = (3*M)/(P + 3*k);- e+ i( g( D: L6 K
7 W7 J3 q; ^" O) W9 y/ ?8 I
Q = (3*k*M + 3*k*P + 7*M*P)/(M*P + 3*k*M + 3*k*P);" M! t; t, i5 u
s = (7*M + 6*k)/(P + 3*k);$ j1 R* n6 ^' |8 J. l
a = E^((M*P*t + 3*k*M*t + 3*k*P*t)/(3*k*d + P*d))
7 P0 a1 J# K0 k4 d4 D, I8 D) x- [) V
G = (z - h)/Sqrt[(r^2 + z^2 - 2*h*z + h^2)] - (z - h)/
0 E6 E6 \1 ^2 K, J% a+ R
Sqrt[(r^2 + h^2 + 2*h*z + z^2)];
, z% [3 N, T+ @6 A: H* H. i& v
H = (z^3 - h^3 - 3*z^2*h + 3*h^2*z)/Power[z^2 - 2*h*z + h^2 + r^2, (* [! {1 D" q) M7 \/ V
3)^-1];
2 Z& j+ N# O7 R& v( m
U = (z^3 + 2*h*z^2 + z*h^2)/Power[z^2 + 2*h*z + h^2 + r^2, (3)^-1];
& v4 z9 F r- p. M* T- \) e
2 ]* ^+ j2 I" i6 Y
J = (18*h^2*z^3 + 18*h^3*z^2 + 6*h^4*z + 6*h*z^4)/Power[2 d; u* ^& j0 y7 b
z^2 + 2*h*z + h^2 + r^2, (5)^-1];1 x. y. b$ I; ^% Q
M = 1 - h^3/Power[h^2 + w^2, (3)^-1];5 l; {6 M: _! |* f
L[r_, z_] =
% D, t- J4 s* S& s: k
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 - ; Y- Z( W; l7 @* X
4*q*M*a*U*Q + 4*q*M*s*a*U + 4*q*M*s;
2 A: p& T5 k3 O6 \) Z) {8 f* M
ContourPlot[c[r, z], {r, 0, 0.0001}, {z, 0.025, 0.035}]

复制代码