clear all
X = [95,118, 85,154,102,96]';
mu = 100.0;
sigma = 10.0;
L = prod(normpdf(X,mu,sigma))


mu = [85:.3:154];
sigma = [0:.5:70];
L = zeros(length(mu),length(sigma));
for i = 1:length(mu)
    for j = 1:length(sigma)
        L(i,j) = prod(normpdf(X,mu(i),sigma(j)));
    end
end

surf(mu,sigma, L')
xlabel('\sigma')
ylabel('\mu')