ssor.txt

VB常用数值算法集 内含有解线性代数方程组

Sub SSOR(A(), N, B(), X(), EPS, OM, II)
    IMAX = 200
    For I = 1 To N
        R = 1# / A(I, I)
        B(I) = B(I) * R
        For J = 1 To N
            A(I, J) = A(I, J) * R
        Next J
    Next I
    For II = 1 To IMAX
        RX = 0#
        For I = 1 To N
            R = B(I)
            For J = 1 To N
                R = R - A(I, J) * X(J)
            Next J
            If Abs(R) > RX Then RX = Abs(R)
            X(I) = X(I) + OM * R
        Next I
        If OM * RX <= EPS Then Exit Sub
    Next II
    Print "  Too many iterations"
End Sub