We present numerical simulation algorithm of magnetic liquid droplet suspended in external field. Boundary integral method is used for calculation of magnetic field on the droplet surface which allows to recalculate surface velocity with a boundary integral. This repository contains numerical code which allows to run droplet simulation for constant field, slow rotating and fast rotating field. You can run these simulations either with interface julia mdrop.jl -h or putting parameters in configuration script config.jl by hand and execute julia -L config.jl FastField.jl.