【正文】 間%n—區(qū)間份數(shù)x=linspace(a,b,n+1)。for k=1:n q=q+F(x(k))+4*F((x(k)+x(k+1))/2)+F(x(k+1))。用辛普森公式有：F=inline(39。)。for k=100:100:5000b=binesimpson(F,1/k,1,k)。k,b]。k=0。T(0+1,0+1)=h/2*(subs(f,x,a)+subs(f,x,b))。h=h/2。 for jj=1:2:2^k1 s=s+subs(f,x,a+jj*h)。 for jj=1:k T(kjj+1,jj+1)=4^jj/(4^jj1)*T(kjj+1+1,jj1+1)1/(4^jj1)*T(kjj+1,jj1+1)。q=T(0+1,k+1)。 T(1:jj1,jj)=zeros(jj1,1)。sqrt(x).*log(x)39。 s=quad(F,0,1,1e4)s = 2.【第二題】計算二重積分.（1）若區(qū)域，試分別用復合辛普森公式及Gauss求積公式（取n=4）求積分。q=0。endq=q*(ba)/n/6。F=exp(x.*y)。x39。F=subs(F,x0,(ba)/2*t+(a+b)/2)。q=(ba)/2*sum(A.*subs(F,t,x))。F=exp(x.*y)。F=exp(x.*y)。m=size(A)。u=zeros(1,n)。 for jj=k+1:n if abs(A(jj,k))abs(A(ik,k)) ik=jj。det=0。無解39。return。A(k,jj)=A(ik,jj)。 end v=b(k)。b(ik)=v。 end for ii=k+1:n A(ii,k)=A(ii,k)/A(k,k)。 end b(ii)=b(ii)A(ii,k)*b(k)。det=0。無解39。return。for ii=n1:1:1 s=0。 end b(ii)=(b(ii)s)/A(ii,ii)。function [A,b,det]=liezhuyuan(A,b)det=1。n=m(1)。for k=1:n1 ik=k。 end end if A(ik,k)==0 x=0。disp(39。)。 end if ik~=k for jj=k:n u(jj)=A(k,jj)。A(ik,jj)=u(jj)。b(k)=b(ik)。det=det。 for jj=k+1:n A(ii,jj)=A(ii,jj)A(ii,k)*A(k,jj)。 endendif A(n,n)==0 x=0。disp(39。)。endb(n)=b(n)/A(n,n)。 for jj=ii+1:n s=s+A(ii,jj)*b(jj)。enddet=A(n,n)*det。n=m(1)。tt=s。if det(A)==0 disp(39。)。endfor r=1:n for ii=r:n A(ii,r)=A(ii,r)A(ii,r:n)*A(r:n,r)。A(r,:)=A(ir,:)。 end A(r,r)=s(r)。 u=0。 end A(r,ii)=A(r,ii)u。 if ii~=u v=b(u)。b(ii)=v。 for k=1:ii1 u=u+A(ii,k)*b(k)。endb(n)=b(n)/A(n,n)。 for k=ii+1:n u=u+A(ii,k)*b(k)。end第六章1.【第一題】給出線性方程組，其中系數(shù)矩陣為Hilbert矩陣，假設，.若取n=6,8,10，分別用Jacobi迭代法及SOR迭代（=1,）求解，比較計算結果。 norm(A*xb)]。p=h(2)。 for t=1:p if t~=s w=w+A(s,t)*x(t)。 end x=z39。k x39。endif flag Lend順便編制Gauss ：function x=gauss_seidel_iteract(A,b,x,n,flag)%Gauss Seidel迭代求線性方程Ax=b的解%x迭代初值%n迭代次數(shù)%flag控制過程是否顯示L=[0 x39。h=size(A)。for k=1:n for s=1:p w=0。 end for t=s+1:p w=w+A(s,t)*x(t)。 end L=[L。 norm(A*xb)]。 norm(A*xb)]。p=h(2)。 for t=1:p w=w+A(s,t)*x(t)。 end L=[L。 norm(A*xb)]。H=hilb(k)。x=*ones(k,1)。Jacobi迭代結果39。J=jacobiiteract(H,b,x,10,1)for v=1:.25:disp(39。)。Gauss Seigel迭代結果39。G=SOR(H,b,x,v,10,1)endend執(zhí)行 xiti6_1結果略第七章1.【第一題】求下列方程的實根：（1）；（2）.要求（1）設計一種不動點迭代法，要使迭代序列收斂，然后用Steffensen加速迭代，計算到為止。y1=x.^23*x+2。 plot(x,y1,x,y2) x=10:.1:6。y2=exp(x)。y1=x.^23*x+2。plot(x,y1,x,y2) x=:.00005:。y2=exp(x)。grid on x0=::。 f1=sqrt(3*x2+exp(x))。 f3=log(x^23*x+2)。 f5=(exp(x)2)/x+3。 y2=subs(diff(f2),x,x0)。 y4=subs(diff(f4),x,x0)。 subplot(3,1,1)。plot(x0,y2) subplot(3,1,3)。plot(x0,y4) subplot(2,1,2)。根據(jù)作圖。n=1。while abs(yx0)tol x0=y。n=n+1。L=[L。endxx=y。 xxxx = nn = 10 L