• API
• FAQ
• Tools
• Archive
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.

Top