Post Reply 
can the prime really not solve this integral?
02-27-2014, 06:47 AM
Post: #1
can the prime really not solve this integral?
I tried solving int(x^3/(e^3-1),x, 0,infinity) (this I ti89 notation, I did not type it like this) on the Prime without success. My nspire does it, even my ti89 does it without a problem. Please tell me this is a user error, otherwise I finally throw the prime out of the window.
Find all posts by this user
Quote this message in a reply
02-27-2014, 07:27 AM
Post: #2
RE: can the prime really not solve this integral?
Where is the problem? It is a simple cubic function and CAS hasn't any problem.

Massimo
from 34c to prime
Find all posts by this user
Quote this message in a reply
02-27-2014, 08:03 AM
Post: #3
RE: can the prime really not solve this integral?yy
This is what I got:

[Image: oswy78.jpg]

First without, then with pressing shift before entering.
Find all posts by this user
Quote this message in a reply
02-27-2014, 08:08 AM
Post: #4
RE: can the prime really not solve this integral?
Did you try +infinity ?

By the way, Emulator and CK permit to do screen shoots

Patrice
“Everything should be made as simple as possible, but no simpler.” Albert Einstein
Find all posts by this user
Quote this message in a reply
02-27-2014, 08:19 AM
Post: #5
RE: can the prime really not solve this integral?
I'm typing this on my phone while trying to sleep Smile
No change with the plus sign.


After more hacking on buttons I found out it works if you place an "approx" in front of the integral. I'm probably going to have nightmares of having to use the Prime on an exam.
Find all posts by this user
Quote this message in a reply
02-27-2014, 08:45 AM (This post was last modified: 02-27-2014 08:46 AM by eried.)
Post: #6
RE: can the prime really not solve this integral?
(02-27-2014 08:19 AM)DeucesAx Wrote:  I'm typing this on my phone while trying to sleep Smile
No change with the plus sign.


After more hacking on buttons I found out it works if you place an "approx" in front of the integral. I'm probably going to have nightmares of having to use the Prime on an exam.

hahaha... thanks the spaggetti monster that I have my trusty 50g as backup---- *won't turn on due inactive time* (every 2-3 days someone writes me thru my website telling me that his 50g does not turn on anymore after long inactivity periods)

My website: ried.cl
Visit this user's website Find all posts by this user
Quote this message in a reply
02-27-2014, 10:33 AM
Post: #7
RE: can the prime really not solve this integral?
(02-27-2014 08:19 AM)DeucesAx Wrote:  I'm typing this on my phone while trying to sleep Smile
No change with the plus sign.


After more hacking on buttons I found out it works if you place an "approx" in front of the integral. I'm probably going to have nightmares of having to use the Prime on an exam.

Did you mean to use e^x (instead of e^3, as in the op)?
Find all posts by this user
Quote this message in a reply
02-27-2014, 11:24 AM
Post: #8
RE: can the prime really not solve this integral?
(02-27-2014 06:47 AM)DeucesAx Wrote:  I tried solving int(x^3/(e^3-1),x, 0,infinity) (this I ti89 notation, I did not type it like this) on the Prime without success. My nspire does it, even my ti89 does it without a problem. Please tell me this is a user error, otherwise I finally throw the prime out of the window.

Might be user error.

With CAS exact unchecked get 6.49 for the definite integral. It does seem to hang in function app if you try to compute signed area for say x=0-100.

