Find all functions $f\colon\Bbb R_+\to\Bbb R_+$ that satisfy the functional equation $$(x+y)\, f\bigl(f(x)\, y\bigr)=x^2 \,f\bigl(f(x)+f(y)\bigr)$$ for all $x,y\in\mathbb{R}_+$. Here, $\mathbb{R}_+$ is the set of positive real numbers.
Attempt: There is no continuous solution $f$ such that $\lim_{x\to 0^+}f(x)$ exists. If such a solution exists, then let $L=\lim_{x\to 0^+}f(x)$. If $L>0$, then by taking $x\to 0^+$, we get $$y\,f(Ly)=0.$$ Therefore $f(Ly)=0$ for all $y>0$. This is a contradiction because the codomain of $f$ is $\mathbb{R}_+$. Therefore, $L=0$.
Now, taking $y\to 0^+$, we get $$0=x^2\,f\big(f(x)\big)$$ for all $x>0$. This means $f\big(f(x)\big)=0$. This is again a contradiction.
How do we solve the functional equation if the continuity and limit assumptions are dropped?
There is no such function $f$.
Assume $f\colon \Bbb R_+\to\Bbb R_+$ satisfies the functional equation $$(x+y)f\bigl(f(x)y\bigr)=x^2f\bigl(f(x)+f(y)\bigr) $$ for all $ x,y>0$.
Lemma. $f$ is injective.
Proof. Assume $f(x_1)=f(x_2)=a$ for some $x_1\ne x_2$. Then with $A=f(a+f(y))$, $B=-f(ay)$, $C=-yf(ay))$, we know that $x_1,x_2$ are distinct roots of $ AX^2+BX+C$. As $A\ne 0$, this implies $$x_1x_2=\frac CA=-y\cdot\frac{f(ay)}{f(a+f(y))}=y\frac BA=-y(x_1+x_2) $$ for all $y$, which is absurd. $\square$
Now whenever $x+y=x^2$, i.e., for $x>1$ and $y=x^2-x$, we can cancel $x^2$ on both sides and find $$f\bigl(f(x)(x^2-x)\bigr)=f\bigl(f(x)+f(x^2-x)\bigr) $$ and by the lemma, $$f(x)(x^2-x)=f(x)+f(x^2-x)$$ and ultimately $$\tag1f(x^2-x)=(x^2-x-1)f(x)$$ for all $x>1$. However, for any $x\in\left(1,\frac{\sqrt5+1}2\right)$, this leads to a contradiction. For example with $x=\frac32$, equation $(1)$ tells us that $$0<f\left(\tfrac34\right)=-\tfrac14 f\left(\tfrac32\right)<0,$$ contradiction.