	function y=bar(r)
	V0=1; a=100; L=5;
	y1=1./(1+exp(-a*r));y2=1./(1+exp(a*(r-L)));
	y=V0*y1.*y2;
	end