Algorithm code for solving unconstrained unitary quantum brachistochrone problems (Mathematica notebook)

We provide a basic version of the algorithm code introduced in our manuscript arXiv:1906.02934. To compare the driving time of different Hamiltonians, we use a uniform energetic constraint, as prescribed in Phys. Rev. Lett. 120, 060409.Other energetic constraint can be considered, such as the standard deviation of the Hamiltonian. The renormalization used is optional, and can be entirely removed. In this case, each Hamiltonian in the sequence will drive the system in unit time, along paths with different length in the space of states.