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

курсовая работа

Текст програми

program leverie;

uses wincrt;

type matr=array[1..10,1..10] of real;

mas=array[1..10] of real;

var n,i,j,k,l,g,v:integer;

p,s:mas;

a,b,c:matr;

stiykist:real;

procedure mnogmatr(a,b:matr;n:integer; var c:matr);

begin

for i:=1 to n do

for j:=1 to n do

c[i,j]:=0;

for k:=1 to n do

for i:=1 to n do

for j:=1 to n do

begin

c[k,i]:= c[k,i]+a[j,i]*b[k,j];

end;

end;

begin

writeln(vvedit rozmirnict matrici n=);

readln(n);

writeln(vvedit koeficienti matrici);

for i:=1 to n do

for j:=1 to n do

begin

write(a[,i,,,j,]=);

readln(a[i,j]);

end;

for i:=1 to n do

for j:=1 to n do

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

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.

Делись добром ;)