Курсовая работа: Метод простой итерации для решения систем линейных алгебраических уравнений

flops(0);

t1=clock;

while any(flag)

for i=1:size_system(1)

temp(i)=b(i);

for ii=1:size_system(1)

if (i~=ii)

temp(i)=temp(i)-a(i,ii)*result(number_iter+1,ii);

else temp(i)=temp(i)-(a(i,ii)-1)*result(number_iter+1,ii);

end;

end;

e(i)=abs(result(number_iter+1,i)-temp(i));

if e(i)<=edop

flag(i)=0;

else flag(i)=1;

end;

end;

edop1=[edop1;e];

result=[result;temp];

number_iter=number_iter+1;

end;

t2=clock;

number_oper=flops;

time=etime(t2,t1);

res=result';

v=size(res);

fprintf('total','\nРезультаты итерационного процесса, реализованного вторым способом\n');

for i=1:size_system(1)

fprintf('total','\nX%g равен:',i);