function [fun df d2f]  = penal(x)

if x<0
    fun=x^2;
    df=2*x;
    d2f=2;;
else
    fun=0;
    df=0;
    d2f=0;
end;

end