kill(all); g[i,j]:='integrate((j*sinh(x)^i)/(i*sinh(x)),x);

M:genmatrix(g, 3,3);