Saving initial energy
function for determining time step based on the nodal locations