We present a numerical method which computes the motion of complex solid/liquid boundaries in crystal growth. The model we solve includes physical effects such as crystalline anisotropy, surface tension, molecular kinetics and undercooling. The method is based on two ideas. First, the equations of motion are recast as a single history-dependent boundary integral equation on the solid/liquid bou...