A Numerical Method for Singular Two Point Boundary Value Problems