Question: Hi currently i do my master thesis and i need help for code fractional system and solve it by homotopy perturbation method

D^γ S(t)=(b-ev)S-α_1 (1-ev)SI_1-α_2 (1-ev)SI_2-dS,
D^γ I_1 (t)=α_1 (1-ev)SI_1-(β_1+d_1+d)I_1,                                           
D^γ I_2 (t)=α_2 (1-ev)SI_2-(β_2+d_2+d) I_2                  , 
D^γ R(t)=β_1 I_1+β_2 I_2+evS-dR 
Please Wait...