A=input('Введите матрицу A='); b=input('Введите вектор b='); eps=input('Укажите точность вычислений eps='); kmax=input('Укажите предельное количество итераций kmax='); n=length(b); fprintf('\n Матрица коэффициентов СЛАУ (Матрица A) \n'); for i=1:n fprintf('%6.2f',A(i,:)); fprintf('\n'); end fprintf('\n Вектор правых частей СЛАУ (Вектор b) \n') fprintf('%6.2f \n',b); x=zeros(n,1); for k=1:kmax z=zeros(n,1); for i=1:n s(i)=b(i); for j=1:n if(i~=j) s(i)=s(i)-A(i,j)*x(j); end end s(i)=s(i)/A(i,i); z(i)=z(i)+abs(x(i)-s(i)); x(i)=s(i); end if(max(z)<eps) x=x1; break end end fprintf('\n Решение СЛАУ по методу Зейделя \n') fprintf('%6.2f \n',x); fprintf('\n Число итераций k=%3d \n',k);
|