function feval = newton_func_penalty(theta)

global penalty;

if theta < 0
    feval = penalty * (-theta);
    return;
else;
feval = - 2 * theta^2 + 1/3 * theta^3;
return;
end;