Find this Determinant

I guess it's easier to go for a reduction formula. I proceed along the generalization mentioned in comment:

Call the determinant $M = M_n(x_1,x_2,\cdots, x_n)$

$$\displaystyle \begin{align}M &= \left|\begin{matrix}\dfrac{1}{x_1+x_1} & \cdots & \dfrac{1}{x_1+x_n}\\ \dfrac{1}{x_2+x_1} & \cdots & \dfrac{1}{x_2+x_n}\\ \cdot & \cdot &\cdot \\ \dfrac{1}{x_n+x_1} & \cdots & \dfrac{1}{x_n+x_n}\end{matrix}\right| \\& \text{subtract R1 from each row} \\&= \left|\begin{matrix}\dfrac{1}{x_1+x_1} & \cdots & \dfrac{1}{x_1+x_n}\\ \dfrac{x_1 - x_2}{(x_2+x_1)(x_1+x_1)} & \cdots & \dfrac{x_1 - x_2}{(x_2+x_n)(x_1+x_n)}\\ \cdot & \cdot &\cdot \\ \dfrac{x_1 - x_n}{(x_n+x_1)(x_1+x_1)} & \cdots & \dfrac{x_1 - x_n}{(x_n+x_n)(x_1+x_n)}\end{matrix}\right| \\&= \prod\limits_{j=2}^{n} (x_1 - x_j)\times \prod\limits_{j=1}^{n}\frac{1}{(x_1+x_j)} \times \left|\begin{matrix}1 & \cdots & 1\\ \dfrac{1}{(x_2+x_1)} & \cdots & \dfrac{1}{(x_2+x_n)}\\ \cdot & \cdot &\cdot \\ \dfrac{1}{(x_n+x_1)} & \cdots & \dfrac{1}{(x_n+x_n)}\end{matrix}\right|\end{align}$$

Again subtracting $\textrm{C1}$ from each column,

$$\begin{align}\left|\begin{matrix}1 & \cdots & 1\\ \dfrac{1}{(x_2+x_1)} & \cdots & \dfrac{1}{(x_2+x_n)}\\ \cdot & \cdot &\cdot \\ \dfrac{1}{(x_n+x_1)} & \cdots & \dfrac{1}{(x_n+x_n)}\end{matrix}\right| &= \left|\begin{matrix}1 & \cdots & 0\\ \dfrac{1}{(x_2+x_1)} & \cdots & \dfrac{x_1 - x_n}{(x_2+x_n)(x_2+x_1)}\\ \cdot & \cdot &\cdot \\ \dfrac{1}{(x_n+x_1)} & \cdots & \dfrac{x_1 - x_n}{(x_n+x_n)(x_n+x_1)}\end{matrix}\right|\\&= \prod\limits_{j=2}^{n} (x_1-x_j) \times \prod\limits_{j=2}^n \frac{1}{x_j+x_1} \times \left|\begin{matrix}1 & \cdots & 0\\ 1 & \cdots & \dfrac{1}{(x_2+x_n)}\\ \cdot & \cdot &\cdot \\ 1 & \cdots & \dfrac{1}{(x_n+x_n)}\end{matrix}\right|\\&= \prod\limits_{j=2}^{n} (x_1-x_j) \times \prod\limits_{j=2}^n \frac{1}{x_j+x_1} \times M_{n-1}(x_2,x_3,\cdots ,x_n)\end{align}$$

Hence, $$M_n(x_1,x_2,\cdots, x_n) = \frac{\prod\limits_{1 \le j < i \le n} (x_i-x_j)^2}{\prod\limits_{1 \le i,j \le n} (x_i+x_j)}$$

We consider a general case. Let $$ |A|=\left|\begin{array}{}{\dfrac1{a_1+b_1}\cdots\dfrac1{a_1+b_n}\\ \vdots\hspace{20 mm} \vdots \\\dfrac1{a_n+b_1}\cdots\dfrac1{a_n+b_n}} \end{array}\right| $$ In your case, just set $a_i=i,b_j=j$.

