I found the following equation form an answer written for a question.
$$2013 y^2 -xy -4026 x=0$$
But I'm confused that can I really learn how to find the positive integer solutions for $x,y$ with having only the middle school mathematics knowledge. How do I find the positive integer solutions for this equation?
Your equation is equivalent to the following $$x=\frac{2013y^2}{y+4026}$$ so you can restrict yourself on looking for $y$ such that this fraction is an integer (note that, for positive $y$, automatically $x$ is positive). Call $y+4026=t$ (so that $y=t-4026$). Recall that, if $y \ge 1$, then $t \ge 4027$. Then $$x=\frac{2013(t-4026)^2}{t} = 2013 t - 2^2 \cdot 2013^2+ \frac{2^2 \cdot 2013^3}{t}$$ The RHS is an integer if and only if $t$ divides $2^2 \cdot 2013^3$. Factor that number: $$2^2 \cdot 2013^3=2^2 \cdot 3^3 \cdot 11^3 \cdot 61^3$$ So, now you have to find all divisors $t$ of this number satisfying $t \ge 4027$ (there are a lot of them). For every such value of $t$, you can solve $y= t-4026$, and then $x = \frac{2013y^2}{t}$.