function B = oshift(A,m)
% OSHIFT(A,m): Downward circular shift of the row indices
% of matrix A, by m positions.  Same as (P^m)*A, where P
% was defined in class.
[r c] = size(A);
B = A(1+mod((-m:-m+r-1),r),:);
return
