看楼主的意思，每次求得的解应该是相同的大小的一维数组，因为不知道具体的解的维度，大概写一下，应该像下面这样吧：for i=1:20 x(:,i)=fsolve(@energe,x0,[],Qxi(i),Qyi(i),cx(i),cy(i));end注意要把i＝i+1去掉，for循环会自动加的。。。如果不对的话请把具体维数告诉我，才能写出程序～～
