I am attempting a plot where:

x axis: x_1

and y axis: f(x_1)

where x_k is determined by:

x_1 = x_0 - f(x_0)/f'(x_0) (Newton's Method)

So far my program is:

NEWTPATH(t,n)

BEGIN

LOCAL k,w,y;

RECT();

FOR k FROM 1 TO n DO

t-(F1(t)/der(F1(X),X=t))- > w;

F1(w)- > y;

PIXON(w,y);

w- > t;

END;

FREEZE;

END;

The "der" is the derivative funtion.

The program either freezes at the input screen or gives me an "out of bounds" error. I have F1(X) already stored before running the program. Any ideas? Thanks!

Eddie

*Edited: 4 Nov 2012, 11:43 a.m. *