Mesh analysis for Chebyshev method

[MatthewPeterKelly]

So far, only the trapezoid and hermiteSimpson methods support mesh analysis. The key idea here is that you can compute an estimate of the error in the dynamics along the trajectory and then use this to construct a new (more accurate) mesh for the next iteration. Look up hp-adaptive meshing if your curious, or send me an email.