lab41.m

利用龙格库塔法（等步长

clear all
global N hx
t0=0;
tend=1;

N=10;
hx=1/N;
x=hx:hx:1-hx;
u0=zeros(1,N);
u0=u0';
    result=[u0'];
    time=[t0];
    u=u0;
    t=t0;
    M=200;
    for k=1:M
        h=(tend-t0)/M;
        u=u+h*fpar1(t,u);
        t=t+h;
        result=[result;u'];
        time=[time;t];
    end
      for i=1:size(time)
        if time(i)<=1&&time(i)>=0
            A(i)=1;
        else
            A(i)=0;
        end
    end



resultp=[A' result];
xp=0:hx:1
mesh(xp,time,resultp)