7
$\begingroup$

I'm starting to become more functional in Mathcad, so I wanted to take the opportunity to look into known drawbacks, mathematically speaking, of CAS/computer math systems. For instance, one of the exercises we did in school involved "breaking" TI Derive by asking for the limit of a discontinuous function (at the point of discontinuity, of course). Another example would be p. 13 of this;

http://people.virginia.edu/~jcb6t/Mathcad/Johns_Tutorial_on_Everyday_Mathcad.pdf

So basically, I'm hoping folks will contribute known problems with these kinds of apps, whether from their experience with a particular package or generalized problems inherent in computer algebra.

I myself am most concerned with Mathcad and Mathematica, but it would be awesome for people to talk about CAS/computer math in general and drawbacks of various systems (matlab et al) in particular.

  • 0
    I removed the tag computability.2012-02-26
  • 1
    I imagine this question is a good candidate for being community wiki!2012-02-26
  • 1
    IMO the scope of this question is too general. Limitations of numerical packages such as MATLAB are way different than limitations of symbolic packages such as Mathematica.2012-02-26
  • 0
    @harald - Why thank you, kind sir! *blush*2012-02-26
  • 0
    Here's a relevant collection of issues: http://www.math.niu.edu/~rusin/known-math/99/calc_errors It's a bit old and some of the bugs mentioned may have been fixed by now, but nevertheless it never hurts to be cautious when doing mathematics with technology.2012-02-26
  • 0
    Mathoverflow: http://mathoverflow.net/questions/11517/2012-04-30

2 Answers 2

8

There are many other problems: branch cuts, zero-recognition, keeping track of domains of validity, etc. A good place to find such information is to browse the web pages of leading researchers, and conference proceedings (ISSAC,SYMSAC,Sigsam,Eurosam, etc). For example, see Richard Fateman's papers, e.g. his 33 page critique of Mathematica, and Why Computer Algebra Systems Can't Solve Simple Equations and Branch Cuts in Computer Algebra, etc.

  • 1
    re Mathematica... In my distant past as a philosophy major, I had occasion to review (okay, review some reviews) of "A new kind of science", Wolfram's doorstopper book in which he basically claims to have discovered recursion and the math behind Penrose tiles (name escapes me atm). Nice to see a "professional" confirm my opinion of the dude. I'm using Mathematica basically because it's got a relatively well-behaved connection to Mathcad and I found a portable version on teh interwebz.2012-02-26
  • 3
    My personal impression? Wolfram has a lot in common with the late unlamented Mr. Jobs - his ability to speak marketese overshadows his deficiencies as a scientist/engineer.2012-02-26
3

I only tried this on Mathematica 4, so I don't know if the same problem happens in more modern versions. If you ask Mathematica about $$ \int_0^\infty (1-e^{-x})^2/x^2 dx, $$ it tells you that the integral does not converge. But if you split the integral, it gives you

$$ \int_0^1 (1-e^{-x})^2/x^2 dx\simeq 0.645751, $$

$$ \int_1^\infty (1-e^{-x})^2/x^2 dx\simeq 0.740543 $$

  • 0
    Mathematica 6.0 doesn't complain about divergence, but delivers another strange result: `Integrate[(1-Exp[-x])^2 / x^2, {x,0,Infinity}]` returns `(4 I) Pi + Log[4]`. The real part $\ln 4$ agrees with the numerical values, but I wonder where that nonzero imaginary part came from! Got to try this with a newer version at work tomorrow...2012-02-26
  • 0
    Alpha can do it http://www.wolframalpha.com/input/?i=int%20%5B%2F%2Fmath%3A((1-e%5E(-x))%5E2)%2Fx%5E2%20dx%2F%2F%5D%20from%20%5B%2F%2Fmath%3Ax%2F%2F%5D%3D%5B%2F%2Fmath%3A0%2F%2F%5D%20to%20%5B%2F%2Fmath%3Ainfinity%2F%2F%5D2012-02-26
  • 0
    Mathematica 8.0 returns `Log[4]`.2012-02-27