function y=rho(x)

y = x.^3 - x.^2;