Matlab气液相界面,MATLAB,气液相平衡程序，求帮忙改一下。

function obj=lsqfun3(t12,t21,t13,t31,t23,t32)

%t12,t21,t13,t31,t23,t32，是所需要关联的NRTL方程参数；

%a1,a2,a3,是组分1,2,3，液相平衡数据；

%y1,y2,y3,是组分1,2,3，气相平衡数据；

%ps1，ps2，ps3，是组分1,2,3的饱和蒸汽压；

%p==101.35;

%obj为

p=101.35;

y1=[];

y2=[];

y3=[];

a1=[];

a2=[];

a3=[];

ps1=[];

ps2=[];

ps3=[];

k1=(a2.^2*t21*exp(-0.3*t21)^2+a3.^2*t31*exp(-0.3*t31)*exp(-0.3*t31)+a2.*a3*exp(-0.3*t21)*exp(-0.3*t31)*(t21+t31))./((a1+a2.*exp(-0.3*t21)+a3*exp(-0.3*t31)).^2);

k2=(a2.^2*t12*exp(-0.3*t12)+a2.*a3*exp(-0.3*t12)*exp(-0.3*t32)*(t12-t32))./((a1*exp(-0.3*t12)+a2+a3*exp(-0.3*t32)).^2);

k3=(a3.^2*t13*exp(-0.3*t13)+a2.*a3*exp(-0.3*t13)*exp(-0.3*t23)*(t13-t23))./((a