title: Show confidence limits and prediction limits in scatter plot
body: <p>I have two arrays of data for height and weight:</p>
  <pre><code>import numpy as np, matplotlib.pyplot as plt
  
  heights = np.array([50,52,53,54,58,60,62,64,66,67,68,70,72,74,76,55,50,45,65])
  weights = np.array([25,50,55,75,80,85,50,65,85,55,45,45,50,75,95,65,50,40,45])
  
  plt.plot(heights,weights,'bo')
  plt.show()
  </code></pre>
  <p>How can I produce a plot similar to the following?</p>
  <p><a href="http://www.sas.com/en_us/software/analytics/stat.html#m=screenshot6" rel="nofollow noreferrer">http://www.sas.com/en_us/software/analytics/stat.html#m=screenshot6</a></p>
  <p><a href="https://i.stack.imgur.com/z7SP1.jpg" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/z7SP1.jpg" alt="enter image description here" /></a></p>