错误使用e04wc，输出参数太多

Uqzalq

如下为NAG library里面e04wc的程序，是别人文章里给的，用来求某函数最小值的，但是我运行时会说错误使用e04wc，输出参数太多，球球大佬怎么搞？

- ^! f3 o8 Q4 |4 ]3 J# H
function e04wc_example
  {+ w0 H: V7 `
5 O+ c% Q; }' l' f' Y6 Cfprintf('e04wc example results\n\n');
$ |0 J& u! t& ~- i. Oa  = [1, 1, 1, 1];
0 j/ z( w4 O) F9 {2 q6 Hbl = [1, 1, 1, 1, -1e25, -1e25,   25];
bu = [5, 5, 5, 5,    20,    40, 1e25];
istate = zeros(7, 1, 'int64');
% Q: }2 ^! G' ?: J% t" jccon   = zeros(2,1);
cjac   = zeros(2,4);
clamda = zeros(7,1);
hess   = zeros(4,4);
, U7 P& q1 K1 ?x      = [1; 5; 5; 1];
5 @; e7 [( D3 o  h0 C9 {[iw,rw,ifail]=e04wc;
[majits, istate, ccon, cjac, clamda, objf, grad, hess, x, ...
                                          iw, rw, user, ifail] = ...
8 J/ j% G! e1 f; B$ x  e04wd(...
5 {( ^2 S( y  ~* z        a, bl, bu, @confun, @objfun, istate, ccon, cjac, clamda, ...
        hess, x, iw, rw);
fprintf('Final objective value    :   %8.1f\n',objf);
* A2 K4 a) r4 ]" dfprintf('Optimal values for x     :  ');
fprintf(' %9.2f',x);
' L- c3 M4 w5 P) dfprintf('\nGradients at x           :  ');
fprintf(' %9.2f',grad);
fprintf('\nConstraint functions at x:  ');
! b( L3 o8 O2 Z% b$ Dfprintf(' %9.2f',ccon);
: O) e( V+ g! k6 x/ Cfprintf('\nNumber of major iterations  = %5d\n',majits);

  Z1 N: x' `& ?2 t: p5 Zfunction [mode, ccon, cjac, user] = ...
    confun(mode, ncnln, n, ldcj, needc, x, cjac, nstate, user)
! z6 R) H$ q1 ]! l4 s" A  ccon = zeros(ncnln, 1);
  if (nstate == 1)
%  first call to confun.  set all jacobian elements to zero.
%  note that this will only work when 'derivative level = 3'
%  (the default; see section 11.2).
    cjac = zeros(ncnln, n);
  end
  if (needc(1) > 0)
- P' e* z& X5 C4 L7 C# g  z: ?    if (mode == 0  ||  mode == 2)
      ccon(1) = x(1)^2 + x(2)^2 + x(3)^2 + x(4)^2;
& c  q+ \' y$ t& T9 [    end
5 T$ T: b% Z, v5 J  \9 {4 h    if (mode == 1  ||  mode == 2)
& b% h# n/ J0 K      cjac(1,1) = 2*x(1);
      cjac(1,2) = 2*x(2);
      cjac(1,3) = 2*x(3);
      cjac(1,4) = 2*x(4);
2 k( n7 _8 m' T2 _( w6 z" D& C    end
" U+ g- S# i& ^/ r  v- a0 ~  end
5 r8 X% ]7 G1 H* J5 N. @  R  if (needc(2) > 0)
# R" _# W; M/ c3 |' V2 \: Z; g5 E    if (mode == 0  ||  mode == 2)
+ f; F+ ^# l5 c5 z1 x: Q      ccon(2) = x(1)*x(2)*x(3)*x(4);
    end
  b  P( Z' Y2 `) d2 G! v( s    if (mode == 1  ||  mode == 2)
, I- g; j9 V& K5 I      cjac(2,1) = x(2)*x(3)*x(4);
      cjac(2,2) = x(1)*x(3)*x(4);
      cjac(2,3) = x(1)*x(2)*x(4);
. Y3 `1 |+ Z. ?$ w+ a7 G      cjac(2,4) = x(1)*x(2)*x(3);
    end
4 B) e- y+ D0 W" `; ]) |/ Q& F& p  end

0 v& J5 [  [! `( M4 ~, Vfunction [mode, objf, grad, user] = objfun(mode, n, x, grad, nstate, user)
  if (mode == 0 || mode == 2)
    objf = x(1)*x(4)*(x(1)+x(2)+x(3)) + x(3);
8 U% ^, H; ?/ C) a7 m5 `  end
  if (mode == 1 || mode == 2)
    grad(1) = x(4)*(2*x(1)+x(2)+x(3));
, {/ O( i5 I" Z: J    grad(2) = x(1)*x(4);
. t( {& r% V0 a  L- I- a! X  }    grad(3) = x(1)*x(4) + 1;
& R+ c5 A% {9 B7 j! T2 @# D9 r    grad(4) = x(1)*(x(1)+x(2)+x(3));
  end
' W" d' }( ^5 ~2 |' x' Q/ C+ n

yin123

来了，来了，专业赚积分顶起来

regngfpcb

哪里有问题啊