FilesExpand file tree

master

/

chap18.py

Copy path

Blame

More file actions

Blame

More file actions

Latest commit

 

History

History

History

30 lines (20 loc) · 658 Bytes

master

/

chap18.py

Copy path

Top

File metadata and controls

• Code
• Blame

30 lines (20 loc) · 658 Bytes

Raw

Copy raw file

Download raw file

Open symbols panel

Edit and raw actions

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

from modsim import *

def make_system(params, data):

G0, k1, k2, k3 = params

t_0 = data.index[0]

t_end = data.index[-1]

Gb = data.glucose[t_0]

Ib = data.insulin[t_0]

I = interpolate(data.insulin)

init = State(G=G0, X=0)

return System(init=init, params=params,

Gb=Gb, Ib=Ib, I=I,

t_0=t_0, t_end=t_end, dt=2)

from modsim import *

def slope_func(t, state, system):

G, X = state

G0, k1, k2, k3 = system.params

I, Ib, Gb = system.I, system.Ib, system.Gb

dGdt = -k1 * (G - Gb) - X*G

dXdt = k3 * (I(t) - Ib) - k2 * X

return dGdt, dXdt