In these lecture notes, we describe the current state-of-the-art for numerical simulations of large-scale structure and galaxy formation. Numerical simulations play a central role in the preparation and the exploitation of large-scale galaxy surveys, in which galaxies are the fundamental observational objects. We first describe basic methods for collisionless N-body dynamics that enable us to model dark matter accurately by solving the Vlasov-Poisson equations. We then discuss simple methods to populate dark matter halos with galaxies, such as Halo and Sub-halo Abundance Matching techniques and baryonification techniques for capturing baryonic effects on the matter distribution. We finally describe how to model the gas component by solving the Euler-Poisson equations, focusing on the foundational assumptions behind these equations, namely local thermodynamical equilibrium, and the nature of the truncation errors of the numerical scheme, namely numerical diffusion. We show a few examples of simulations of a Milky-Way-like halo without cooling, with cooling and with star formation. We finally describe different subgrid prescriptions recently developed to model star formation, supernovae feedback and active galactic nuclei and how they impact cosmological simulations.