a,E0,R,c,t = var('a,E0,R,c,t')      #Declaring all as a variable
assume(a>0)                     #We need to specify that a is non-negative
I_t = E0/R + c*exp(-a*t)       #Equation for I_t
limit(I_t,t=Infinity)