#Métodos Numéricos - 2023. IMERL, Fing
#Código de ejemplo para el ejercicio 5a del práctico 4

#Parte a
l1=1;
l2=1.5;
gamma=3*pi/5;
f=@(a) l2.*cos(pi-gamma-a)./(sin(pi-gamma-a)).^2-l1.*cos(a)./(sin(a).^2);

#Hallo alfa usando el método de bisección
#Alfa debe estar entre 0 y pi/2
#Hay problemas si sin(alfa)=0 y si sin(pi-gamma-alfa)=0.
k=3;
a_=[0.1];
b_=[pi-gamma-0.1];
kmax=20;
a=a_(end);
b=b_(end);
m=(a+b)/2;

X=linspace(0.1,pi-gamma-0.1,100);
Y=f(X);
plot(X,Y)

#Usando método de bisección
while (kmax>k) & (b-a>eps)
    f_a = f(a)
    f_b = f(b)
    f_m=f(m)

    if f_a*f_m>0
      a=m;
      a_=[a_;a];
    else
      b=m;
      b_=[b_;b];
    endif

    m=(a+b)/2
    f_m=f(m);
    k=k+1;
endwhile

alfa=m
L=l2/sin(pi-gamma-alfa)+l1/sin(alfa)