Add constraints to ExIntegrator; test against FEM-style solution. Correct frequencies are now generated.