function s=mec4bard(C,q,s)
% posição de mecanismo de 4 barras deslizante via Newton-Raphson
% notação S. Doughty

st=s;
eps=1e-6; % condição de convergência
f=pos4bard(C,q,st);
while norm(f) > eps
  J=jaco4bard(C,st);
  dp=J\(-f);
  st(1)=st(1)+dp(1);
  st(2)=st(2)+dp(2);
  f=pos4bard(C,q,st);
end;
s=st;
end


