1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
|
===================================
Test of polynomial routines
===================================
Real polynomials
Roots, r = [-0.283885 1.364462 0.906001]
Polynomial, p = [1.000000 -1.986579 0.591655 0.350939]
r = roots(p) = [1.364462+0.000000i 0.906001+0.000000i -0.283885+0.000000i]
Roots, r = [1.249419 1.557830 0.564970 -0.596145 -1.197318 -1.267564 -0.300385]
Polynomial, p = [1.000000 -0.010807 -3.896433 -0.597658 4.298162 1.153145 -1.021590 -0.298857]
r = roots(p) = [1.557830+0.000000i 1.249419+0.000000i -1.267564+0.000000i -1.197318+0.000000i 0.564970+0.000000i -0.596145+0.000000i -0.300385+0.000000i]
x = [-1.074602 2.364758 2.365075 1.228044 0.234136 -0.732016 -2.105534 -0.923457 -1.892274]
polyval(p, x) = [-0.087985 165.384526 165.592625 0.078881 -0.424166 -0.085999 -68.059907 -0.154228 -23.858406]
Complex polynomials
Roots, r = [-0.469741-0.614789i -0.035692+0.565461i 0.727250-0.636742i]
Polynomial, p = [1.000000+0.000000i -0.221817+0.686070i -0.034581+0.042280i -0.109854+0.409246i]
r = roots(p) = [0.727250-0.636742i -0.469741-0.614789i -0.035692+0.565461i]
Roots, r = [-0.160972-1.299157i 0.036193+1.019169i -1.046929+0.760191i -0.280465+0.361084i 0.780099+0.643584i 0.794379+1.032883i -0.190616-0.710921i]
Polynomial, p = [1.000000+0.000000i 0.068312-1.806834i -0.158305-0.930377i -1.202161-3.461733i -3.908780+0.529733i -1.009816+0.365843i -1.368117+1.443013i 0.265664+0.718430i]
r = roots(p) = [-0.160972-1.299157i -1.046929+0.760191i -0.190616-0.710921i 0.794379+1.032883i 0.780099+0.643584i 0.036193+1.019169i -0.280465+0.361084i]
x = [1.323672+0.349763i 0.510483+0.174148i 0.475202-0.708159i -0.127081+0.388325i -0.572027-1.624474i -0.320883+0.183600i -1.023713+0.214321i -0.983530+0.261137i 1.133963-0.225782i]
polyval(p, x) = [10.010566-7.090964i -1.057338+0.508034i 6.313101+3.020747i -0.067120+0.186664i 33.098080-42.921040i 0.345802-0.000105i -1.501636-4.154270i -1.503561-3.174843i -17.655443-0.481204i]
Chebyshev polynomial
x = [1.523584 -1.251292 0.658525 -0.859404 0.075397 -0.172018 -1.200352 -1.822961]
cheb(10, x) = [9312.529843 520.885927 -0.617377 0.608798 -0.728487 0.157326 253.608153 88254.183938]
cheb(15, x) = [1270915.019652 -16812.346582 0.977463 0.195158 -0.905276 0.521343 -5711.593426 -37078141.121001]
|
