wolffd@0: function plot_axis_thru_origin()
wolffd@0: 
wolffd@0: lnx=line(get(gca,'xlim'),[0 0]); lny=line([0 0],get(gca,'ylim'));