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

clc
clear all
close all

R=[12 4 10 7];
teta=[0 0 45*pi/180 135*pi/180];
[teta(3), teta(4)]=m4bar(R,teta);
[teta(3) teta(4)]