wine_data;
c1_ind=find(wine_ins(:,1)==1);
wine_ins_c1=wine_ins(c1_ind,2:end);

clf;
subplot(4,1,1);
%hh=plot(wine_ins_c1(:,1),wine_ins_c1(:,1)*0,'.'); set(hh,'MarkerSize',[12]);
hold on;
%set(hh,'XLim',[10 20]);
sigma2 = 1;
subplot(4,1,2);
xx=5:0.01:20;
hold on;
cum_p=xx*0+1;
cum_p_log=xx*0;
for ii=1:size(wine_ins_c1,1),
   subplot(4,1,2);
   p1=exp(-(xx-wine_ins_c1(ii,1)).^2/(2*sigma2));
   p1=p1/sum(p1);
   plot(xx,p1);
   subplot(4,1,1);
   hh=plot(wine_ins_c1(ii,1),0,'.');  set(hh,'MarkerSize',[12]);set(gca,'XLim',[5 20]);
   subplot(4,1,3);
   %hold on;
   cum_p=cum_p.*p1;
   plot(xx,cum_p);
   subplot(4,1,4);
   %hold on;
   cum_p_log = cum_p_log + log(p1);
   plot(xx,cum_p_log);
   pause;
end



