using EcoEvoSim
config = EcoEvoConfig(
ecoDyn = unstructuredModel() do i, n, z, aux, S
n[i] * (1 - z[i][1]^2 - sum(n[j] for j in 1:S))
end,
mutationGenerator = generateMutant(
invaderPopsize = 0.001,
variance = 0.002^2
),
integrationParams = IntegrationParams(maxTime = 1e10),
extThreshold = 0.003
)
initCommunity = Community([1.0], [0.2])
evoHistory = evolve(initCommunity, config, 1500)