Simulation of elastic flows subject to constraints (e.g. length and area) and numerical computation of the energy profile curve