function Q1(n)

f = @(x) 2 * x / ( (1+x)^(2/3) + (1-x^2)^(1/3) + (1-x)^(2/3) ) - 10^(-10);
df = @(x) (1/3) * (1+x)^(-2/3) + (1/3) * (1-x)^(-2/3);
g = @(x) x - f(x) / df(x);

P0 = 1.5 * 10^(-10);

for i = 1:n
	P = g(P0);
	vpa(g(P0), 30)
	P0 = P;
end

end