w=30, s =50,s*exp(-w^2/(2*s^2))/sqrt(2*pi) + (w/2)*erfc(-w/(s*sqrt(2)))