Quadratic non-polynomial spline approach to the solution of a system of second-order boundary-value problems

A quadratic non-polynomial spline functions based method is developed to find approximations solution to a system of second-order boundary-value problems associated with obstacle, unilateral, and contact problems. The present approach has less computational cost and gives better approximations than those produced by other collocation, finite-difference and spline methods. Convergence analysis of the method is discussed. A numerical example is given to illustrate practical usefulness of the new method.