how to plot bifurcation diagram in matlab

a=7.2*10^-6;

Gm=2.5*10^-2

Sp=2.75*10^-5;

yon=6*10^-2;

b=4.7;

son=4.5*10^-1;

B=10^-4;

vm=0:0.5:2.5;

gp=sinh(vm);

A=10^-10;

soff=1.3*10^-2;

beta=500;

yoff=1.3*10^-2

v=sawtooth(2*pi*1*t)

i=v.*(y*Gm+(1-y)*a.*exp(b*sqrt(abs(v))))

p=i.*v

dydt = B*sinh(v/son).*exp(-y.^2/yon^2).*exp(p/Sp)*heaviside(v)+...

A*sinh(v/soff).*exp(-yoff^2./y.^2).*exp(1/(1+beta*p))*heaviside(-v)

i want to bifurcation diagram for i,v

Answer 1

[Pranav Verma]

Hi Nune,

Please see the below answers on the same lines which discuss the steps to plot the bifurcation diagrams:

https://www.mathworks.com/matlabcentral/answers/319767-how-to-plot-a-bifurcation-diagram-for-differential-equation-is-there-any-formula-to-plot-the-bifurc

https://www.mathworks.com/matlabcentral/answers/20382-bifurcation-diagram

Thanks