偏微分方程数值解迎风格式代码

差分格式是数值计算方法中微分以及偏微分导数的一种离散化方法，即用相邻两个或者多个数值点的差分取代偏微分方程中导数或者偏导数的一种算法。 选择差分格式是离散化偏微分方程的第一步。本代码是迎风格式的代码。具体参加李荣华的数值解书目。

clear
h=0.1;
t=0.05;
a=2/h;
b=1/t;
r=t/h;
u=zeros（a+1b+1）;
for i=1:a+1
    if （（i-1）*h>0.35）&&（（i-1）*h<0.65）
        u（i1）=1;
    else
        u（i1）=0;
    end
end
for k=1:b
    u（1k