We model the flow past a sphere in a viscous, incompressible low Reynolds number stratified fluid. The flow has uniform velocity and linear stratification. The governing set of equations, Navier-Stokes, describing the above flow were numerically solved by direct integration using boundary-fitted coordinates and a Pressure Correction Method with finite-difference approximations. In order to solve the pressure linear system with high accuracy we use the BiCGstab method. Numerical experiments are performed for different flow conditions by varying the Reynolds number. Results are presented for density and velocity fields. Simulated sphere drag is also compared to the classical drag curve cited in the literature.