diff --git a/Plots.m b/Plots.m
deleted file mode 100644
index 2118f09..0000000
--- a/Plots.m
+++ /dev/null
@@ -1,24 +0,0 @@
-% mex command is given by: 
-% mex CXXFLAGS="\$CXXFLAGS -std=c++11 -O3" Thalamus_mex.cpp Thalamic_Column.cpp
-
-function Plots(T)
-
-if nargin == 0
-    Con     	= [ 0.01;       % g_h
-                    0.05];      % g_LK
-
-    T       	= 30;  		% duration of the simulation
-end
-[Vt, Vr, ah] = Thalamus_mex(T, Con, var_stim);
-
-L        = max(size(Vt));
-timeaxis = linspace(0,T,L);
-
-% Plot the data
-figure(1)
-subplot(311), plot(timeaxis,Vt)
-title('Thalamic relay membrane voltage'), xlabel('time in s'), ylabel('Vt [mV]')
-subplot(312), plot(timeaxis,Vr)
-title('Thalamic reticular membrane voltage'), xlabel('time in s'), ylabel('Vr [mV]')
-subplot(313), plot(timeaxis,ah)
-title('Thalamic relay I_h activation'), xlabel('time in s'), ylabel('m_h')
diff --git a/Thalamic_Column.h b/Thalamic_Column.h
index fba5c27..52b3830 100644
--- a/Thalamic_Column.h
+++ b/Thalamic_Column.h
@@ -59,8 +59,7 @@ class Thalamic_Column {
 
 	/* Constructor for simulation */
 	Thalamic_Column(double* Par)
-	: 	g_LK		(Par[1]),	g_h 		(Par[0]),
-		N_tr		(Par[2]),	N_rt		(Par[3]),	N_rr		(Par[4])
+	: 	g_LK		(Par[1]),	g_h 		(Par[0])
 	{set_RNG();}
 
 	/* Iterate one time step through SRK4 */