First commit. The 1D calculations work. The 2D part is not implemented.
add comments