ctPoly function

Plots uncertainty bands with shading