fpar.m

利用龙格库塔法（等步长

function rhs=fpar(t,u);
global N hx
if t<=1&&t>=0 
    b=1;
else
    b=0;
end
rhs=zeros(N,1);
rhs(1)=(b-2*u(1)+u(2))/(hx*hx);
rhs(2:N-1)=(u(1:N-2)-2*u(2:N-1)+u(3:N))/(hx*hx);
rhs(N)=(2*u(N-1)-2*u(N))/(hx*hx);