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

Uqzalq

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

' X8 j7 p6 |3 e' w
function e04wc_example

fprintf('e04wc example results\n\n');
5 B, A* p- Z+ Oa  = [1, 1, 1, 1];
" _% b) s3 x* H: c1 ]- vbl = [1, 1, 1, 1, -1e25, -1e25,   25];
: w: Z1 l" E9 l1 nbu = [5, 5, 5, 5,    20,    40, 1e25];
istate = zeros(7, 1, 'int64');
ccon   = zeros(2,1);
cjac   = zeros(2,4);
* \% L7 n/ A4 W: v7 y" oclamda = zeros(7,1);
hess   = zeros(4,4);
x      = [1; 5; 5; 1];
" h0 c: }' B% O. o7 k# I: g8 D[iw,rw,ifail]=e04wc;
; T# ]% d- Q. g4 |4 e[majits, istate, ccon, cjac, clamda, objf, grad, hess, x, ...
                                          iw, rw, user, ifail] = ...
  e04wd(...
9 I, p7 z( F$ G3 {, R        a, bl, bu, @confun, @objfun, istate, ccon, cjac, clamda, ...
% e& Z2 u4 O' _* D" A, L        hess, x, iw, rw);
fprintf('Final objective value    :   %8.1f\n',objf);
1 \* p8 _7 H, u' Hfprintf('Optimal values for x     :  ');
2 Y5 u* |; I3 {4 |/ Dfprintf(' %9.2f',x);
fprintf('\nGradients at x           :  ');
5 v9 \5 c2 `% S: w7 m' gfprintf(' %9.2f',grad);
fprintf('\nConstraint functions at x:  ');
fprintf(' %9.2f',ccon);
! `1 G6 x3 t+ [0 L7 i5 dfprintf('\nNumber of major iterations  = %5d\n',majits);

function [mode, ccon, cjac, user] = ...
3 O: p9 Y; k; T- l& f    confun(mode, ncnln, n, ldcj, needc, x, cjac, nstate, user)
: W( x& H7 Y1 h) F# Y' f2 F# o1 m  ccon = zeros(ncnln, 1);
  if (nstate == 1)
%  first call to confun.  set all jacobian elements to zero.
; y0 d( M$ L: i6 _! ^3 P; L%  note that this will only work when 'derivative level = 3'
; W/ @& R! E7 v- z* f%  (the default; see section 11.2).
    cjac = zeros(ncnln, n);
  end
  if (needc(1) > 0)
    if (mode == 0  ||  mode == 2)
      ccon(1) = x(1)^2 + x(2)^2 + x(3)^2 + x(4)^2;
    end
    if (mode == 1  ||  mode == 2)
& I: @; m! v+ u0 S: e      cjac(1,1) = 2*x(1);
      cjac(1,2) = 2*x(2);
6 K( |4 J* x) c      cjac(1,3) = 2*x(3);
' A  @) `9 O( r( y' f" d( P      cjac(1,4) = 2*x(4);
+ }+ D: ]5 `2 V; s2 i; L    end
  end
  if (needc(2) > 0)
    if (mode == 0  ||  mode == 2)
      ccon(2) = x(1)*x(2)*x(3)*x(4);
) p! f$ f7 l  r5 T" m! P2 c    end
    if (mode == 1  ||  mode == 2)
. _* q" s$ w- m' j7 ^& j3 a; J      cjac(2,1) = x(2)*x(3)*x(4);
      cjac(2,2) = x(1)*x(3)*x(4);
$ E: F0 k. c& l( m; Z% U      cjac(2,3) = x(1)*x(2)*x(4);
      cjac(2,4) = x(1)*x(2)*x(3);
: v# E  c+ `* `3 d/ y5 Z    end
! N: `2 N  j. e- c7 y# |2 D  end
8 w! B% N8 N8 ?
function [mode, objf, grad, user] = objfun(mode, n, x, grad, nstate, user)
! d: Z* O+ U/ j# Z! K* |9 s5 s& M$ ]  if (mode == 0 || mode == 2)
1 _) Z. V" _: H7 x    objf = x(1)*x(4)*(x(1)+x(2)+x(3)) + x(3);
- r" }5 y( q' l0 f9 M  end
  if (mode == 1 || mode == 2)
    grad(1) = x(4)*(2*x(1)+x(2)+x(3));
    grad(2) = x(1)*x(4);
    grad(3) = x(1)*x(4) + 1;
$ Z- Y  `8 j: W2 w9 L& m9 J9 k" ^$ h    grad(4) = x(1)*(x(1)+x(2)+x(3));
/ l, v! @  X" f- H+ d  end

regngfpcb

哪里有问题啊

yin123

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