BlindSimu <- function(X1=rexp(100, rate=1.0),C1=rexp(100, rate=1.0),
                  X2=rexp(100, rate=1.0),C2=rexp(100, rate=1.0)) {

          Z1<-pmin(X1,C1)
          D1<-as.numeric(C1>X1)
          Z2<-pmin(X2,C2)
          D2<-as.numeric(C2>X2)

          temp1<-KMNAatT(Z1, D1, 0.5)
          temp2<-KMNAatT(Z2, D2, 0.5)
          Tvar<- (temp1$std.err)^2 + (temp2$std.err)^2 

          temp3<-KMNAatT(c(Z1,Z2), c(D1,D2), 0.5)
          Bvar<- 4*(temp3$std.err)^2

          list(Tvariance=Tvar,Bvariance=Bvar,
               TnelsonVar=(temp1$NelsonAalenVar + temp2$NelsonAalenVar),
               BnelsonVar= 4*(temp3$NelsonAalenVar), 
               Best=temp3$surv,BGW=temp3$std.err) 
}