function y=my_func(x)
y = x.^3+3.*x.^2-5.*x+2;