broke out code into single and double pen programs separately. The refactoring continues
found one typo in the double pen rk4 formulae. No change in wrong data
removed old r4k implementation from double_pen.cpp (it was previously commented out)
major refactor of single_pen.cpp to aid the introduction of stability tests
began refactoring double_pen.cpp for stability analysis