预订演示
前页 后页

poldiv

Polynomial Division

SYNOPSIS:

int maxpol, na, nb, nc;
double a[], b[], c[];

nc = na + nb
polini( MAXPOL );
i = poldiv( a, na, b, nb, c );


DESCRIPTION:

i = poldiv( a, na, b, nb, c );  c = b / a,  nc = MAXPOL

returns i = the degree of the first nonzero coefficient of a.
The computed quotient c must be divided by x^i. 
An error message is printed if a is identically zero.


a, b, c are polynomials of degree na, nb, nc respectively.
The degree of a polynomial cannot exceed a run-time value MAXPOL.
An operation that attempts to use or generate a polynomial of higher degree might produce a result that suffers truncation at degree MAXPOL.
The value of MAXPOL is set by calling the function
polini( MAXPOL );

Each polynomial is represented by an array containing its coefficients, together with a separately declared integer equal to the degree of the polynomial. 
The coefficients appear in ascending order; that is,

                                        2                      na
a(x)  =  a[0]  +  a[1] * x  +  a[2] * x   +  ...  +  a[na] * x  .