In[]:=
ClearAll["`*"];Remove["`*"]
In[]:=
b=0.03408;d=0.00061;d1=0.0046;ω=0.0095;α=0.0001;γ=0.1896;r=0.05;a=0.46885;
In[]:=
v1=rS[t];
In[]:=
eqns={CaputoD[S[t],{t,0.95}]==b-α(1-a)S[t]iI[t]-αaS[t]iI[t]-av1-dS[t],CaputoD[iI[t],{t,0.95}]==α(1-a)S[t]iI[t]+αaS[t]iI[t]-γiI[t]-(d+d1)iI[t],CaputoD[R[t],{t,0/95}]==av1+γiI[t]-dR[t],S[0]==10000,iI[0]==1,R[0]==0};
In[]:=
sol=DSolve[eqns,{S[t],iI[t],R[t]},{t,0,12}]
Out[]=
DSolve[{CaputoD[S[t],{t,0.95}]0.03408-0.0240525S[t]-0.0001iI[t]S[t],CaputoD[iI[t],{t,0.95}]-0.19481iI[t]+0.0001iI[t]S[t],R[t]0.1896iI[t]-0.00061R[t]+0.0234425S[t],S[0]10000,iI[0]1,R[0]0},{S[t],iI[t],R[t]},{t,0,12}]
Plot[Evaluate[{S[t],iI[t],R[t]}/.sol],{t,0,15}]