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

Uqzalq

如下为NAG library里面e04wc的程序，是别人文章里给的，用来求某函数最小值的，但是我运行时会说错误使用e04wc，输出参数太多，球球大佬怎么搞？
& Z2 E( s4 e! z; ^& y  C
, v; e% j& S, u( V
function e04wc_example

: T* r8 N7 A$ c5 F  ^4 `- Rfprintf('e04wc example results\n\n');
a  = [1, 1, 1, 1];
0 y4 T$ e* h2 p3 l9 ~, Tbl = [1, 1, 1, 1, -1e25, -1e25,   25];
bu = [5, 5, 5, 5,    20,    40, 1e25];
4 P" F' R/ N% R! O* tistate = zeros(7, 1, 'int64');
ccon   = zeros(2,1);
cjac   = zeros(2,4);
' a2 [% @- S! Rclamda = zeros(7,1);
hess   = zeros(4,4);
x      = [1; 5; 5; 1];
' H. w- U+ M8 K2 _9 N& [+ G  t[iw,rw,ifail]=e04wc;
[majits, istate, ccon, cjac, clamda, objf, grad, hess, x, ...
$ X6 M% D$ x" ^                                          iw, rw, user, ifail] = ...
/ H" X1 b/ T8 D3 C9 K  m  e04wd(...
6 I6 l& x! Z/ P0 d4 _        a, bl, bu, @confun, @objfun, istate, ccon, cjac, clamda, ...
        hess, x, iw, rw);
fprintf('Final objective value    :   %8.1f\n',objf);
fprintf('Optimal values for x     :  ');
: ^. z# {, u% ]$ a1 Rfprintf(' %9.2f',x);
fprintf('\nGradients at x           :  ');
fprintf(' %9.2f',grad);
  F! u2 z0 I# R* Yfprintf('\nConstraint functions at x:  ');
; N0 y6 [# d8 C, Dfprintf(' %9.2f',ccon);
fprintf('\nNumber of major iterations  = %5d\n',majits);
3 f* m# C' O" m
function [mode, ccon, cjac, user] = ...
    confun(mode, ncnln, n, ldcj, needc, x, cjac, nstate, user)
  ccon = zeros(ncnln, 1);
  if (nstate == 1)
%  first call to confun.  set all jacobian elements to zero.
8 k9 ?6 v% F/ n%  note that this will only work when 'derivative level = 3'
0 {7 x* P& @9 R1 b%  (the default; see section 11.2).
& Z/ c! q/ y6 ?1 O. Z& ^' R    cjac = zeros(ncnln, n);
0 X* s: d8 L) N) {# d, R  end
  if (needc(1) > 0)
    if (mode == 0  ||  mode == 2)
% A- n( U0 t) A      ccon(1) = x(1)^2 + x(2)^2 + x(3)^2 + x(4)^2;
7 T4 G3 H+ T7 i: n    end
    if (mode == 1  ||  mode == 2)
" X7 [; w. U8 R      cjac(1,1) = 2*x(1);
      cjac(1,2) = 2*x(2);
      cjac(1,3) = 2*x(3);
      cjac(1,4) = 2*x(4);
9 e0 d, \! W8 S( H    end
% @. [. R* S- E5 M8 k  end
$ m! o- Z; N; W9 y# z* j& t  z4 e  if (needc(2) > 0)
    if (mode == 0  ||  mode == 2)
( ]# I1 g- P: M0 G0 I      ccon(2) = x(1)*x(2)*x(3)*x(4);
    end
    if (mode == 1  ||  mode == 2)
/ s* `' \- b7 |' \$ y- U9 f      cjac(2,1) = x(2)*x(3)*x(4);
      cjac(2,2) = x(1)*x(3)*x(4);
' l- y- T! @8 v- t" x) @4 ?0 S1 C      cjac(2,3) = x(1)*x(2)*x(4);
! L7 y3 O% n7 U( ?$ k: m      cjac(2,4) = x(1)*x(2)*x(3);
    end
  end

function [mode, objf, grad, user] = objfun(mode, n, x, grad, nstate, user)
$ u9 b! f2 ~& r  if (mode == 0 || mode == 2)
    objf = x(1)*x(4)*(x(1)+x(2)+x(3)) + x(3);
9 N; z0 f5 J* ^4 J0 _+ V$ k# N  end
9 p4 W1 U' ^" a& @; I7 [  if (mode == 1 || mode == 2)
2 k1 ^9 O9 q* C/ b    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));
  end
! ?* R" Q2 ^% s, _8 {3 [

yin123

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

regngfpcb

哪里有问题啊