kill(all)$

van=facebook.com/LearningCalculus$

f(x,n,p):=log(x)^n/x^(p+1);

F(x,n,p):=''(integrate(log(x)^n/x^(p+1),x))
         = -(log^n(x)*(p*log(x))^(-n)*gamma_incomplete(n+1, p*log(x)))/p + C;

h[n,p]:=integrate(f(x,n,p),x,1,inf);

M:genmatrix(h,5,10);

factor(M);