Home License -- for personal use only.  Not for government,
	academic, research, commercial, or other organizational use.

12-Jan-2026 11:03:20

wdk_test():
  MATLAB/Octave version: 9.11.0.2358333 (R2021b) Update 7
  wdk() uses the Weierstass-Durand-Kerner algorithm
  to compute all roots of a polynomial.

wdk_test01():
  Try wdk() on a degree four polynomial with a pair of complex roots.

  Polynomial to be analyzed: = 
     +1.000000 * x^4
    -26.000000 * x^2
    -75.000000 * x
    -56.000000

  Iter  Change  ||Poly(root)||

     1       28.158   1.0533e+07
     2       17.182   2.7767e+06
     3       12.511   6.5607e+05
     4        8.852   1.4844e+05
     5        6.055        33066
     6       4.1884       7330.6
     7        3.157       1652.9
     8       2.6402       415.79
     9        1.736       115.72
    10      0.77883       10.237
    11      0.31326      0.69565
    12     0.041028    0.0040484
    13   0.00023829   4.1498e-08
    14   1.1505e-09   2.2737e-13
    15   7.1089e-16   2.2737e-13

  Computed roots:
   6.2749 + 0.0000i  -1.2749 + 0.0000i  -2.5000 + 0.8660i  -2.5000 - 0.8660i

  Exact roots:
   6.2749 + 0.0000i  -1.2749 + 0.0000i  -2.5000 + 0.8660i  -2.5000 - 0.8660i


wdk_test():
  Normal end of execution.
12-Jan-2026 11:03:20