摘要： (代码在后面) m=80,n=100; MATLAB代码: clc clear M=80;%x N=80;%T L=1;T=1; h=L/M;%空间步长 tao=T/N;%时间步长 x=0:h:L; t=0:tao:T; r=L*(tao/h)^2 f=@(x,t) 0; %% 构造边值矩阵 u=z 阅读全文