post=m_disc(m,prior,data)

function post=m_disc(m,prior,data)
%
% M_DISC Posterior distribution for a normal mean with discrete models.
%	POST = M_DISC(P,PRIOR,DATA) returns a vector of posterior probabilities.
%	M is the vector of values of the mean, PRIOR is the corresponding
%	vector of prior probabilities and DATA is the vector of data 
%	(sample mean, sample size, and population standard deviation).

xbar=data(1); n=data(2); sigma=data(3);

like=-.5*n/sigma^2*(m-xbar).^2;
like=exp(like-max(like));

product=like.*prior;
post=product/sum(product);