Buenas, ya pude resolver el problema, estaba poniendo np.arctan2(y,x); ya que pensaba que era lo mismo poner arctan(ys-y/xs-x)=arctan(y/x) ya que xs=ys=0; pero al ser np.arctan2(ys-y,xs-x) = np.arctan2(0-y,0-x) tenia que poner el signo negativo, ahora con eso ya me salio, era ese el error que me daba mal todo, ahora lo cambie y me da lo que m e deberia dar