10-Aug-2023 21:46:41

disk01_positive_rule_test():
  MATLAB/Octave version 5.2.0
  Test disk01_positive_rule().

disk01_positive_rule():
  Determine a quadrature rule for the unit positive disk.

  Seeking a rule using 3 points
  The degree of precision will be 2
  The number of degrees of freedom is 9
  The number of monomial integrals to match is 6

  Initial guess for W, X, Y values:

        0.261799        0.482963        0.353553
        0.261799        0.353553         0.12941
        0.261799         0.12941        0.482963

  FSOLVE solves for W, X, Y values:

        0.241185        0.572424        0.621705
        0.303028        0.621705        0.028522
        0.241185        0.028522        0.572424

  Residuals for WXY:

     1.11022e-16
    -1.12758e-07
    -1.12758e-07
     1.09836e-06
     1.09883e-06
     1.09836e-06

  Monomial       Estimate           Exact       Error

  x^0*y^0        0.785398        0.785398     1.11022e-16
  x^1*y^0        0.333333        0.333333     1.12758e-07
  x^0*y^1        0.333333        0.333333     1.12758e-07
  x^2*y^0        0.196351         0.19635     1.09836e-06
  x^1*y^1        0.125001           0.125     1.09883e-06
  x^0*y^2        0.196351         0.19635     1.09836e-06

disk01_positive_rule_test():
  Normal end of execution.

10-Aug-2023 21:46:41