Let $Sp(2n):=\{A\in\mathbb{R}^{2n\times 2n}\mid A^tA_0A=A_0\}$ be the group of symplectomorphisms from $(\mathbb{R}^{2n},\omega_0)$ to itself, where:
\begin{gather} A_0 := \begin{bmatrix}{} 0 & I\\ -I & 0\\ \end{bmatrix}\in\mathbb{R}^{2n\times 2n} \end{gather}
represents the standard symplectic form $\omega_0$ with respect to the canonical basis of $\mathbb{R}^{2n}$. Prove that $Sp(2n)$ is an embedded submanifold of $GL(2n)$ and has dimension $2n^2+n$.
I know the essential idea is to look at the map: \begin{align*} f:GL(2n)&\to \text{Sympl}(2n)\\ A & \mapsto A^tA_0A \end{align*}
where $\text{Sympl}(2n):=\{A\in\mathbb{R}^{2n\times 2n}\mid A=-A^t\text{ and }\det A\neq 0\}$, which is the submanifold of symplectic forms and has dimension $\frac{(2n)^2-2n}{2}$.
Considering that $f$ is a submersion, we have $Sp(2n)=f^{-1}(A_0)$ is an embedded submanifold with dimension $\dim (GL(n))-\dim (\text{Sympl}(2n))=(2n)^2-\left(\frac{(2n)^2-2n}{2}\right)=2n^2+n$ by the Regular Value Theorem.
My question is really basic: how do I prove $f$ is a submersion? I've tried to calculate $df_A(M)$ by taking the curve $\alpha(t)=A+tM$, so that: \begin{align*} df_A(M)&=(f\circ \alpha)'(0)\\ &=(A^tA_0A+tA^tA_0M+tM^tA_0A+t^2M^tA_0M)'(0)\\ &=A^tA_0M+M^tA_0A \end{align*} but I don't know how to prove this is surjective, and it seems complicated. Is there some trick to it? Or is there a better way?