i = n+1
h = i
k = 0
z = 2

for j in range(1,n):
    while i <= h + n - k-1:
        if i == h:
            if i == node_number-1:
                    y[i] = 0
                    theta[i] = 0
                    M[i] = Me
                    nu[i] = PM(M[i], gamma)
                    mu[i] = m.asin(1 / M[i]) * RTOD
                    mI= np.tan(m.radians(mu[i - n + k] - theta[i - n + k] + mu[i]) / 2)
            else:
                    y[i] = 0
                    theta[i] = 0
                    nu[i] = theta[i - n + k] + nu[i - n + k] + (1 / (np.sqrt((M[i - n + k]**2) - 1) - mp.cot(m.radians(theta[i - n + k]))))
                    M[i] = InversePrandtlMeyer(nu[i])
                    mu[i] = m.asin(1 / M[i]) *RTOD
                    mI= np.tan(m.radians(mu[i - n + k] - theta[i - n + k] + mu[i]) / 2)
                    print(mI)
        i = i + 1
    k = k + 1
    h = i

2.184370031112393
1.7380673494238976
1.5419781887088249
1.4039585635698797
1.2989038769190393
1.214861123121338
1.1452669347438043
1.0861568995207378
1.0349662574998555
0.98994725077706
0.9498593079129821
0.9137924304933154
0.8810607737523871
0.8511355525702466
0.8236011144240196
0.7981252463827491
0.7744385459289631
0.752319743774369
0.7315850413556333
0.7120802207977384
0.6936747098613351
0.6762570513386782
0.6597313984083366
0.6440147708625135
0.6290348834011446
0.6147284094475984
0.6010395803489821
0.5879190455792305
0.5753229380400084