% method of class @signal
% 
%   INPUT VALUES:
%  
%   RETURN VALUE:
%
% 
% (c) 2003-2008, University of Cambridge, Medical Research Council 
% Maintained by Tom Walters (tcw24@cam.ac.uk), written by Stefan Bleeck (stefan@bleeck.de)
% http://www.pdn.cam.ac.uk/cnbh/aim2006
% $Date: 2008-06-10 18:00:16 +0100 (Tue, 10 Jun 2008) $
% $Revision: 585 $

function sig=secondderivation(sig)
% calculates the second derivation by calculation of the differences


val=getvalues(sig);

nr=length(val);


nval=val;
nval(1)=0;
nval(2)=0;
for i=3:nr

    if i==370
        a=0;
    end
    nval(i)=(val(i-2)-2*val(i-1)+val(i))/2;
end

sig=setvalues(sig,nval);
