cuberoot <- function(x)
   return(x^(1/3))