graphics_toolkit("gnuplot");

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

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

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

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

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

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

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

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


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

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

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

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

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

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

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

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

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