数值分析上机答案computer中科院.docx

源程序如下：function [x_jac,x_gauss,x_sor,A,b]=diedai(n,w,e,N,Ain,bin)if nargin3 e=1e-5;N=1000;endif nargin4 N=1000;endif nargin2 error(2?êy2?1?);endA=hilb(n);b=A*ones(n,1);if nargin4 if nargin==5 error(??ê?è?b); else A=Ain; b=bin; endend %%????±èμü′údiedai_n=0;x=zeros(n,2);while max(abs(b-A*x(:,2)))e diedai_n=diedai_n+1; if diedai_nN disp(????±èμü′ú?′′?μ????è); break; end x(:,1)=x(:,2); for i=1:n if i1 x(i,2)=(b(i)-A(i,1:(i-1))*x(1:(i-1),1)-A(i,(i+1):n)*x((i+1):n,1))/A(i,i); else x(i,2)=(b(i)-A(i,(i+1):n)*x((i+1):n,1))/A(i,i); end end end x_jac.x=x(:,2);x_jac.n=diedai_n-1;x_jac.error=max(abs(b-A