//Données géométriques

clear
clf
H=109;R=81. //h= pas de la vis%
p=0.0001; //précision du calcul


function r=f(teta_31,teta_21)
   r=360/(2*%pi)*atan((H+R*sin(teta_21*2*%pi/360))/(R*cos(teta_21*2*%pi/360)))-teta_31
   endfunction


for i=-25:210
            teta_21=i
        teta_31_min=-360;teta_31_max=360
    while
        
        
        
        
        
        
        
    end
   
   
plot(teta_21,teta_31,"*r")
xtitle("loi entrée-sympact","teta_21 en °","teta_31 en °")


end






