kill(all);

f[i,j]:=sin(i*x)*cos(j*x);

genmatrix(f,3,3);

integrate(%,x);