First time here? Check out the FAQ!

Ask Your Question
0

integration ends up with hypergeometric function can not be done by sage

asked 11 years ago

chenming gravatar image

updated 11 years ago

kcrisman gravatar image

we know the integration of

1/(a^b+1)

with respect to a is going to get

2F1(1,1/a,1+1/a,-a^x)

where 2f1 is hypogeometric function. However, once we to that in sage, it is not possible to get results

sage: integrate(1/(x^b+1),x)  
integrate(1/(x^b + 1), x)

hope this can be fixed up latter or have some alternative way to work around.

Preview: (hide)

2 Answers

Sort by » oldest newest most voted
0

answered 11 years ago

Shashank gravatar image

It seems sage's default algorithm can't do it. You can use mathematica's algorithm to do it in sage. You don't need mathematica installed on you computer. Just type the following command.

b=var('b')
integral(1/(x^b+1),x,algorithm='mathematica_free')
Preview: (hide)
link
1

answered 11 years ago

kcrisman gravatar image

Correct - see Trac 2516 for the current status.

See this Bessel function doc for an example that might help you by using Maxima more directly in your own context.

sage: m = maxima(bessel_J(2, x))
sage: m.integrate(x)
hypergeometric([3/2],[5/2,3],-x^2/4)*x^3/24

Also related - Trac 9908.

Preview: (hide)
link

Comments

Though apparently Maxima can't do this particular integral either.

kcrisman gravatar imagekcrisman ( 10 years ago )

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

Stats

Asked: 11 years ago

Seen: 716 times

Last updated: Oct 16 '13