d11r12.frm

方程组的矩阵解法

VERSION 5.00
Begin VB.Form Form1 
   Caption         =   "Form1"
   ClientHeight    =   5820
   ClientLeft      =   60
   ClientTop       =   345
   ClientWidth     =   6435
   LinkTopic       =   "Form1"
   ScaleHeight     =   5820
   ScaleWidth      =   6435
   StartUpPosition =   3  'Windows Default
   Begin VB.CommandButton Command1 
      Caption         =   "Command1"
      Height          =   375
      Left            =   4680
      TabIndex        =   0
      Top             =   5160
      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 D11R12
    'Driver for routine PROBKS
    Dim Text(50) As String
    Print
    Print Tab(5); "Probability func. Kolmogorov-Smirnov statistic"
    Print
    Print Tab(5); "Lambda:    Value:     Graph:"
    NPTS = 20
    EPS = 0.1
    SCALE1 = 40#
    For I = 1 To NPTS
        ALAM = I * EPS
        Value = PROBKS(ALAM)
        Text(1) = "*"
        For J = 1 To 50
            If J <= CInt(SCALE1 * Value) Then
                Text(J) = "*"
            Else
                Text(J) = ""
            End If
        Next J
        Print Tab(5); Format$(ALAM, "#.#0");
        Print Tab(13); Format$(Value, ".###.000");
        Print Tab(25);
        For J = 1 To 50
            Print Text(J);
        Next J
    Next I
End Sub
Function PROBKS(ALAM)
    EPS1 = 0.001
    EPS2 = 0.00000001
    A2 = -2# * ALAM ^ 2
    FAC = 2#
    PROBKS = 0#
    TERMBF = 0#
    For J = 1 To 100
        TERM = FAC * Exp(A2 * J ^ 2)
        PROBKS = PROBKS + TERM
        AAA = Abs(TERM)
        If AAA < EPS1 * TERMBF Or AAA < EPS2 * PROBKS Then
            Exit Function
        End If
        FAC = -FAC
        TERMBF = Abs(TERM)
    Next J
    PROBKS = 1#
End Function