“…The first step is to use modified Boussinesq's approximation to convert the higher order non-linear PDEs [ [1] , [2] , [3] , [4] , [5] ] into second order simultaneous ODEs. With the help of the similarity technique, it is again transformed into an initial value problem described by Ali et al [ 46 , 47 ] and Uddin and Nasrin [ 48 ]. Using numerical methods, we successfully solve coupled non-linear ODEs [ [24] , [25] , [26] , [27] ] with boundary conditions [ 28 , 29 ] applying the shooting technique.…”