integrate((π^4 x y (a+9 c+8 b cos((π z)/l)+18 c cos((2 π z)/l)) sin^2((π z)/l))/l^4)dz