Fourier Series as 2 loops: TwoLoop.m
Replace inner loop with  u'*v: OneLoop.m
Replace loops with (u*v')*a: NoLoop.m