Daniel@0: d = 2; M = 3; Q = 4; T = 5; Sigma = 10;
Daniel@0: N = sample_discrete(normalize(ones(1,M)), 1, Q);
Daniel@0: data = randn(d,T);
Daniel@0: mu = randn(d,M,Q);
Daniel@0: 
Daniel@0: [BM, B2M] = parzen(data, mu, Sigma, N);
Daniel@0: [B, B2] = parzenC(data, mu, Sigma, N);
Daniel@0: 
Daniel@0: approxeq(B,BM)
Daniel@0: approxeq(B2,B2M)