% zweite Version mit Matrizendarstellung
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); 
dfh = cos(x); % zum testen des Ergebnisses

Dm=DiffMat(1,h,DIM);
dfm = Dm*fh';

plot(x,dfm,'mo',x,dfh,'k-');
grid on

