It is known for $f$ a smooth function and vector fields $X$, the covariant derivative obeys product rule
$$\nabla_Y(fX) = (Yf)\nabla_YX+f\nabla_Y X$$
I just did this calculation and i keep getting $\nabla_Y(fX) = X \nabla_Yf +f\nabla_Y X$.
Basically, $\nabla_Y(fx) = (Y(fX^1), \dots, Y(fX^n))$.
Now $Y(fX^1) = <grad(fX^1), Y> = <X^1gradf,Y>+ <fgrad X^1, Y>$.
So putting it back I should get $\nabla_Y(fX) = X \nabla_Yf +f\nabla_Y X$.
But my first term does not match