set.seed(1000);
nobs <- 100000;
sdm <- 1.0;
sde <- 0.3;

mx1  <- rnorm(n=nobs,mean=0.0,sd=sdm);
mx2  <- (mx1+rnorm(n=nobs,mean=0.0,sd=sdm))/sqrt(2.0);

ex1  <- rnorm(n=nobs,mean=0.0,sd=sde);
ex2  <- rnorm(n=nobs,mean=0.0,sd=sde);

dflt <- 1.2*(1.0-mx1) + 0.2*(mx2-0.5)^2;

x1 <- plogis(mx1+ex1);
x2 <- plogis(mx2+ex2);
y  <- (dflt>2.0);

out <- cbind(x1,x2,y);
write(t(out),file="eval.dat",ncol=ncol(out));