Can any one help me to plot bifurcation diagram for the system function

 function xdot = fzr(t, x)
r=5;k=830;b=0.011;e=0.1;c=0.4;d1=0.5;d=0.5;g=0.0092;h=0.09;v=1.58;m=0;v1=0;
xdot = zeros(3,1);    
xdot(1)= r*x(1)*(1-(x(1)+x(2))/k)-c*x(1)*x(2)-b*x(1)*x(3)+m+v1*x(2);
xdot(2)= c*x(1)*x(2)-e*x(2)*x(3)-d1*x(2)-v1*x(2);
xdot(3)= g*x(1)*x(3)+h*x(2)*x(3)-d*x(3)-o*((x(1)+x(2))/(v+x(1)+x(2)))*x(3); 
with o as paprameter