Kevin Ryde on Wed, 21 Jan 2015 23:20:00 +0100
|
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
Re: polynomial partial fractions
|
- To: pari-users@pari.math.u-bordeaux.fr
- Subject: Re: polynomial partial fractions
- From: Kevin Ryde <user42_kevin@yahoo.com.au>
- Date: Thu, 22 Jan 2015 09:18:35 +1100
- Delivery-date: Wed, 21 Jan 2015 23:20:00 +0100
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=yahoo.com.au; s=s2048; t=1421878792; bh=35ILkKkM/Pm2MBMuzArgpSziVvYzJyasZSk7NahWirs=; h=From:To:Subject:References:Date:From:Subject; b=HwU1JqAYJZq3HpkLHsR6t0Qsmbn4Sbuy/Vk6x+bOrLEHfS7uAH+LwA8oT8VS/xxp8jryNeC0DFYH9LoD+vahv8x1lu1CWhYxAmmtXMa4TBKemYPsXWGL/ikinVsBZSHoaqTQUSLNhxGptKi/YaM5/gjhOJE/aiduuWnUkJsGI03hXIbk8G7tKuonNLdsOHx31Lm6naG2MDTc6JkGMSAzPuWsui9JyoZKbuzbLs2ecB7UQs9SFdyhFTmH5JAUwMV4Yz3O/t2o9Kv45RVNBgKgHN/T8DWSRiQbTCf6ksCl+lwm0tm9RvqISvLajSjAOZ4TDDPpsHuzEL5XM2zG324BIw==
- References: <87oarafuxi.fsf@blah.blah> <20141211094419.GB7418@math.u-bordeaux1.fr> <878ui7x5j6.fsf@blah.blah> <20141216145632.GB13201@math.u-bordeaux1.fr> <87fvcbxkcv.fsf@blah.blah>
- User-agent: Gnus/5.13 (Gnus v5.13) Emacs/24.4 (gnu/linux)
I got to these few lines. Limited to exact coefficients due to an
add-back self test but maybe that could be loosened if there would be a
use for inexact fracs.
http://user42.tuxfamily.org/pari-pol-pfrac/index.html
When using it with some generating functions I had the urge to print
polynomials with powers ascending and the denominator starting 1 (and
fractions elsewhere as necessary). Eg.
print( (8*x + 3)/(-48*x^2 + 24) )
=>
(1/8 + 1/3*x)/(1 - 2*x^2)
I tried some rough bits for this, but does an option already exist and
is hiding?