import sympy as sp
import numpy as np
from IPython.display import display

sp.init_printing(pretty_print=True)

sp.var("x")

f1=(x-1)**2
f2=(x-2)
ff=f1*f2

print("f(x)= ")
display(ff)

ab_d = [[-2,-1],[-2,0],[-2,1],[-2,2],[-1,0]]
for ab in ab_d:
    hh=x**2+ab[0]*x+ab[1]
    r1=sp.div(hh**7,ff)
    print("a=",ab[0]," b=",ab[1]," のとき")
    print("h(x)=")
    display(hh)
    print("h1(x)=")
    display(r1[1])
    h2=sp.div(r1[1]**7,ff)
    print("h2(x)=")
    display(h2[1])