Курсовая работа: Стійкість системи лінійних алгебраїчних рівнянь

for i:=1 to n do

s[1]:=s[1]+a[i,i];

p[1]:=s[1];

for l:=2 to n do

begin

for i:=1 to n do

for j:=1 to n do

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

mnogmatr(a,b,n,c);

for j:=1 to n do

s[l]:=s[l]+c[j,j];

for k:=1 to l do

p[l]:=p[l]-p[k]*s[l-k];

p[l]:=(s[l]+p[l])/l;

end;

writeln('haracteristichniy mnogochlen');

g:=n; v:=0;

repeat

write(p[g]:2:3,'*l^',v,'+');

g:=g-1; v:=v+1;

until g=0;

writeln('-l^',v,'=0');

for i:=1 to n do

stiykist:=stiykist+abs(p[i]);

if stiykist>1 then writeln('sistema ne stiyka') else writeln('sistema stiyka');

end.

Приклад

Перевірити на стійкість систему Аx=B

-0.77 -0.44 0.21 -0.18 -1,24