#!/usr/bin/env python

coding: utf-8

import time

if name == "main": l = 0.1 m = 10 a = 1.39e-5 dat = [0.18, 18] p = 100 alf = 1163 tf = 300 lam = 50 tma = 1800 dx = l/m m1 = m+1 fo = a+dat[0]/dx*dx

if fo > 0.5:
    print("fo : ", fo)
else:
    tim = 0
    tn = [0 for i in range(m+2)]
    print("tim : ", 0)
    for i in range(m+1):
        print("T"+str(i)+" : ", tn[i])

    cn = 0
    while tim <= tma:
        cn += 1
        tim += dat[0]
        tn[m+1] = tn[m-1]
        for i in range(m+1):
            tn[i] = (1-2*fo)*tn[i]+fo*(tn[i+1]+tn[i-1])
        tn[0] = (tn[1]+alf*dx*tf/lam)/(1+alf*dx/lam)
        if cn > p:
            print("tim : ", 0)
            for i in range(m+1):
                print("T"+str(i)+" : ", tn[i])
            cn =