octave_v = load ("octave_run/v.dat");
neuron_v = load ("neuron_run/v.dat");

subplot(2,1,1);
h=plot (octave_v(1:1.01e3,1),octave_v(1:1.01e3,2),"1");
set(h(1),"linewidth",2);
title ("Octave simulation");
ylabel("Membrane potential [mV]");


% add an offset of 12ms to the NEURON results
subplot(2,1,2);
h=plot(neuron_v(1:2.46e4,1),neuron_v(1:2.46e4,2),"3");
set(h(1),"linewidth",2);
title ("NEURON simulation");
xlabel("Time [ms]");
ylabel("Membrane potential [mV]");