Kevin Ryde on Sat, 21 Jul 2018 10:42:36 +0200
|
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
generating function solution to poly
|
- To: pari-users@pari.math.u-bordeaux.fr
- Subject: generating function solution to poly
- From: Kevin Ryde <user42_kevin@yahoo.com.au>
- Date: Sat, 21 Jul 2018 18:38:30 +1000
- Delivery-date: Sat, 21 Jul 2018 10:42:36 +0200
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=yahoo.com.au; s=s2048; t=1532162553; bh=N0a2TssUS1AZVmhMJwGxDxH7moQ6kkv4cSaM5j9KrpI=; h=From:To:Subject:Date:From:Subject; b=JrZ3Ebt1Anik0+oCIFGd8sF4G6bS/L3/10aHgwKpZeonasu4+59Dcp3mFsq0SGXpoVsdDQhD4uklijdo2XsGenP4iJUajVGO6m9k+DgrMTZdcZ624ibnHg38hG5/l6wyl0NXvsCd2V5DrRCi9qwplI9ycvx3JOy4E3zW2TbTfjp0ZTJOnomT89iK+jYDZJqfAL1iIi2nlkBhjoohrKtq8L/HfgCRA1yNZVYYc30bxjvD/pweeaOa5xXYSa8+iyOVC+WvzFo7c/tspYsHKBqSUsneRHskqbO2p8it2ZYdfJBV8jfjAst6QYF0ztWwV3BXVIc8EEWieRgTDffvcFxzCA==
- Organization: Bah Humbug
- User-agent: Gnus/5.13 (Gnus v5.13) Emacs/24.5 (gnu/linux)
I have a generating function (and more terms too)
g = x^2 + x^3 + x^4 + 3*x^5 + 6*x^6 + 12*x^7 + 29*x^8 + 67*x^9 + O(x^10);
which satisfies a cubic
(1+x)*g^3 - 2*g^2 + (1-x+2*x^2)*g - x^2 == 0
Is there an easy or good way to have gp solve that for series g?
I know how to work upwards to get g term by term (after deciding lowest
should be 0), but maybe gp already has it. I wondered only for interest
or generality though, since this one comes from recurrences which are
easy to calculate. I attempted serreverse() without joy (change
variables to solve in x, but I may have confused myself).