Numerical representation

[Robert Kern]

On 3/4/11 5:49 PM, Santoso Wijaya wrote:
> Have you taken a look at numpy? [1] It was written for exactly this kind of usage.
>
> ~/santa
>
> [1] http://numpy.scipy.org/

While numpy does provide arrays much like MATLAB's, it won't help at all for the 
precision issues the OP is encountering (and honestly doesn't help much when 
implementing an RK integrator, except as input and output data structures). 
Python, numpy, and MATLAB all use double-precision floating point numbers by 
default.

-- 
Robert Kern

"I have come to believe that the whole world is an enigma, a harmless enigma
  that is made terrible by our own mad attempt to interpret it as though it had
  an underlying truth."
   -- Umberto Eco