//Données géométriques

clear
clf
a=595;b=460;c=460;h=3;r=20/70;L0=443 //h= pas de la vis;L0=lambda0
p=0.0001; //précision du calcul


function r=f(teta_m,teta_30)
   r=(1/r)*(360/h)*(sqrt(a^2+b^2+c^2+2*a*c*cos(teta_30*2*%pi/360)-2*b*c*sin(teta_30*2*%pi/360))-L0)-teta_m
   endfunction


for i=74:106
        teta_30=i
        teta_m_min=0;teta_m_max=120000
    while
        
        
        
        
        
        
        
    end
   
   
plot(teta_30,teta_m,"*r")
xtitle("loi entrée-pilote electrique","teta_30 en °","teta_m en °")


end






