Matlab-codes

n=10000;

tdata=randn(n,1)*10;

plotDensityComp(tdata,tdata+randn(n,1))