There are a few problems that have to be addressed for the successful implementation of the algorithm described in the preceding section. The most important are the boundary values, the existence of a double solution (7-10), and the complications of finding the solution in the vicinity of the cusp points.