graphics_toolkit("gnuplot");

CaHVA1       = load ("neuron_vclamp/CaHVA.dat");
CaHVA2       = load ("octave_vclamp/CaHVA.dat");

KA1       = load ("neuron_vclamp/KA.dat");
KA2       = load ("octave_vclamp/KA.dat");

Kir1      = load ("neuron_vclamp/Kir.dat");
Kir2      = load ("octave_vclamp/Kir.dat");

KM1      = load ("neuron_vclamp/KM.dat");
KM2      = load ("octave_vclamp/KM.dat");

KV1      = load ("neuron_vclamp/KV.dat");
KV2      = load ("octave_vclamp/KV.dat");

Na1      = load ("neuron_vclamp/Na.dat");
Na2      = load ("octave_vclamp/Na.dat");

Nar1      = load ("neuron_vclamp/Nar.dat");
Nar2      = load ("octave_vclamp/Nar.dat");

pNa1      = load ("neuron_vclamp/pNa.dat");
pNa2      = load ("octave_vclamp/pNa.dat");


subplot(3,3,1);
plot(CaHVA1(:,1),CaHVA1(:,2),CaHVA2(:,1),CaHVA2(:,2),'linewidth',2);
title ("CaHVA current");

subplot(3,3,2);
plot(KA1(:,1),KA1(:,2),KA2(:,1),KA2(:,2),'linewidth',2);
title ("KA current");

subplot(3,3,3);
plot(Kir1(:,1),Kir1(:,2),Kir2(:,1),Kir2(:,2),'linewidth',2);
title ("Kir current");

subplot(3,3,4);
plot(KM1(:,1),KM1(:,2),KM2(:,1),KM2(:,2),'linewidth',2);
title ("KM current");

subplot(3,3,5);
plot(KV1(:,1),KV1(:,2),KV2(:,1),KV2(:,2),'linewidth',2);
title ("KV current");

subplot(3,3,6);
plot(Na1(:,1),Na1(:,2),Na2(:,1),Na2(:,2),'linewidth',2);
title ("Na current");

subplot(3,3,7);
plot(Nar1(:,1),Nar1(:,2),Nar2(:,1),Nar2(:,2),'linewidth',2);
title ("Nar current");

subplot(3,3,8);
plot(pNa1(:,1),pNa1(:,2),pNa2(:,1),pNa2(:,2),'linewidth',2);
title ("pNa current");

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