迭代方程
Private Sub Command1_Click()
Dim TB1(76), TB(76), TA1(76), TA(76) As Single
Dim pa, ca, lma, pb, cb, lmb As Single
pa = 2680#: ca = 1340#: lma = 1.93: pb = 0.00152: cb = 950: lmb = 0.00165
Dim tl, ts, L, Ta0, Tb0, Tf As Single
Ta0 = 1500: Tb0 = 20
Dim dx, dt, ma, mb, hb As Single
Dim b, q, p, a, st, r, dr, ri, ab, aa  As Single
Dim fa, fb, f0 As Single
Dim i As Integer
Dim time As Single
Dim f1, f2 As Single
dr = 1: dt = 0.5: r = 75
time = Val(Text1.Text)
f1 = pa * lma * dt / (ca * dt * (dr * dr))
f2 = pb * lmb * dt / (cb * dt * (dr * dr))
st = p * dt / 60
q = 0#
For i = 1 To 60
TB(i) = Tb0
Next i
For i = 61 To 75
TA(i) = Ta0
Next i
Print TA(i); TB(i)
 Print CStr(time) & "分" & "time' tempreture:"
For p = 1 To 10 Step dt
ri = dr * p
   For i = 2 To 59
     TB(1) = 20
     ab = pb * lmb / (cb * dt)
     TB1(i) = (f2 / 2) * ((2 + dr / ri) * TB(i - 1) + (2 / f2 - 4) * TB(i) + (2 - dr / ri) * TB(i + 1))
  Next i
Print Tab(20); Format$(i, TB(i))
  If i = 60 Then
     ab = pb * lmb / (cb * dt)
     TB1(60) = ab * dt * ((TB(59) - TA(61)) / (2 * ri * dr) + (TB(59) - 2 * TB(60) + TA(61)) / (dr * dr)) + TB(60)
   Print TB(60)
     End If
   If i = 61 Then
   TA(61) = 780
     aa = pa * lma / (ca * dt)
     TA1(61) = aa * dt * ((TB(60) - TA(62)) / (2 * ri * dr) + (TB(60) - 2 * TA(61) + TA(61)) / (dr * dr)) + TA(61)
    Print Tab(20); i = 61; Format$(TA(61))
     End If
   For i = 62 To 74
     TA(61) = 1480
     aa = pa * lma / (ca * dt)
     TA1(i) = f1 / 2 * ((2 + dr / ri) * TA(i - 1) + (2 / f1 - 4) * TA(i) + (2 - dr / ri) * TA(i + 1))
  Print TA(i)
    
   Next i
Print Tab(20); Format$(TA(i))
   
   TA1(75) = 4 * (TA(74) + (1 / f1 - 4) * TA(75))
 Print Format$(TA(75), "###.####")
 Next p
 Print TA(i)