Calculations of the behavior of liquid crystal droplets in magnetic fields applied along the axis of cylindrical symmetry have been carried out using relaxation methods on the finite difference approximation to the Euler–Lagrange equations derived from the elastic strain energy. The nature of the solutions and how they change with field is discussed. The application of these methods to problems in magnetism is suggested.