fork download
copy 
  1. program ideone;
  2. Uses Crt;
  3. const n=4;
  4. begin
  5.  
  6.  
  7. type mass1=array[1..n,1..n] of real;
  8.      mass2=array[1..n] of real;
  9. var  a,s:mass1;b,y,x:mass2;i,j,k:integer;
  10.      f:boolean;
  11.  
  12. procedure vvod;
  13. var i,j:integer;
  14. begin
  15. a[1,1]:=4;a[1,2]:=2;a[1,3]:=6;a[1,4]:=8;
  16. a[2,1]:=2;a[2,2]:=5;a[2,3]:=5;a[2,4]:=10;
  17. a[3,1]:=6;a[3,2]:=5;a[3,3]:=14;a[3,4]:=17;
  18. a[4,1]:=8;a[4,2]:=10;a[4,3]:=17;a[4,4]:=30;
  19. b[1]:=-80;
  20. b[2]:=156;
  21. b[3]:=-304;
  22. b[4]:=221;
  23. writeln('Variant 10:');
  24. for i:=1 to n do begin
  25. for j:=1 to n do
  26. write(a[i,j]:2:0,' ');
  27. writeln('| ',b[i]:2:0);end;
  28. writeln;
  29. end;
  30.  
  31.  
  32. procedure proverka(var f:boolean);
  33. var i,j:integer;
  34. begin
  35. f:=true;
  36. for i:=1  to n do
  37. for j:=1 to n do
  38. if a[i,j]<>a[j,i] then f:=false;
  39. end;
  40.  
  41. procedure matr_S;
  42. var i,j,k:integer;s1:real;
  43. begin
  44. for i:=1 to n do begin
  45. s1:=0;
  46. for k:=1 to (i-1) do
  47. s1:=s1+s[k,i]*s[k,i];
  48. s[i,i]:=sqrt(a[i,i]-s1);
  49. for j:=(i+1) to n do begin
  50. s1:=0;
  51. for k:=1 to (i-1) do
  52. s1:=s1+s[k,i]*s[k,j];
  53. s[i,j]:=(a[i,j]-s1)/s[i,i];
  54. end;end;
  55. writeln('Matrica S:');
  56. for i:=1 to n do begin
  57. for j:=1 to n do
  58. write(s[i,j]:2:0,' ');writeln; end;
  59. end;
  60.  
  61. procedure matr_Y;
  62. var i,k:integer;s1:real;
  63. begin
  64. for i:=1 to n do begin
  65. s1:=0;
  66. for k:=1 to (i-1) do
  67. s1:=s1+s[k,i]*y[k];
  68. y[i]:=(b[i]-s1)/s[i,i];end;
  69. writeln;
  70. writeln('Vspomogatelniy vektor Y:');
  71. for i:=1 to n do begin
  72. write(y[i]:2:0,' ');writeln;end;
  73. end;
  74.  
  75. procedure matr_X;
  76. var i,k:integer;s1:real;
  77. begin
  78. for i:=n downto 1 do begin
  79. s1:=0;
  80. for k:=(i+1) to n do
  81. s1:=s1+s[i,k]*x[k];
  82. x[i]:=(y[i]-s1)/s[i,i];end;
  83. writeln;
  84. writeln('Otvet (vektor X):');
  85. for i:=1 to n do begin
  86. write(x[i]:2:0,' ');writeln;end;
  87. end;
  88.  
  89. begin
  90. clrscr;
  91.  
  92. vvod;
  93. proverka(f);
  94. if f=true then begin
  95. matr_S;matr_y;matr_x;end
  96. else writeln('metod ne primenim dlya dannoy matrici!');
  97. end.
Compilation error #stdin compilation error #stdout 0s 0KB
comments (?)
stdin
copy 
Standard input is empty
compilation info 
prog.pas:2: error: module/unit interface `Crt' could not be imported
stdout
copy 
Standard output is empty
https://ideone.com/dXa4Hs
language:
Pascal (gpc 20070904)
created:
6 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!