function [X,Y,delta]=artificialData(n,d)
X=[0.5+1.5.*randn(n*0.2,d);0.6+1.5.*randn(n*0.2,d);0.7+1.5.*randn(n*0.2,d);0.8+1.5.*randn(n*0.2,d);0.9+1.5.*randn(n*0.2,d);];
X(:,d*0.05+1:end)=1.5.*randn(n,d-d*0.05);
l=find(rand(1,n)<0.4);
X(l,d*0.05+1:d*0.1)=2+randn(length(l),d*0.05);
l=find(rand(1,n)<0.5);
X(l,d*0.1+1:d*0.2)=1+randn(length(l),d*0.1);
l=find(rand(1,n)<0.7);
X(l,d*0.2+1:d*0.3)=0.5+randn(length(l),d*0.1);
stime=sum(X(:,1:n*0.05),2)+0.01.*randn(n,1);
ctime=n*0.1+randn(n,1);
delta=stime<ctime;
stime(delta==0)=ctime(delta==0);
Y=stime;