The aim of a present article is to investigate the laminar unsteady two‐dimensional boundary layer flow of a nanofluid with Stefan blowing and slip effect. First, governing boundary layer equations are converted in the ordinary form of the differential equations (ODEs) from partial differential equations using appropriate coordinate transformations. The obtained ODEs are then solved by applying a shooting method with the Runge‐Kutta fourth order method by implementation of the Maple software. The influences of different controlling dimensionless parameters over the dimensionless velocity, temperature, concentration, friction factor, local heat as well as mass transfer have been discussed and represented by plots. It is found that there exist dual solutions for the different applied nanofluid parameters along with the blowing parameter. The results reveal that by increasing the values of the Brownian motion (Nb), thermophoresis (Nt) and blowing parameters (fw), the skin friction increases (decreases) in the first (second) solution.