function HW5_1(n)
 
f = @(x) exp(x) - x - 1;
df = @(x) exp(x) - 1;
ddf = @(x) exp(x);
%g = @(x) x - f(x) / df(x); 
g = @(x) x - f(x) / df(x) - ddf(x) * ( f(x)^2 ) / ( 2 * ( df(x)^3 ) );
P_0 = 1; P = 0;
for i = 1:n
	P_1 = g(P_0);
	fprintf('n = %3.0f %16.8f\n', i, abs(P_1-P) / abs(P_0-P) );
	P_0 = P_1;
end
 
 end