Eukaryotic cell motility involves complex interactions of signalling molecules, cytoskeleton, cell membrane, and mechanics interacting in space and time. Collectively, these components are used by the cell to interpret and respond to external stimuli, leading to polarization, protrusion, adhesion formation, and myosin-facilitated retraction. When these processes are choreographed correctly, shape change and motility results. A wealth of experimental data have identified numerous molecular constituents involved in these processes, but the complexity of their interactions and spatial organization make this a challenging problem to understand. This has motivated theoretical and computational approaches with simplified caricatures of cell structure and behaviour, each aiming to gain better understanding of certain kinds of cells and/or repertoire of behaviour. Reaction-diffusion (RD) equations as well as equations of viscoelastic flows have been used to describe the motility machinery. In this review, we describe some of the recent computational models for cell motility, concentrating on simulations of cell shape changes (mainly in two but also three dimensions). The problem is challenging not only due to the difficulty of abstracting and simplifying biological complexity but also because computing RD or fluid flow equations in deforming regions, known as a free-boundary problem, is an extremely challenging problem in applied mathematics. Here we describe the distinct approaches, comparing their strengths and weaknesses, and the kinds of biological questions that they have been able to address.