import numpy as np
import scipy.stats as stats
import matplotlib.pyplot as plt

SEED_POINTS = 5
JITTER_SAMPLES = 3

xs,ys = np.random.random((2,SEED_POINTS))
plt.scatter(xs, ys, c='b')

# create jittered data for x and y coords
xs_jit = stats.norm.rvs(loc=xs, scale=0.05, size=(JITTER_SAMPLES,SEED_POINTS))
ys_jit = stats.norm.rvs(loc=ys, scale=0.05, size=(JITTER_SAMPLES,SEED_POINTS))

plt.scatter(xs_jit, ys_jit, c='r')

plt.show()