intfun = @(t) t.*cos(t);
y = integral(intfun,0,pi/2)
