function y=sigm(u, Vd, g0, q)
%%funkcja do sigmoidy do ers/erd
if u<=Vd
y=g0*exp(q*(u-Vd));
else
y=g0*(2-exp(q*(Vd-u)));    
end