clear all
clc

E=input('de o valor do deslocamento do cursor em relaçao ao eixo de manivelas E ');
r=input('de o valor da relaçao L/R ');
S=input('de o valor do curso do cursor S ');

ef=10^-8;                 % tolerancia para o metodo de Newton Raphson
R=input('de o valor inicial do raio da manivela ');
f=1;                      %funçao f cuja norma e maior que ef

% coordenadas secundarias

while norm(f)>ef
    f=(((r+1)*R)^2-E^2)^(0.5)-(((r-1)*R)^2-E^2)^(0.5)-S;
    J=(r+1)^2*R*(((r+1)*R)^2-E^2)^(-0.5)-(r-1)^2*R*(((r-1)*R)^2-E^2)^(-0.5);
    R=R-inv(J)*f;
end
R
L=r*R;

fi1=asin(E/(L-R));
fi2=asin(E/(L+R));

alfa=pi+fi1-fi2;
alfa_graus=alfa*180/pi
beta=pi-fi1+fi2;
beta_graus=beta*180/pi

disp('razao temporal alfa/beta')
rt=alfa/beta
