d11r5.txt

常用的数值算法的VB程序

Private Sub Command1_Click()
    'PROGRAM D11R5
    'Driver for routine AMOEBA
    NP = 3
    MP = 4
    FTOL = 0.000001
    Dim P(4, 3), X(3), Y(4)
    P(1, 1) = 0#: P(1, 2) = 0#: P(1, 3) = 0#
    P(2, 1) = 1#: P(2, 2) = 0#: P(2, 3) = 0#
    P(3, 1) = 0#: P(3, 2) = 1#: P(3, 3) = 0#
    P(4, 1) = 0#: P(4, 2) = 0#: P(4, 3) = 1#
    NDIM = NP
    For I = 1 To MP
        For J = 1 To NP
            X(J) = P(I, J)
        Next J
        Y(I) = FAMOEB(X)
    Next I
    Call AMOEBA(P(), Y(), MP, NP, NDIM, FTOL, ITER)
    Print
    Print Tab(5); "Iterations: ", Format$(ITER, "##")
    Print Tab(5); "Vertices of final 3-D simplex and"
    Print Tab(5); "function values at the vertices:"
    Print Tab(5); "I        X(I)       Y(I)       Z(I)      FUNCTION"
    For I = 1 To MP
        Print Tab(5); Format$(I, "#");
        For J = 1 To NP
            Print Tab(J * 11); Format$(P(I, J), ".#####0");
        Next J
        Print Tab(46); Format$(Y(I), ".#####0")
    Next I
    Print
    Print Tab(5); "True minimum is at ( 0.5, 0.6, 0.7)"
End Sub
Function FAMOEB(X())
    AAA = (X(1) - 0.5) ^ 2 + (X(2) - 0.6) ^ 2 + (X(3) - 0.7) ^ 2
    FAMOEB = 0.6 - BESSJ0(AAA)
End Function