matlab问题 求大家指点

菠萝饭

想要优化一个函数，程序如下：
* Z  I: m1 Y  t5 i) e. O9 @%给定初值
cm=3;
cp=5;
l=20;
h=0.05;
p0=1/h;
q0=2*l/h;
nn=0;
( ^+ w$ X% A+ v- i/ mcc=cm*l*l
' n) v2 T; T; K! F; L- H. H) t" T$ \%主程序
. Z: I7 }4 ~1 t7 Jfor p=1:h:p0
+ o8 l4 R% o/ [- u8 [    for q=1:h:q0
        a=p*h;
        b=q*h;
. w" I  b4 i+ O3 _2 ~+ n' z4 _        for n=1:12*l/b) %行向量
7 a: A+ j' D9 E! D- n" T- }            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;
            else
/ T) l& U) a' x                for k=1:12*l/b)
                    if (c(k)-c(k-1)<=0 && c(k+1)-c(k)>=0)
                        nn=k;
/ ]" @# k3 D' O  C: }0 J! Y7 ~2 o5 c                    end
$ L" i0 o* [3 e& B- l                end
: i6 E% x5 S8 G: ]& W5 O            end
7 f# J: M$ V3 d" v* L        end
8 t8 m3 u2 f; ^5 R: Q        if(c(nn)<cc)
% j/ i+ w' X" f$ d5 C            n=nn;
3 N( _$ r1 I/ C$ e            aa=a;
            bb=b;
; N3 D; T0 c' m4 A* B6 ^           cc=c(nn)
        end
    end
' o" x, ?2 Y5 }/ O& wend
, C, ~" v' q' b$ \, @fprinft('%d',n)
9 j8 a8 D- v8 Y; G! c" v/ k+ Bfprintf('%d',cc)
fprintf('%d',aa)
fprintf('%d',bb)

+ {) N1 h, |, U
. W6 F) y0 W+ D+ z/ d运行不出来，错误提示：
cc =

  T4 l4 F& }4 {! o$ |        1200
$ l1 t! s7 i. }' w9 A
??? Attempted to access c(2); index out of bounds because numel(c)=1.


1 b, {6 }) s) ~+ {% p& ~这是什么问题呀......