options=optimset('TolFun',1e-6,'Display','Final');
xx=fsolve('A1_Kval_PR',[.29 .71 0.665 0.335],options);
x = xx(1:2)
y = xx(3:4)

K = y./x
alpha_12 = K(1)/K(2)