ini <- function(x){
  x*(x<1/2)+(1-x)*(x>=1/2)
}