# how can I draw bifurcation diagram for given code

38 views (last 30 days)
mallela ankamma rao on 31 Aug 2022
Answered: Abderrahim. B on 31 Aug 2022
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

Abderrahim. B on 31 Aug 2022
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 .
Hope this helps