share/cobyla/ex/ex2.mac

   1 /*
   2  * 2D unit circle
   3  *
   4  * f = x1*x2
   5  * 1-x1^2-x2^2 >= 0
   6  *
   7  * True solution: x1=1/sqrt(2), x2=-1/sqrt(2)
   8  */
   9
  10 fmin_cobyla(x1*x2, [x1, x2], [1,1], constraints = [1>=x1^2+x2^2], iprint=1);
  11
  12 /* Output
  13
  14    Normal return from subroutine COBYLA
  15
  16    NFVALS =   66   F =-5.000000E-01    MAXCV = 1.999956E-12
  17    X = 7.071058E-01  -7.071077E-01
  18 [[x1 = .7071058493484819,x2 = -.7071077130247994],-.4999999999992633,
  19          66, 0];
  20
  21 */
  22
  23 /* True output
  24
  25        Output from test problem 2 (2D unit circle calculation)
  26
  27    Normal return from subroutine COBYLA
  28
  29    NFVALS =   66   F =-5.000000E-01    MAXCV = 1.999956E-12
  30    X = 7.071058E-01  -7.071077E-01
  31
  32 */
  33
  34