matlab问题 求大家指点

菠萝饭

想要优化一个函数，程序如下：
) {3 t) }. d$ W* |& ~( n%给定初值
cm=3;
cp=5;
% U9 ?( G. @" Il=20;
h=0.05;
p0=1/h;
q0=2*l/h;
nn=0;
1 N0 l6 P, w: G) _' T) [cc=cm*l*l
%主程序
; q# m5 ~. c  @, B, n1 `for p=1:h:p0
    for q=1:h:q0
* J9 w/ w5 ~  @+ ^        a=p*h;
' _  o( x& u+ X. T+ }* M9 a6 [" g        b=q*h;
) P% J( U' h- _/ c. H/ w4 _        for n=1:12*l/b) %行向量
1 @( W6 Q+ ^, a7 O( z5 K1 Z$ P            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;
& q4 p# ~2 b3 O/ M            if c(2)>=c(1)
; c+ m6 H7 {. J9 }% O                nn=1;
1 p' K) ^" y& }1 O& X: h6 `& }            elseif c(2*l/b)-c(2*l/b-1)<=0
  ?3 h- b7 ?4 B' j  Y                nn=2*l/b;
            else
  ]  W' y% K" k9 ~                for k=1:12*l/b)
                    if (c(k)-c(k-1)<=0 && c(k+1)-c(k)>=0)
  d# b! ?% R( f  R                        nn=k;
                    end
                end
            end
        end
        if(c(nn)<cc)
& H+ R- _! g" D2 E, P            n=nn;
" C; `3 v- |* O8 G7 e6 r            aa=a;
/ }$ H. m7 m& G: V+ |$ q! g4 G            bb=b;
           cc=c(nn)
        end
    end
, v- s0 d8 O# ?* E6 o" p) ]end
( O3 l5 U3 ~3 R# p7 ffprinft('%d',n)
& |* w& |* B  b0 b+ pfprintf('%d',cc)
fprintf('%d',aa)
, S6 m9 B& g4 A% O( A! gfprintf('%d',bb)
& \* q+ l( ?% B: C, _
5 q2 r$ ~7 v" O! V
运行不出来，错误提示：
4 A0 r) o% ?. k* i6 [cc =

  d$ L5 R2 z7 _, _0 m. |- f, [0 n4 }        1200

??? Attempted to access c(2); index out of bounds because numel(c)=1.
" @$ D+ e3 z5 M! ]% u$ n$ p( s8 ~

" i) v" X( I9 c2 N% p) f8 q这是什么问题呀......