how can I draw bifurcation diagram for given code

Good afternoon sir

can you tell me anyone how to draw bifurcation diagram like this

N=1390000000;
zetaa = 0.4;
zetai = 0.3;
zetag = 0.3;
theta = 0.2;
omega = 0.5;
beta = 0.3311;
lamdai = 0.0046;
lamdau = 0.5129;
gammai = 0.0046;
gammau = 0.5129;
gammag = 0.1318;
mui = 0.0046;
muu = 0.4852;
mug = 0.0128;
mu = 0.0014;
pi=2000
Ki = lamdai+gammai+mui+mu; Ku = lamdau+gammau+muu+mu; Kg = gammag+mug+mu; Km = omega+mu;

R0 = (beta*zetag/Km)+(beta*zetai*theta*omega/(Ki*Km))+...
    ((beta*zetag*lamdai*theta*omega*Ku+beta*zetag*lamdau*(1-theta)*omega*Ki)/(Ki*Kg*Ku*Km));


A1 =(pi-mu*N/R0)/Km;  
I1 =(theta*omega)/Ki; 
U1=((1-theta)*omega*A1)/Ku;
G1 =(lamdai*theta*omega*Ku+lamdau*(1-theta)*omega*Ki)*A1/(Ki*Kg*Ku);
R1 =(((lamdai*theta*omega)/mu*Ki)+((1-theta)*omega/Ku*mu)+((gammag*lamdai*theta*omega)/(mu*Kg*Ki))...
    +(gammag*lamdau*(1-theta)*omega)/(mu*Kg*Ku))*A1;
V1 = beta*(zetaa*A1+zetai*I1+zetag*G1)*(1/N);          
plot(R0,V1)
xlabel('R0');
ylabel('V1');
box on
hold on

Answer 1

[Abderrahim. B]

Hi!

MATLAB does not have a function to plot this kind of diagrams, but it has capabilities that you can use to write your own specific codes.

Consider using function the below file echanges that you may find useful .

1. Plotting bifurcation diagram.
2. Chaotic generators demo

Hope this helps