Minor updates.
- Changes Pk_int_orig[:] to Pk_int_orig.copy() as previous wasn't
copying as expected.
- Small update in calculate() and main().