By multiplying $i^{th}$ row with $\prod\limits_{j=1}^{n}(a_i+b_j)$ each, we have $$ |A|=\dfrac1{\prod\limits_{i,j=1}^{n}(a_i+b_j)}\left|\begin{array}{}{\prod\limits_{j=2}^{n}(a_1+b_j)\cdots\prod\limits_{j=1}^{n-1}(a_1+b_j)\\ \vdots\hspace{20 mm} \vdots \\\prod\limits_{j=2}^{n}(a_n+b_j)\cdots\prod\limits_{j=1}^{n-1}(a_n+b_j)} \end{array}\right| \hspace{20 mm}(1) $$ Since each $(i,j)$ element in $(1)$ can be expanded into polynomial, i.e $$ \prod\limits_{k=1,k\ne j}^{n}(a_i+b_k)=a_i^{n-1}+\sum\limits_{k=1,k\ne j}^{n}{b_k}a_i^{n-2}+\cdots+\prod\limits_{k=1,k\ne j}^{n}{b_k}=\pmatrix{1&a_i&\cdots &a_i^{n-1}}\pmatrix{\prod\limits_{k=1,k\ne j}^{n}{b_k}\\ \vdots \\ \sum\limits_{k=1,k\ne j}^{n}{b_k}\\1} $$ So \begin{align} \left|\begin{array}{}{\prod\limits_{j=2}^{n}(a_1+b_j)\cdots\prod\limits_{j=1}^{n-1}(a_1+b_j)\\ \vdots\hspace{20 mm} \vdots \\\prod\limits_{j=2}^{n}(a_n+b_j)\cdots\prod\limits_{j=1}^{n-1}(a_n+b_j)} \end{array}\right|&=\left|\begin{array}{}{1 \hspace{4 mm} a_1 \hspace{4 mm}\cdots\hspace{4 mm}a_1^{n-1}\\ \vdots\hspace{20 mm} \vdots \\ 1 \hspace{4 mm} a_n \hspace{4 mm}\cdots\hspace{4 mm}a_n^{n-1}} \end{array}\right|\left|\begin{array}{}{\prod\limits_{k=1,k\ne 1}^{n}{b_k} \hspace{4 mm} \prod\limits_{k=1,k\ne 2}^{n}{b_k} \hspace{4 mm}\cdots\hspace{4 mm}\prod\limits_{k=1,k\ne n}^{n}{b_k}\\ \vdots\hspace{40 mm} \vdots \\ \sum\limits_{k=1,k\ne 1}^{n}{b_k} \hspace{4 mm} \sum\limits_{k=1,k\ne 2}^{n}{b_k} \hspace{4 mm}\cdots\hspace{4 mm}\sum\limits_{k=1,k\ne n}^{n}{b_k}\\1\hspace{20 mm}1\hspace{7 mm}\cdots\hspace{7 mm} 1} \end{array}\right| \\ &=\left|\begin{array}{}{1 \hspace{4 mm} a_1 \hspace{4 mm}\cdots\hspace{4 mm}a_1^{n-1}\\ \vdots\hspace{20 mm} \vdots \\ 1 \hspace{4 mm} a_n \hspace{4 mm}\cdots\hspace{4 mm}a_n^{n-1}} \end{array}\right|f(b_1\cdots b_n) \\ &=R_1(b_1\cdots b_n)\tag{2} \end{align} where $f(b_1\cdots b_n)$ is a polynomial of $b_1\cdots b_n$.

Now by multiplying $j^{th}$ column with $\prod\limits_{i=1}^{n}(a_i+b_j)$ each, we have $$ |A|=\dfrac1{\prod\limits_{i,j=1}^{n}(a_i+b_j)}\left|\begin{array}{}{\prod\limits_{i=2}^{n}(a_i+b_1)\cdots\prod\limits_{i=1}^{n-1}(a_i+b_1)\\ \vdots\hspace{20 mm} \vdots \\\prod\limits_{i=2}^{n}(a_i+b_n)\cdots\prod\limits_{i=1}^{n-1}(a_i+b_n)} \end{array}\right| \hspace{20 mm}(3) $$

