请各位MATLAB高手，不吝赐教，想用matlab编写logit模型的效用函数，然后用极大似然估计求解效用函数的系数值。其中，mydata都是输入数据，X为想求解的各个变量的系数，我实在是头疼，不懂啊！
mydata=xlsread('D:\wo1.xls');
a=mydata(:,6);
b=mydata(:,7);
c=mydata(:,8);
d=mydata(:,9);
e=mydata(:,10);
f=mydata(:,11);
g=mydata(:,12);
h=mydata(:,13);

X0=[1,1,1,1,1,1,1,1,1,1,1]
options=optimset('largescale','off','display','iter','tolx',1e-5,'tolfun',1e-5);
[X,fval,exitflag,output,grad,hessian]=fminunc('fun1',X0,options);

function fun1=myfun(x)
f= x(1)+x(2)*a+x(3)*b+x(4)*c+x(5)*d+x(6)*e+x(7)*(ln(x(8))+exp(x(9)*f)+exp(x(10))*g+exp(x(11))*h);
fun1=-sum(log(f));