where in the code I have squared the inside of the square root then applied another square root in order to avoid complex values since depending on how the solution converges the inside may be negative. The domain is given by $x \in [a,b]$. So, the free parameters here are $u_i(a)$ and $z_i(a)$ which we can input "guess" values to start the calculation.
