% Programa de solucion de la funcion
%   Ue, este programa es llamado 
%   por el  programa ssbbp

function f=ssubb(Ue)
global Phi1 Phi2 Yave m1 m2 T R L C E

% Funcion a trabajar para resolverla
%      numericamente
Ut=(1-Ue);
num1=m2*(Phi1^Ut)-m1*(Phi2^Ut);
den=m2*(1-(Phi1^Ut))-m1*(1-(Phi2^Ut));
fact1=(E*Ue*T);
num2=Phi2^Ut-Phi1^Ut;
fact2=E;
f=((num1/den)*fact1+(num2/den)*fact2)/L+(E*Ue*T)/(2*L)-Yave;
end
