//Données géométriques

clear
clf
a=100;b=150;c=460;l=280;c=20;d=350
p=0.0001; //précision du calcul


function r=f(A_teta_31,B_teta_31,C_teta_31,teta_21)
   r=-360/(2*%pi)*acos(-C_teta_31/(sqrt(A_teta_31^2+B_teta_31^2)))+360/(2*%pi)*atan(B_teta_31/A_teta_31)-teta_21
      endfunction


for i=-145:0
        teta_31=i
        A_teta_31=-2*b*c+2*a*d+2*l*c*sin(teta_31*2*%pi/360)-2*l*d*cos(teta_31*2*%pi/360)
        B_teta_31=2*a*c+2*b*d-2*l*d*sin(teta_31*2*%pi/360)-2*l*c*cos(teta_31*2*%pi/360)
        C_teta_31=a^2+b^2+c^2+d^2-2*b*l*sin(teta_31*2*%pi/360)-2*a*l*cos(teta_31*2*%pi/360)
        teta_21_min=-360;teta_21_max=10000
               
          while
              
              
              
              
              
              
              
          end        
    
 plot(teta_31,teta_21,"*r")
xtitle("loi entrée-sortie portail","teta_31 en °","teta_21 en °")


end






