function Taout = findTa(Ha,Ta0,Fa)
Ta1 = Ta0;
tol = 1.e-6;
while(1)
    Ta2 = Ta1 + 0.01;
    Tanew = Ta1 - 0.01*fHa(Ha,Ta1,Fa)/(fHa(Ha,Ta2,Fa)-fHa(Ha,Ta1,Fa));
    if abs((Tanew-Ta1)/Tanew) < tol ; break; end
    Ta1 = Tanew;
end
Taout = Tanew;

