2.3.33 Gradient and directional derivative

condition: Find the derivative \( d z / d l \) of the function \( z=x^{2} y^{2}-x y^{3}-3 y+1 \) at the point \( M(2; 1) \) in the direction from this point to the origin.