Ask Your Question
0

How to define an integral to not be evaluated

asked 4 years ago

tomcho gravatar image

I'm trying to write the following script:

var('alpha beta x')

num = x**(alpha-1) * (1-x)**(beta-1)
denom = integral(x**(alpha-1) * (1-x)**(beta-1), x)

gammadist = num / denom

diff(gammadist, x)

So it is mich easier if the integral on denom does not get evaluated. However, I cannot find a way to do that on Sage and Sage always tries to evaluate it (without success). Is there a way to tell Sage to keep it as an integral?

Thanks

Preview: (hide)

1 Answer

Sort by » oldest newest most voted
2

answered 4 years ago

rburing gravatar image

updated 4 years ago

See the documentation of integral:

To prevent automatic evaluation use the hold argument.

Indeed,

sage: var('alpha beta x')
sage: denom = integral(x**(alpha-1) * (1-x)**(beta-1), x, hold=True); denom
integrate(x^(alpha - 1)*(-x + 1)^(beta - 1), x)
sage: diff(denom,x)
x^(alpha - 1)*(-x + 1)^(beta - 1)
Preview: (hide)
link

Your Answer

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

Add Answer

Question Tools

1 follower

Stats

Asked: 4 years ago

Seen: 287 times

Last updated: Jul 06 '20