Add Courant-Snyder and Normal Form calculations

[awegsche]

Add the following functionality:

• Courant-Snyder and complex Courant-Snyder from x and p
• transform C-S to normal form and back to arbitrary order