powfunc <- function(x,a,b) a * x^b
ipowfunc <- function(y,a,b) (y/a)^1/b