share/hompack/hompack_ex1.mac

   1 /*
   2  * This is the example for fortran/mainp.f which comes from the file fortran/innhp.dat
   3  */
   4 eqs:[-0.98e-3*x1^2 + 0.978e6*x2^2-0.98e1*x1*x2-0.235e3*x1+0.889e5*x2-1,
   5  -0.1e-1*x1^2 - 0.984*x2^2 -0.297e2*x1*x2+0.987e-2*x1-0.124*x2-0.25];
   6
   7 hompack_polsys(eqs, [x1, x2], epsbig = 0.1e-3, epssml = 0.1e-13, iflg1 = 11, sspar = [0,0,0,0,1,0,0,0]);
   8
   9 /*
  10  * Output:
  11 [0, [[x1 = 2342.338519591274 - 1.086941277045357e-11 %i,
  12 x2 = 3.83101739296796e-15 %i - 0.7883448240941404],
  13 [x1 = 1.684969554988813 %i + 0.01614785792343519,
  14 x2 = 0.004428029939736611 %i + 2.679947396144589e-4],
  15 [x1 = 0.01614785792343615 - 1.684969554988814 %i,
  16 x2 = 2.679947396144615e-4 - 0.004428029939736612 %i],
  17 [x1 = 0.09089212296153959 - 1.018053148133034e-15 %i,
  18 x2 = 9.436995207960155e-18 %i - 0.09114970981974996]], [1, 1, 1, 1],
  19 [1.0, 0.9999999999999996, 1.0, 1.0], [10.05533507452728, 1.721129364831411,
  20 2.023295971193018, 4.163263622709848], [53, 37, 35, 46]]
  21
  22  Compare with the Fortran output:
  23
  24   PATH NUMBER =    1
  25
  26   FINAL VALUES FOR PATH
  27
  28   ARCLEN = 0.100553350729765D+02
  29   NFE =   53
  30   IFLG2 =    1
  31  REAL SOLUTION
  32    FINITE SOLUTION
  33   LAMBDA = 0.100000000000000D+01
  34  X    = 0.234233851959128D+04-0.488863540872078D-11
  35  X    =-0.788344824094143D+00 0.180075471927527D-14
  36
  37   PATH NUMBER =    2
  38
  39   FINAL VALUES FOR PATH
  40
  41   ARCLEN = 0.172112936480935D+01
  42   NFE =   37
  43   IFLG2 =    1
  44  COMPLEX SOLUTION
  45    FINITE SOLUTION
  46   LAMBDA = 0.100000000000000D+01
  47  X    = 0.161478579234351D-01 0.168496955498881D+01
  48  X    = 0.267994739614459D-03 0.442802993973661D-02
  49
  50   PATH NUMBER =    3
  51
  52   FINAL VALUES FOR PATH
  53
  54   ARCLEN = 0.202329597129656D+01
  55   NFE =   35
  56   IFLG2 =    1
  57  COMPLEX SOLUTION
  58    FINITE SOLUTION
  59   LAMBDA = 0.100000000000000D+01
  60  X    = 0.161478579234362D-01-0.168496955498881D+01
  61  X    = 0.267994739614462D-03-0.442802993973661D-02
  62
  63   PATH NUMBER =    4
  64
  65   FINAL VALUES FOR PATH
  66
  67   ARCLEN = 0.416326362289792D+01
  68   NFE =   46
  69   IFLG2 =    1
  70  REAL SOLUTION
  71    FINITE SOLUTION
  72   LAMBDA = 0.100000000000000D+01
  73  X    = 0.908921229615396D-01-0.929821875294837D-15
  74  X    =-0.911497098197500D-01-0.424664784358207D-16
  75
  76 */