factor 1+cos(x)+cos(y)+cos(x-y)