Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

PDF / GAUSS

.pdf
Скачиваний:
42
Добавлен:
07.01.2014
Размер:
29.04 Кб
Скачать

1: 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:

Соседние файлы в папке PDF