d12r5.txt

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

Private Sub Command1_Click()
    'PROGRAM D12R5
    'Driver for routine COSFT
    EPS = 0.001
    NP = 16
    NPP2 = NP + 2
    Width1 = 30#
    PI = 3.14159
    Dim DATA(16), SIZE1(16)
    Print
    Print Tab(5); "Period of sinusoid in channels (2-";
    Print Format$(NP, "##"); "; or 0 to stop)"
    PER = 4
    If PER <= 0# Then Exit Sub
    For I = 1 To NP
        DATA(I) = Cos(2# * PI * (I - 1) / PER)
    Next I
    Call COSFT(DATA(), NP, 1)
    BIG = -10000000000#
    SMALL = 10000000000#
    For I = 1 To NP
        If DATA(I) < SMALL Then SMALL = DATA(I)
        If DATA(I) > BIG Then BIG = DATA(I)
    Next I
    SCAL1 = Width1 / (BIG - SMALL)
    For I = 1 To NP
        NLIM = Int(SCAL1 * (DATA(I) - SMALL) + EPS)
        Print Tab(5); Format$(I, "##");
        For J = 1 To NLIM + 1
            Print "*";
        Next J
    Next I
    Call COSFT(DATA(), NP, -1)
    BIG = -10000000000#
    SMALL = 10000000000#
    For I = 1 To NP
        If DATA(I) < SMALL Then SMALL = DATA(I)
        If DATA(I) > BIG Then BIG = DATA(I)
    Next I
    SCAL1 = Width1 / (BIG - SMALL)
    For I = 1 To NP
        NLIM = Int(SCAL1 * (DATA(I) - SMALL) + EPS)
        Print Tab(5); Format$(I, "##");
        For J = 1 To NLIM + 1
            Print "*";
        Next J
    Next I
End Sub