Find suitable intrinsic growth rates

Usage

findr(
  alpha,
  k0,
  theta = 1,
  lambda0 = 0,
  interval = 0.01,
  sigma = 0.01,
  maxit = 1000
)

alpha: Numeric matrix. n_species x n_species interaction matrix
k0: Numeric. Total carrying capacity for basal species combined
theta: Numeric. Shape parameter for a Dirichlet distribution that determines relative equilibrium densities of basal species
lambda0: Numeric. Initial lambda value for the exponential decay of equilibrium densities with trophic position
interval: Numeric. Increment of lambda value.
sigma: Numeric. Degree of noise added to equilibrium densities
maxit: Integer. Maximum number of iterations to find a suitable lambda