shell.txt

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

Sub SHELL(N, ARR())
    ALN2I = 1.442695: TINY = 0.00001
    LOGNB2 = Int(Log(CSng(N)) * ALN2I + TINY)
    M = N
    For NN = 1 To LOGNB2
        M = Int(M / 2)
        K = N - M
        For J = 1 To K
            I = J
            Do
              DONE% = -1
              L = I + M
              If ARR(L) < ARR(I) Then
                T = ARR(I)
                ARR(I) = ARR(L)
                ARR(L) = T
                I = I - M
                If I >= 1 Then DONE% = 0
              End If
            Loop While Not DONE%
        Next J
    Next NN
End Sub