matlab问题 求大家指点

菠萝饭

想要优化一个函数，程序如下：
1 y2 k1 N- N9 W, M3 U%给定初值
cm=3;
cp=5;
l=20;
5 A( I: A' o9 K" g- G" S: B2 Ph=0.05;
p0=1/h;
q0=2*l/h;
- M; A% }) h7 L; ]  ?) }0 N5 w- tnn=0;
cc=cm*l*l
%主程序
) w7 o4 O6 `+ N. a7 r+ F6 ufor p=1:h:p0
    for q=1:h:q0
  X5 a, F% K1 ]  c9 x3 m$ Q        a=p*h;
        b=q*h;
) f" u/ H. L/ f# x        for n=1:12*l/b) %行向量
% j) i  Q3 r: B% D            c(n)=cm*b*b/4*((1-a^n)/(1-a)+2*n/(1-a)-2*(1-a^n)/(1-2*a+a*a))+cm*(a^n*(l-n*b/2)*(l-n*b/2)+(1-a^n)*b*(l-n*b/2)/(1-a))+n*cp;
            if c(2)>=c(1)
                nn=1;
            elseif c(2*l/b)-c(2*l/b-1)<=0
                nn=2*l/b;
1 M3 s' L/ z  Q( j/ R            else
  ]! n6 X" K! ?$ C8 G                for k=1:12*l/b)
                    if (c(k)-c(k-1)<=0 && c(k+1)-c(k)>=0)
+ u, w/ ?) i. S' q8 L# r! i                        nn=k;
                    end
                end
            end
& _1 N: u/ B5 v; T# T+ A$ J        end
        if(c(nn)<cc)
            n=nn;
            aa=a;
            bb=b;
. e2 {& {3 G. r- P! Z           cc=c(nn)
+ r, o3 B+ b& d- X8 j        end
    end
) d. y8 C8 B3 h7 V" ^2 fend
fprinft('%d',n)
  [1 Q4 {( v) X3 z4 [- ]fprintf('%d',cc)
( N" C1 @5 n4 f: T7 m  x2 n+ @fprintf('%d',aa)
fprintf('%d',bb)

; c' i9 L1 I' ?8 l$ e
- r" A5 Q" K* f9 ~运行不出来，错误提示：
1 X' c# w: _/ i- n: Pcc =

  H3 E: L% T* ]$ C. O7 z        1200
# s  D4 n$ R+ j: ], q  z/ K% v
??? Attempted to access c(2); index out of bounds because numel(c)=1.

# a. s" B9 ~2 S# G/ H& I8 _1 G+ w
这是什么问题呀......