for i=1:rows(para.labels)
      pdf(:, i) = mvnpdf(pts, para.mu(i, :), para.sig(:, :, para.labels(i)));
    endfor