第一个
c1=0;
c2=0;
n=6;
m=11;
a=1;
b=0.2;
c=1;
x=0:.2:1%（没有也可以）


h=a/(n-1);
k=b/(m-1);
r=c^2*k/h^2;
s=1-2*r;
U=zeros(n,m);
U(1,1:m)=c1;
U(n,1:m)=c2;
U(2:n-1,1)=f(h:h:(n-2)*h); %feval(f,h:h(n-2)*h)';
for j=2:m
    for i=2:n-1
        U(i,j)=s*U(i,j-1)+r*(U(i-1,j-1)+U(i+1,j-1));
    end
end
U=U'


第二个
function U=forwdif2(r,c1,c2,a,b,c,n,m)
x=0:0.2:1;
h=a/(n-1);
k=b/(m-1);
r=c^2*k/h^2;
s=1-2*r;
U=zeros(n,m);
U(1,1:m)=c1;
U(n,1:m)=c2;
U(2:n-1,1)=r(h:h:(n-2)*h)%feval('r',h:h:(n-2)*h);
for j=2:m
    for i=2:n-1
        U(i,j)=s*U(i,j-1)+r*(U(i-1,j-1)+U(i+1,j-1));
    end
end
U=U'


function Y=r(x)
Y=4.*x-4.*x.^2