Initial commit

Discretization seems to work, but system very sensible to parameter
changes.
2DO: Implement easy switch between cyclic/periodic and reflecting
boundary conditions.