title: Alternative for scipy.stats.norm.pdf?
body: <p>Does anyone know of an alternative for scipy.stats.norm.pdf()? I'm hosting my python site on Google App Engine and Google doesn't support SciPy.</p>
  
  <p>I've tried this function, but that didn't return the same results as scipy:</p>
  
  <pre><code>def normpdf(x, mu, sigma):
      u = (x-mu)/abs(sigma)
      y = (1/(sqrt(2*pi)*abs(sigma)))*exp(-u*u/2)
      return y
  </code></pre>
  
  <p>For example:</p>
  
  <pre><code>print scipy.stats.norm.pdf(20, 20, 10)
  print normpdf(20, 20, 10)
  
  print scipy.stats.norm.pdf(15, 20, 10)
  print normpdf(15, 20, 10)
  
  print scipy.stats.norm.pdf(10, 20, 10)
  print normpdf(10, 20, 10)
  </code></pre>
  
  <p>Returns these values:</p>
  
  <pre><code>0.0398942280401
  0.0398942280401
  
  0.0352065326764
  0.0146762663174
  
  0.0241970724519
  0.0146762663174
  </code></pre>