kill(all);

f[i,j]:=%e^((i-1)*x)/(%e^((j-1)*x)+%e^((i-1)*x));

M:genmatrix(f,4,4)$

'integrate(M,x);

integrate(M,x);