share/cobyla/ex/ex6.mac

   1 /*
   2  * Equation 9.1.15 in Fletcher's book Practical Methods of
   3  * Optimization
   4  *
   5  * f = -x1-x2
   6  * x2-x1^2 >= 0
   7  * 1-x1^2-x2^2 >= 0
   8  *
   9  * True solution: x1=1/sqrt(2), x2=1/sqrt(2)
  10  */
  11
  12 fmin_cobyla(-x1-x2, [x1, x2], [1, 1], constraints = [x2-x1^2 >= 0, 1-x1^2-x2^2 >= 0], iprint=1);
  13
  14 /* Output
  15
  16    Normal return from subroutine COBYLA
  17
  18    NFVALS =   63   F =-1.414214E+00    MAXCV = 2.000067E-12
  19    X = 7.071064E-01   7.071072E-01
  20 [[x1 = .7071063653426506,x2 = .7071071970316142],-1.414213562374265,
  21          63, 0]
  22 */
  23
  24 /* True output
  25
  26        Output from test problem 6 (Equation (9.1.15) in Fletcher)
  27
  28    Normal return from subroutine COBYLA
  29
  30    NFVALS =   63   F =-1.414214E+00    MAXCV = 2.000067E-12
  31    X = 7.071064E-01   7.071072E-01
  32
  33 */
  34
  35