Duffing Attractor

The Duffing attractor is a 2-dimensional strange attractor produced by the following equation:

xn+1 =xn+yndt yn+1 =yn+ gcos(wdt) -dyn +axn -bxn3