Duffing's equation, x+dxx+x =ycos(wt), exhibits chaotic behavior for its parameters in certain ranges. Simulate this equation for (8, y, w) = (0.05, 0.4, 1.3) and plot the phase portrait (i.e., plot trajectories on the x vs. x plane from several initial conditions). You may use a solver such as ode45 in MATLAB or odeint in SciPy.