  function P3()
% Checks out MyCubeRoot...

clc
N = 100000;
V = linspace(1/8,1,N);
z = zeros(1,N);
tic
for i=1:N
    z(i) = MyCubeRoot(V(i));
end
t1 = toc;
nBadError = length(find(abs(z-V.^(1/3))>(10^-14)*V.^(1/3)));
fprintf('Time for 100000 cube roots = %5.3f\n\n',t1)
fprintf('Number insufficiently accurate cube roots = %1d\n\n',nBadError)

