Spline-based numerical treatments of Bratu-type equations