Power series computations

pm5 has asked for the wisdom of the Perl Monks concerning the following question:

Hi. I want to do some ordinary algebraic operations (addition and multiplication) on power series over real or complex numbers. There is no real world purpose except to give myself several examples to play with. In fact, I'll probably just need exp(x) and some trigonometric functions to begin with. A quick CPAN search didn't shed a light, though. Is there any modules that might be helpful?

Comment on Power series computations

Replies are listed 'Best First'.
Re: Power series computations by McDarren (Abbot) on Feb 14, 2006 at 03:22 UTC
Maybe you aren't using appropriate search terms. I tried "math", "algebraic" and "trigonometric" and got dozens of hits. For example: Math::Series Math::Calculus::Expression Math::Trig Math::Algebra::Symbols etc, etc... Cheers, Darren :)	[reply]
Re: Power series computations by rhesa (Vicar) on Feb 14, 2006 at 03:26 UTC
You might want to check out PARI and the wonderful Math::Pari perl wrapper, which lives on CPAN and optionally installs PARI for you. Here's a quick example run in the gp calculator: ? taylor(exp(y),y) %366 = 1 + y + 1/2y^2 + 1/6y^3 + 1/24y^4 + 1/120y^5 + 1/720y^6 + +1/5040y^7 + 1/40320y^8 + 1/362880y^9 + 1/3628800y^10 + 1/39916800 +y^11 + 1/479001600y^12 + 1/6227020800y^13 + 1/87178291200y^14 + 1 +/1307674368000y^15 + O(y^16) ? taylor(sin(y),y) %374 = y - 1/6y^3 + 1/120y^5 - 1/5040y^7 + 1/362880y^9 - 1/3991680 +0y^11 + 1/6227020800y^13 - 1/1307674368000y^15 + O(y^16) ? taylor(exp(y),y)+taylor(sin(y),y) %375 = 1 + 2y + 1/2y^2 + 1/24y^4 + 1/60y^5 + 1/720y^6 + 1/40320y +^8 + 1/181440y^9 + 1/3628800y^10 + 1/479001600y^12 + 1/3113510400* +y^13 + 1/87178291200*y^14 + O(y^16) [download]	[reply] [d/l]
Re: Power series computations by spiritway (Vicar) on Feb 14, 2006 at 03:32 UTC
You've got atan2, cos, exp, log, sin, and sqrt, all in the core distribution. To perform exponentiation, use x ** y. If you don't want to roll your own complex numbers, you can check out Math::Complex. For the missing trig functions, you can easily make functions (e.g., `tan(x)=sin(x)/cos(x)`, etc.) to do it. For more complete modules, you might want to tryMath::Taylor or Math::Calculus::TaylorSeries.	[reply] [d/l]