For indefinite solution I'm not sure.
Find all posts by this user
Quote this message in a reply
02-27-2014, 12:21 PM
Post: #9
RE: can the prime really not solve this integral?
You get a numeric answer with:
int(x^3/(exp(x)-1),x,0,inf)
then shift-enter.
x^3/(exp(x)-1) does not have an antiderivative than you can express with elementary function (special functions required, polylogs here). There is probably a trick than can give you the exact answer (pi^4/15) for the definite integral, any idea? On a voyage 200, you don't get the exact value by the way.
Find all posts by this user
Quote this message in a reply
02-27-2014, 12:37 PM
Post: #10
RE: can the prime really not solve this integral?
(02-27-2014 08:19 AM)DeucesAx Wrote:  I'm probably going to have nightmares of having to use the Prime on an exam.
I already mentioned that. I was initially thinking of offering the Prime to my daughter until I realized I had no intention to sabotage her studies.
As I write this, I have not yet been able to do the integration on your example...
Find all posts by this user
Quote this message in a reply
02-27-2014, 12:39 PM
Post: #11
RE: can the prime really not solve this integral?
(02-27-2014 12:21 PM)parisse Wrote:  You get a numeric answer with:
int(x^3/(exp(x)-1),x,0,inf)
then shift-enter.
x^3/(exp(x)-1) does not have an antiderivative than you can express with elementary function (special functions required, polylogs here). There is probably a trick than can give you the exact answer (pi^4/15) for the definite integral, any idea? On a voyage 200, you don't get the exact value by the way.

Thanks - I keep forgetting to use shift enter.

Hitting "a b/c" seems to convert 6.49 to the fraction equivalent to pi^4/15.
Find all posts by this user
Quote this message in a reply
02-27-2014, 02:11 PM
Post: #12
RE: can the prime really not solve this integral?
(02-27-2014 12:21 PM)parisse Wrote:  You get a numeric answer with:
int(x^3/(exp(x)-1),x,0,inf)
then shift-enter.
x^3/(exp(x)-1) does not have an antiderivative than you can express with elementary function (special functions required, polylogs here). There is probably a trick than can give you the exact answer (pi^4/15) for the definite integral, any idea? On a voyage 200, you don't get the exact value by the way.

Possibly a dumb question, but why is the result dependent on the Rad/Deg setting when solving this? If the calc is in Degrees the result (on my calculator) is 372.07.... Other integrals involving exp(x) (for instance int(exp(x^2),x,0,.5)) seems not to be influenced by the Deg/Rad setting. Am I missing something?

Cheers, Terje
Find all posts by this user
Quote this message in a reply
02-27-2014, 02:14 PM
Post: #13
RE: can the prime really not solve this integral?
(02-27-2014 02:11 PM)Terje Vallestad Wrote:  
(02-27-2014 12:21 PM)parisse Wrote:  You get a numeric answer with:
int(x^3/(exp(x)-1),x,0,inf)
then shift-enter.
x^3/(exp(x)-1) does not have an antiderivative than you can express with elementary function (special functions required, polylogs here). There is probably a trick than can give you the exact answer (pi^4/15) for the definite integral, any idea? On a voyage 200, you don't get the exact value by the way.

Possibly a dumb question, but why is the result dependent on the Rad/Deg setting when solving this? If the calc is in Degrees the result (on my calculator) is 372.07.... Other integrals involving exp(x) (for instance int(exp(x^2),x,0,.5)) seems not to be influenced by the Deg/Rad setting. Am I missing something?

Cheers, Terje

Do you also have complex mode on? Since \( e^{i\theta} = cos\theta + i \sin \theta \) the angle mode would affect the result.

Graph 3D | QPI | SolveSys
Find all posts by this user
Quote this message in a reply
02-27-2014, 02:29 PM
Post: #14
RE: can the prime really not solve this integral?
(02-27-2014 02:14 PM)Han Wrote:  Do you also have complex mode on? Since \( e^{i\theta} = cos\theta + i \sin \theta \) the angle mode would affect the result.

Yes I thought of that, but it does not seem to have any effect on the result whether complex mode is ticked or not, except for getting a warning message about complex answers when it is set.

Cheers, Terje
Find all posts by this user
Quote this message in a reply
02-27-2014, 03:47 PM
Post: #15
RE: can the prime really not solve this integral?
(02-27-2014 02:11 PM)Terje Vallestad Wrote:  Possibly a dumb question, but why is the result dependent on the Rad/Deg setting when solving this? If the calc is in Degrees the result (on my calculator) is 372.07.... Other integrals involving exp(x) (for instance int(exp(x^2),x,0,.5)) seems not to be influenced by the Deg/Rad setting. Am I missing something?
Cheers, Terje
Had the same question... so I took my 15C and tried both. As a result I always get pi^4/15 deg or rad which makes much more sense.
Find all posts by this user
Quote this message in a reply
02-27-2014, 05:44 PM
Post: #16
RE: can the prime really not solve this integral?
Hi,

