package("polmat")

math

 

The polynomial spectral factorisation of the polynomial a(s) will be computed.

a:=polmat(-s^6-s^4+s^2+1, s)

math

 

[x,n,eps]:=polmat::spf(a, s):

math

math

math

math

math

math

math

 

float(expr(x))

math

 

The polynomial spectral factorisation of the polynomial b(s) will be computed.

b:=polmat(-d^3+5*d^2-16*d+44)

math

 

[x,n,eps]:=polmat::spf(b, d):

 

float(expr(x)), n, eps

math