function [xx] = readakernel_dat(fname)

fid = fopen(fname,'rb');
NN = fscanf(fid,'%d',2);
xx = fscanf(fid,'%f', [NN(1), NN(2)]);
fclose(fid);
image(xx);
h=title(fname);
set(h,'Interpreter','none');
colormap(gray(256));
axis('image');

return

