title: truncated normal distribution with scipy in python
body: <p>I am trying to use a truncated normal distribution with <code>scipy</code> in Python3. I want to do something simple: plot the pdf of a truncated normal centered at 0.5 and ranging from 0 to 1. I have the following code line</p>
  
  <pre><code>from scipy import truncnorm
  import matplotlib.pyplot as plt
  plt.plot([truncnorm.pdf(p,0,1, loc=0.5) for p in np.arange(0,1.1,0.1)])
  </code></pre>
  
  <p>However, this does not give the nice bell-shaped probability distribution function I would expect. Rather, it equals 0 before 0.5, and I cannot figure out why. Any advice on this?</p>