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