function y=sigma(x)

y = (9*x.^3 + 19* x.^2 - 5.*x + 1)/24;