An elegant way to do this is via function handles:
x = linspace(0, 11, 100);
y = @(x) .0012*x.^1.5;
yforArea = @(q) max(y(x)) - y(q);
area = integral(yforArea, 0, max(x))
Or just via trapezoidal integration directly:
x = linspace(0, 11, 100);
y = .0012*x.^1.5;
z = trapz(x, max(y)-y)
