h=(10^-1).^[1:8];
x=pi;
r=5e-6.*rand(1,16);
fminush1=cos(x-h)+r(1:8);
fplush1=cos(x+h)+r(9:16)
fder=(fplush-fminush)./(2*h);
plot(log(h),log(abs(-sin(x)-fder)),'r')
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!