第一张图是求解模型所用的方程组，我想得到在那个矩形区域内的温度分布图（ContourPlot）。
第二张图是我在mma上的代码，但是它提示我说方程个数大于了未知量数，删掉1个条件又出现错误，求同志们帮忙。
代码在最下方（可复制）。
谢谢。
a = 0.5*10^-3; b = 20*10^-3; \[Rho] =
11.3*10^3; c = 130; \[Lambda] = 300; i = 10; R = 0.01;
H = 15;
fun = \[Rho]*c*D[T[r, z, t], t] - \[Lambda]*
Laplacian[T[r, z, t], {r, z}] - (i^2 R)/(\[Pi]*a^2*b);
ic = T[r, z, 0] == 300.15;
bc1 = fun == NeumannValue[H*(T[r, z, t] - T[r, z, 0]), z == 0];
bc2 = fun == NeumannValue[-H*(T[r, z, t] - T[r, z, 0]), r == a];
bc3 =
fun == NeumannValue[-H*(T[r, z, t] - T[r, z, 0]), z == b];
sol = NDSolveValue[{bc2, bc1, bc3, ic},
T[r, z, t], {r, 0, a}, {z, 0, b}]
Clear["Global`*"]

你对NeumannValue的使用方法有误，NeumannValue必须要是方程的一部分，具体请参看自带帮助，注意“更多信息”部分也要看。