from cmath import inf
import numpy as np
imax = 1572
for i in range(imax + 1):
result = np.tan(i / imax * np.pi / 2)
print(i, result)
def get_i(l):
return 2 / np.pi * np.arctan(l) * imax
print(get_i(0.0))
print(get_i(0.0005))
print(get_i(2.500))
print(get_i(2.5057922722284847))
print(get_i(inf))