PDF / GAUSS
.pdf1: program GAUSS; 2:
3:const
4:maxn = 10;
5:
6:type
7:matrix = array [1..maxn,1..maxn] of real;
8:vector = array [1..maxn] of real;
9:
10:var
11:ns,i,j: integer;
12:as: matrix;
13:xs,bs: vector;
15:procedure gauss_method(N: integer; var A: matrix; var B,X: vector );
16:var
17:m:real;
18:i,j,k,l:integer;
19:begin
20:for k:=1 to n-1 do
21:for i:=k+1 to n do begin
22:m:=a[i,k]/a[k,k];
23:b[i]:=b[i]-m*b[k];
24:for j:=k to n do a[i,j]:=a[i,j]-m*a[k,j];
25:end;
26:for j:=n downto 1 do begin
27:m:=b[j];
28:for l:=n downto j+1 do m:=m-a[j,l]*x[l];
29:x[j]:=m/a[j,j];
30:end;
31:end;
32:
33:begin
34:WriteLn('-- Решение СЛАУ методом Гаусса --');
35:Write(' Порядок СЛАУ : ');
36:ReadLn(ns);
37:for i:=1 to ns do
38:for j:=1 to ns do
39:begin
40:Write(' Введите : A(',i,',',j,') = ');
41:ReadLn(as[i,j]);
42:end;
43:for j:=1 to ns do
44:begin
45:Write(' Введите : B(',j,') = ');
46:ReadLn(bs[j]);
47:end;
48:gauss_method(ns,as,bs,xs);
49:for j:=1 to ns do
50:WriteLn(' Решение : X(',j,') = ',xs[j]:3:0);
51:ReadLn;
52:end.
53:
54:
55:
56: