a function rfNewtonR that implements the Newton-Raphson method. Your function should accept the following inputs

x = computed root estimate

function x = rfNewtonR(fun,x,TolX,itmax)

% Find the root of f(x) = 0 using the Newton-Raphson method.

% Start the iteration.