If a 1D random walk, along the x-axis has has a closed for expression (Due to the central limit theorem)
$P(N)$ = $\frac{1}{\sqrt{2\pi N}}*e^{\frac{-n^2}{2N}}$, where $N=2n$.
This is a normal distribution with variance $N$, and standard deviation $N^2$.
I used this expression to generate a closed form for a 2D symmetric random walk on a lattice with unit step length.
I think that the x-component $P_x(n)$ will be the same as the 1D probability closed form, except with variance $\frac{N}{2}$ instead of $N$. This makes the x and y components of the equation,
$P_x(N)$ = $\frac{1}{\sqrt{\pi N}}*e^{\frac{-x^2}{N}}$,
$P_y(N)$ = $\frac{1}{\sqrt{\pi N}}*e^{\frac{-y^2}{N}}$,
and hence $P(r)=P_x(n)*P_y(n)$, which gives an expression
$P(e)$ = $\frac{1}{{\pi N}}*e^{\frac{-r^2}{N}}$, where $r=[x,y]$ <-- r is a vector with components
Is this logic correct, because when I try small values of N, the expression has a very large error. Or, what is the correction that I have to make for the expression to make sense?