Finding the solution of a linear boundary problem is assisted by the linear structure of the equation and the use of two special initial value problems. Suppose that u(t) is the unique solution to the I.V.P. with . Furthermore, suppose that v(t) is the unique solution to the I.V.P. with . Then the linear combination . is a solution to with . Program (Linear Shooting Method) : To approximate the solution of the boundary value problem with over the interval [a,b] by using the Runge-Kutta method of order n=4. The method involves solving a two systems of equations over . First solve with , and . Then solve with , and . Finally, the desired solution x(t) is the linear combination . The subroutine Runge2D will be used to construct the two solutions , and .
Introductory Methods of Numerical Analysis