% erste Version 
clear all
close all

I   = [0 2*pi];
DIM = 40;

h = (I(2)-I(1))/(DIM-1);
x = I(1):h:I(2);

fh  = sin(x); % vorher GPS(x); Dazu kennen wir aber keine exakte Ableitung
dfh = cos(x); % zum testen des Ergebnisses

dfr(1)     = 0;
dfr(2:DIM) = (fh(2:DIM)-fh(1:DIM-1))/h;

dfv(1:DIM-1) = (fh(2:DIM)-fh(1:DIM-1))/h;
dfv(DIM)     = 0;

dfm = (dfr+dfv)/2;

plot(x,dfr,'bo',x,dfv,'ro',x,dfm,'mo',x,dfh,'k-');
grid on


