reservoirs in series and Jordan’s form

reservoirs in series and Jordan’s form

# Identical reservoirs in series

# example of Jordan form

# x[i]: specie in reservoir i [gr]

# V= volume [lt]

# b= volume flux [lt/seg]

# sigma = deposition coeff

# f = entry concentration [gr/lt]

# First reservoir...

x1'= -sigma*x1-b*x1/V+f*(1+cos(omega*t))*b

init x[1]=1

# and the other reservoirs...

%[2..4]

x[j]' = b*x[j-1]/V-sigma*x[j]-b*x[j]/V

init x[j]=0

# Parameters

param V=1, sigma=0.1

param b=1, f=0, omega=2

@ TOTAL=10

done

more about Jordan’s form:

more about compartimental modelling: