% programm numdiff
% tabelliert numerische Ableitungen von sin(x) bei x=1
% als Funktion der Schrittweite
%
help numdiff
%
h = 10.^(-[1:16]); % zu testende Schrittweiten
x = 1;
dex = cos(x); % exakt
d1 = (sin(x+h)-sin(x))./h;        % asym. Formel
d2 = (sin(x+h)-sin(x-h))./(h+h);  %  sym. Formel
y = [ h ; (d1-dex)/dex ; (d2-dex)/dex ];
%
fprintf('   h       asymm.     symm. \n\n')
fprintf(' %5.0e %10.1e %10.1e \n',y)