Сципи решает уравнения trigonemtri c с ограничениями и множественными корнями

Question

Я пытаюсь добавить ограничения для альфа между 0 и 90 градусами при решении уравнений, как в коде ниже. Буду признателен за любую помощь.

Кроме того, есть ли способ сохранить несколько решений?

from scipy.optimize import fsolve
from scipy.interpolate import spline
import math
import numpy as np

%pylab
from matplotlib import pyplot as plt
import matplotlib.dates as mdates
%matplotlib inline

x = 151
y = 988

def equations(p):
    alphaAngle, betaAngle = p
    f1 = 800*math.cos(math.radians(alphaAngle)) + 600*math.cos(math.radians(betaAngle)) - x
    f2 = 800*math.sin(math.radians(alphaAngle)) + 600*math.sin(math.radians(betaAngle)) - y
    return (f1,f2)

alphaAngle,betaAngle = fsolve(equations, (0,0))
print((alphaAngle,betaAngle))

Результаты:

Использование бэкэнда matplotlib: Qt5Agg Заполнение интерактивного пространства имен от numpy и matplotlib (44.4179365293905, 134.4809022164253)