学习笔记21 有限差分法求解热传导方程的代码

距离9月份全国大学生数学建模竞赛仅剩10天，但目前的准备工作进展较为缓慢。在有限差分法这一部分内容上，我可能已经停滞了一个星期，不过值得庆幸的是，现在并非毫无思路，对算法原理已有一定掌握，因此在阅读相关代码时也变得更加顺畅了。
接下来的准备安排是，在完成2018A题之后将正式进入冲刺阶段。在强化训练期间，我们学习了多种模型，并研读了一些论文资料，目前已具备一定的知识储备，包括评价模型、预测模型以及微分方程模型等，这使得我们有能力应对部分问题。计划在9月4日至6日期间选择一个赛题进行练习，以发现并弥补自身的不足之处。随后再进一步明确自身存在的薄弱环节。
在此不再赘述，以下为MATLAB代码（同样来源于B站视频内容，已将其整理下来，并与学习笔记20所涉及的内容属于同一套体系）。

    clc,clear
    a=1;
    dx=0.02;
    x=0:dx:1;
    dt=0.0001;t=0:dt:1;
    u=zeros(length(x),length(t));
    %行数和x的个数一样
    u(:,1)=sin(pi*x);%x的函数，初始函数，第一列
    m1=0+0.1*sin(t);
    m2=0-0.1*sin(t);
    A=-2*eye(length(x))+diag(ones(1,length(x)-1)