MATLAB非线性方程组问题.M文件有了,问题是工作区的代码怎么写啊,最好能输出到EXCEL中.function f=o

MATLAB非线性方程组问题
.M文件有了,问题是工作区的代码怎么写啊,
最好能输出到EXCEL中.
function f=oilyy(x,yy)
a=300;
m=5;
I=800;
e=0.5;
n=293;
h=0.01;
M=0.9;
B=x(1);
P=x(2);
c=x(3);
g=x(4);
z=x(5);
r=x(6);
A=x(7);
C=x(8);
o=x(9);
D=x(10);
E=x(11);
F=x(12);
f(1)=5.67*0.05*pi*0.158*((B/100)^4-(0.0552*(a^1.5)/100)^4)/2-0.05*pi*0.158*(c-B)/(2*h)+(0.54*0.0259*pi*(B-a)*(451100*(B-a))^0.25)/2;
f(2)=I*m*M*0.9*0.02-(5.67*(P/100)^4-F)/1.17-5.67*0.85*0.158*((P/100)^4-(0.0552*(a^1.5)/100)^4)-0.0259*0.54*(P-a)*(451100*(P-a))^0.25;
f(3)=I*m*M*0.9*0.91*0.02*0.1/0.158-0.05*pi*0.158*(c-B)/(2*h)-(0.158/2)*pi*e*(5.67*(c/100)^4-E)/(1-e);
f(4)=I*m*M*0.9*0.02*0.91*(0.58/0.158+0.1*0.9/0.158)+5.67*pi*0.047*((z/100)^4-(g/100)^4)/20.143-(5.67*(g/100)^4-D)/3.043;
f(5)=I*m*M*0.9*0.91*0.93*(0.058/0.158+0.1*0.9/0.158)-5.67*pi*0.047*((z/100)^4-(g/100)^4)/20.143-2*pi*1.1*(z-r)/(log(47/45));
f(6)=2*pi*236*(r-A)/(log(45/43))-2*2*pi*398*(A-C)/(log(8/6));
f(7)=2*2*pi*398*450*(A-C)/(log(8/6))-2000*900*yy*pi*(0.006^2)*(o-n)/4;
f(8)=2000*900*yy*pi*(0.006^2)*(o-n)/4-450*0.6*pi*2*0.023*(12000*yy)^0.8*(6.9)^0.4*(C-(o+n)/2);
f(9)=D-0.85*5.67*(g/100)^4-0.15*(0.612*E+0.388*F);
f(10)=E-e*5.67*(c/100)^4-(1-e)*(0.449*D+0.551*F);
f(11)=F-0.85*5.67*(P/100)^4-0.15*(0.447*D+0.553*E);
f(12)=1.1*(z-r)/(log(47/45))-236*(r-A)/(log(45/43));
end