Dim i, j
Dim dt, dx
Dim T(20), TP(20)

Open "c:\classes\mee410-00\hw473t\hw473t.dat" For Output As #1

'set finite diff parameters
dt = 0.1 's
dx = 0.003 'm

'set properties
ks = 5 'W/mK
ps = 2200 'kg/m3
cs = 745 'J/kgK

'set init temps
For i = 1 To 20
   T(i) = 20
Next i

Picture1.Cls

For i = 1 To 20000
   
   'solve for "next" temp
   TP(1) = T(1) + (dt / (dx ^ 2 * ps * cs)) * (ks * (T(6) - T(1)) + ks * (T(2) - T(1)) + hs * dx * (Ta - T(1)))
   TP(2) = T(2) + (dt / (dx ^ 2 * ps * cs)) * (ks * (T(7) - T(2)) + ks * (T(1) - T(2)) + ks * (T(3) - T(2)) + hs * dx * (Ta - T(2)))
   
   
   'pass "next" temps back for use in next iteration
   For j = 1 To 20
      T(j) = TP(j)
   Next j
   
   'print results
   If i = 200 Or i = 500 Or i = 1000 Or i = 2000 Or i = 3000 Or i = 4000 Or i = 5000 Or i = 10000 Or i = 15000 Or i = 20000 Then
      Picture1.Print i * dt; "  ";
      'Print #1, i * dt; "  ";
      For j = 1 To 20
         Picture1.Print Format(T(j), "0.0"); "  ";
         
         Print #1, Format(T(j), "0.00"); "  ";
         If j = 5 Or j = 10 Or j = 15 Or j = 20 Then
            Print #1,
         End If
      Next j
      
      Picture1.Print
      Print #1,
   End If

Next i

Close #1