function f = myfun(x)

f = 3.06^2 * x.^2 .* (1+x) - (1-x).^2 .* (3+x);