The current analysis deals with radiative aspects of magnetohydrodynamic boundary layer flow with heat mass transfer features on electrically conductive Williamson nanofluid by a stretching surface. The impact of variable thickness and thermal conductivity characteristics in view of melting heat flow are examined. The mathematical formulation of Williamson nanofluid flow is based on boundary layer theory pioneered by Prandtl. The boundary layer nanofluid flow idea yields a constitutive flow laws of partial differential equations (PDEs) are made dimensionless and then reduce to ordinary nonlinear differential equations (ODEs) versus transformation technique. A built-in numerical algorithm bvp4c in Mathematica software is employed for nonlinear systems computation. Considerable features of dimensionless parameters are reviewed via graphical description. A comparison with another homotopic approach (HAM) as a limiting case and an excellent agreement perceived.