We present the development of pyLLE, a freely accessible Lugiato-Lefever equation solver programmed in Python and Julia and optimized for the simulation of microresonator frequency combs. Examples illustrating its operation and performance are presented.