For given
p(x) = PROD[i=1,m]{SUM[j=2,n+2]{(A(i,j)*x^(j2))^A(i,1)}}
we shall get
p(x) = SUM[s=1,N+1]{p(s)^(N+1s)}
For example
If
p(x) = (x4)^5 * (3x^67x^3+5x+2)^2 * (x^3+8)^3 * x^2
or
A = [ 5 4 1 0 0 0 0 0
2 2 5 0 7 0 0 3
3 8 0 0 1 0 0 0
1 0 0 1 0 0 0 0 ]
then from
p = polyget(A)
we get
p = [ 9 180 1440 5586 .... 7864320 209715 0 0 ]
or
p(x) = 9x^28180x^27+1440x^265586x^25+ ... 7864320x^32097152x^2.
This routine is mainly to be used for creating test polynomials to
(a) determine the polynomial GCD of a pair of polynomials,
(b) find the roots with muliplicities of a given polynomial.
References in MATLAB Central:
(1) "GCD of polynomials,"
File ID 20859, 12 Apr 2009
(2) "Factorization of a polynomial with multiple roots,"
File ID: 20867, 27 Jul 2008
(3) "Multipleroots polynomial solved by partial fraction expansion,"
File ID: 22375, 10 Dec 2008
F C Chang 04/25/09
