function y = df(x)
y = 3*x^2 + 1;