Проверка при трех неизвестных

import numpy as np
from scipy.optimize import fsolve
 
XArray = []#Массивы подходящих значений
YArray = []
ZArray = []
q1 = 0#Л
q2 = 0
q3 = 0
L = 2000
# Функция, возвращающая значения уравнений
def equations():
    q1 = 0
    q2 = 0
    q3 = 0
    for x in range(0,1000,100):
      for y in range(0,1000,100):
        for z in range(0,1000,100):
          eq1 = (x - q1)**2 - y**2 - z**2 - L ** 2
          eq2 = (x - q2*np.cos(120))**2 + (y - q2*np.sin(120))**2 - L**2 - z
          eq3 = (x - q3*np.cos(240))**2 + (y - q3*np.sin(270))**2 - L**2 - z
          if eq1 > -100 and eq1 < 100 and eq2 > -100 and eq2 < 100 and eq3 > -100 and eq3<100:
            print(eq1, eq2, eq3)#XArray.append(x)
            #XArray.append(x)
            #YArray.append(y)
            #ZArray.append(z)
# Решение уравнений
equations()
print("Массивы подходящих значений:")
print(XArray)
print(YArray)
print(ZArray)

Cos1 = math.radians(np.cos(120))
Cos2 = math.radians(np.cos(240))
Sin1 = math.radians(np.sin(120))
Sin2 = math.radians(np.sin(240))
eq1 = (x - q1)**2 - y**2 - z**2 - L ** 2
eq2 = (x - q2*Cos1)**2 + (y - q2*Sin1)**2 - L**2 - z ** 2
eq3 = (x - q3*Cos2)**2 + (y - q3*Sin2)**2 - L**2 - z ** 2

from math import sin,cos,radians
 
print(sin(radians(30))) # вывод 0.49999 = 0.5