Int( 0 to pi) Int( 0 to 1/2) x cos(xy) cos^2(pi x) dxdy
=Int( 0 to pi) Int( 0 to 1/2) x cos(xy) (1+ cos (pi x))/2 dxdy
=(1/2) Int( 0 to pi) Int( 0 to 1/2) x cos(xy) (1+ cos( 2 pi x)) dxdy
=(1/2) Int( 0 to pi) Int( 0 to 1/2) x cos(xy) + x cos(xy) cos (2 pi x) dxdy
=(1/2) Int( 0 to pi) Int( 0 to 1/2) x cos(xy) + x cos(xy) cos (2 pi x) dxdy
=(1/4) Int( 0 to pi) Int( 0 to 1/2) 2 x cos(xy) + x.[2 cos(xy) cos (2 pi x)] dxdy
=(1/4) Int( 0 to pi) Int( 0 to 1/2) 2 x cos(xy) + x.[ cos(2pi x + yx) + cos(2 pi x - yx) ] dxdy
=(1/4) Int( 0 to pi) Int( 0 to 1/2) 2 x cos(xy) + xcos(2pi x + yx) + xcos(2 pi x - yx) dxdy
Integration by parts
=(1/4) Int( 0 to pi) 2x sin(yx) - 2(1/y^2) cos(yx) + x[ sin(2 pi x+ yx)/(2pi +y) + sin(2 pi x- yx)/(2pi -y) ] - [ -cos (2 pi x+ yx)/(2pi +y)^2 - cos(2 pi x- yx)/(2pi -y)^2 ] (0 to pi) dy