expec = integrate('x^3', 'x', 0, 1);
disp(expec, "The expectation is")