求助老哥们出现这个错误怎么解决啊

主要每次重启matlab,第一次运行都是可以出正常结果的，(如图),但是第二次再运行就出错了，这一般是什么原因啊

用theta格式解一维抛物问题，显示错误出现在这里

从你给出的图片代码，不连贯，一下子看不出问题，可以将代码发我看一下吗？

老哥这样可以吗，还是复制代码发给你

function U = Theta(a,b,M,N,phi,alpha,beta)
h = 1 / M; tau = 1 / N; r = a * tau / h / h;
M1 = zeros(M-1, M-1); M2 = zeros(M-1, M-1); M3 = zeros(M-1, 1); U = zeros(N+1, M-1);
for i=1:M-1 M1(i,i) = 1 + 2*b*r; M2(i,i) = 1 - 1-2*(1-b)*r; U(1,i) = phi(h*i);
if i-1 >= 1 M1(i, i-1) = - b*r; M2(i, i-1) = (1-b)*r; end
if i+1 <= M-1 M1(i, i+1) = -b*r; M2(i, i+1) = (1-b)*r; end end
for k = 1:N M3(1) = (1-b)*r* alpha(tau*(k-1)) + b*r*alpha(tau*k) ; M3(end) = (1-b)*r* beta(tau*(k-1)) + b*r*beta(tau*k) ;
temp = inv(M1) * M2 * U(k,:)' + inv(M1) * M3 ; U(k+1,:) = temp'; endend

a = 1;b = 0.4;M = 40;N = 40;
syms x t phi(x) = exp(x); alpha(t) = exp(t); beta(t) = exp(1+t);
U = Theta(a,b,M,N,phi,alpha,beta);for j = 0.2:0.2:0.8 temp = U(end, round(j * M)); fprintf("输出点（%0.1f,1.0)的数值解为：%f\n",j, temp)end

你的程序没有问题啊。能正常运行（运行环境matlab2022a）。

我可以解决，丝我