% 3 point differentiator

function cc=drpdiff(cc,sf)

cof=[-1 0 1];               
cof=cof.*(sf/2);             
temp=cc'*cof;                
temp=temp';
aantal=length(temp);

temp2=[temp(1,1:aantal-2);temp(2,2:aantal-1);temp(3,3:aantal)];   
e=sum(temp2);                
cc=[e(1) e e(aantal-2)];