f=@(x) sqrt(1-x.^2);
x0=sqrt(2)/2

difc=@(x,h) (f(x+h)-f(x-h))/(2*h)

for i=1:15
derivadaap=difc(x0,10^(-i));
error=(-1)-derivadaap;
fprintf('Para h=10^(-%2d), el error es %d.\n',i,error)
end