Calculating (and plotting) prediction interval from a function in r with ggplot2

Question

I would like to compare prediction intervals from a theoretical model to experimental data I have generated. Below is an example of the kind of data and function I want to v