function [xi] = newton(x0, error, baller, ballerprime) 

format long;

j = 100000;

while (error < abs(j))
    xi = x0 - (baller(x0)/ballerprime(x0));
    j = abs(xi - x0);
    x0 = xi;
    %disp(xi);
end