towards alternative implicit integration
towards lagraneg multipliers
somewhat working method
just works with chentanez integration method and energy term from bergou