Convergence issue scaling square root

Ideally I want users to be able to give a value for c_0, apply the scaled BC, solve for theta, and scale again to show users c.