y=(2*x**4-14*x**3+51*x**2-89*x-82)/((x-1)*(x*x-4*x+5)^2)
z=y.partial_fraction()
z.integral(x)

Hirst-10.2 (last edited 2010-12-18 10:17:52 by chandrahas.kolla)