已知某一400mm普通机床，电机功率为7.5kw，主轴的内径d-=48mm，作用在主轴外伸端处的载荷F=15000N，许用挠度[y]=0.124mm，需用切应力[t]=240MPa，允许转角[sita]=0.01rad。主轴材料的密度rou=7.8x10^-6kg/mm^3，主轴材料的弹性模量为E=2x10^5MPa,剪切模量为G=8x10^5MPa。主轴的计算转速为50r/min，经过计算主轴的最大输出扭矩为T=1432Nm。
 根据简历的主轴优化模型，结合Matlab优化工具箱对主轴结构进行优化设计。
 
下面是代码：
 (1）编写目标函数的m文件：objfun.m,返回x处的函数值f
 function f=objfun(x)
 f=6.126*(10e-006)*(x(2)+x(3))(x(1)^2-48^2);
 (2)因设计约束含非线性约束，需编写一个描述非线性约束的nonlinconstr.m文件：
 function [g,ceq] =nonlinconstr(x)
 g(1)=0.5093*x(2)^2*(x(2)+x(3))/(x(1)^4-49^4)-0.124;
 g(2)=7162.5/(x(1)^3-48^3)-240;
 g(3)=0.5093*x(1)*x(2)/(x(1)^4-48^4)-0.01;
 ceq=[];
 (3)给定变量的初值，并调用优化函数：
 x0=[140,130,550];
 A=[0,0,-1;0,0,1;-1,0,0;1,0,0;0,-1,0;0,1,0]
 b=[-200;650;-50;180;-80;150]
 Ib=[0,0,0]
 options=optimset('display','iter','latgescale'.'off');
 [x,fval,exitflag,ouput]=fmincon('objfun',x0,A,b,[],[],Ib,[],'nonlinconstr',options)
 
但是输入进去有错误，求高手指出错误，十万火急啊！！！！！！！