graphics_toolkit("gnuplot");

CaHVA0       = load ("parts/original_vclamp/CaHVA.dat");
CaHVA1       = load ("parts/neuron_vclamp/CaHVA.dat");

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

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

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

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

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

NaR0      = load ("parts/original_vclamp/NaR.dat");
NaR1      = load ("parts/neuron_vclamp/NaR.dat");

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


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

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

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

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

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

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

subplot(3,3,7);
plot(NaR0(:,1),NaR0(:,2),NaR1(:,1),NaR1(:,2),'linewidth',2);
title ("NaR current");

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

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