function q = twoCornersOneBall2DnonConstant(X1,X2,varargin) % two corners and one ball with non constant contrast seti = varargin{end}; R = seti.rCD/2; q = twoCornersOneBall2D(X1,X2,seti); fac = seti.k; c = 0.2; q = (1+c*sin(fac*X1)).*q; q = double(q); end