clear?all;
hold?off;

x（1）=0.02;y（1）=0.01;z（1）=0.03;
m0=-1/7;m1=2/7;a=10;b=15;
T=1e-2;
N=5e4;
for?i=1:N
????i
????t（i）=i;
????f（i）=m1*x（i）+0.5*（m0-m1）.*（abs（x（i）+1）-abs（x（i）-1））;
????x（i+1）=x（i）+T*a*y（i）-T*a*f（i）;
????y（i+1）=T*x（i）+（1-T）*y（i）+T*z（i）;
????z（i+1）=-T*b*y（i）+z（i）;
end

n1=N/2;n2=N;
figure（1）
plot（t（n1:n2）x（n1:n2）‘b-‘）;
xla