#infint(y_,x_):= integral(y,x,0,oo);

infint(t_^a_/(e^t_-x_),t_):= gamma(a+1)*polylog(a+1,x)/x;
infint(x_^n_/(1-e^x_),x_):= -gamma(n+1)*zeta(n+1);
infint(x_^n_/(e^x_-1),x_):= gamma(n+1)*zeta(n+1);
infint(x_^n_/(e^x_+1),x_):= gamma(n+1)*eta(n+1);
infint(sin(b_*x_)*exp(a_*x_),x_):= b/((-a)^2+b^2);
infint(cos(b_*x_)*exp(a_*x_),x_):= -a/((-a)^2+b^2);
infint(sin(x_)*exp(a_*x_),x_):= 1/((-a)^2+1);
infint(cos(x_)*exp(a_*x_),x_):= -a/((-a)^2+1);


infint(x_^n_*exp(a_*x_),x_):= n!/(-a)^(1+n);
infint(x_^n_*exp(x_),x_):= n!/(-1)^(1+n);
infint(log(x_)*exp(x_),x_):= psi(1);

infint(x_*sin(b_*x_)*exp(a_*x_),x_):= -2a*b/((-a)^2+b^2)^2;
infint(x_*cos(b_*x_)*exp(a_*x_),x_):= ((-a)^2-b^2)/((-a)^2+b^2)^2;
infint(x_*sin(x_)*exp(a_*x_),x_):= -2a/((-a)^2+1)^2;
infint(x_*cos(x_)*exp(a_*x_),x_):= ((-a)^2-1)/((-a)^2+1)^2;
infint(exp(-x_)*sin(x_)/x_,x_):= pi/4;

infint(y_):= infint(y,x);