function?chua_memristor
clear?all
figure（1）
%This?program?originally?computed?by?S.?Sabarathinam?
%based?on?‘Memristor?based?Chaotic?Circuits-M.?Bharathwaj‘
%don‘t?forgot?to?rate?my?files...
%To?vary?the?value?of?‘r‘?you?can?get?the?bifurcation?scenario.
???????display（‘S.?SABARATHINAM>>>please?rate?my?files..‘）;
????????????r=2000.0;?c1=6.8e-9;	c2=68.0e-9;	bind=18.0e-3;
????????????alpha?=?-0.667D-03;beta=0.029D-03;?xx=bind*c2;tau=3.499e-05;
??????????cn1=tau/c1*r;cn2=tau/c2*r;?cn3=tau/bind;
?%?????????????????????
????????????tspan?=?0:0.01:500;??
????????????x10?=?0.1;?x20?=?-0.2;?x30?=?-0.003;
????????????x40?=?0.05;?

????????????y0?=?[x10;?x20;?x30?;?x40?];
???????????op=odeset（‘absto