d4r29.frm

来自「常用数值算法集哈」· FRM 代码 · 共 74 行

VERSION 5.00
Begin VB.Form Form1 
   Caption         =   "Form1"
   ClientHeight    =   2700
   ClientLeft      =   60
   ClientTop       =   345
   ClientWidth     =   3915
   LinkTopic       =   "Form1"
   ScaleHeight     =   2700
   ScaleWidth      =   3915
   StartUpPosition =   3  'Windows Default
   Begin VB.CommandButton Command1 
      Caption         =   "Command1"
      Height          =   375
      Left            =   2160
      TabIndex        =   0
      Top             =   2040
      Width           =   1335
   End
End
Attribute VB_Name = "Form1"
Attribute VB_GlobalNameSpace = False
Attribute VB_Creatable = False
Attribute VB_PredeclaredId = True
Attribute VB_Exposed = False
Private Sub Command1_Click()
    'PROGRAM D4R29
    'Driver for routine EX
    Dim ACTUAL(5), X(5)
    ACTUAL(1) = 0.5597736
    ACTUAL(2) = 0.2193839
    ACTUAL(3) = 0.1000196
    ACTUAL(4) = 0.01304838
    ACTUAL(5) = 0.001148299
    X(1) = 0.5
    X(2) = 1#
    X(3) = 1.5
    X(4) = 3#
    X(5) = 5#
    Print
    Print "    X       Actual         EX(X)"
    For I = 1 To 5
        Print Tab(4); Format$(X(I), "#0.0#");
        Print Tab(10); Format$(ACTUAL(I), ".0000000E+00");
        Print Tab(25); Format$(EX(X(I)), ".0000000E+00")
    Next I
End Sub
Function EX(X)
    R1 = 0.00107857
    R2 = -0.00976004
    R3 = 0.05519968
    R4 = -0.24991055
    R5 = 0.99999193
    R6 = -0.57721566
    T1 = 8.5733287401
    T2 = 18.059016973
    T3 = 8.6347608925
    T4 = 0.2677737343
    S1 = 9.5733223454
    S2 = 25.6329561486
    S3 = 21.0996530827
    S4 = 3.9584969228
    If X < 1# Then
        AAA = (R1 * X + R2) * X + R3
        EX = ((AAA * X + R4) * X + R5) * X + R6 - Log(X)
    Else
        Y# = (((X + T1) * X + T2) * X + T3) * X + T4
        W# = (((X + S1) * X + S2) * X + S3) * X + S4
        EX = (Y# / W#) / (Exp(X) * X)
    End If
End Function