--- /dev/null
+function[phiInit,rhoInit,piInit,gamInit] = basicInitParameters(n,p,m,k)
+
+ phiInit = zeros(p,m,k);
+
+ piInit = (1.0/k) * ones(1,k);
+
+ rhoInit = zeros(m,m,k);
+ for r=1:k
+ rhoInit(:,:,r) = eye(m,m);
+ end
+
+ gamInit = 0.1 * ones(n,k);
+ R = random('unid',k,n,1);
+ for i=1:n
+ gamInit(i,R(i)) = 0.9;
+ end
+ gamInit = gamInit / (sum(gamInit(1,:)));
+
+end