Analytical Series Expansion for Efficient Gradient Evaluation in Multi-Qubit Optimal Control