def METHOD_NAME():
    global energies
    system.integrator.run(int_steps)
    visualizer.update()
    energies = system.analysis.energy()
    plot.set_xdata(np.append(plot.get_xdata(), system.time))
    plot.set_ydata(np.append(plot.get_ydata(), energies['total']))