d4r29.frm

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

VERSION 5.00
Begin VB.Form Form1 
   Caption         =   "Form1"
   ClientHeight    =   3480
   ClientLeft      =   60
   ClientTop       =   345
   ClientWidth     =   4290
   LinkTopic       =   "Form1"
   ScaleHeight     =   3480
   ScaleWidth      =   4290
   StartUpPosition =   3  'Windows Default
   Begin VB.CommandButton Command1 
      Caption         =   "Command1"
      Height          =   375
      Left            =   2520
      TabIndex        =   0
      Top             =   2760
      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
    Open "D:\VB常用数值算法集\DATA\FNCVAL.DAT" For Input As #1
    Do
      Line Input #1, Text$
    Loop While Text$ <> "Exponential Integral Ex"
    Line Input #1, NVAL$
    NVAL = Val(NVAL$)
    Print Text$
    Print
    Print "    X       Actual         EX(X)"
    For I = 1 To NVAL
        Line Input #1, DUM$
        X = Val(Mid(DUM$, 1, 7))
        ACTUAL = Val(Mid(DUM$, 8))
        Print Tab(4); Format$(X, "#0.0#");
        Print Tab(10); Format$(ACTUAL, ".0000000E+00");
        Print Tab(25); Format$(EX(X), ".0000000E+00")
    Next I
    Close #1
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