%Newton.m takes a hardwired function and a hardwired starting point and
%tries to find the zero around that starting point

f=inline('816x*^3-3835*x^2+6000*x-3125');
fprime=inline('2448*x^2-7670*x^+6000');

xprev=1;