[求助]如何使用矩阵编此matlab程序

mathtype

原程序如下，我想把温度T用向量500到1000来表示，应如何改动原方程<BR>function xFsolve<BR>% xFsolve.m<BR>% 求解热力学平衡计算<BR><BR>clear all<BR>clc<BR><BR>n10=0.2;n20=1-n10;<BR>T=700;             %体系温度，K<BR>K1=exp(7.9305*log(T)-5.6347e-03*T+5.8622e-07*T^2-2.7848e+04/T-19.03390284);<BR>K2=exp(0.248617*log(T)-1.4693e-03*T+2.28e-07*T^2-4.1539e+04/T+4.213066);<BR>K3=exp(6.833053*log(T)-3.7606e-03*T+3.48e-07.*T^2-7.276675e+03/T-33.8969);<BR>x0 = [n10 n20 0 0 0 0]';<BR>x1 =fsolve(@NonlinEqs,x0,[],K1,K2,K3,n10,n20)<BR><BR><BR>% ------------------------------------------------------------------<BR>function f = NonlinEqs(x,K1,K2,K3,n10,n20)<BR>P=1;<BR>xtotal=x(1)+x(2)+x(3)+x(4)+x(5)+x(6);<BR>f(1) = x(1)+x(2)+x(3)+x(6)-1;<BR>f(2) = 2*x(1)+x(4)+x(5)-2*n10;<BR>f(3)=2*x(2)+x(3)+x(5)-2*n20;<BR>f(4)=x(3)^2*x(4)^2-K1*x(1)*x(2)/(P/xtotal)^2;<BR>f(5)=x(3)*x(5)-K2*x(2)*x(4);<BR>f(6)=x(6)*x(4)^2-K3*x(1);<BR><BR>

suffer

这个好像必须要用循环才行