%optimF2F3
TF2b=ToutF2tot(end);
TnewF2=TF2b;
uF3b=uF3;

uF3a=uF3b-0.1;
uF3=uF3a;
sim('al');
TF2a=ToutF2tot(end);

uF3c=uF3b+0.1;
uF3=uF3c;
sim('al');
TF2c=ToutF2tot(end);

linTF2=[TF2a TF2b TF2c];
linuF2=[uF3a uF3b uF3c];

while linTF2(2)>max(linTF2(1),linTF2(3))+0.01
[TminF2,indexF2]=min(linTF2);
uF3=(linuF2(indexF2)+linuF2(2))/2;
sim('al');
if indexF2==3
    if    ToutF2tot(end)>linTF2(2)
          linTF2(1)=linTF2(2);
          linuF2(1)=linuF2(2);
          linTF2(2)=ToutF2tot(end);
          linuF2(2)=uF3;
    else  linTF2(3)=ToutF2tot(end);
          linuF2(3)=uF3;
    end
else
    if    ToutF2tot(end)>linTF2(2)
          linTF2(3)=linTF2(2);
          linuF2(3)=linuF2(2);
          linTF2(2)=ToutF2tot(end);
          linuF2(2)=uF3;
    else  linTF2(1)=ToutF2tot(end);
          linuF2(1)=uF3;
    end    
end
end
%update
ToldF2=TnewF2;
uF3new=linuF2(2);
TnewF2=linTF2(2);
uF3=uF3new;