import numpy as np
import matplotlib.pyplot as plt
plt.style.use('ff.mplstyle')

alpha = np.arange(0, 180, 15) * 3.14 / 180  # рад
J =  [34.932, 35.605, 38.659, 43.843, 52.709, 62.682, 67.729, 66.948, 55.613, 45.787, 38.473, 36.451] # кг*м^2
Jt = [46.578, 47.585, 50.706, 56.572, 63.896, 72.126, 76.043, 72.126, 63.896, 56.172, 50.706, 47.585] # кг*м^2

fig, ax = plt.subplots(subplot_kw={'projection': 'polar'})

ax.plot(alpha, J, 'o', color='C0', label='эксперимент')
ax.plot(alpha + 3.14, J, 'o', color='C0')
ax.plot(alpha, Jt, 'o', color='C1', label='расчёт')
ax.plot(alpha + 3.14, Jt, 'o', color='C1')

ax.set_rlim(0, 80)
ax.set_rticks([20, 40, 60])  # Меньше делений
ax.set_rlabel_position(-90)  # Угол для подписей к делениям
ax.grid(True)
ax.legend(loc=(0.8,0.94))

ax.set_title(r'Сечение эллипсоида инерции, кг$\cdot$м$^2$')
fig.savefig('example_8_polar.pdf', bbox_inches='tight')