matlab问题 求大家指点

菠萝饭

想要优化一个函数，程序如下：
%给定初值
cm=3;
cp=5;
+ v9 g/ A* l; l$ S. h/ O/ i: c" }l=20;
( z* c& m+ ~3 ^! s+ ph=0.05;
p0=1/h;
q0=2*l/h;
/ T  K, l8 o& R# W  }, wnn=0;
cc=cm*l*l
%主程序
' o( O$ Z- E# L/ wfor p=1:h:p0
    for q=1:h:q0
        a=p*h;
        b=q*h;
        for n=1:12*l/b) %行向量
  B5 Q& d1 R& \+ ]            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;
/ |- b* m: C& Z            if c(2)>=c(1)
                nn=1;
9 \7 s( ^2 e; v8 I( e            elseif c(2*l/b)-c(2*l/b-1)<=0
: T& U& d3 ]  f5 W( {7 X5 \7 d                nn=2*l/b;
' s! M: s! u& h+ e5 {' a. q            else
3 n* B+ X' U3 V* m                for k=1:12*l/b)
                    if (c(k)-c(k-1)<=0 && c(k+1)-c(k)>=0)
  O/ I( w. K9 j& P, ]* d2 G                        nn=k;
% G6 Y3 n) E& x2 Q% Y                    end
                end
            end
# m, h2 d# _: N' g. t5 Y- _        end
        if(c(nn)<cc)
7 Y4 @( {+ z. m( V/ V            n=nn;
            aa=a;
            bb=b;
. P# G) l4 i. @. r9 L* V& v0 Z1 z" y: R           cc=c(nn)
! V8 M/ z9 j8 F) B. `3 d/ I        end
    end
end
$ a' g5 ^# }6 ]! o  x* S8 ifprinft('%d',n)
" \8 D% m. q  M9 V4 Vfprintf('%d',cc)
fprintf('%d',aa)
, x6 [( c6 m  Efprintf('%d',bb)
1 H9 j. r5 N7 A5 G- P* f# G  j
5 M& x' m$ i+ B( {% M) ^
运行不出来，错误提示：
cc =
+ V8 L- y& L' V9 A
+ O/ p$ m8 }9 B5 b% y: x" }2 j        1200

??? Attempted to access c(2); index out of bounds because numel(c)=1.
8 v* q( m! G5 u+ U7 V

4 [+ w7 C+ t) T这是什么问题呀......