Generating differential operators for the basis functions in the variational LCAO‐type methods