The novelty of this research is surveying the effects of CuO nanoparticles in jeffery-hamel flow with the high levelized of magnetic field in converging and diverging channels using quasi-linearization method (QLM) and also the meshless method which based on radial basis functions. In order to evaluate the convergence analysis of the method, error estimations are made by a residual function denoted. Also, the ability of the present method is shown by comparing it with the numerical method to solve this problem, which is in good agreement. Effects of multivariable parameters is analyzed on magnetic field, nanoparticles volume fraction and angle of converging and diverging channels. The obtained results show that at angles or Reynolds number of greater in divergent channels, backflow occurs so the high levelized of magnetic field is eliminated this phenomenon. Increasing viscous force is due to velocity reduction close to the walls of converging and diverging channels.