syms x a;
f=x.*exp(-a*x);int(f)
syms b x c;
f=1/sqrt(x^2+b*x+c);
int(f)
f=x*tan(x);
int(f)
syms y;
f=sin(y*x);
int(f^2,y)
f=exp(-3*x)*cos(x);
int(f,0,inf)
f=sqrt(x);
int(pi*f^2,0,1)
f=x^2;
int(pi*f^2,0,1)
x=linspace(0,pi,50);y=cos(pi*x);a=trapz(x,y);100*abs((int(cos(pi*r),r,0,pi)-a)/int(cos(pi*r),r,0,pi))
x=linspace(-2.2,2.2,1000);y=exp(-x.^2);a=trapz(x,y);100*abs((int(exp(-r.^2),r,-2.2,2.2)-a)/int(exp(-r.^2),r,-2.2,2.2))