n_t=maple('rsolve({n(t+1)-n(t)=(1300-16*n(t))/10000,n(0)=0.13},n(t))')
n_t =
-2028/25*(624/625)^t+325/4