For a given discretised interval:
- Discretise lagrangian
- Form differential equation system
- Solve differential equation system
Then loop for each interval
- Quadrature computation
- Discretise lagrangian
- Mostly done, test
- Form differential equation system
- Solve differential equation system