Status: New
Owner: ----
Labels: Type-Defect Priority-Medium

New issue 3586 by [email protected]: Exponentiation Integration involving polynomial term with mult/added coefficients
http://code.google.com/p/sympy/issues/detail?id=3586

The following commands did not complete the integration process:
     import sympy as sy
     alpha1, alpha2 = sy.symbols('alpha1 alpha2', positive=True)
ExpV1a = sy.integrate(sy.exp(2*alpha1**2*x**2 + (4*alpha1*alpha2 - 2)*x + 2*alpha2**2), x) ExpV1b = sy.integrate(sy.exp(2*alpha1**2*x**2 + 4*alpha1*alpha2*x + 2*alpha2**2)*sy.exp(-2*x), x) ExpV1c = sy.integrate(sy.exp(2*alpha1**2*x**2 + 4*alpha1*alpha2*x + 2*alpha2**2)*sy.exp(-2*x), x)
     ExpV1d = sy.integrate(sy.exp(2*(alpha1*x + alpha2)**2)*sy.exp(-2*x), x)

The following command variation performed the equivalent operation correctly:
     alpha3 = sy.symbols('alpha3', positive=True)
ExpV2a = sy.integrate(sy.exp(2*alpha1**2*x**2 + alpha3*x + 2*alpha2**2), x)
     ExpV2a = ExpV2a.subs(alpha3, (4*alpha1*alpha2 - 2))

[email protected] commented on an earlier communication regarding this issue that the "risch_integrate() won't have a hope of computing this until the day that erf() extensions are added to the code ...".


--
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.

Reply via email to