Sampling the Banana-shaped distribution

[junpenglao]

What about below:

with pm.Model() as m:
    mu = pm.Normal('mu', 0., 10.)
    x = pm.Normal('x', mu, 1., observed=np.array([2., 1.4]))
    trace = pm.sample()