//Données géométriques

clear
clf
a=88;b=125; 
p=0.0001; //précision du calcul


function r=f(teta_20,teta_10)
   r=360/(2*%pi)*atan(a*sin(teta_10*2*%pi/360)/(-b+a*cos(teta_10*2*%pi/360)))-teta_20
   endfunction


for i=-45:45
            teta_10=i
        teta_20_min=-50;teta_20_max=50
    while
        
        
        
        
        
        
    end
   
   
plot(teta_10,teta_20,"*r")
xtitle("loi entrée-sortie capsuleuse","teta_10 en °","teta_20 en °")


end






