Re: Numerical routines in maxima



Greetings, and thanks so much for this example.  It works!  Maxima is
great ...

Take care,

Virgil <virgilinux@yahoo.com> writes:

> --- Camm Maguire <camm@enhanced.com> wrote:
> > 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) 
> 
> z(c):=NEWTON(%E^d-d-c,d,1.0,0.01) fails because, it seems to me, you
> cannot define a "generic" function in terms of the newton function
> applied with an unknown parameter.
> 
> I believe I did something similar to what you want by using:
> i) create_list
> ii) openplot_curves
> 
> It would go something like:
> 
> load("newton");
> cpts:create_list(j/10,j,12,23);
> pts:create_list([c,NEWTON(%E^d-d-c,d,1.0,0.001)],c,cpts);
> openplot_curves([pts]);
> 
> See also http://www.ma.utexas.edu/users/wfs/netmath/pmaxima.html
> 
> Good luck!
> 
> __Virgil
> 
> 
> 
> 
> __________________________________________________
> Do You Yahoo!?
> Get personalized email addresses from Yahoo! Mail
> http://personal.mail.yahoo.com/
> 
> 

-- 
Camm Maguire			     			camm@enhanced.com
==========================================================================
"The earth is but one country, and mankind its citizens."  --  Baha'u'llah