Symbolic calculation in development of algorithms: split-step methods for the Gross–Pitaevskii equation