function y=sig7(x)
%Vd = 0 - upper part sigmoid, Vd = 10 -lower part of sigmoid
Vd =7; 
q=0.34;
g = 25;
if x<=Vd y=g*exp(q*(x-Vd)); else y=g*(2-exp(q*(Vd-x)));
end;