function s= zerofind(beginningpt,tolerance)
x=beginningpt;
diff=100;
while abs(diff) > tolerance)
    f=funct(x);
    fp=functp(x);
    y=x-(f/fp);
    diff=x-y;
    x=y;
end
s=x;