n=10000;
x=linspace(0,1.5,(n+1));
h=(1.5-0)/n
fx=0.84885406+31.51924706*x-137.66731262*x.^2+240.55831238*x.^3-171.45245361*x.^4+41.95066071*x.^5;
I=0;
for ii=1:n
I=I+(h/2)*(fx(ii)+fx(ii+1));
end
I