integrate 1/(s*sqrt(2*pi)) * e^-((y - x)^2/(2*s^2)) * 1/(r*sqrt(2*pi)) * e^-((x-0)^2/(2*r^2)) dx from -inf to + inf where s=2, r=15