 F = @(x)(sin(5*x)+20)/10-2;
 
 root = newton_iteration(F, 1.5, 0.001)
 y = F(root)
 
 X = linspace(-pi, pi, 100);
 Y = F(X);
 plot(X,Y)
 hold on
 plot(root,y,'o')