[mathjiang]

我只是把你贴出来的程序精简了一下：
% from lanfans.com
function f = shufcn(x)
if size(x,2)<2
error('输入的矩阵至少含2列')
end
i =1:5;
for j =1:size(x,1)
temp1 = i.*cos((i+1)*x(j,1)+i);
temp2 = i.*cos((i+1)*x(j,2)+i);
f(j) =sum(temp1)*sum(temp2);
end