inlagrn2mainmodule.bas

几种常用的内插数值计算方法

Attribute VB_Name = "MainModule"
Sub Main()
    Dim m As Integer, n As Integer, i As Double, j As Double
    Dim u As Double, v As Double, w As Double
    Dim str As String
    
    ' x方向结点数
    n = 11
    
    ' y方向结点数
    m = 11
    
    ' 分配内存
    ReDim x(n) As Double, y(m) As Double, z(n, m) As Double
        
    ' 结点的值与结点的函数值
    For i = 1 To n
        x(i) = 0.1 * (i - 1)
    Next i
    
    For j = 1 To m
        y(j) = x(j)
    Next j
    
    For i = 1 To n
        For j = 1 To m
            z(i, j) = Exp(-(x(i) - y(j)))
        Next j
    Next i
    
    ' 待插值点的值
    u = 0.35
    v = 0.65
    ' 插值
    w = INLagrn2(n, m, x, y, z, u, v)
    str = str & "z(" & u & ", " & v & ") = " & w & Chr(13)
    
    ' 待插值点的值
    u = 0.45
    v = 0.55
    ' 插值
    w = INLagrn2(n, m, x, y, z, u, v)
    str = str & "z(" & u & ", " & v & ") = " & w
        
    MsgBox str
    
End Sub