source("stat.s")

 plot(x0,1-0.6*pnorm(x0,0.5,0.1)-0.4*pnorm(x0,1,0.1),type="l")
 statS("sim3/ICM/surv1.1")
#M.C.stderr= 0.0126732609225946 
 title(sub="NPMLE, maxMCse=0.0127")