polcof.txt

本书共有数值计算中常用的VisualBasic子过程近200个

Sub POLCOF(XA(), YA(), N, COF())
    Dim X(15), Y(15)
    For J = 1 To N
        X(J) = XA(J)
        Y(J) = YA(J)
    Next J
    For J = 1 To N
        Call POLINT(X, Y, N + 1 - J, 0#, COF(J), DY)
        XMIN = 1E+38
        K = 0
        For I = 1 To N + 1 - J
            If Abs(X(I)) < XMIN Then
                XMIN = Abs(X(I))
                K = I
            End If
            If X(I) <> 0# Then Y(I) = (Y(I) - COF(J)) / X(I)
        Next I
        If K < N + 1 - J Then
            For I = K + 1 To N + 1 - J
                Y(I - 1) = Y(I)
                X(I - 1) = X(I)
            Next I
        End If
    Next J
End Sub