MadonnaModeløu {Lytic Phage Model - Chemostat culture with Refuge} {CRISPR - CN has phage-induced resistance} METHOD EULER STARTTIME = 0 STOPTIME= 1200 DT = 0.001 DTOUT = 4 init N= 1e8 {Sensitive non-CRSPR} init C =1e2 {Sensitive CRSPR} init CR=1e-4 {RESISTANT CRSPR } init Nx=10{Resistant N-C} init P = 1e4 {Free phage 1} init R=0.065 {Resource} Vn=1 {Max growth sensitive N} vc = .9 {Max growth C} vcr=.9 {Maximum growrth CR} vx=0.7{Maximum growth rate resistant} cc=50 {Input resource} w=0.2 {flow rate} e=5E-7 {Conversion efficiency} k=0.25 {Monod constant} dn=5e-9 {Adsorption rate P on N} dc= 5e-9 {Adsorption rate P on C} dcrmax =5e-9 {Maxumum Adsorption rate P on CN} dcrmin = 1e-14 {Minimum adsorption rate} bn = 50 {Burst P on N} bc = 50 {Burst P on C} bcr =50 {Burst P on CN} srefn= 1e2 {Refuge density of bacteria} srefc = 1e2 srefcr =1e2 m = 1e-4 {CRISPR conversion C--> CR} nm =1e-4 {CR --> C} Rmax=100 MOI = P/CR nn=2 kk=1e2 z=0 dcr=dcrmin +z*(dcrmax*MOI^nn)/(kk+MOI^nn) d/dt (R) = w*(cc-R) - (vn*N +vc*C +vcr*CR+vx*Nx)* (R/(r+k))*e d/dt (N) = (vn*R/(R+k))*N - refn*N*dn*P - w*N d/dt (C) = (vc*R/(R+k))*C - refc*C*dc*P*(1-m) +CR*nm - w*C d/dt (CR) = (vcr*R/(R+k))*CR- refcr*dcr*CR*P + refc*C*dc*P*m -nm*CR- w*CR d/dt (P) = (P* (refn*dn*N*bn+refc*dc*C*bc +refcr*dcr*CR*bcr -dcrmax*CR) - w*P)*PREF d/dt (Nx) =Nx*(vx*R/(R+k)) -w*Nx refn = IF N > srefn THEN 1 ELSE 0 refc= IF C> srefc THEN 1 ELSE 0 refcr= IF CR > srefcr THEN 1 ELSE 0 NT=N+C+CR +Nx PC=(C+CR)/NT Pref= IF P < 1 THEN 0 ELSE 1 ø Genevaþaþ5ö ½Ýéñ  ùøö]µUý@Éùøa@ø eulerööþöúü©ñÒMbP?úü©ñÒMbP?úü©ñÒMbP?öööö $bf^ZH & ,6jö. µy ¹ ‘ñ¡5‰e™Ý¡± !$©&±(¹*Á, .©É