最后几句改了下:
fprintf('最后时刻数值解与精确解分别为:\n');disp(u);
fprintf('差分法得到的结果与正确结果的最大相对误差:');
disp([num2str(E*100) '%']); %百分号前要有空格
plot(x0,u(:,1),'r',x0,u(:,2),'b-');
legend('数值解','精确解')
xlabel('x');ylabel('u(x,t)') %中间要分号隔开
title('最后时刻热传导问题数值解与精确解比较')
调整之后还有错,前面的t没有定义,自己改改吧。
祝顺利!
本文如未解决您的问题请添加抖音号:51dongshi(抖音搜索懂视),直接咨询即可。