#!/usr/bin/env python

#This program generates a Gaussian random no with 0 mean and unit variance

#Importing numy, scipy, mpmath and pyplot
import numpy as np
import mpmath as mp
import scipy 
import scipy.stats as sp
import matplotlib.pyplot as plt
import subprocess



simlen = 1e5 #No of samples

n = np.random.normal(0,1,simlen)#Random vector

mean = np.sum(n)/simlen #Mean value

print mean

var = np.sum(np.square(n - mean*np.ones((1,simlen))))/simlen

print var
