## 1 Introduction  
Flow through non-parallel walls is an important area of research due to its many practical, industrial, physical and biological applications, such as flows through rivers and canals. Blood also exhibits such a flow, where the arteries and capillaries are linked with each other. Pioneering work in this area was done by [PERSON] and [PERSON] [1,2]. Plenty of studies are available in this area inspecting the various properties of flow by considering various external forces, such as external magnetic fields and the other forces. One can easily find enough literature on this topic in [3,4] and the references therein. The study of stretching or shrinking sheets has been an area of modern research due to the practical applications in the real world. These applications include glass fiber production, engineering melt spinning, the manufacturing of rubber sheets, _etc._ The work presented by [PERSON] [5] was the first study on stretching flows. Since then, several researchers from all over the world have used the said problem to formulate and present the various properties of flows over the stretching/shrinking surfaces that can be seen in [6, 7, 8] and some of the references therein. They used heat and mass transfer effects and analyzed the various properties of flows over stretching/shrinking surfaces. Traditional fluids, like water, kerosene, ethylene glycol, _etc._, in general, are not good conductors. To cope with this problem and to enhance the thermal properties of various fluids, [PERSON] [9,10] presented a model that used the term nanofluids. The idea is to add an appropriate amount of nanoparticles to the traditional fluids. The experimental results have proven that the thermal properties of traditional fluids can be enhanced appreciably by using the said technique. After this groundbreaking idea, many models were presented by various researchers that had also incorporated the Brownian motion and thermophoresis effects. [PERSON] [11], [PERSON] [12], [PERSON] and [PERSON] [13] and [PERSON] [14] are the major contributors who have presented some useful and efficient models for nanofluids. In recent times, many authors have been inclined to work in the field of nanofluids. Using [PERSON]'s model, [PERSON] _et al._[15] extended the work on nanofluids to a stretching surface. Since then, many studies have been presented using the models discussed earlier for different geometries and various types of nanoparticles. One can find enough literature on these topics [16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35]. The complex nature of the equations in daily life problems makes exact solutions least likely. Therefore, many analytical, as well as numerical techniques are developed to solve the problems representing the physical nature. Some of the analytical techniques are [PERSON]'s decomposition method (ADM), the homotopy perturbation method (HPM), the variational iteration method (VIM), the variation of parameters method (VPM), the homotopy analysis method (HAM), _etc._ In this article, VIM is successfully applied to the problem, and the efficiency and accuracy of VIM can be seen in [36,37] and the references therein. In a recent study, [PERSON] [38] extended the traditional Jeffery-Hamel flow problem to a case incorporating stretching/shrinking walls. He also studied the heat transfer effects on the flow. In the study under consideration, we examined a further extended form of the same problem dealing with the nanofluid flow under the influence of magnetohydrodynamic forces. Due to the very abstract nature of the governing equations and the expected non-availability of the exact solutions, an analytical approximation technique called the variational iteration method (VIM) has been used to approximate the solution. A comparison of analytical results is made with the numerical ones obtained by employing the Runge-Kutta-Fehlberg technique and is presented in the form of tables.