Six compartments: differential equation model



A simple pharmokinetic compartmental model. The solution is a constant.

   model
   {
   solution[1:n.grid, 1:dim] <- ode.solution(init[1:dim], grid[1:n.grid], D(C[1:dim], t), 0, tol)

   D(C[1], t) <- PER1 * C[7] - R * kT1 * C[1]
D(C[2], t) <- PER2 * C[8] - R * kT2 * C[2] - CLR * C[8] / V2
D(C[3], t) <- PER3 * C[8] - R * kT3 * C[3]
D(C[4], t) <- (QHEP * C[8] + R * kT3 * V3 * C[3] -
         R * kT4 * V4 * C[4] * (1 + CLHEP / (R * Q4 - CLHEP))) / V4
D(C[5], t) <- PER5 * C[8] - R * kT5 * C[5]
D(C[6], t) <- PER6 * C[8] - R * kT6 * C[6]
D(C[7], t) <- (R * (kT2 * V2 * C[2] + kT4 * V4 * C[4] + kT5 * V5 * C[5] +
         kT6 * V6 * C[6]) - Q1 * C[7]) / VVEN
D(C[8], t) <- (R * kT1 * V1 * C[1] - Q1 * C[8]) / VART

PER1 <- Q1 / V1
PER2 <- Q2 / V2
PER3 <- Q3 / V3
PER5 <- Q5 / V5
PER6 <- Q6 / V6

Q4 <- Q3 + QHEP
}

Data
list(
n.grid = 73, dim = 8, tol = 1.0E-3,
R = 1, Q1 = 51.2,
Q2 = 14.48, Q3 = 16, QHEP = 4, Q5 = 16.32, Q6 = 0.4,
V1 = 1, V2 = 4.8, V3 = 13, V4 = 10.5, V5 = 180.2, V6 = 10, VVEN = 10, VART = 6.67,
kT1 = 10.119, kT2 = 0.596, kT3 = 0.223, kT4 = 0.134, kT5 = 0.033, kT6 = 0.004,
CLR = 0, CLHEP = 10,
init = c(0, 0, 0, 0, 0, 0, 25, 0),
grid = c(0, 10, 20, 30, 40, 50, 60, 70, 80, 90,
100, 110, 120, 130, 140, 150, 160, 170, 180, 190,
200, 210, 220, 230, 240, 250, 260, 270, 280, 280,
300, 310, 320, 330, 340, 350, 360, 370, 380, 390,
400, 410, 420, 430, 440, 450, 460, 470, 480, 490,
500, 510, 520, 530, 540, 550, 560, 570, 580, 590,
600, 610, 620, 630, 640, 650, 660, 670, 680, 690,
700, 710, 720)
)
( click to open )


Results

      
[sixcomp1]

On this page