fminbnd的优化问题

羽逐云

function f=H(n,a1)  
, ~/ [6 o/ n, m5 Y0 C8 Pc=23.43/180*pi;
/ [/ C, t: J" {- C: _a0=32.05*pi/180;
& w2 l7 Z. A! R5 @0 k6 g4 Xh=zeros(n,1);
% L9 k- S, {6 N. a6 B& M%a1=zeros(1,n);
    t=0;
    for k = 1: n
, s! K3 `$ k& L        b=asin(sin(c)*cos(pi*k/n));
5 o: k0 R' s% n0 N' i& U& l        A=sqrt(1-(sec(a0))^2*(sin(b))^2);
        B=sin(b)*acos(tan(a0)*tan(b));
        e=atan(B/A);
9 D8 G" h7 v! a( E& ^/ m        h(k)=sqrt(A^2+B^2)*cos(a0+e-a1);
* L' K# k  v+ \0 m3 v- Q" ~        %T(k)=k/n;
        %a1(k)=a1(k)*180/pi;
        t=t+h(k);
    end
    f=189.8325-t;
& ]- @0 b" l! V; `
) r0 _$ U, Z9 r! o4 s& f* C然后用下面程序调用，怎么老出错呢，希望高手能解答一下
5 i6 [) c0 v, v5 q7 ~, I1 v
n=182;
. u9 m* A7 [3 d: G) I* v9 Ta1=fminbnd(@H,0,pi/2,[],n)
不知道哪里出错了，谢谢帮助了

吖信

报错？还是结果错了？我刚试试没报错啊~