function q1=cinv_pendulo(x0,y0)
%cinemática inversa del péndulo
q1=atan(y0/x0);
end