% Square wave bursting

close all;
clear all;
x0 = [-0.3 0 -0.06];
[t, phi] = ode45('morrisLecar',[0 350],x0);

plot3(phi(:,1),phi(:,3),phi(:,2));
xlabel('V');
ylabel('u');
zlabel('w');
figure
plot(t,phi(:,1));
xlabel('t');
ylabel('V');