matlab问题 求大家指点

菠萝饭

想要优化一个函数，程序如下：
/ x4 k$ I" _/ k; O) k%给定初值
8 Z* H0 k7 |, [1 Mcm=3;
cp=5;
5 G$ y* e- ^* v  }8 `l=20;
3 c% i7 G* G3 th=0.05;
p0=1/h;
9 b6 b8 X! h, O8 {q0=2*l/h;
# |; |- _) S2 l6 @9 Xnn=0;
( U4 K1 v) g2 u6 Ycc=cm*l*l
%主程序
& B" \3 a: r2 ~3 U4 A: u5 N2 ifor p=1:h:p0
( H  ^; B# Z+ F    for q=1:h:q0
        a=p*h;
+ J+ O8 f* Y/ x        b=q*h;
2 d3 U5 P6 |; J1 s        for n=1:12*l/b) %行向量
            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)
6 ?8 S6 U  b7 A( F; m2 s3 r$ T                nn=1;
            elseif c(2*l/b)-c(2*l/b-1)<=0
  ^; j4 W% P. {* ?- P* s                nn=2*l/b;
            else
                for k=1:12*l/b)
5 X. x6 y2 `6 h6 _                    if (c(k)-c(k-1)<=0 && c(k+1)-c(k)>=0)
                        nn=k;
5 |1 ?+ w0 h& N1 ^                    end
! L7 i6 ?+ v5 _2 E& t                end
            end
. d3 N7 B* ~+ ]        end
& y6 ]- }' E* R6 ]) N/ ~2 e  r) I        if(c(nn)<cc)
            n=nn;
- _# u$ O) D8 `: K& e            aa=a;
$ H2 B/ O( j5 X, g! d' X            bb=b;
           cc=c(nn)
3 F& f. Q7 H, E( j* _        end
/ f" u9 w& f% w    end
% a; T  B- r  P; [" l$ Vend
fprinft('%d',n)
- G, h4 p/ [2 W/ _fprintf('%d',cc)
8 f1 C" x. R8 h2 `0 ?fprintf('%d',aa)
" L5 {! _) U4 \fprintf('%d',bb)
+ [3 T1 q: |- V
) {( D. F7 s) i) H/ l: v
运行不出来，错误提示：
4 W; H0 g, B* X& D3 d* M, Vcc =
4 F" t1 L& L2 \* `8 \
( s( x5 Q0 W5 [( _! K        1200
5 K2 M( z3 g7 O& S5 }
) C  @7 m$ f+ O2 b/ y% }1 Q??? Attempted to access c(2); index out of bounds because numel(c)=1.

# e9 D& N0 R+ S) G& I
这是什么问题呀......