--- 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/