SHARE
TWEET

Untitled

a guest Dec 15th, 2018 81 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. function dxdt = sys(t,X,Q)
  2.     [be, delta, d1, k2, lambda1, Kb] = Q
  3.     d2 = 0.01
  4.     k1 = 8e-7
  5.     lambda2 = 31.98
  6.     eps = .1
  7.     de = 0.25
  8.     Nt = 100
  9.     f = 0.34
  10.     m1 = m2 = 1e-5
  11.     lambda_e = 1
  12.     delta_e = 0.1
  13.     p1 = p2 = 1
  14.     Kd = 1
  15.     T1 = X(1)
  16.     T2 = X(2)
  17.     T1s = X(3)
  18.     T2s = X(4)
  19.     V = X(5)
  20.     E = X(6)
  21.    
  22.     dT1 = lambda1 - d1*T1 - (1-eps)*k1*V*T1
  23.     dT2 = lambda2 - d2*T2 - (1-f*eps)*k2*V*T2
  24.     dT1s = (1-eps)*k1*V*T1 - delta*T1s - m1*E*T1s
  25.     dT2s = (1 - f*eps)*k2*V*T2 - delta*T2s - m2*E*T2s
  26.     dV = Nt*delta*(T1s + T2s) - c*V -((1-eps)*p1*k1*T1 + (1-f*eps)*p2*k2*T2)*V
  27.     dE = lambda_e + E*(be*(T1s + T2s))/(T1s + T2s + Kb) - E*(de*(T1s + T2s))/(T1s + T2s + Kd) - delta_e*E
  28.    
  29.     dxdt = [dT1,dT2,dT1s,dT2s,dV,dE]
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top