matthiasm@8: function addmax(varargin)
matthiasm@8: if length(varargin)>0 && strcmp(varargin,'rows')
matthiasm@8:     dimension = 2;
matthiasm@8: else
matthiasm@8:     dimension = 1;
matthiasm@8: end
matthiasm@8: 
matthiasm@8: chil=get(gca,'children');
matthiasm@8: z = get(chil(1),'CData');
matthiasm@8: x = get(chil(1),'XData');
matthiasm@8: [maxind, maxind] = max(z,[],dimension);
matthiasm@8: 
matthiasm@8: if length(x) == 2
matthiasm@8:     x = x(1):x(2);
matthiasm@8: end
matthiasm@8: 
matthiasm@8: hold on
matthiasm@8: m = plot(x, maxind, '.k','Erasemode','xor')
matthiasm@8: hold off