## CAMBAM Summer School 2022
## Mathematical Oncology
## A selection of tumour growth models

# Exponential
x1' = r*x1

# Logistic
x2' = r*x2*(1-x2/K)

# Allee
x3' = r*x3*(1-x3/K)*(x2 - a)

# Gompertz 
x4' = r*x4*log(K/x4)

# von Bertalanffy
x5' = r*x5^(2.0/3.0) - mu*x5

# Surface
x6' = r*x6*(x6 + K)^(-1.0/3)

init x1=0.10
init x2=0.10
init x3=0.10
init x4=0.10
init x5=0.10
init x6=0.10

params K=1.8
params a=0.05
params mu=0.1
params r=2.0

@ total=4
@ bound=100000
@ Xlo=0,Xhi=4

done