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.