//Données géométriques

clear
clf
a=595;b=460;c=460;
p=0.0001; //précision du calcul


function r=f(lambda,teta_30)
   r=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))-lambda
   endfunction


for i=74:106
        teta_30=i
        lambda_min=400;lambda_max=800
    while
        
        
        
        
        
        
        
    end
   
   
plot(teta_30,lambda,"*r")
xtitle("loi entrée-pilote hydraulique","teta_30 en °","lambda en mm")


end






