Cosmology and massive neutrinos have made good friends with one another in the contemporary era. The large-scale structure surveys with increasing precision show promise to provide the tightest constraint on the sum of neutrino masses. The push towards smaller scales and lower redshifts in the observational data means the modelling of non-linear clustering of matter becomes unavoidable. The N-body simulation remains the only reliable technique for such a task. In this talk I will discuss the inclusion of massive neutrinos in N-body simulations via several methods. Each method comes with its own question of cost-benefit analysis and suitability for a particular context. The computational resource constraint forces us to confront these questions and be appropriately prepared for the upcoming data from a wide range of observables.