dml = load ("Izhikevich03_ML.data");
doct = load ("Izhikevich03_Octave.data");

subplot(3,2,1);
plot (dml(1:24000,1),dml(1:24000,2),doct(1:2400,1),doct(1:2400,2));
title("regular spiking (RS)");
subplot(3,2,2);
plot (dml(1*24000+8:2*24000,1),dml(1*24000+8:2*24000,2),doct(1*2400+109:2*2400,1),doct(1*2400+109:2*2400,2));
title("intrinsically bursting (IB)");
subplot(3,2,3);
plot (dml(2*24000+24:3*24000,1),dml(2*24000+24:3*24000,2),doct(2*2400+221:3*2400+200,1),doct(2*2400+221:3*2400+200,2));
title("chattering (CH)");
subplot(3,2,4);
plot (dml(3*24000+48:4*24000,1),dml(3*24000+48:4*24000,2),doct(3*2400+349:4*2400+300,1),doct(3*2400+349:4*2400+300,2));
title("fast spiking (FS)");
subplot(3,2,5);
plot (dml(4*24000+80:5*24000,1),dml(4*24000+80:5*24000,2),doct(4*2400+484:5*2400+400,1),doct(4*2400+484:5*2400+400,2));
title("low-threshold spiking (LTS)");

set (findall (gcf, "-property", "fontsize"), "fontsize", 12);
r = 3;
c = 3;
dx = 1/(10*r);
dy = 1/(10*c);
hax = findobj (gcf, "type", "axes");
for n = 1:numel(hax)
  set (hax(n), "position", get (hax(n), "position") + [dx,dy,-2*dx,-2*dy]);
endfor