function y=StDe(x)

y = (abs(x).^2).*(abs(x)<1) + (2*abs(x)-1).*(abs(x)>=1) ;

