Algorithms for partial fraction decomposition and rational function integration

• Authors: Ellis Horowitz
• Venue: SYMSAC '71 Proceedings of the second ACM symposium on Symbolic and algebraic manipulation
• Year: 1971

Algorithms for symbolic partial fraction decomposition and indefinite integration of rational functions are described. Two types of partial fraction decomposition are investigated, square-free and complete square-free. A method is derived, based on the solution of a linear system, which produces the square-free decomposition of any rational function, say A/B. The computing time is shown to be O(n4(1n nf)2) where deg(A) 4(1n nf)2). A thorough analysis is then made of the classical method for rational function integration, due to Hermite. It is shown that the most efficient implementation of this method has a computing time of O(k3n5(1n nc) 2), where c is a number closely related to f and k is the number of square-free factors of B. A new method is then presented which avoids entirely the use of partial fraction decomposition and instead relies on the solution of an easily obtainable linear system. Theoretical analysis shows that the computing time for this method is O(n5 (in nf) 2) and extensive testing substantiates its superiority over Hermite's method.