kill(all);

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

genmatrix(f,3,3);

integrate(%,x);