import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
PC =[]
x = np.arange(-10,11,1)for i in x:#[1,2]:
y =list(np.arange(-np.sqrt(100-i**2),np.sqrt(100-i**2)+1))for it in y:
PC.append([i,it])#print(PC)#print(np.array(PC).shape, PC[-1])
PC_a =[]for i,j in np.array(PC):#ii# print(ii)# i,j = iiprint(i,j,-np.sqrt(100-i**2-j**2), np.sqrt(100-i**2-j**2))if i**2+ j**2>=100:continue
z =list(np.arange(-np.sqrt(100-i**2-j**2), np.sqrt(100-i**2-j**2)+1))for it in z:
PC_a.append([i,j,it])print(PC_a)print(np.array(PC_a).shape)
PC_a = np.array(PC_a