function frq=cinematica_pendulo(l,q1,beta1)
     frq=[l*cos(q1); l*sin(q1); beta1];
end