Statistics for A new characteristic numerical approach with evolutionary residual error analysis to nonlinear boundary value problems occurring in heat and mass transfer via combinatoric Mittag-Leffler polynomial