Skip to content

Instantly share code, notes, and snippets.

@muscovitebob
muscovitebob / CompartSysFunc.m
Last active March 19, 2018 18:03
MATLAB four compartment pharmacological model with regular mass intake
function s = CompartSysFunc(t, x)
ka = 5;
kc = 2;
kp = 2;
ke = 0.4;
s1 = -ka.*x(1);
s2 = ka.*x(1) -(kc + ke).*x(2) + kp.*x(3);
s3 = kc.*x(2) -kp.*x(3);
s4 = ke.*x(2);
s = transpose([s1 s2 s3 s4]);
@muscovitebob
muscovitebob / keybase.md
Created March 19, 2018 19:27
keybase.md

Keybase proof

I hereby claim:

  • I am muscovitebob on github.
  • I am muscovitebob (https://keybase.io/muscovitebob) on keybase.
  • I have a public key ASBTIVsKlpB4SmpX6ozRrtWfcjgASZKYYCTir-hVmy2hmgo

To claim this, I am signing this object:

@muscovitebob
muscovitebob / Goldbeter95Model.py
Last active May 17, 2018 11:51
Reproducing results from Goldbeter, Albert. "A model for circadian oscillations in the Drosophila period protein (PER)." Proc. R. Soc. Lond. B 261.1362 (1995): 319-324.
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint
def PERModel(y, t, vd = 0.95):
#here y0 is Pn , y1 is M, y2 is P0, y3 is P1, y4 is P2
vs, vm = 0.76, 0.65
Km, Ki, Kd, = 0.5, 1, 0.2
ks = 0.38