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