[Dev] Constrain endpoint controls

[jbcaillau]

to avoid singularities (check num scheme for dynamics)

• for euler + midpoint add constraint uN = uN-1 (as uN is possibly free)
• for euler_n add constraint u0 = u1 (as u0 is possibly free)

NB. nothing to do for trapeze

check #94 (comment)