GPLv3

Pari/GP pol-pfrac

pol-pfrac.gp splits a Pari/GP polynomial fraction into partial fractions which sum to the original. See the comments in pol-pfrac.gp for details.

read("pol-pfrac.gp");
f = 2/(x^2-1);
v = pol_partial_fractions(f);
\\ gives  v == [ 1/(x-1), -1/(x+1) ]
\\ with   vecsum(v) == f

pol-pfrac.gp is free software (free as in freedom), published under the terms of the GNU General Public License (v3 or higher). Download version 2 here,

pol-pfrac.gp (16k, and sig)
pol-pfrac-2.tar.gz (26k, and sig)

The tar file has some self-tests but just pol-pfrac.gp is enough to run. The sig files are Gnu PG ascii armoured signatures generated from my key.


This page Copyright 2015, 2016, 2017 Kevin Ryde, except for the GPLv3 logo which is Copyright Free Software Foundation and used here in accordance with its terms.

(Back to the sitemap, or the Pari/GP section there)