水文地质数值法.docVIP

一维显式有限差分法： 程序代码： Private Sub Command1_Click() Dim l As Integer, T1 As Integer, u As Single, x As Integer, t As Single Dim k As Integer, i As Integer Dim h(10, 10) As Single Dim s As String i = 0 k = 0 l = 10000 T1 = 1000 u = 0.001 x = 1000 For k = 0 To 10 Step 1 h(0, k) = 10 Next k t = InputBox(请输入时间步长, 一维显式有限差分法, 0) For i = 1 To 9 Step 1 For k = 0 To 9 Step 1 h(i, k + 1) = (T1 * (h(i - 1, k) - h(i, k)) - T1 * (h(i, k) - h(i + 1, k))) * t / (u * x * x) + h(i, k) Next k Next i Print t=; Print t For k = 0 To 10 Step 1 For i = 0 To 10 Step 1 s = s Format$(h(i, k)