for k in range(class_count):
        norm = scipy.stats.norm(class_mean[k], class_std[k])
        class_gaussian[:, k] = class_prob[k] * norm.pdf(level)