function A=fun1(t,h,T)
%上機實習一顯式格式;
%T是時間，t迭代的時間步長。h空間步長。
r=t^2/h^2;
D=100;E=5000;F=3000;L=1;
N=T/t+2;
J=L/h+2;
A(2,J)=0;                         %另第一,二行為零
for i=3:N
    A(3,1)=0;
    A(3,end)=0;
    for k=2:J-1
        A(3,k)=r*D*A(3-1,k+1)+(2-2*r*D)*A(3-1,k)+r*D*A(3-1,k-1)-A(3-2,k)+t^2*(E*(k-1)*h*(L-(k-1)*h)+F);
    end
   A(1,:)=A(2,:);
   A(2,:)=A(3,:);
end