>> x=-10:0.32:10;
>> M=randint(1,10,[1 10]);%%產(chǎn)生a,b,c...n,n=10
>> yy=zeros(10,length(x));
for i=1:10
      yy(i,:)=M(i).^x;
   end%%%%得到a^x b^x c^x n^x,yy的第一行對(duì)應(yīng)a^x,第二行對(duì)應(yīng)b^x,第三行對(duì)應(yīng)c^x...
>> y=sum(yy).^(1./x);%%得到最終的y,sum是按列求和就是(a^x+b^x+c^x+.+n^x)
>> plot(x,y)%%x和y是一一對(duì)應(yīng)的關(guān)系,所以找到y(tǒng)最大時(shí)的列號(hào),那么x的取值就是該列號(hào)對(duì)應(yīng)的元素.
>> [m,n]=find(y==max(y))%%畫出的圖像是單峰值,是最大值,不是最小值.
m =
     1
n =
    33%%最大值的列號(hào)
>> x=x(n)
x =
    0.2400%%y最大時(shí)對(duì)應(yīng)的x