Numerical routines in maxima



Hi - I still haven't figured out myself how to plot discrete data points
with Maxima, which I think is what you're trying to do.  The problem lies
with the plotting, not Newton.  Sorry...         Dave Holmgren

> -----Original Message-----
> From:	Camm Maguire [SMTP:camm at enhanced]
> Sent:	Tuesday, July 17, 2001 2:03 PM
> To:	David Holmgren
> Cc:	'maxima@www.ma.utexas.edu'
> Subject:	Re: [Maxima] Numerical routines in maxima
> 
> Greetings, and thanks so much for your helpful reply!  Now do you know
> how to get something like this to work?
> 
> z(c):=NEWTON(%E^d-d-c,d,1.0,0.01);
> plot2d(z(c),[c,1.2,2.3]);
> 
> 
> I get the following error:
> 
> MACSYMA was unable to evaluate the predicate:
> ABS(c - 1.718281828459045) < 0.01
> #0: NEWTON(exp=%E^d-d-c,var=d,x0=1.0,eps=0.01)
> #1: zz(c=c)
> #2: z(c=1.5)
>  -- an error.  Quitting.  To debug this try DEBUGMODE(TRUE);)
> (C242) 
> 
> David Holmgren <DavidHolmgren@smarttech.com> writes:
> 
> > Hi Camm,
> > 
> >  Certainly!  At the beginning of a session, type numer:true - this will
> let
> > you do nearly everything in floating point.  There is a routine in the
> share
> > directory for Newton's method for your nonlinear equation.  I've done a
> few
> > numerical routines by essentially transcribing published Fortran or C
> codes.
> 
> I've a little library built up myself.  Have you interfaced into
> maxima?  If so, how?
> 
> Many thanks!
> 
> > 
> >  Dave Holmgren
> > 
> > Dr. David E. Holmgren,
> > Imaging Scientitst,
> > SMART Technologies, Inc.
> > Calgary, AB, Canada
> > http://www.smarttech.com
> > [403]-235-1452, ext. 251
> > DavidHolmgren@smarttech.com
> > 
> > > -----Original Message-----
> > > From:	Camm Maguire [SMTP:camm at enhanced]
> > > Sent:	Monday, July 16, 2001 6:48 PM
> > > To:	maxima@www.ma.utexas.edu
> > > Subject:	[Maxima] Numerical routines in maxima
> > > 
> > > 
> > > Greetings!  I've been looking at a number of the numerical routines
> > > listed by the info pages as accessible under maxima, but I can't get
> > > any of them to work.  Like zsolve, for example.  Are these routines
> > > still available?
> > > 
> > > If not, is there any other convenient way to numerically, rather than
> > > symbolically, process a given stage of a maxima session?  Can one for
> > > example approximate the root of d=%E^d-1.5?
> > > 
> > > Is there any interest in providing some numerical functions, say from
> > > the gsl library, under maxima?  If so, is there a howto or simple
> > > example? I could probably add what I need.
> > > 
> > > Take care,
> > > 
> > > -- 
> > > Camm Maguire
> camm@enhanced.com
> > >
> ==========================================================================
> > > "The earth is but one country, and mankind its citizens."  --
> Baha'u'llah
> > > _______________________________________________
> > > Maxima mailing list
> > > Maxima@www.math.utexas.edu
> > > http://www.math.utexas.edu/mailman/listinfo/maxima
> > 
> > 
> 
> -- 
> Camm Maguire			     			camm@enhanced.com
> ==========================================================================
> "The earth is but one country, and mankind its citizens."  --  Baha'u'llah