I want to show that the sum of integer squares from $i=1$ to $n$ is $\frac{n(n+1)(2n+1)}{6}$ I've watched some videos and read other posts about it but haven't been able to find anything that makes it click. I know the rate of the difference of consecutive partial sums for $S(n)$ shows $S$ is a cubic function. I tried making use of this using a system of equations with $S(n)=An^3+Bn^2+Cn+D$ at $n=0, 1, 2, 3$, and I ended up getting $S(n)=-\cfrac{643}{24}n^3-\cfrac{131}{8}n^2+\cfrac{137}{12}n$ but I'm unsure how to factor this
Is there a simple approach to this that doesn't involve some weird collapsing sum?

A beautiful proof without words that I recently learned from a friend of mine (credits to KK):
Some words of explanation: we have a function $f$ with a constant gradient defined over a domain which is an equilateral triangle. When we consider $g=f+f_{\omega}+f_{\omega^2}$, where $f_\omega$ and $f_{\omega^2}$ are the functions defined over the domain rotated by $120^\circ$ and $240^\circ$, we have that $g$ has a null gradient, hence it is constant.