¿Cómo dividir un polinomio en otro usando NumPy en Python?

En este artículo, haremos un programa NumPy para dividir un polinomio en otro. Se dan dos polinomios como entrada y el resultado es el cociente y el resto de la división.

• El polinomio p(x) = C3 x2 + C2 x + C1  se representa en NumPy como: ( C1, C2, C3 ) { los coeficientes (constantes)}.
• Tomemos dos polinomios p(x) y g(x) y luego divídalos para obtener el cociente q(x) = p(x) // g(x) y el resto r(x) = p(x) % g(x) como resultado.

If p(x) = A3 x2 + A2 x + A1
and 
g(x) = B3 x2 + B2 x + B1 

then result is 
q(x) = p(x) // g(x) and r(x) = p(x) % g(x)

and the output is coefficientes of remainder and
 the coefficientes of quotient.

Esto se puede calcular utilizando el método polydiv() . Este

numpy.polydiv(p1, p2)

A continuación se muestra la implementación con algunos ejemplos:

Ejemplo 1 :  

# importing package
import numpy
  
# define the polynomials
# p(x) = 5(x**2) + (-2)x +5
px = (5, -2, 5)
  
# g(x) = x +2
gx = (2, 1, 0)
  
# divide the polynomials
qx, rx = numpy.polynomial.polynomial.polydiv(px, gx)
  
# print the result
# quotiient
print(qx)
  
# remainder
print(rx)

Producción :

[-12.   5.]
[ 29.]

Ejemplo 2: 

# importing package
import numpy
  
# define the polynomials
# p(x) = (x**2) + 3x + 2
px = (1,3,2)
  
# g(x) = x + 1
gx = (1,1,0)
  
# divide the polynomials
qx,rx = numpy.polynomial.polynomial.polydiv(px,gx)
  
# print the result
# quotiient
print(qx)
  
# remainder
print(rx)

Producción :

[ 1.  2.]
[ 0.]