Предыдущая страница реферата | 1  2  3  4  5  6  7  8  9  10  11 | Следующая страница реферата

For i := k + 1 to n do begin q := a[i, k] / a[k, k];

For j := 1 to n do

If j = k then a[i, j] := 0 else a[i, j] := a[i, j] - q * a[k, j]; b[i] := b[i] - q * b[k]; end;

end;

{ Вычисляем решение } x[n] := b[n] / a[n, n];

For i := n - 1 downto 1 do begin t := 0;

For j := 1 to n-i do t := t + a[i, i + j] * x[i + j]; x[i] := (1 / a[i, i]) * (b[i] - t); end;

Gauss := true;
End;

Var n, i: Integer; a: Matrix ; b, x: Vector;
Begin

ClrScr;

Writeln('Программа решения систем линейных уравнений по методу
Гаусса');

Writeln;

Writeln('Введите порядок матрицы системы (макс. 10)');

Repeat

Write('>');

Read(n);

Until (n > 0) and (n 4

Введите расширенную матрицу системы
A 1 2 3 4 b

1 3.2 5.4 4.2 2.2 2.6
2 2.1 3.2 3.1 1.1 4.8
3 1.2 0.4 -0.8 -0.8 3.6
4 4.7 10.4 9.7 9.7 -8.4

Результат вычислений по методу Гаусса x1 = 5.0000000000E+00 x2 = -4.0000000000E+00 x3 = 3.0000000000E+00 x4 = -2.0000000000E+00

2.2 Программа решения систем линейных уравнений по методу Зейделя

2.2.1. Постановка задачи. Требуется решить систему линейных алгебраических уравнений с вещественными коэффициентами вида

a11x1 + a12x2 + … + a1nxn = b1 ,

a21x2 + a22x2 + … + a2nxn = b2 ,

. . . . . . . . . . . . . an1x1 + an2x2 + … + annxn = bn

Рекомендуем скачать другие рефераты по теме: скачать доклад на тему, отчет по практике.





Предыдущая страница реферата | 1  2  3  4  5  6  7  8  9  10  11 | Следующая страница реферата