Finds all function $f:\mathbb{R}^*\to\mathbb{R}^*$ such that $$\forall x,y\in\mathbb{R}^*,y\neq-x^2\qquad f\left(x^2+y\right)=f(x)^2+\frac{f(xy)}{f(x)}$$ where $\mathbb R^*=\mathbb R\setminus\{0\}$.
Let $P(x,y)$ denote the functional equation.
What I found is:
- $f(-1)=-1$ (we can deduce it from $P(1,1),P(-1,1),P(-1,-2)$) and then conclude that $f(1)=1$
- $f(x+1)=f(x)+1$ (from $P(1,x)$)
- $f\left(x^2\right)=f(x)^2$
- $f(nx)=nf(x)$ for all integer $n$ and real $x$ (because $f(x)^2+\frac{f(nx)}{f(x)}=f\left(x^2+n\right)=f\left(x^2\right)+n$)