Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- program project1;
- uses Crt;
- Type
- Matrix = array[1..20,1..20] of real;
- Vec = array [1..20] of real;
- function det(A: Matrix; N: Integer): real;
- var i,j,k:Integer;
- t,q,r:Real;
- begin
- r:=1;
- for k:=1 to n do
- begin
- q:=0;
- for i:=k to n do
- begin
- t:=a[i,k];
- if abs(t) > abs(q) then
- begin
- q:=t;
- j:= i;
- end
- end;
- if j<>k then
- begin
- r:=-r;
- for i:=k to n do
- begin
- t:=a[j,i];
- a[j,i]:=a[k,i];
- a[k,i]:=t
- end;
- end;
- for i:= k+1 to n do
- begin
- t:= a[i,k]/q;
- for j:= k+1 to n do
- a[i,j]:= a[i,j]-t*a[k,j]
- end;
- R:=R*a[k,k]
- end;
- det:=r;
- end;
- var
- a:matrix;
- b:vec;
- c:array[1..20]of matrix;
- n,i,j,k:integer;
- z,z1,x,temp:real;
- begin
- write('N=');
- readln(n);
- writeln('Koeficienti');
- for i:=1 to n do
- for j:=1 to n do
- read(a[i,j]);
- readln;
- writeln('Chleni');
- for j:=1 to n do
- read(b[j]);
- readln;
- write(a[i,j]:7:2);
- write(b[i]:9:2);
- writeln;
- z := Det(a,n);
- for i:=1 to n do
- begin
- for k:=1 to n do
- begin
- for j:=1 to n do
- c[i][k,j]:=a[k,j];
- c[i][k,i]:=b[k];
- write(c[i][k,j]:7:2);
- writeln;
- end;
- writeln('Det= ',Det(c[i],n):7:2);
- z1:= Det(a,n) ;
- x:=z1/z;
- writeln('x(', i, ') = ', x:3:0);
- for j:=1 to n do
- begin
- temp := b[j];
- b[j] := a[i, j];
- a[i, j] := temp;
- end;
- end;
- repeat until keypressed;
- end.
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement