The present study concerns investigating the two-dimensional Magnetohydrodynamics (MHD) boundary layer flow of Williamson nanofluid over a non-linear stretching sheet. The focus of this study is based on the global influence of the non-Newtonian Williamson fluid parameter (λ) rather than the local one that exists in the literature for linear and non-linear stretching cases. The mathematical model of the problem is based on the law of conservation of mass, momentum, and energy. The derived partial differential equations are transformed into ordinary differential equations by applying an appropriate similarity transformation. The subsequent equations are solved numerically by using the Shooting method. The physical quantities Skin friction coefficient, as well as the Sherwood and Nusselt numbers are computed locally. To validate the implemented shooting method, a comparison is made with the results obtained by Matlab function bvp4c, and good agreement is found. The Prandtl number, Pr, has an increasing impact of 25.14% on the wall temperature gradient. The impact of various physical parameters are presented through graphs and tables.