a2=28
ctheta3=cos(1:1:360)
stheta3=sin(1:1:360)
a3=39
stheta2=sin(1:1:360)
ctheta2=cos(1:1:360)
s=(a2*ctheta2)+(a3*ctheta3)
e=-((a2*stheta2)+(a3*stheta3))
subplot(3,2,1),plot(s,'c'),xlabel('Posición'),ylabel('Tiempo'),title('Grafica de posición s','FontSize',12)
subplot(3,2,2),plot(e,'r'),xlabel('Posición de e'),ylabel('Tiempo'),title('Grafica de posición e','FontSize',12)
DCt3=diff(ctheta3)
DSt3=diff(stheta3)
DCt2=diff(ctheta2)
DSt2=diff(stheta3)
W3=(a2*DCt2)./(a3*ctheta3)
subplot(3,2,3),plot(W3,'g'),xlabel('Velocidad de theta3'),ylabel('Tiempo'),title('Grafica de velocidad angular de theta3','FontSize',12)