ct=0;
for theta = 0:.01:5
	ct = ct+1;
	pf = 1./(sqrt(2*pi)*sqrt(.25));
	post(ct) = pf*exp((-1/2)*((1-theta)^2./sqrt(.25))) * pf*exp((-1/2)*((2-theta)^2./sqrt(.25))) * pf*exp((-1/2)*((3-theta)^2./sqrt(.25))) * pf*exp((-1/2)*((4-theta)^2./sqrt(.25)));
end

theta = 0:.01:5;
return
ct=0;
for t = 0:.01:15
	ct = ct+1;
	x = t:.01:20;
	tmp = pf1*exp((-1/2)*(x.^2));
	pc2g1(ct) = sum(tmp)./100;
	pc1g2(ct) = .5;
end

th = 0:.01:15;