Expand each $(i,j)$ element in $(3)$ into polynomial of $b_j$, we have $$ \prod\limits_{k=1,k\ne i}^{n}(a_k+b_j)=b_j^{n-1}+\sum\limits_{k=1,k\ne i}^{n}{a_k}b_j^{n-2}+\cdots+\prod\limits_{k=1,k\ne i}^{n}{a_k}=\pmatrix{1&b_j&\cdots &b_j^{n-1}}\pmatrix{\prod\limits_{k=1,k\ne i}^{n}{a_k}\\ \vdots \\ \sum\limits_{k=1,k\ne i}^{n}{a_k}\\1} $$ So \begin{align} \left|\begin{array}{}{\prod\limits_{i=2}^{n}(a_i+b_1)\cdots\prod\limits_{i=1}^{n-1}(a_i+b_1)\\ \vdots\hspace{20 mm} \vdots \\\prod\limits_{i=2}^{n}(a_i+b_n)\cdots\prod\limits_{i=1}^{n-1}(a_i+b_n)} \end{array}\right| &=\left|\begin{array}{}{1 \hspace{4 mm} b_1 \hspace{4 mm}\cdots\hspace{4 mm}b_1^{n-1}\\ \vdots\hspace{20 mm} \vdots \\ 1 \hspace{4 mm} b_n \hspace{4 mm}\cdots\hspace{4 mm}b_n^{n-1}} \end{array}\right|\left|\begin{array}{}{\prod\limits_{k=1,k\ne 1}^{n}{a_k} \hspace{4 mm} \prod\limits_{k=1,k\ne 2}^{n}{a_k} \hspace{4 mm}\cdots\hspace{4 mm}\prod\limits_{k=1,k\ne n}^{n}{a_k}\\ \vdots\hspace{40 mm} \vdots \\ \sum\limits_{k=1,k\ne 1}^{n}{a_k} \hspace{4 mm} \sum\limits_{k=1,k\ne 2}^{n}{a_k} \hspace{4 mm}\cdots\hspace{4 mm}\sum\limits_{k=1,k\ne n}^{n}{a_k}\\1\hspace{20 mm}1\hspace{7 mm}\cdots\hspace{7 mm} 1} \end{array}\right| \\ &=\left|\begin{array}{}{1 \hspace{4 mm} b_1 \hspace{4 mm}\cdots\hspace{4 mm}b_1^{n-1}\\ \vdots\hspace{20 mm} \vdots \\ 1 \hspace{4 mm} b_n \hspace{4 mm}\cdots\hspace{4 mm}b_n^{n-1}} \end{array}\right|f(a_1\cdots a_n) \\ &=R_2(b_1\cdots b_n)\tag{4} \end{align}

Now consider $R_1-R_2$ as polynomial of $b_1\cdots b_n$. Since the number of roots of $R_1-R_2$ is more than its degree, $R_1-R_2$ vanishes and the coefficients of $b_1\cdots b_n$ in $R_1$ must be equal to those of $R_2$. By comparing $(2)$ and $(4)$, we have $$ f(b_1\cdots b_n)=\left|\begin{array}{}{1 \hspace{4 mm} b_1 \hspace{4 mm}\cdots\hspace{4 mm}b_1^{n-1}\\ \vdots\hspace{20 mm} \vdots \\ 1 \hspace{4 mm} b_n \hspace{4 mm}\cdots\hspace{4 mm}b_n^{n-1}} \end{array}\right| $$ So from $(1)$ and $(2)$, we have \begin{align} |A|&=\dfrac1{\prod\limits_{i,j=1}^{n}(a_i+b_j)}\left|\begin{array}{}{1 \hspace{4 mm} a_1 \hspace{4 mm}\cdots\hspace{4 mm}a_1^{n-1}\\ \vdots\hspace{20 mm} \vdots \\ 1 \hspace{4 mm} a_n \hspace{4 mm}\cdots\hspace{4 mm}a_n^{n-1}} \end{array}\right|\left|\begin{array}{}{1 \hspace{4 mm} b_1 \hspace{4 mm}\cdots\hspace{4 mm}b_1^{n-1}\\ \vdots\hspace{20 mm} \vdots \\ 1 \hspace{4 mm} b_n \hspace{4 mm}\cdots\hspace{4 mm}b_n^{n-1}} \end{array}\right| \\ &=\dfrac{\prod\limits_{1\le i<j\le n} (a_j-a_i)\prod\limits_{1\le i<j\le n} (b_j-b_i)}{\prod\limits_{i,j=1}^{n}(a_i+b_j)} \end{align} In your case $$ |A|=\dfrac{\prod\limits_{1\le i<j\le n} (j-i)^2}{\prod\limits_{i,j=1}^{n}(i+j)} $$