The Museum of HP Calculators

HP Forum Archive 21

[ Return to Index | Top of Index ]

Hard or wanted CAS stuff
Message #1 Posted by Oliver Unter Ecker on 15 Apr 2012, 7:19 a.m.

I'm testing my CAS and am looking for some "hard" problems.

Have you run into CAS use cases that were just a little too much for your 50g (or your TI CAS whatever calc--for those of us who crossed into enemy territory, lately ;-))?

Something like symbolically integrating

'sin(ln(x))+cos(ln(x))' 
which will not yield a result on the 50g.

Apart from those, what are the top three things you always wished the 50g had (in the CAS department)?

Thank you!

      
Re: Hard or wanted CAS stuff
Message #2 Posted by peacecalc on 15 Apr 2012, 8:23 a.m.,
in response to message #1 by Oliver Unter Ecker

Hello Oliver,

I think this example isn't a very good choice, because with a simple substitution the integral can be performed. It's a lack, that the CAS (of the HP 50g) that doesn't recognizes this. But if you are working in the STEP BY STEP mode you get the information "Invalid unary operator" which indicates a necessary substitution, this and the following integration can easily performed with this handy CAS.

Sincerely peacecalc

            
Re: Hard or wanted CAS stuff
Message #3 Posted by Oliver Unter Ecker on 15 Apr 2012, 8:41 a.m.,
in response to message #2 by peacecalc

Hi peacecalc,

Thanks for your note. I had not noticed that this integral can be done in step-by-step mode, nor the simple substitution. This may well be not a good example. It was just meant to show a function on which the CAS chokes. I'm looking for cases like that (and even more so cases that don't work at all), to see how my other CAS performs.

My biggest curiosity is about what's felt to be missing in the 50g CAS (so I can look into patching those holes).

Cheers.

Edited: 15 Apr 2012, 8:43 a.m.

                  
Re: Hard or wanted CAS stuff
Message #4 Posted by peacecalc on 15 Apr 2012, 9:53 a.m.,
in response to message #3 by Oliver Unter Ecker

Hello Oliver,

my problem with your question is, that I'm not as good in algebra like the HP 50g ;-P. That means, when the calc isn't able to perform for instance an symbolic integral you have more possibilities:

a) there doesn't exist a closed form of the integral b) it is a problem which can be solved with a little changement (substituion or partial fraction and so on) c) there exists a solution, but in form of a sum with an infinite number of terms (like the solution of the mathematic pendicular without the approximation for small angles). d) there exists a closed solution but the CAS doesn't find it. and so on...

My solution for this is: IF the CAS failed THEN I keyed in the wrong input. Because there exist a) (see above), it is not always the right decision, but as a working hypothesis it leads you correctly (I'm hoping) through the mathematical fog.

Sincerely peacecalc

      
Re: Hard or wanted CAS stuff
Message #5 Posted by Marcus von Cube, Germany on 15 Apr 2012, 10:47 a.m.,
in response to message #1 by Oliver Unter Ecker

I tried on my early version TI-92 which does not find a solution for the integral. Derive 6 on the PC has no problems returning a solution. I doubt that Derive is really the ancestor of the TI-92 CAS because there are way to many differences.

The Nspire CAS PC software (1.3) does compute the integral symbolically. Time to test my more recent TI-92plus or Voyage 200.

Edit: My TI-89 with OS 2.09 can handle the integral. As a general note: AFAIK TI-CAS cannot solve step by step so you cannot see what's going wrong if it cannot find a solution.

Edited: 15 Apr 2012, 11:05 a.m.

      
Re: Hard or wanted CAS stuff
Message #6 Posted by Gerson W. Barbosa on 15 Apr 2012, 1:33 p.m.,
in response to message #1 by Oliver Unter Ecker

What about symbolically integrating Kahan's integral? I don't know whether TI CAS can do it. W|A can:

integrate (sqrt(x)/(x - 1) - 1/ln(x)) from x = 0 to x = 1

Edited: 15 Apr 2012, 1:38 p.m.

            
Re: Hard or wanted CAS stuff
Message #7 Posted by George Litauszky on 15 Apr 2012, 2:33 p.m.,
in response to message #6 by Gerson W. Barbosa

Hi Gerson,

Quote:
I don't know whether TI CAS can do it.

My TI 92 Plus can.

     /1
     |
 I = | (sqrt(x)/(x-1) - 1/ln(x))dx = .036489973974, time= 31 sec
     /0

George
                  
Re: Hard or wanted CAS stuff
Message #8 Posted by Gerson W. Barbosa on 15 Apr 2012, 4:09 p.m.,
in response to message #7 by George Litauszky

Hi George,

Quote:
My TI 92 Plus can.

It cannot find a closed form solution however. By the way, the last significant digits in the TI-92 Plus answer are wrong (39 instead of 85).

TI-92 Plus: .0364899739739
 HP-15C LE: .03648997377    
    Actual: .03648997397857652.. (2 - ln(4) - EulerGamma)

Cheers,

Gerson.

            
Re: Hard or wanted CAS stuff
Message #9 Posted by Oliver Unter Ecker on 15 Apr 2012, 4:09 p.m.,
in response to message #6 by Gerson W. Barbosa

That's a good one, Gerson. It works, and leaves me with li() in the result expression, which I know I have to implement now. Thank you for that!

      
Re: Hard or wanted CAS stuff
Message #10 Posted by Mike (Stgt) on 16 Apr 2012, 4:32 a.m.,
in response to message #1 by Oliver Unter Ecker

I did not try it yet with CAS' help, but even with a pencil and a pile of paper you should know the way how to solve it: the cumulative distribution function (integrating the Gaussian distribution) in terms of elementary functions. According to this wiki it is not possible, but IIRC there is simple trick.

At least, CAS could help to go the 'long' way: integrate the Taylor's expansion of the Gaussian distribution. I did it once w/o CAS to make an HP41 program for the cumulative distribution function (like Q on the HP32E).

Ciao.....Mike


[ Return to Index | Top of Index ]

Go back to the main exhibit hall