clc
clear all
close all

fcn = @(x) log(2*x)
dfcn = @(x) 1/x

x_0 = 10;
niter = 50;

for i = 1:niter
     %% NEWTON
     x_1 = x_0 - 0.1*fcn(x_0)/dfcn(x_0)
     x_0 = x_1;
     plot(i,x_1,'*k');
     hold on;
end