1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
===================================
Test of polynomial routines
===================================
Real polynomials
Roots, r = [0.764873 -1.13073 0.452336]
Polynomial, p = [1 -0.0864783 -1.03036 0.39121]
r = roots(p) = [-1.13073+0i 0.764873+0i 0.452336+0i]
Roots, r = [-0.451817 1.36468 0.355365 -0.556704 1.00977 0.251912 -0.902459]
Polynomial, p = [1 -1.07075 -1.62686 1.27376 0.821026 -0.322887 -0.0949017 0.0280021]
r = roots(p) = [1.36468+0i 1.00977+0i -0.902459+0i -0.556704+0i -0.451817+0i 0.355365+0i 0.251912+0i]
x = [-0.946328 0.403829 -0.129499 0.28002 -0.439783 -0.0494359 -0.235643 -0.268905 -1.53793]
polyval(p, x) = [-0.0595982 0.00460212 0.0335058 -0.00121379 0.000936527 0.0318135 0.0265788 0.0226503 -16.9714]
Complex polynomials
Roots, r = [1.18706+0.721909i 0.887142+0.707331i -0.79531+0.260976i]
Polynomial, p = [1+0i -1.2789-1.69022i -1.48017+0.884713i 0.817695+1.03555i]
r = roots(p) = [-0.79531+0.260976i 1.18706+0.721909i 0.887142+0.707331i]
Roots, r = [0.432532-0.0565171i -0.0718804+0.727733i -0.975402-0.377971i 1.07181-0.0820541i 0.841014-0.0749022i 0.0729376-0.609413i -0.297478+0.210593i]
Polynomial, p = [1+0i -1.07353+0.262532i -0.453137-0.315718i 0.439976+0.459462i -0.224018-0.643275i 0.26218+0.288645i 0.0730715-0.0609737i -0.0642264+0.0215646i]
r = roots(p) = [-0.975402-0.377971i -0.0718804+0.727733i -0.297478+0.210593i 0.0729376-0.609413i 1.07181-0.0820541i 0.841014-0.0749022i 0.432532-0.0565171i]
x = [-0.20257+0.104277i -0.829977+0.0806219i 0.369791-0.393648i -0.117384+0.314781i 0.255673-0.524515i 0.343174+0.178492i 0.704202-0.0609598i -0.284082+0.668861i -0.108213+0.131127i]
polyval(p, x) = [-0.0437532+0.0357859i 0.359059+1.1608i -0.0374979-0.0709598i -0.0774064+0.000394322i -0.0923904-0.0431629i -0.0122133+0.0558953i 0.0213529-0.00735109i -0.332684+0.133051i -0.0574082+0.0252281i]
|
