/* Unprimed + Malaria Model 2 1. Use the population results of Model 1 in Unprimed+MalariaModel2. 2. Duration should be 13 days (312 hours). 3. Use population results of Model 2 in Unprimed+MalariaModel3. 4. Duration should be 3 days (72 hours). Expected Results The TH1 population should rise by about 5.6% to 30%. The TH2 population should show a significant decrease by 76% to 78%. */ recruit_TH1_rate= 0.01; //T helper cell type 1 recruit rate recruit_TH2_rate= 0.01; //T helper cell type 2 recruit rate divide_TH1_rate=0.006; //T helper cell type 1 division rate (Divide every 7 days(168 hours) divide_TH2_rate=0.006; //T helper cell type 2 division rate (Divide every 7 days(168 hours) big_r = 1000000000; //passive rate th_rate = 0.0001; //For TH components on left side of the equation //Naive cells can come into contact with APCs and recruit signals of TH cells. Naive =(recruit_TH1,big_r).NaiveBounceTH1 + (recruit_TH2,big_r).NaiveBounceTH2; NaiveBounceTH1 = (recruit_new_TH1,recruit_TH1_rate).TH1; NaiveBounceTH2 = (recruit_new_TH2,recruit_TH2_rate).TH2; //TH1 and TH2 used here for numbers as their behaviour is carried out by BehaviourTH1/TH2 TH1 =(recruit_and_divide,th_rate).TH1; TH2 =(recruit_and_divide,th_rate).TH2; //T helper cell division and recruitment components BehaviourTH1 = (recruit_TH1,recruit_TH1_rate).BehaviourTH1 + (divide_TH1,divide_TH1_rate).BehaviourTH1; BehaviourTH2 = (recruit_TH2,recruit_TH2_rate).BehaviourTH2 + (divide_TH2,divide_TH2_rate).BehaviourTH2; //Should be equal number of Dormant to Naive + PrecursorDivisionComponent Dormant = (produce_new_TH1,big_r).BehaviourTH1 + (produce_new_TH2,big_r).BehaviourTH2 + (recruit_new_TH1,big_r).BehaviourTH1 + (recruit_new_TH2,big_r).BehaviourTH2; //Components for division PrecursorDivisionComponent = (divide_TH1,big_r).PrecursorDivisionComponentTH1 + (divide_TH2,big_r).PrecursorDivisionComponentTH2; PrecursorDivisionComponentTH1 = (produce_new_TH1,divide_TH1_rate).TH1; PrecursorDivisionComponentTH2 = (produce_new_TH2,divide_TH2_rate).TH2; //Populations taken from Model 1 ((Naive[15000] || PrecursorDivisionComponent[750] || TH1[580] || TH2[580]) (BehaviourTH1[580] || BehaviourTH2[580] || Dormant[15750]))