Comment #3 on issue 2713 by [email protected]: meijerint should use partial
fractions
http://code.google.com/p/sympy/issues/detail?id=2713
I started working on a branch where you can do the following:
In [4]: f = ((p**3 - 6*p + 1)/(p**4 + 4*p**3 + 3*p**2))
In [5]: f
Out[5]:
3
p - 6⋅p + 1
────────────────
4 3 2
p + 4⋅p + 3⋅p
In [6]: apart(f*exp(p), p)
Out[6]:
⎛ 4 3 22 1 ⎞ p
⎜───────── + ───── - ─── + ────⎟⋅ℯ
⎜9⋅(p + 3) p + 1 9⋅p 2⎟
⎝ 3⋅p ⎠
Note that apart(f*exp(p)) (no generators) still will give an exception,
treating f*exp(p) as bivariate rational function. It works this way because:
In [7]: apart(f.subs(p, exp(p)))
Out[7]:
-p -2⋅p
22⋅ℯ ℯ 4 3
- ────── + ───── + ────────── + ──────
9 3 ⎛ p ⎞ p
9⋅⎝ℯ + 3⎠ ℯ + 1
--
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sympy-issues?hl=en.
