Given a system of second order difference equations, how do I use Matlab to calculate numerical simulations to determine the asymptotic behavior of the solution?

Consider the system of second order difference equations n(t+1)=(r*n(t)/1+n(t-1)^5)+0.001. For each of the following values of r, use numerical simulations in Matlab to determine the asymptotic behavior of the solution if n(0)=n(1)=1.

a. r=1.1

b. r=1.3

c. r=1.55

d. r=2.0