How to change function to get correct bifurcation behavior?
1 次查看(过去 30 天)
显示 更早的评论
I have a system of ODE's that was modelled and now I am trying to vary the parameters for an ODE of interest (dDSB). The equation is: dDSB:=alpha/H1*dD-k1*c0*(1-h1) where: dD=.001, H1=3000, k1=100, c0=25/3000, h1=(25+2)/3000. Solving for alpha I get 0.0000003333333333*alpha - 991/1200. I then tried to graph the plot:
clf; hold on;
for alpha=0.1:.1:25
Yfix=roots([3.3*10^7*alpha -991/1200]);
I=imag(Yfix)==0;
Yfix=Yfix(I);
Lfix=2*Yfix;
plot(repmat(alpha,size(Yfix)), Yfix, 'o', 'LineWidth',5);
end
hold off
How can I get the graph to bifurcate?
0 个评论
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Ordinary Differential Equations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!