This problem deals with the effects of double diffusive, mixed convective flow of an incompressible viscous fluid through a vertical heated plate embedded in a non-Darcy porous medium under the influence of variable fluid properties numerically. The governing equations are modeled for the double diffusive boundary layer flow to understand the behaviour of velocity, temperature and concentration for variable fluid properties namely permeability, porosity, thermal conductivity and solutal diffusivity of the physical modal. Using a suitable similarity transformation, the highly nonlinear coupled PDE's are reduced into a set of coupled ordinary differential equations. By applying the Shooting technique with the help of Runge-Kutta-Fehlberg integral scheme and Newton-Raphson method, numerical computations have been carried out and are illustrated through the graphs to analyze the behaviour of velocity, temperature, concentration, Skin friction, Nusselt number and Sherwood numbers for the non-dimensional parameters of the physical system for both uniform permeability and variable permeability. To validate our numerical scheme, present results are compared with previously published work in the absence of few non-dimensional parameters and found to be in good agreement up to six decimal places of accuracy. © 2018, Global Digital Central. All rights reserved.