5a.R is R code which generates noisy data at 100 points on [0,1] from a given curve with specified noise standard deviation (ss), replicated 10 times, generates a cubic spline curve fit with the gcv estimate of lambda, generates the "bayesian confidence intervals", generates a table result.txt with the smoothing parameter, the gcvscore, the ratio of the estimated to the true standard deviation, counts the number of "true" data points inside the confidence interval, and averages the counts.