“…Due to the flexibility and ability, the VIM has emerged as a promising tool to solve non-linear or complex systems. Consequently, this method was successfully implimented to find the approximate solution of Hilfer fractional advection diffusion equations with power law initial condition [2], of reaction diffusion system with fast reversible reaction [1], of wave-like and heat-like equations in large domains [9], of delay differential equations [10], of nonlinear singular boundary value problems [16], of the time fractional Fornberg-Whitham equation [25], of Burger's flow with fractional derivatives [28], of non-isothermal reaction diffusion model equations in a spherical catalyst [26], of Optical bright and dark soliton solutions for coupled nonlinear Schrdinger (CNLS) equations [27], to Bratu-like equation arising in electrospinningand [12], to free vibration analysis of a tapered beam mounted on two-degree of freedom subsystems [8] and to others [7,20,21,22]. The paper is organized as follows: The continuous model problem is defined in Section 2.…”