复习市公开课特等奖市赛课微课一等奖PPT课件.pptx

第0章绪论第一章插值方法第二章数值积分第三章常微分方程差分方法第四章方程求根迭代法第六章线性方程组直接法考试安排上机作业functionA=gauss(c,b)
k=length(b);
c=[cb(:)];
forii=1:k
[temp,j]=max(abs(c(ii:k,ii)));
j=j+ii-1;
if~(j==ii)
temp=c(ii,ii:k+1);
c(ii,ii:k+1)=c(j,ii:k+1);
c(j,ii:k+1)=temp;
end
ifc(ii,ii)==0
disp('对角线上元素为零，程序退出！');
pause;
return;
end
forn=ii+1:k
c(n,ii+1:k+1)=c(n,ii+1:k+1)-c(ii,ii+1:k+1).*c(n,ii)/c(ii,ii);
end
end>>c=[234;432;121]
c=
234
432
121
>>b=[484]'
b=
4
8
4
>>x1=inv(c)*b
x1=
1.00000000000000
2.00000000000000
-1.00000000000000functiony=simpson(x,h)
N=length(x);
ifrem(N,2)==1
n=(N-1)/2;
Flag=1;
else
n=(N-2)/2;
Flag=0;
end
y=h/3*(x(1)+4*sum(x(2:2:2*n))+2*sum(x(3:2:2*n-1))+x(2*n+1));
ifFlag==0
y=y+h/2*(x(N-1)+x(N));
end
return;>>h=0.01;
>>t=0:h:1;
>>x=t.^2;
>>y=simpson(x,h)

y=
0.33333333333333

>>x=t.^3;
>>y=simpson(x,h)

y=
0.25000000000000