function pval = rasp ( x )
n = length ( x );
pval = mod ( [0:n-1], 4 );
