用matlab写一段完整的Hindmarsh_Rose模型的代码

我可以回答这个问题。以下是一个简单的Hindmarsh_Rose模型的Matlab代码：

clear all;
a=1;
b=3;
c=1;
d=5;
r=0.001;
x=-1.6;
y=1;
z=0;
t=0;
dt=0.01;
tmax=1000;
n=round(tmax/dt);
for i=1:n
x1=x;
y1=y;
z1=z;
x=x+dt*(y-a*x^3+b*x^2-z+r);
y=y+dt*(-x-c*y+d*z);
z=z+dt*(r*(x-y-z));
t=t+dt;
if mod(i,100)==0
plot3(x,y,z,'.','MarkerSize',5);
hold on;
end
end
xlabel('x');
ylabel('y');
zlabel('z');
title('Hindmarsh-Rose Model');
grid on;
hold off;