I m not getting the exact bifurcation plot for matlab code. my bifurcation plot for the particular system is not showing the divergence to nonlinearity from linear region. what should be the error in my codes?

% mlcmainfile

 clf

 clc

 global f 

 for f=0:0.01:50

   initialx(1)=0.8;initialx(2)=0.04;

 [t,x]=ode45('mlceq',[0 450],[initialx(1),initialx(2)]);

 plot(f,x(440:450,1),'r');

 hold on

 end

 xlabel('f');ylabel('x')

 axis ([0 50 -1.2 1.2])

--------------------------------------------------------------------------------------------------------------------------------------

 % filename:mlceq(t,x)

 function s=mlceq(t,x)

 w=0.7;

 global f

 xdot=1.65* mlcg(x(1))-0.55*x(1)+x(2);

 ydot=-x(2)-x(1)+f*sin(w*t);

 s=[xdot;ydot];

 % % xdot =1.65*G(x)-0.55*x+y;
 % % ydot =-y-x+fsin(w*t);
-------------------------------------------------------------------
 % Filename: mlcg(x)

 function s = mlcg(x) 

 c=0.52;

 if (x(1))>c

          s=c;

 else

     if (x(1))<-c

           s=-c;

    else
           s=x(1);

     end

 end

 end