eulsum.txt

来自「VB常用数值算法集 内含有解线性代数方程组」· 文本 代码 · 共 22 行

Sub EULSUM(SUM, TERM, JTERM, WKSP())
    Static NTERM As Integer
    If JTERM = 1 Then
        NTERM = 1
        WKSP(1) = TERM
        SUM = 0.5 * TERM
    Else
        TMP = WKSP(1)
        WKSP(1) = TERM
        For J = 1 To NTERM
            DUM = WKSP(J + 1)
            WKSP(J + 1) = 0.5 * (WKSP(J) + TMP)
            TMP = DUM
        Next J
        If Abs(WKSP(NTERM + 1)) <= Abs(WKSP(NTERM)) Then
            SUM = SUM + 0.5 * WKSP(NTERM + 1)
            NTERM = NTERM + 1
        Else
            SUM = SUM + WKSP(NTERM + 1)
        End If
    End If
End Sub