[Maxima-bugs] [ maxima-Bugs-1588623 ] romberg with % or %o argument

[SourceForge.net <no...@so...>]

Bugs item #1588623, was opened at 2006-11-01 08:42
Message generated for change (Comment added) made by willisbl
You can respond by visiting: 
https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1588623&group_id=4933

Please note that this message will contain a full copy of the comment thread,
including the initial issue submission, for this request,
not just the latest update.
Category: Lisp Core - Integration
Group: None
Status: Open
Resolution: None
Priority: 5
Private: No
Submitted By: Barton Willis (willisbl)
Assigned to: Nobody/Anonymous (nobody)
Summary: romberg with % or %o argument

Initial Comment:
(%i1) sqrt(1+x^3);
(%o1) sqrt(x^3+1)
(%i2) romberg(%, x,0,1);
Maxima encountered a Lisp error:

Also, the command romberg(%o1,x,0,1) gives a Lisp
error. This works OK:

(%i4) romberg(sqrt(1+x^3), x,0,1);
(%o4) 1.111447999508385


(%i1) build_info();
Maxima version: 5.10.0
Maxima build date: 19:9 9/21/2006
host type: i686-pc-mingw32
lisp-implementation-type: GNU Common Lisp (GCL)
lisp-implementation-version: GCL 2.6.7

----------------------------------------------------------------------

>Comment By: Barton Willis (willisbl)
Date: 2006-11-08 06:00

Message:
Logged In: YES 
user_id=895922

I agree. It would be good to have a function for big float
numerical integration.
Not that romberg works with big floats:

(%i2) romberg(sin(x^2),x,0,1);
(%o2) 0.31026884083166872
(%i3) rombergtol : 1.0b-12;
(%o3) 1.0b-12
(%i4) romberg(sin(x^2),x,0,1);
Maxima encountered a Lisp error:
 Error in MACSYMA-TOP-LEVEL [or a callee]: ((BIGFLOAT SIMP
56) ...




----------------------------------------------------------------------

Comment By: Robert Dodier (robert_dodier)
Date: 2006-11-05 10:14

Message:
Logged In: YES 
user_id=501686

I don't recommend trying to fix romberg. The quadpack
functions are more extensive implementations of the same
basic idea (extrapolation from simple quadrature ruls);
quadpack subsumes romberg and goes much farther. I'd like to
cut out romberg entirely.

----------------------------------------------------------------------

You can respond by visiting: 
https://sourceforge.net/tracker/?func=detail&atid=104933&aid=1588623&group_id=4933