options=optimset('TolFun',1e-10,'Display','Iter');
xx0 = [.5 .5 .5 .5 250 .5 .5 .5 .5 230]';

xx=fsolve('A2_P2_f',xx0,options);
x3 = xx(1:2);
y2 = xx(3:4);
T2 = xx(5)-273.15;
x4 = xx(6:7);
y5 = xx(8:9);
T4 = xx(10)-273.15;

Stage1_xyT = [x3 y2 T2*ones(2,1)]
Stage2_xyT = [x4 y5 T4*ones(2,1)]
