function y = ID(x)
y = (1-tanh(3*x))/2;
