Create cdif.py

Differentiation matrices of arbitrary order and stencil width. Periodic boundary conditions implemented using the circulant property of the diagonals.