多元非线性规划Matlab,matlab多元非线性规划，程序有点小问题，急求修改

A=[];B=[];

LB=zeros(23,1);UB=ones(23,1);%下限都为0上限为1

Aeq=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1];Beq=[1];%线性等式矩阵A为系数B为右端常数

x0=zeros(23,1);

[x,fval]=fmincon('myfun',x0,A,B,Aeq,Beq,LB,UB)

function Y=myfun1(x)

c1=[0.4404 0 0.1180 0 0 1.0000 1.0000 0 1.0000 0.5476 1.0000 1.0000 0 1.0000 1.0000 1.0000 1.0000 0.7048 1.0000 1.0000 1.0000 1.0000 0];

c2=[0 0.0062 0 0.4370 0.9310 0.1845 0 0.8640 0.4438 0.5476 1.0000 1.0000 1.0000 0 0 0.5476 0.5476 0.7048 0.5476 0 0.5476 0 0.8047];

c3=[0.8450 0.7967 0.7998 0.2257 1.0000 0 0.6592 1.0000 0.4438 0 0.2945 0 1.0000 0 0 0 0 1.0000 0 0 0 0 1.0000];

c4=[1.0000 1.0000 1.000 1.0000 0.6049 0 0.6592 0.8640 0 1.0000 0 0.5476 0 0 0 0.5476 0.5476 0 0.5476 0 0.5476 0 0.4603];

c=[c1;c2;c3;c4];

M1=[-1 -1 -1 0 0 0 -1 0 0 0 0 0 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1];

M2=[0 0 0 -1 -1 -1 0 -1 -1 -1 -1 -1 0 0 0 0 0 0 0 0 0 0 0 ];%定义权重

i=0;

r=0;

while i<4

i=i+1;

j=0;

while j<23

j=j+1;

r11=sqrt(sumsqr(c(i,j)*x(j)+M1(j)*x(j)));

end

k=0;

while k<23

k=k+1;

r12=sqrt(sumsqr(c(i,k)*x(k)+M2(k)*x(k)));

end

r=r11-r12+r;

end

Y=r;