d13r1.txt

常用的数值算法的VB程序

Private Sub Command1_Click()
    'PROGRAM D13R1
    'Driver for routine MOMENT
    PI = 3.14159265
    NPTS = 10000
    NBIN = 100
    NDAT = NPTS + NBIN
    Dim DATA(100100)
    I = 1
    For J = 1 To NBIN
        X = PI * J / NBIN
        NLIM = CInt(Sin(X) * PI / 2# * NPTS / NBIN)
        For K = 1 To NLIM
            DATA(I) = X
            I = I + 1
        Next K
    Next J
    Print
    Print Tab(5); "Moments of a sinusoidal distribution"
    Call MOMENT(DATA(), NPTS, AVE, ADEV, SDEV, VAR, SKEW, CURT)
    Print
    Print Tab(30); "Calculated    Expected"
    Print Tab(5); "Mean : "; Tab(32); Format$(AVE, "#.###0");
    Print Tab(44); Format$(PI / 2#, ".####")
    Print Tab(5); "Average Deviation : ";
    Print Tab(32); Format$(ADEV, ".###0"), Tab(44); 0.5708
    Print Tab(5); "Standard Deviation : ";
    Print Tab(32); Format$(SDEV, ".###0"), Tab(44); 0.6837
    Print Tab(5); "Variance : ";
    Print Tab(32); Format$(VAR, ".###0"), Tab(44); 0.4674
    Print Tab(5); "Skewness : ";
    Print Tab(32); Format$(SKEW, "##.###0"), Tab(44); 0#
    Print Tab(5); "Kurtosis : ";
    Print Tab(32); Format$(CURT, "##.###0"), Tab(44); -0.8062
End Sub