//Données géométriques

clear
clf
a=80;b=80;c=70;h=4; //h= pas de la vis%

//calculs intermédiaires
w=atan(a/c)
y=a^2+b^2+c^2
z=2*b*sqrt(c^2+a^2)
L0=sqrt(a^2+(b+c)^2)//lambda_0%;
p=0.0001; //précision du calcul


function r=f(x,teta_32)
   r=360/(2*%pi)*(acos(((L0-h*(teta_32/360))^2-y)/z)-w)-x
endfunction


for i=1:100:8000
        teta_32=i
        teta_51min=-0;teta_51max=100
    while
        
        
        
        
        
        
        
        
    end
   
   
plot(teta_32,teta_51,"*r")
xtitle("loi entrée-sortie Maxpid","teta_32 en °","teta_51 en °")

end






