Fitting gravitational microlensing models to the data is a challenging task and will be more challenging when the upcoming NASA flagship mission - Nancy Grace Roman Telescope (formerly WFIRST) - begins its observing program. We present MulensModel package, which aims at allowing newcomers to conduct microlensing research. MulensModel allows several higher-order effects to be modeled (like binary-lens events with source limb-darkening). MulensModel is written in an object-oriented Python 3 and is already able to model simulated Roman Telescope data. Full documentation of the code is provided and the code is continuously developed.