samer@3: % integrate - integrate sequence
samer@3: %
samer@3: % integrate ::
samer@3: %    seq [[N1,N2]]    ~'sequence to integrate', 
samer@3: %    D:natural        ~'dimension to integrate along',
samer@3: %    [[M1,M2]]        ~'initial value',
samer@3: %    ( [[M1,M2]] -> [[M1,M2]]) ~'function to filter initial value each block'
samer@3: % -> seq [[N,M]].
samer@3: %
samer@3: % cumulative sum taken along Dth dimension of the array
samer@3: 
samer@3: function y=integrate(x,n,a,ff)
samer@3: 	if nargin<4, ff=@(t)t; end
samer@23: 	y=cumsum(map(ff,x),a,n);
samer@3: end