The equation is a consequence of the symmetry of the reflection of the light about the normal:

We know from the law of reflection that the angle between R
and N
is the same as the one between L
and N
. Therefore it's also clear that:
R•N = L•N (assume all are unit vectors)
From symmetry, it's also evident that
R - (R•N)N = - (L - (L•N)N)
Now we can replace (R•N)N
with (L•N)N
from the previous equation. So we have:
R - (L•N)N = - (L - (L•N)N)
Adding (L•N)N
to both sides gives:
R = - (L - (L•N)N) + (L•N)N
= - L + (L•N)N + (L•N)N
= 2(L•N)N - L