A finite difference continuation method for computing energy levels of Bose-Einstein condensates