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);