matlab问题 求大家指点

菠萝饭

想要优化一个函数，程序如下：
) _4 g. H0 Q+ A5 v' D%给定初值
cm=3;
; d( W+ {1 k4 Tcp=5;
l=20;
: D: T. r$ Y" O+ B3 bh=0.05;
/ x& W& H! |' p6 {, g8 ap0=1/h;
, z* S8 T% H- ^! ~5 \q0=2*l/h;
nn=0;
cc=cm*l*l
  F+ H; j1 O2 X/ h9 Y) Z6 X%主程序
4 p* l: B* H/ O) B! M* K8 zfor p=1:h:p0
/ d6 L2 A! U/ y" X& e    for q=1:h:q0
0 h/ x. W1 n) n2 t# E% m# _) C        a=p*h;
8 ^4 i3 a  l1 N4 e  g& |        b=q*h;
        for n=1:12*l/b) %行向量
' m/ J+ R- m8 @0 j            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;
; y# `* u/ v5 V9 S! {; T3 q3 e4 a- w            elseif c(2*l/b)-c(2*l/b-1)<=0
                nn=2*l/b;
+ k7 b0 o* A" t5 S  ^+ B4 o' x            else
                for k=1:12*l/b)
' Y# o% K, s1 t. ^) a                    if (c(k)-c(k-1)<=0 && c(k+1)-c(k)>=0)
2 V( y" D: u: m' H/ T& J* b                        nn=k;
                    end
                end
            end
" h7 |0 Z  h5 |5 G  |: S5 f        end
        if(c(nn)<cc)
7 u- O# F3 y/ e2 h; V. U) M" J7 R2 c            n=nn;
            aa=a;
: p' A3 K; t( T6 k2 n' e& M- Z( |            bb=b;
5 H& p, i' d5 L0 ~7 h4 m5 \& A. P4 _% v           cc=c(nn)
: i2 v9 n: E! @# @6 H' P        end
    end
" E7 E/ l' l5 e. a7 r8 Jend
2 ~1 U/ I: C' R6 [fprinft('%d',n)
2 F* x, k6 _3 M9 V- q1 Z) D& Afprintf('%d',cc)
fprintf('%d',aa)
* u7 l& ]" `9 U+ A" `6 }# Afprintf('%d',bb)

( T5 o# p" {$ z* _
运行不出来，错误提示：
cc =

: k' h( o; @% m) @% A) x2 A        1200
1 I1 T0 E' f2 h* l* e
5 q1 g- Q$ W$ D9 a7 K??? Attempted to access c(2); index out of bounds because numel(c)=1.


这是什么问题呀......