1、精確解：p=1,5,10; x=0:1/19:1;for i=1:1:3 for j=1:1:20 y(i,j)=(exp(p(1,i)*19*x(1,j)-1)/(exp(p(1,i)*19)-1); end plot(x,y(i,:); hold on;end由題對中心差分、一階迎風(fēng)、混合格式進(jìn)行模塊編程：他們之間可以通用，只需更改ae關(guān)于p的函數(shù)即可：程序如下:(1)中心差分p=1,5,10;for i=1:1:3 ae=1-0.5*p(1,i); aw=p(1,i)+ae; ap=ae+aw;for i=1:1:18 for j=1:1:20 a(i,j)=0; endendfor i

2、=1:1:18 j=i; a(i,j)=aw; a(i,j+1)=-ap; a(i,j+2)=ae;endfor i=1:1:17 n=i+1; for m=i:-1:1 b(1,1)=a(m,n); a(m,n)=-a(i+1,n)/a(i+1,n)*b(1,1)+a(m,n); a(m,n+1)=-a(i+1,n+1)/a(i+1,n)*b(1,1)+a(m,n+1); a(m,n+2)=-a(i+1,n+2)/a(i+1,n)*b(1,1)+a(m,n+2); endendF(1)=0;F(20)=1;F(19)=(-a(1,20)*F(20)-a(1,1)*F(1)/a(1,19);f

3、or i=2:1:18 F(i)=(-a(i,20)*F(20)-a(i,19)*F(19)/a(i,i);endx=0:1/19:1;y(1,:)=F;plot(x,y);hold onend（2）一階迎風(fēng)更改ae表達(dá)式：if p(1,i)2 ae=0;else if (p(1,i)0) ae=1-p(1,i)/2;else ae=-p(1,i);end（4）quick格式 用matlab編quick格式：程序如下：p1=1,5,10;for i=1:1:3p=p1(1,i);F(1,1)=0;F(1,20)=1;for i=2:1:19 F(1,i)=0;endF1(1,1)=0;F1(1

4、,20)=1;for i=2:1:19 F1(1,i)=0;endwhile (F1(1,2)=0)|max0.00000001for i=1:1:20 F(1,i)=F1(1,i);endA=1;B=1/(2+p);C=(1+p)/(2+p);D(1,2)=p*(6*F(1,2)-3*F(1,1)-3*F(1,3)/8)/(2+p);for i=3:1:19 D(1,i)=p*(5*F(1,i)-F(1,i-1)-F(1,i-2)-3*F(1,i+1)/8)/(2+p);endP(1,2)=B;Q(1,2)=C*F(1,1)+D(1,2);for i=3:1:19 P(1,i)=B/(A-C*P(1,i-1); Q(1,i)=(D(1,i)+C*Q(1,i-1)/(A-C*P(1,i-1);endfor i=19:-1:2 F1(1,i)=P(1,i)*F1(1,i+1)+Q(1,i);endmax=abs(F1(1,2)-F(1,2);for i=3:1:19 if abs(F1(1,i)-F(1,i)max