neuron_K1    = load ("neuron_vclamp/K1.dat");
neuron_CaBK  = load ("neuron_vclamp/CaBK.dat");
neuron_CaP   = load ("neuron_vclamp/CaP.dat");
neuron_Narsg = load ("neuron_vclamp/Narsg.dat");

octave_K1    = load ("octave_vclamp/K1.dat");
octave_CaBK  = load ("octave_vclamp/CaBK.dat");
octave_CaP   = load ("octave_vclamp/CaP.dat");
octave_Narsg = load ("octave_vclamp/Narsg.dat");


f=subplot(2,2,1);
h=plot(neuron_K1(:,1),neuron_K1(:,2),octave_K1(:,1),octave_K1(:,2),"r");
set(h(1),"linewidth",2);
legend("NEURON","Octave");
title ("K fast current");
ylabel("Current [nA]");
set(f,"ylim",[0 0.20]);

f=subplot(2,2,2);
h=plot(neuron_CaBK(:,1),neuron_CaBK(:,2),octave_CaBK(:,1),octave_CaBK(:,2),"r");
set(h(1),"linewidth",2);
title ("CaBK current");
set(f,"ylim",[0 0.20]);

f=subplot(2,2,3);
h=plot(neuron_CaP(:,1),neuron_CaP(:,2),octave_CaP(:,1),octave_CaP(:,2),"r");
set(h(1),"linewidth",2);
title ("CaP current");
xlabel("Time [ms]");
ylabel("Current [nA]");
set(f,"ylim",[-0.0002 0 ]);

g=subplot(2,2,4);
h=plot(neuron_Narsg(:,1),neuron_Narsg(:,2),octave_Narsg(:,1),octave_Narsg(:,2),"r");
set(h(1),"linewidth",2);
title ("Narsg current");
xlabel("Time [ms]");

print ("NEURON_and_Octave_vclamp.tex", "-depslatexstandalone");