1 2 3 4 5 6 7 8 9 # Converted from onednormal.m - standard normal PDF import numpy as np from misc import save_ascii npts = 100 x = np.linspace(-4.0, 4.0, npts) y = np.exp(-0.5*x**2) / np.sqrt(2*np.pi) table = np.column_stack([x, y]) save_ascii('onednormal.dat', table)