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");


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

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

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

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

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

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

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

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