mixt_model

rand('seed',sed)
[l,c]=size(m)
Pacc=P(1);
for i=1:c
     t=Pacc(i-1)+P(i);
     Pacc=[Pacc t];
end
X=[];
y=[];
for i=1:N
     t=rand;
     ind=sum(t>Pacc)+1
    X=[X; mvnrd(m(:,ind)',S(:,:,ind),1)];
    y=[y ind];
end
X=X'