Исправление ошибок в программе

import matplotlib.pyplot as plt
import string 
 
def button1_click():
    x = [1,2,3,4,5]
    y = [10,20,15,25,30]
    
    plt.figure(1)
    plt.plot(x, y)
    plt.title('Chart 1')
    
    plt.figure(2)
    plt.plot(x, y)
    plt.title('Chart 2')
    
    memo1 = "Memo 1 content"
    memo2 = "Memo 2 content"
    memo3 = "Memo 3 content"
    
    print(memo1)
    print(memo2)
    print(memo3)
 
plt.show()
# инициализация переменных
dT = [[[0.0 for _ in range(30)] for _ in range(16)] for _ in range(16)]
T = [[[0.0 for _ in range(30)] for _ in range(16)] for _ in range(16)]
time = 0.0
tmin = 0.0
dx = dy = dz = a1 = a2 = lam1 = lam2 = Lx = Ly = Lz = 0.0
u1 = lapl = lapl1 = sum = raz = raz1 = tz = [[0.0 for _ in range(16)] for _ in range(16)]
u = at1 = at2 = at3 = ptau = Tp = Kf = fx = fy = 0.0
i = j = x = y = z = tau = tt = 0
k1 = k3 = t1 = t3 = ta1 = ta3 = e1 = e12 = e2 = e4 = n1 = n2 = n4 = f1 = f3 = w = G1 = G3 = df = dm = dtau = wc1 = wc3 = mg1 = mg3 = dd = dw = wcm1 = wcm3 = k = ee4 = dww = G = fr = Mor = Mo = fp = dfaz = am1 = am2 = Pee4 = Pe4 = Pee2 = pe2 = ee2 = wvv = mra = a3 = ind2 = ind4 = indsr = indd = wk2 = wk4 = Snag = temp = 0.0
 
# данные
dtau = 0.01
a2 = 0.0000514  # воздух
lam2 = 0.0393  # воздух
Lx = 1
Ly = 0.8
Lz = 0.32
dx = Lx / 16
dy = Ly / 16
dz = Lz / 29
temp = 10
 
# Регулятор
E1 = 1.640613
n1 = 34.562029
E2 = 4.790245
n2 = 1000000000000.
E4 = 0.000367
n4 = 63.555858
raz = [[0 for _ in range(16)] for _ in range(16)]
V = [[0 for _ in range(16)] for _ in range(16)]
Sum = [[0 for _ in range(16)] for _ in range(16)]
dT = [[[0 for _ in range(30)] for _ in range(16)] for _ in range(16)]
T = [[[0 for _ in range(30)] for _ in range(16)] for _ in range(16)]
u = 1
at3 = 0.0
at2 = 0.0
time = 0
tau = 3000
 
for tt in range(1, tau+1):
    for i in range(1, 61):
 
        # Тепловые процессы, протекающие в печи
 
        # воздух
        for x in range(1, 16):
            for y in range(1, 16):
                for z in range(29, 1, -1):
                    dT[[x][y][z]] = a2 * dtau * ((T[[x-1][y][z]] - 2*T[[x][y][z]] + T[[x+1][y][z]]) / (dx*dx) +
                                                (T[[x][y-1][z]] - 2*T[[x][y][z]] + T[[x][y+1][z]]) / (dy*dy) +
                                                (T[[x][y][z-1]] - 2*T[[x][y][z]] + T[[x][y][z+1]]) / (dz*dz))
 
        for x in range(1, 16):
            for y in range(1, 16):
                for z in range(29, 1, -1):
                    T[x][y][z] = T[x][y][z] + dT[x][y][z]
for x in range(2, 16):
    for y in range(2, 16):
        raz1[x][y] = temp - T[x][y][17]
        
        for x in range(2, 16):
            for y in range(2, 16):
                lapl1[x][y] = (raz1[x-1][y] - 2*raz1[x][y] + raz1[x+1][y]) / (dx*dx) + \
                             (raz1[x][y-1] - 2*raz1[x][y] + raz1[x][y+1]) / (dy*dy)
        
        for x in range(2, 16):
            for y in range(2, 16):
                raz[x][y] = E1 * (((n1-1)/n1) * raz1[x][y] - (1/n1) * lapl1[x][y])
        
        for x in range(2, 16):
            for y in range(2, 16):
                lapl[x][y] = (raz[x-1][y] - 2*raz[x][y] + raz[x+1][y]) / (dx*dx) + \
                             (raz[x][y-1] - 2*raz[x][y] + raz[x][y+1]) / (dy*dy)
        
        for x in range(2, 16):
            for y in range(2, 16):
                Sum[x][y] = Sum[x][y] + E4 * (((n4-1)/n4) * raz[x][y] - (1/n4) * lapl[x][y]) * dtau
                u1[x][y] = (raz[x][y] + E2 * ((((n2-1)/n2) * raz[x][y] - (1/n2) * lapl[x][y] - V[x][y]) / dTau) + Sum[x][y])             
                T[x][y][30] = u1[x][y] * dz / lam2 + T[x][y][29]
                V[x][y] = raz[x][y]             
# граничные условия
for x in range(2, 16):
    for y in range(2, 16):
        for z in range(29, 1, -1):
            T[0][y][z] = 0  # сторона S2
            T[15][y][z] = 0  # сторона S3
            T[x][0][z] = 0  # сторона S4
            T[x][15][z] = 0  # сторона S5
            T[x][y][0] = T[x][y][1]
while time < tau:
    time += dtau
    
    am2 = T[8][8][17]
    am1 = u1[8][8]  # T[3][8][17]
 
    wvv = time / 60
    if wvv % 1 == 0:
        Series1.append((wvv, am1))
        Series2.append((wvv, am2))
        memo1.append(str(wvv))
        memo2.append(str(am1))
        memo3.append(str(am2))