how can I draw bifurcation diagram for given code
22 views (last 30 days)
Show older comments
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
0 Comments
Answers (1)
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
0 Comments
See Also
Categories
Find more on Get Started with MATLAB in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!