neuron_K1    = load ("neuron_vclamp/K1.dat");
neuron_K2    = load ("neuron_vclamp/K2.dat");
neuron_K3    = load ("neuron_vclamp/K3.dat");
neuron_CaBK  = load ("neuron_vclamp/CaBK.dat");
neuron_CaP   = load ("neuron_vclamp/CaP.dat");
neuron_Ih    = load ("neuron_vclamp/Ih.dat");
neuron_Leak  = load ("neuron_vclamp/Leak.dat");
neuron_Narsg = load ("neuron_vclamp/Narsg.dat");

octave_K1    = load ("octave_vclamp/K1.dat");
octave_K2    = load ("octave_vclamp/K2.dat");
octave_K3    = load ("octave_vclamp/K3.dat");
octave_CaBK  = load ("octave_vclamp/CaBK.dat");
octave_CaP   = load ("octave_vclamp/CaP.dat");
octave_Ih    = load ("octave_vclamp/Ih.dat");
octave_Leak  = load ("octave_vclamp/Leak.dat");
octave_Narsg = load ("octave_vclamp/Narsg.dat");


subplot(3,3,1);
plot(neuron_K1(:,1),neuron_K1(:,2),octave_K1(:,1),octave_K1(:,2));
legend("NEURON","Octave");
title ("K fast current");

subplot(3,3,2);
plot(neuron_K2(:,1),neuron_K2(:,2),octave_K2(:,1),octave_K2(:,2));
legend("NEURON","Octave");
title ("K mid current");

subplot(3,3,3);
plot(neuron_K3(:,1),neuron_K3(:,2),octave_K3(:,1),octave_K3(:,2));
legend("NEURON","Octave");
title ("K slow current");

subplot(3,3,4);
plot(neuron_CaBK(:,1),neuron_CaBK(:,2),octave_CaBK(:,1),octave_CaBK(:,2));
legend("NEURON","Octave");
title ("CaBK current");

subplot(3,3,5);
plot(neuron_CaP(:,1),neuron_CaP(:,2),octave_CaP(:,1),octave_CaP(:,2));
legend("NEURON","Octave");
title ("CaP current");

subplot(3,3,6);
plot(neuron_Ih(:,1),neuron_Ih(:,2),octave_Ih(:,1),octave_Ih(:,2));
legend("NEURON","Octave");
title ("Ih current");

subplot(3,3,7);
plot(neuron_Leak(:,1),neuron_Leak(:,2),octave_Leak(:,1),octave_Leak(:,2));
legend("NEURON","Octave");
title ("Leak current");

subplot(3,3,8);
plot(neuron_Narsg(:,1),neuron_Narsg(:,2),octave_Narsg(:,1),octave_Narsg(:,2));
legend("NEURON","Octave");
title ("Narsg current");

print ("Khaliq03vclamp.eps", "-depsc");

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