import matplotlib.pyplot as plt
import math

x_def = [pow(math.pi, 2) / 2]
y_def = [2 * math.pi]
h = -0.0001
i = 0
while 1:
    if x_def[i] < math.e:
        break
    x_def.append(x_def[i] + h)
    y_def.append((pow(math.tan(math.cos(x_def[i] * y_def[i] * y_def[i])), 2) + 1) / (math.sin(x_def[i] / y_def[i])))
    i += 1
print(y_def[len(y_def) - 1])
plt.plot(x_def, y_def)
plt.show()


решение: y(e) = -6.660331226646006