Stabilized multistep methods for index 2 Euler-Lagrange DAEs

We consider multistep discretizations, stabilized by -blocking, for Euler-Lagrange DAEs of index 2. Thus we may use nonstiff multistep methods with an appropriate stabilizing difference correction applied to the Lagrangian multiplier term. We show that orderp =k + 1 can be achieved for the differential variables with orderp =k for the Lagrangian multiplier fork-step difference corrected BDF methods as well as for low orderk-step Adams-Moulton methods. This approach is related to the recently proposed half-explicit Runge-Kutta methods.