I am just wondering, how did you integrate up to infinity on your HP-15C ?

I just curious, I have a weak practice of this calculator especially with the solver and the integration ... only use it with defined numeric intervals !
Find all posts by this user
Quote this message in a reply
02-27-2014, 07:08 PM
Post: #17
RE: can the prime really not solve this integral?
We almost never test the CAS in degree mode, therefore more bugs are expected there. The reason here is that for numeric integration with an infinite boundary you must make a change of variable, here we set tan(y)=x for y in 0..pi/2, dx=(1+tan(y)^2)*dy ... except that in degree mode you must add a pi/180 factor.
Find all posts by this user
Quote this message in a reply
02-27-2014, 09:11 PM
Post: #18
RE: can the prime really not solve this integral?
(02-27-2014 12:21 PM)parisse Wrote:  You get a numeric answer with:
int(x^3/(exp(x)-1),x,0,inf)
then shift-enter.
x^3/(exp(x)-1) does not have an antiderivative than you can express with elementary function (special functions required, polylogs here). There is probably a trick than can give you the exact answer (pi^4/15) for the definite integral, any idea? On a voyage 200, you don't get the exact value by the way.

No, no, you don't use primitives for this, come on. You could use contour integration/Residue Theorem, but it's a bit tricky for this. See for instance:

http://en.wikipedia.org/wiki/Stefan_bolt...w#Appendix

The easy way is:

http://math.stackexchange.com/questions/...u3-eu-1-du

That integral appears when you integrate the Planck distribution to all frequencies in order to recover the Stefan-Boltzmann law, also in the Debye theory of specific heat. The functions defined by these integrals are called Debye functions. See Abramowitz-Stegun §27.1

http://people.math.sfu.ca/~cbm/aands/page_998.htm

(Now that I'm thinking about it, it shouldn't be too difficult to implement a large class of definite integrals by the Residue theorem. That would be come in handy.)
Find all posts by this user
Quote this message in a reply
02-28-2014, 01:51 AM
Post: #19
RE: can the prime really not solve this integral?
(02-27-2014 11:24 AM)CR Haeger Wrote:  
(02-27-2014 06:47 AM)DeucesAx Wrote:  I tried solving int(x^3/(e^3-1),x, 0,infinity) (this I ti89 notation, I did not type it like this) on the Prime without success. My nspire does it, even my ti89 does it without a problem. Please tell me this is a user error, otherwise I finally throw the prime out of the window.

Might be user error.

With CAS exact unchecked get 6.49 for the definite integral. It does seem to hang in function app if you try to compute signed area for say x=0-100.

For indefinite solution I'm not sure.

Are you performing a signed area computation for X^3/(e^X-1), from 0 to 100, in the Function app Plot view? I just tried this on the emulator and on the device. Neither hanged for me — both quickly brought up a numerical approximation.

Can you reliably reproduce the hang?
Find all posts by this user
Quote this message in a reply
02-28-2014, 04:52 AM
Post: #20
RE: can the prime really not solve this integral?
(02-27-2014 07:08 PM)parisse Wrote:  We almost never test the CAS in degree mode, therefore more bugs are expected there. The reason here is that for numeric integration with an infinite boundary you must make a change of variable, here we set tan(y)=x for y in 0..pi/2, dx=(1+tan(y)^2)*dy ... except that in degree mode you must add a pi/180 factor.

Interestingly enough, both HP50g, and TI-92 Plus, give a numerically correct approximation to the exact answer (pi^4/15) independent of the degree/radiant setting :-)
Quote this message in a reply
Post Reply 




User(s) browsing this thread: 1 Guest(s)