function [y]=std_cols(x)
% standardise columns of x
%
 [m,n]=size(x);
 y=x-repmat(mean(x),m,1);
 y=y./repmat(std(y),m,1);
%