#!/usr/bin/env python3
# -*- coding: utf-8 -*-

"""
Scatter multivariate normal (version 2)
"""

import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import axes3d

# Build datas ###############

num_points = 1000
mu = np.array([0., 0., 0.])
cov = np.array([[1., 0., 0.],
                [0., 1., 0.],
                [0., 0., 1.]])

x = np.random.multivariate_normal(mu, cov, num_points)

# Plot data #################

fig = plt.figure()
ax = axes3d.Axes3D(fig)
ax.scatter(x[:,0], x[:,1], x[:,2], color='b')

# SAVE FILES ######################
plt.savefig("demo4b_mplot3d_scatter_multivariate_normal.png")

plt.show()

