function [d]=df33(x)
if x>=0
    d=3*x^2+1;
else
    d=-3*x^2+1;
end