Given an inverse square law $\ddot{\vec{r}}=-\frac{\mu}{r^2}\hat{r}$, I define the Angular momentum per unit mass as $\vec{H}=\vec{r}\times\dot{\vec{r}}$. Showing it's constant is strightfoward. Then I defined the Laplace–Runge–Lenz vector per unit mass
$\vec{c}=\dot{\vec{r}}\times\vec{H}-\mu\hat{r}$
I think (yet not sure) this is the normal one scaled down by $m^2$. I now try to take the derivative and get
$\frac{d\vec{c}}{dt}=\ddot{\vec{r}}\times\vec{H}+\dot{\vec{r}}\times\dot{\vec{H}}-\mu\dot{\hat{r}}=\ddot{\vec{r}}\times\vec{H}-\mu\dot{\hat{r}}=-\frac{\mu}{r^2}\hat{r}\times\vec{H}-\mu\dot{\hat{r}}$
I'm not sure how to get to $0$. I thought of substituting $\vec{H}$'s definition and using a triple cross identity but that wont yield anything.