I get a situation where you need to keep track of two numbers for every operation. Maybe I should be doing my power series as a vector (or some kind of object) sum to do both at the same time with overridden +-*/^ operators. This is getting complicated and messy since I don't know how to create a new type of Maxima object for this.
------------Original Message------------
From: Barton Willis <willisb at unk.edu>
To: "Barton Willis" <willisb at unk.edu>
Cc: "Maxima List" <maxima at math.utexas.edu>, maxima-bounces at math.utexas.edu
Date: Tue, May-6-2008 3:45 PM
Subject: Re: [Maxima] Big float Round off errors
A Maxima version:
(%i1) bigfloatepsilon () := bfloat(2^-?fpprec);
Barton
-----maxima-bounces at math.utexas.edu wrote: -----
>(%i2) :lisp(defun $bigfloatepsilon () ($bfloat (div 1 (expt 2 fpprec)))));
>$BIGFLOATEPSILON
>
_______________________________________________
Maxima mailing list
Maxima at math.utexas.edu
http://www.math.utexas.edu/mailman/listinfo/maxima