function feval = newton_func(theta)

feval = - 2 * theta^2 + 1/3 * theta^3;