function y = cubicPoly(x)
y = x.^3 + x.^2 + x + 1; 
end