f=@(x) (3+2*x-x.^2).*(x.^2+3);
g=@(x) (x.^2+5).*exp(-x);
intf=quad(f,-1,3)
xx=linspace(-1,3);
fxx=feval(f,xx);
gxx=feval(g,xx);
plot(xx,fxx,xx,gxx)
legend('f','g');
F=@(x) (3+2*x-x.^2).*(x.^2+3)-(x.^2+5).*exp(-x);
x1=fzero(F,-1);
x2=fzero(F,3);
area=quad(F,x1,x2)
