For the division of univariate polynomials, given the dividend b(x) of degree n and the divisor a(x) of degree m, such that b(x) = q(x)*a(x) + r(x), the quotient q(x) of degree n-m and the remainder r(x) of degree m-1 are then obtained.

This very short code is similar to the MATLAB built-in function 'deconv.m'.

The presented routine [q,r] = poly_div(b,a) can avoid the un-wanted data preceeding the remaider's coefficients applying the built-in routine [q,r] = deconv(b,a). See author's self comment.