R/potence.R

In FractalParameterEstimation: Simulation and Parameter Estimation of Randomized Sierpinski Carpets using the p-p-p-q-Model

Defines functions potence

Documented in potence

potence <-
function(a, b) {
  c = a
  if(b < 0 || a < 0)
    stop("invalid input")
    if( b == 0) return(1)
      while(b > 1) {
        c = c*a
        b = b-1
      }
  return(c)
}

Try the FractalParameterEstimation package in your browser