function out=q1(x, N_theta, N_phi)

r=abs(x);
alpha=angle(x);


theta1=(alpha+acos(r/2));
theta1=round(theta1/2/pi*N_theta)/N_theta*2*pi;
phi1=(alpha-acos(r/2));
phi1=round(phi1/2/pi*N_phi)/N_phi*2*pi;
out=exp(1i*theta1)+exp(1i*phi1);



 