As with initial value problems, there are techniques other than shooting for solving two-point boundary problems. The most commonly used alternative technique is called relaxation. Relaxation methods divide the integration range into a 1-D grid of points. The ODE is represented by finite-difference equations that are solved at each point over the integration domain. The solution is iterated on until the required boundary conditions are met. We won't implement a relaxation method in this chapter, but if you wanted to you probably know how to do it by now. You would define a public , static method that would take an ODE object as one of its input arguments. The body of the method would then implement whatever relaxation technique was desired. |