Курсовая работа: Решение системы линейных уравнений методом Гаусса и Жордана-Гаусса

begin

hold:=b[l,j];

b[l,j]:=b[i,j];

b[i,j]:=hold

end;

hold:=w[l];

w[l]:=w[i];

w[i]:=hold

end;

for j:=i1 to n do

begin

t:=b[j,i]/b[i,i];

for k:=i1 to n do

b[j,k]:=b[j,k]-t*b[i,k];

w[j]:=w[j]-t*w[i]

end

end

end;

if b[n,n]=0.0 then error:=true

else

begin

coef[n]:=w[n]/b[n,n];

i:=n-1;

repeat

sum:=0.0;

for j:=i+1 to n do

sum:=sum+b[i,j]*coef[j];

coef[i]:=(w[i]-sum)/b[i,i];

i:=i-1