An HP48/49 CAS question « Next Oldest | Next Newest »

 ▼ Valentin Albillo Posting Freak Posts: 1,755 Threads: 112 Joined: Jan 2005 10-17-2003, 05:01 AM Hi everybody. Not owning any HP48/49 model, there's no immediate way for me to test this, so I'd appreciate if any of you HP48/49 owners would try this for me and post here the results. The question is: what result does the HP48/49 symbolic algebra system produce for the following indefinite integral: ``` Integral[ sin(x^n) . dx ] ``` where n is a constant ? Specifically, I'm interested to know if it does produce some kind of exact closed result, in terms of known functions (elementary or not). I'm not interested in Taylor-series-expansion based results or any other non-closed approximations (not to mention purely numerical results). In case it can't produce a closed result for general, arbitrary n, I would then be interested to know if it can produce a closed result for some particular values of n, and if so, the specific values of n it can solve. I would expect it to be able to produce a closed result for n = 0 and n = 1 at the very least, but what about other values of n, both integer and non-integer (n = 2, 3, ..., 1/2, 1/3, ...) ? What about negative values of n (n = -1, -2, ..., -1/2, ...) ? Thanks in advance for any results or comments and best regards from V. ▼ Werner Huysegoms Junior Member Posts: 44 Threads: 0 Joined: Jan 1970 10-17-2003, 07:23 AM Hi Valentin. Firstly, the 48 does not dispose of a CAS worth mentioning; And the 49 cannot solve any but the most trivial of these integrals (ie with N=0 or 1). All the rest are simply returned in integral form. Not a great surprise, as for N=2 the result is a Fresnel integral, and for larger n it is expressible as complex incomplete gamma functions, both of which are not builtin to the 49G. Cheers, Werner ▼ Valentin Albillo Posting Freak Posts: 1,755 Threads: 112 Joined: Jan 2005 10-17-2003, 09:26 AM Werner posted: "The 49 cannot solve any but the most trivial of these integrals (ie with N=0 or 1). All the rest are simply returned in integral form." Thanks for replying, Werner, but have you actually tried it on an HP-49, or is your comment just what you think based on your experience ? I would need an actual attempt on the calculator. "Not a great surprise, as for N=2 the result is a Fresnel integral, and for larger n it is expressible as complex incomplete gamma functions, both of which are not builtin to the 49G." But what about non-integer and/or negative values of N ? Have you actually tried some ? Which ? If it cannot deal with values other than N=0 and N=1, it *would* be somewhat of a surprise, at least for me ... Please try it on an actual machine, if you haven't done already so, and tell me what values did you try and what results did it produce. Best regards from V. rolando Junior Member Posts: 1 Threads: 0 Joined: Jan 1970 10-17-2003, 10:11 AM Valentín: I have used an HP-49G emulator downloaded from hpcalc.org and I got the following answers: INTVX(SIN(X^N)) = INT(SIN(EXP(N*LN(X))) RISCH(SIN(X^N),X) = INT(SIN(EXP(N*LN(X))) I did no try the IBP (integration by parts) command because the function to be integrated must be expressed as a product of functions. The ROM in the emulator was 1.18 It seems that this is too much for the HP49G CAS. Rolando hugh Member Posts: 189 Threads: 23 Joined: Jan 1970 10-17-2003, 12:39 PM you're expecting a bit much arent you. these integratals are mostly only in terms of special functions (like complex gamma, sinc, and fresnel) which arent even there. except for the cases n = 0, 1 and 1/2. interstingly the case n = -1, is close to one of my numerical tests: integrate(cos(1/x), 0, 1) numerically or otherwise. so far no calculator has even got close (-0.0844109505) ▼ hugh Member Posts: 189 Threads: 23 Joined: Jan 1970 10-17-2003, 02:41 PM there's some more. 1/m (m = 1,2...) are also representable since int(sin(x^n), dx) can be solved when int(x^n sin(x), dx) is soluble. i'd be interested to know if these come out of the 48/49, using alg48, for example. Werner Huysegoms Junior Member Posts: 44 Threads: 0 Joined: Jan 1970 10-20-2003, 03:30 AM Hi Valentin. I did use a real 49G in the process: n>1 or negative: integral returned n=0,1 OK n=1/2 '2*SIN(SQRT(X))-2*SQRT(X)*COS(SQRT(X))' n=-1/2,n=1/3 integral returned That was in 'real mode'. In complex mode, nothing much changes. For negative n, the integral is split in two parts that are both left unintegrated, eg for n=-1 I get: 'INT(-(i*EXP(i/Xt)/2),Xt,X)+INT(i/(2*EXP(i/Xt)),Xt,X)' The Xt stuff is a change of variable and the 4x's way of saying it can't solve the integral. Cheers, Werner BTW you can play with a 49 yourself using the excellent Emu48 emulator and a 49G ROM, both of which can be found at www.hpcalc.org ▼ Valentin Albillo Posting Freak Posts: 1,755 Threads: 112 Joined: Jan 2005 10-20-2003, 04:33 AM Thanks Werner, it's more or less what I expected and what would be reasonable, by the way. The fact that it manages to solve the n=1/2 case is encouraging. I've always been somewhat amazed by the fact that sin(sqrt(x)) is elementarily integrable while sin(x^2) is not. One would tend to expect exactly the opposite, as the square root function seems more 'irrational', so to speak, than the simple, polynomial, 'square' function :-) Best regards from V.

 Possibly Related Threads... Thread Author Replies Views Last Post HP Prime: CAS taylor Alberto Candel 5 1,841 12-13-2013, 09:45 PM Last Post: Alberto Candel HP Prime CAS curiosity bluesun08 11 2,600 12-10-2013, 01:03 PM Last Post: Han [HP-Prime CAS] "Warning, ^ (Command) Is ambiguous on non square matrices"?? CompSystems 1 1,151 12-07-2013, 07:15 PM Last Post: CompSystems HP Prime: complex numbers in CAS. Alberto Candel 1 983 12-06-2013, 02:36 PM Last Post: parisse [HP Prime] plotfunc() bug in CAS Chris Pem10 2 1,208 12-04-2013, 02:46 PM Last Post: Chris Pem10 HP Prime: Home/CAS? Alasdair McAndrew 11 2,173 11-26-2013, 02:48 PM Last Post: Alberto Candel HP Prime... NOT meant to replace HP48,49,50 ? Chris Pem10 21 3,408 11-18-2013, 03:30 PM Last Post: Chris Smith HP Prime numerical precision in CAS and HOME Javier Goizueta 5 1,258 11-16-2013, 03:51 AM Last Post: Paul Dale Need CAS help Michael de Estrada 6 1,446 11-13-2013, 01:29 AM Last Post: cyrille de Brébisson HP Prime - CAS functions in Spreadsheet App CR Haeger 6 1,408 11-11-2013, 12:37 AM Last Post: Michael de Estrada

Forum Jump: