samer@4: % scancols :; 
samer@4: %    ([[M]], [[N]] -> [[M]])   ~'scannning function',
samer@4: %    [[M]]                     ~'initial value',
samer@4: %    [[N,L]]                   ~'data to scan, sequence of length L'
samer@4: % -> [[M,L]].
samer@4: 
samer@12: function Y=scancols(f,y,X)
samer@4: 	Y=zeros(size(y,1),size(X,2));
samer@12: 	for i=1:size(X,2), y=f(y,X(:,i)); Y(:,i)=y; end
samer@12: end
samer@4: 		
samer@4: