FYI, the ESA (European Space Agency) has a python package, pykep, that implements these types of algorithms for you: https://esa.github.io/pykep/index.html

Check out the pontryagin optimization module, which uses the maximum principle: https://esa.github.io/pykep/documentation/pontryagin.html