lagrangian method