Two collocation-based methods utilizing the novel Bessel polynomials (with positive coefficients) are developed for solving the non-linear Troesch's problem. In the first approach, by expressing the unknown solution and its second derivative in terms of the Bessel matrix form along with some collocation points, the governing equation transforms into a non-linear algebraic matrix equation. In the second approach, the technique of quasi-linearization is first employed to linearize the model problem and, then, the first collocation method is applied to the sequence of linearized equations iteratively. In the latter approach, we require to solve a linear algebraic matrix equation in each iteration. Moreover, the error analysis of the Bessel series solution is established. In the end, numerical simulations and computational results are provided to illustrate the utility and applicability of the presented collocation approaches. Numerical comparisons with some existing available methods are performed to validate our results.