Sigma has to be greater than 0; sigma = 0 is tested in the invalid values. But if y= 0, the boundary behavior should be 0.