Problem: Let $M$ be a sphere of radius $a$ in $\mathbb{R}^3$, defined by $$ x^2 + y^2 + z^2 = a^2. $$ Parametrize the sphere using spherical coordinates $$x = a \sin \phi \cos \theta \\ y = a \sin \phi \sin \theta \\ z = a \cos \phi $$ where $0 \leq \phi < \pi, 0 \leq \theta < 2 \pi$.
Then for each $p \in M, e_1 = \partial/\partial\phi$ and $e_2 = \partial/ \partial \theta$ is a basis for the tangent space $T_p M$ for $p \in M$. Let $N_p$ be the unit outward normal vector at $p$ on the sphere.
Find the matrix of the shape operator of the sphere with respect to the basis $e_1, e_2$.
Attempt: I need to calculate $L(e_j) = - D_{e_j} N$ and then find the matrix $[a^{i}_{j}]$ such that $L(e_j) = \sum a^{i}_{j} e_i. $ Here $L$ denotes the shape operator.
I know that in cartesian coordinates, a normal unit vector is $N = \frac{1}{a} (x, y, z)$. But do I have to write this first in the $(\theta, \phi)$ basis?
I would calculate e.g. $$ D_{\partial_{\phi}} \frac{1}{a} (a \sin \phi \cos \theta, a \sin \phi \sin \theta, a \cos \phi) \\ = ( \cos \phi \cos \theta, \cos \phi \sin \theta, - \sin \phi). $$ But I need to write this as a linear combination of the basis vectors $\partial_{\phi}$ and $\partial_{\theta}$ to find the matrix representation?