vector manipulation, for structure and variable storage

import math

deg2rad = math.pi/180
rad2deg = 1/deg2rad

#Cutting conditions
Adc = 1                    #axial depth of cut [mm]
fpt = 1/1000               #feed per tooth [mm/tooth]
phi_st = 0                 #starting angle [rad]
phi_ex = 180*deg2rad       #exit angle [rad]

#Tool geometry
Nt = 1                     #Number of teeth (flute)
d = 0.5                    #diameter [mm]
re = 2/1000                #edge radius [mm]
Theta_hx = 30*deg2rad       #Helix angle [rad]

####################################################################
K = 5                      #number of angular integration steps
L = 5                      #Number of axial integration steps
del_phi = 2*math.pi/K      #integration angle
phi_p = 2*math.pi/Nt       #cutter pitch angle
dz = Adc/L               #integration height

for i in range(1, K+1, 1):
    phi = phi_st + i * del_phi

    for k in range(1, Nt+1, 1):
        phi1 = phi + (k - 1) * phi_p
        #print(phi1)

        for j in range(1, L+1, 1):
            #a_h = axial height of disk element
            a_h = j * dz
            #print(a_h)
            # Immersion angle due to helix
            phi2 = phi1 - ((2*math.tan(Theta_hx)/d)*a_h)
            #print(phi2)

The problem is that my algorithm is calculating the phi2 values in the wrong way, it is not getting the values of a_h correctly that were calculated in that loop ... I did the bills on paper and realized that it is only using the a_h = 1.0 in the phi2 formula ... where I want a_h = [0.2, 0.4, 0.6, 0.8, 1.0] to match the values of phi1 = [1.2566, 2.5132, 3.7698, 5.0264, 6.2830] so that I have phi2 for each position of the vectors ... sorry but I'm a little lay in programming, can anyone analyze the code ?? Thank you in advance !!

edit: phi2 = [0.79472, 1.5894, 2.3842, 3.1789, 3.9736] and (2 * math.tan (Theta_hx) / d) = 2.3094