unit2.pas

来自「1、说明： 本书中所有的常用数值算法子过程按书中的章数分别放在以C开头」· PAS 代码 · 共 43 行

unit Unit2;

interface
uses
  Windows, Messages, SysUtils, Classes, Graphics, Controls,unit1, Forms, Dialogs;
procedure SSOR(A:matrx2; N:integer; B:array of real;
                      var X:array of real; EPS,OM:real;var II:integer);

implementation
procedure SSOR(A:matrx2; N:integer; B:array of real;
                      var X:array of real; EPS,OM:real;var II:integer);
var
    IMAX,J,I,L:integer; RX,R:real;
begin
    IMAX:=200;
    For I:=1 To N do
    begin
        R:=1 / A[I, I];
        B[I]:=B[I] * R;
        For J:=1 To N do
            A[I, J]:=A[I, J] * R;
    end;
    For L:=1 To IMAX do
    begin
        RX:=0;
        For I:=1 To N do
        begin
            R:=B[I];
            For J:=1 To N do
                R:=R - A[I, J] * X[J];
            If Abs(R) > RX Then RX:=Abs(R);
            X[I]:=X[I] + OM * R;
        end;
        If OM * RX <= EPS Then
        begin
            II:=J;
            Exit;
        end;
    end;
    ShowMessage('  Too many iterations');
end;
end.