Let $a,b \in \mathbb{R}$ such that $0<a<b$ and let $n \in \mathbb{N}$. Consider the function $$f(x_1,...x_n)=\frac{(x_1\cdot... \cdot x_n)}{(a+x_1)(x_1+x_2)...(x_{n-1}+x_n)(x_n+b)}$$ where $x_i \in [a,b]$ and $n$ is fixed. Find the maximum of such function.
I was trying to calculate partial derivatives and find if the maximum occurs in the interior. If we write for symmetry that $x_0=a$ and $x_{n+1}=b$ i got that $\frac{df}{dx_i}= (x_{i-1}x_{i+1} -x_i^2) \cdot S$, where S is nonzero for $i=1,2,...,n$. So maximum in the interior might occure iff $x_{i-1}x_{i+1}-x_i^2=0$ for $i=1,2,...,n$. But then I also have to check boundary.
Are partial derivatives good way to tackle this problem or should I try with some inequalities? Any help will be greatly appreciated.