计算水力学追赶法vb程序

Option Explicit
Option Base 1
Const Nsr As Integer = 200
Private Ct As Single, P As Single, V As Single, DTT As Single
Private IB As Integer, DXX(Nsr) As Single, CNO(Nsr) As Single, ZZ(Nsr) As Single, QQ(Nsr) As Single, PP(Nsr) As Single
Private VV(Nsr) As Single, SS(Nsr) As Single, TT(Nsr) As Single
Private L1 As Integer, L2 As Integer, IC As Integer, N As Integer

Private Sub Command1_Click()
 Dim de As Long
 Dim delay As Long
 Dim Dx As Single
 Dim CN As Single
 Dim MT As Integer
 Dim DT As Integer
 Dim ZU As Single
 Dim i As Integer
 Dim j As Integer
 Dim k As Integer
 Dim x As Single
 Dim Zd As Single
 Dim xx As Single
 Dim yy As Single
 Open App.Path + "\resultz.txt" For Output As #1
 Open App.Path + "\resultq.txt" For Output As #2
 DT = Text2.Text
 Ct = Text3.Text
 CN = Text4.Text
 N = Text5.Text
 For i = 1 To N

 DXX(i) = 1 * 1000
 CNO(i) = CN
 ZZ(i) = 4 + 0.025 * (20 + 1 - i)
 QQ(i) = 0
 Next i
 MT = 60 / DT
 DT = 60 * DT
 DTT = 2 * DT
 L1 = 1
 L2 = N
 IC = CInt(Text1.Text)
 If IC = 1 Then
  IB = 0
  ZU = 4.5
 ElseIf IC = 2 Then
  IB = 1
  ZU = 500
 Else
  IB = 1
 End If
For i = 1 To 49
 For j = 1 To MT
 x = i + j / MT
 Zd = 4 + 1.5 * Sin(3.14159265 * x / 12)
 If IC = 1 Then
 P = ZU
 V = 0
 ElseIf IC = 2 Then
 P = ZU
 V = 0