function pl(mat)
% scatter plot
plot(mat(1:end,1), mat(1:end,2),'.');