% solução numérica sequencial de mecanismo de 4 barras

clc
clear all
close all

R=[12 4 10 7];
teta=[0 0 45*pi/180 135*pi/180];
dt=5*pi/180;

for i=1:72
  [teta(3), teta(4)]=m4bar(R,teta);
  angs(i,:)=[teta(2) teta(3) teta(4)];
  teta(2)=teta(2)+dt;
end

angs=angs*180/pi; % conversão para graus
plot(angs(:,1),angs(:,2),'r-',angs(:,1),angs(:,3)),'b-';
xlabel('Ângulos de entrada (graus)');
ylabel('Ângulos de saída (graus)');
axis([0 360 0 180])
legend('Teta 3','Teta 4')