pan := n**3/3 + 9*n**2/2 - n/3;
