Currently the y_err value just adds to the y value and subtracts to make a shadow region either side. However it would be better if this was a multiplication as the y_err is actually the standard deviation. Then another argument for sigma can be added to allow users to input 1 sigma but plot 2 , 3 or more sigma error