Advertisement
Guest User

Untitled

a guest
Nov 3rd, 2016
214
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.42 KB | None | 0 0
  1. program project1;
  2. uses Crt;
  3. Type
  4. Matrix = array[1..20,1..20] of real;
  5. Vec = array [1..20] of real;
  6. function det(A: Matrix; N: Integer): real;
  7. var i,j,k:Integer;
  8. t,q,r:Real;
  9. begin
  10. r:=1;
  11. for k:=1 to n do
  12. begin
  13. q:=0;
  14. for i:=k to n do
  15. begin
  16. t:=a[i,k];
  17. if abs(t) > abs(q) then
  18. begin
  19. q:=t;
  20. j:= i;
  21. end
  22. end;
  23. if j<>k then
  24. begin
  25. r:=-r;
  26. for i:=k to n do
  27. begin
  28. t:=a[j,i];
  29. a[j,i]:=a[k,i];
  30. a[k,i]:=t
  31. end;
  32. end;
  33. for i:= k+1 to n do
  34. begin
  35. t:= a[i,k]/q;
  36. for j:= k+1 to n do
  37. a[i,j]:= a[i,j]-t*a[k,j]
  38. end;
  39. R:=R*a[k,k]
  40. end;
  41. det:=r;
  42. end;
  43.  
  44.  
  45.  
  46. var
  47. a:matrix;
  48. b:vec;
  49. c:array[1..20]of matrix;
  50. n,i,j,k:integer;
  51. z,z1,x,temp:real;
  52.  
  53.  
  54.  
  55. begin
  56. write('N=');
  57. readln(n);
  58. writeln('Koeficienti');
  59. for i:=1 to n do
  60. for j:=1 to n do
  61. read(a[i,j]);
  62. readln;
  63. writeln('Chleni');
  64. for j:=1 to n do
  65. read(b[j]);
  66. readln;
  67. write(a[i,j]:7:2);
  68. write(b[i]:9:2);
  69. writeln;
  70. z := Det(a,n);
  71.  
  72.  
  73. for i:=1 to n do
  74. begin
  75. for k:=1 to n do
  76. begin
  77. for j:=1 to n do
  78. c[i][k,j]:=a[k,j];
  79. c[i][k,i]:=b[k];
  80. write(c[i][k,j]:7:2);
  81. writeln;
  82. end;
  83.  
  84. writeln('Det= ',Det(c[i],n):7:2);
  85.  
  86. z1:= Det(a,n) ;
  87. x:=z1/z;
  88. writeln('x(', i, ') = ', x:3:0);
  89.  
  90. for j:=1 to n do
  91. begin
  92. temp := b[j];
  93. b[j] := a[i, j];
  94. a[i, j] := temp;
  95. end;
  96. end;
  97. repeat until keypressed;
  98.  
  99. end.
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement