samer@57: function z=plothist_lin(h,m)
samer@57: % plothist_lin - Plot histogram assuming triangular kernel (linear interpolation)
samer@57: %
samer@57: % plothist_lin :: [[N]->natural], dmap(N-1) -> handle.
samer@57: 
samer@57: 	e=edges(m);
samer@57: 	e=[2*e(1)-e(2),e,2*e(end)-e(end-1)];
samer@57: 
samer@57: 	z=ezplot(@pdf,[e(1),e(end)]);
samer@57: 
samer@57: 	function p=pdf(x)
samer@57: %		if ~isscalar(x), error('must be scalar'); end
samer@57: %		j=find(x>e(1:end-2) & x<e(3:end));
samer@57: 		p=0;
samer@57: 		for i=1:length(h) %j
samer@57: 			p=p+triangle(e(i:i+2),h(i));
samer@57: 		end
samer@57: 
samer@57: 		function y=triangle(a,h)
samer@57: 			y=zeros(size(x));
samer@57: 			I1=(x>a(1) & x<=a(2));
samer@57: 			I2=(x>a(2) & x<a(3));
samer@57: 			y(I1)=h*((x(I1)-a(1))/(a(2)-a(1)));
samer@57: 			y(I2)=h*((x(I2)-a(3))/(a(2)-a(3)));
samer@57: 		end
samer@57: 	end
samer@57: end
samer@57: