How to plot ,bifurcation diagram and FFT of lander pol oscillator?

[Sreeshankar P]

function VanderPol() [t,y] = ode23(@vdp1,[0 20],[2; 0]); title('Solution of van der Pol Equation (\mu = 1) with ODE23'); xlabel('Time t'); ylabel('Solution y'); legend('y_1','y_2') figure plot(y(:,1),y(:,2)) title('Phase plane plot') function dydt = vdp1(t,y) dydt = [y(2); mu*(1-y(1)^2)*y(2)-y(1)+A*sin(omega*t)]; end end idis = 2; ivel = 0; A = 8; omega = 3.9; mu = 4.1; [t,y] = VanderPol(idis, ivel, mu, A, omega); figure plot(t,y) grid

[Sreeshankar P]

This code is not giving the bifurcation diagram for the question asked.

I want to plot bifurcation diagram and perform fft on the oscillator given by the equation in the question,

Also the parameters are given