Calculus.df

class Calculus
   function Derivative(channel)
      private derivatives
      
      for (Private.Counter = 1, Counter < NumRows(channel), Counter++)
         private currentTime = GetTime(channel[Counter])
         private previousTime = GetTime(channel[Counter - 1])
         
         derivatives[Counter - 1] = (channel[Counter] - channel[Counter - 1]) / (currentTime - previousTime)
         
         derivatives[Counter - 1].Time = (currentTime + previousTime) / 2 // Set derivative time to midpoint of sampled data.
      endfor
      
      return derivatives
   endfunction
   
   function RiemannTrapezoidal(channel)
      private integral = 0
      
      for (Private.Counter = 1, Counter < NumRows(channel), Counter++)
         integral += (channel[Counter] + channel[Counter - 1]) * Abs(GetTime(channel[Counter]) - GetTime(channel[Counter - 1])) / 2
      endfor
      
      return integral
   endfunction
endclass

global Calculus = new(Calculus)