The smoothed-particles hydrodynamics method (SPH) – in which the continuum is modeled as the arbitrary lattice of interacting particles/interpolation points – is often advantageous in modeling scenarios involving extreme deformations, fluid motion, and fragmentation over traditional mesh-based techniques. This study is focused on developing an SPH solver for simulating the motion of incompressible fluids, including the following tasks: •??Development of the interface with the pre- and post-processing software. •??Implementation of the SPH solver for Navier-Stokes equations in MATLAB. •??Identification of the most efficient technique for density calculation within the SPH solver through comparison of the “Summation Density” and the “Rate Density” algorithms.