Fixed the formula for k-equivariant equation
Nicer harmonic.py