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

Uqzalq

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

function e04wc_example
8 V# q8 k4 s+ d
fprintf('e04wc example results\n\n');
a  = [1, 1, 1, 1];
bl = [1, 1, 1, 1, -1e25, -1e25,   25];
bu = [5, 5, 5, 5,    20,    40, 1e25];
istate = zeros(7, 1, 'int64');
% z% c5 I4 q* i+ Y* |ccon   = zeros(2,1);
0 n  L! ?$ D& ^2 w3 x) n' n" fcjac   = zeros(2,4);
clamda = zeros(7,1);
/ d5 w" b  t  E8 x, Y' l4 ~# |hess   = zeros(4,4);
x      = [1; 5; 5; 1];
[iw,rw,ifail]=e04wc;
9 G* ~; [- }- a# y[majits, istate, ccon, cjac, clamda, objf, grad, hess, x, ...
                                          iw, rw, user, ifail] = ...
* i, i, G6 q) N: o( C, K! ?8 E4 ~  e04wd(...
* m4 ~7 |7 y+ m! N3 @- a8 C' H* A        a, bl, bu, @confun, @objfun, istate, ccon, cjac, clamda, ...
' b3 I- y$ A- g5 ?        hess, x, iw, rw);
fprintf('Final objective value    :   %8.1f\n',objf);
fprintf('Optimal values for x     :  ');
fprintf(' %9.2f',x);
fprintf('\nGradients at x           :  ');
5 h# u. E3 G. y; t1 v- Rfprintf(' %9.2f',grad);
fprintf('\nConstraint functions at x:  ');
fprintf(' %9.2f',ccon);
fprintf('\nNumber of major iterations  = %5d\n',majits);
& L5 n. X& N3 C' A7 h& {
function [mode, ccon, cjac, user] = ...
) f" `! a! K& C. [9 ~    confun(mode, ncnln, n, ldcj, needc, x, cjac, nstate, user)
  ccon = zeros(ncnln, 1);
  if (nstate == 1)
. y' Y$ I1 A- x7 V0 D7 k0 f%  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);
2 P. c9 w- p8 p: S& D  end
+ o/ ?1 H) V- ?  if (needc(1) > 0)
    if (mode == 0  ||  mode == 2)
      ccon(1) = x(1)^2 + x(2)^2 + x(3)^2 + x(4)^2;
: e! N: u2 |7 m% ?) p: m* \    end
9 ?( I% {% l1 n$ J* k    if (mode == 1  ||  mode == 2)
& W9 }8 e* Q- e4 A8 Y% O      cjac(1,1) = 2*x(1);
      cjac(1,2) = 2*x(2);
1 [, f* |9 \& d" ^) w- w3 F( ?      cjac(1,3) = 2*x(3);
6 o2 ~3 y  R: K- }9 Z0 [$ Z+ K      cjac(1,4) = 2*x(4);
9 u0 c) s; f( Z& F  ]7 k  ?1 r    end
  end
+ W2 W6 ]' E! i  if (needc(2) > 0)
    if (mode == 0  ||  mode == 2)
9 [3 d" d: q% F* H# _6 g      ccon(2) = x(1)*x(2)*x(3)*x(4);
    end
- O3 v+ u' U6 Y7 i" M1 d    if (mode == 1  ||  mode == 2)
      cjac(2,1) = x(2)*x(3)*x(4);
      cjac(2,2) = x(1)*x(3)*x(4);
& f; D6 |( i# @( k) N1 x  ^      cjac(2,3) = x(1)*x(2)*x(4);
% H- S' E+ a6 W6 z, U# [& g# J4 d      cjac(2,4) = x(1)*x(2)*x(3);
* }; @& Y5 ]/ T% Q    end
  end
+ S4 p* j9 e- R2 p* i" @" ]
function [mode, objf, grad, user] = objfun(mode, n, x, grad, nstate, user)
( V; Q4 p, G9 O' [! t. F  if (mode == 0 || mode == 2)
    objf = x(1)*x(4)*(x(1)+x(2)+x(3)) + x(3);
! h  a9 Q3 N/ p0 Q: _3 C' `' M- K  end
& B4 p( H2 h1 ~: F* b& Z  if (mode == 1 || mode == 2)
& R' [$ T4 Q/ x. R# n, e; W7 s    grad(1) = x(4)*(2*x(1)+x(2)+x(3));
    grad(2) = x(1)*x(4);
    grad(3) = x(1)*x(4) + 1;
    grad(4) = x(1)*(x(1)+x(2)+x(3));
% d6 {% S. p9 O- [  end

yin123

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

regngfpcb

哪里有问题啊