function [xx] = load_matrix(fname)

fid = fopen(fname,'rb');
NN = fscanf(fid,'%d',2);
xx = fscanf(fid,'%f', [NN(2), NN(1)]);
xx = xx';
fclose(fid);
return

