Given stochastic processes $X_t, Y_t$ of form given below, we want to find $d(X_t Y_t)$. Assume $$X_t=\sigma_tW_t + \mu_t t$$ $$Y_t=\rho_tW_t + \nu_t t$$
I have two different answers depending on how I derive them, and both of them are wrong, but I don't know why.
Method 1:
$X_tY_t=\sigma_t\rho_tW_t^2+(\sigma_t\nu_t+\mu_t\rho_t)tW_t+\mu_t\nu_tt^2$
Therefore $X_tY_t=f(W_t,t).$ Using Taylor expansion: $df(W_t,t)=\frac{\delta f}{\delta W_t}dW_t+(\frac{1}{2}\frac{\delta^2f}{\delta W_t^2}+\frac{\delta f}{\delta t})dt.$
Therefore: $d(X_tY_t)=(2\sigma_t\rho_tW_t+\sigma_t\nu_t+\mu_t\rho_t)dW_t+2(\sigma_t\rho_t+\mu_t\nu_tt)dt$
Method 2:
However, if we use a different approach, by defining instead $g(X_t,Y_t)=X_tY_t$, then a Taylor polynomial gives: $$d(X_t,Y_t)=\frac{\delta g}{\delta X_t}dX_t+\frac{\delta ^2g}{\delta X_t^2}(dX_t)^2+\frac{\delta g}{\delta Y_t}dY_t+\frac{\delta^2 g}{\delta Y_t^2}(dY_t)^2+\frac{\delta^2 g}{\delta X_t\delta Y_t}dX_tdY_t$$ This results in: $$d(X_tY_t)=Y_tdX_t+X_tdY_t+2dX_tdY_t$$ If I were to expand this, it would be quite big and definitely not equal to the result from (1).
Both of these are wrong, because my textbook gives an even different result, namely $$d(X_tY_t)=Y_tdX_t+X_tdY_t+ \sigma_t\rho_tdt$$ I don't know how exactly they got there other than that they used $X_tY_t=\frac{1}{2}((X_t+Y_t)^2 - X_t^2 - Y_t^2)$.
I must be fundamentally misunderstanding something about stochastic calculus (I've only just started with it).
What am I doing wrong, and what misunderstanding do I have?