diff --git a/Cortical_Column.h b/Cortical_Column.h
index 1ea95d0..f4ece2d 100644
--- a/Cortical_Column.h
+++ b/Cortical_Column.h
@@ -151,7 +151,7 @@ class Cortical_Column {
 	const double 	gamma_i		= 58.6E-3;
 
 	/* Axonal flux time constant */
-    const double 	nu			= 60E-3;
+	const double 	nu			= 120E-3;
 
 	/* Conductivities in mS/cm^-2 */
 	/* Leak */
@@ -182,8 +182,8 @@ class Cortical_Column {
 	const double 	N_ei		= 72;
     const double 	N_ie		= 90;
     const double 	N_ii		= 90;
-	const double 	N_te		= 0;
-	const double 	N_ti		= 0;
+	const double 	N_te		= 2.5;
+	const double 	N_ti		= 2.5;
 
 	/* Pointer to thalamic column */
 	Thalamic_Column* Thalamus;
diff --git a/Plots.m b/Plots.m
index 9a113f0..0fedd7e 100644
--- a/Plots.m
+++ b/Plots.m
@@ -19,11 +19,11 @@ function Plots(type)
                            120E-3];     % dphi
                        
     Param_Thalamus      = [0.051;       % g_h
-                           0.02];       % g_LK
+                           0.020];       % g_LK
 end
                         
-Connectivity            = [ 2.5;        % N_et
-                            2.5;        % N_er
+Connectivity            = [ 2.6;        % N_et
+                            2.6;        % N_er
                             5;          % N_te
                             10];        % N_ti   
 
diff --git a/Test_Parameters.m b/Test_Parameters.m
index 14e2b9f..9c2c78b 100644
--- a/Test_Parameters.m
+++ b/Test_Parameters.m
@@ -38,7 +38,7 @@ function Test_Parameters(type)
                            120E-3];     % dphi
                        
     Param_Thalamus      = [0.051;       % g_h
-                           0.02];       % g_LK
+                           0.0205];       % g_LK
                        
     fn_data = '/nfshome/schellen/Documents/MATLAB/TC_model/Data/SO_Average_N3';
     Model_Range_ERP = [-75, -45]; 
@@ -50,7 +50,7 @@ function Test_Parameters(type)
                         
 Connectivity            = [ 2.5;        % N_et
                             2.5;        % N_er
-                            5;          % N_te
+                            15;          % N_te
                             15];        % N_ti   
                         
 load(fn_data);
diff --git a/Test_Stimulation.m b/Test_Stimulation.m
index 4007992..3f3f9d0 100644
--- a/Test_Stimulation.m
+++ b/Test_Stimulation.m
@@ -2,7 +2,7 @@
 
 function Test_Stimulation(type)
 if nargin == 0
-    type = 1;
+    type = 4;
 end
 
 
@@ -35,7 +35,7 @@ function Test_Stimulation(type)
         fn_data         = '/nfshome/schellen/Documents/MATLAB/TC_model/Data/Open_Loop_3Klick.mat';
         Model_xRange    = [-70, -45]; 
         Model_Range_FSP = [-0.25, 1.25]; 
-        Data_xRange     = [-80, 25]; 
+        Data_xRange     = [-75, 25]; 
         Data_Range_FSP  = [-4, 8]; 
         xRange          = -0.5:0.5:4;
     case 2;
@@ -53,7 +53,7 @@ function Test_Stimulation(type)
         Data_Range_FSP  = [-2, 4]; 
         xRange          = -1:0.5:3.5;
     case 4;
-        fn_data         = 'Data/ClosedLoopStim.mat';
+        fn_data         = '/nfshome/schellen/Documents/MATLAB/TC_model/Data/ClosedLoopStim.mat';
         Model_xRange = [-75, -45]; 
         Model_Range_FSP = [-0.25, 1.25]; 
         Data_xRange  = [-80, 50]; 
@@ -69,14 +69,14 @@ function Test_Stimulation(type)
 % 1 == semi-periodic
 % 2 == phase dependend
     
-var_stim    = [ 1;          % mode of stimulation
-                40;         % strength of the stimulus              in Hz (spikes per second)
-                100;       	% duration of the stimulus              in ms
+var_stim    = [ 2;          % mode of stimulation
+                80;         % strength of the stimulus              in Hz (spikes per second)
+                70;       	% duration of the stimulus              in ms
                 5;          % time between stimulation events       in s  (ISI)
                 0;          % range of ISI                          in s  [ISI-range,ISI+range]  
-                3;          % Number of stimuli per event
+                2;          % Number of stimuli per event
                 1050;        % time between stimuli within a event   in ms         
-                450];       % time until stimuli after minimum      in ms
+                400];       % time until stimuli after minimum      in ms
 
 T       	= 300;           % duration of the simulation
 
@@ -88,8 +88,8 @@ function Test_Stimulation(type)
 x_SO        = Marker_Stim;
 
 % Remove those events, that are too close to begin/end
-x_SO        = x_SO(x_SO<(T-2)*Fs); 
-x_SO        = x_SO(x_SO>    2*Fs);
+x_SO        = x_SO(x_SO<(T-xRange(end))*Fs); 
+x_SO        = x_SO(x_SO>  -xRange(1)*Fs);
 
 % Set the variables 
 N_Stim      = length(x_SO);
diff --git a/Thalamic_Column.h b/Thalamic_Column.h
index 9706e77..c600560 100644
--- a/Thalamic_Column.h
+++ b/Thalamic_Column.h
@@ -185,7 +185,7 @@ class Thalamic_Column {
 	const double	g_T_r		= 2.3;
 
 	/* h current */
-	const double	g_h			= 0.07;
+	const double	g_h			= 0.051;
 
 	/* Reversal potentials in mV */
 	/* Synaptic */
@@ -208,7 +208,7 @@ class Thalamic_Column {
 	/* Calcium parameters */
 	const double	alpha_Ca	= -51.8E-6;			/* influx per spike in nmol		*/
 	const double	tau_Ca		= 10;				/* calcium time constant in ms	*/
-	const double	Ca_0		= 2.4E-4;				/* resting concentration 		*/
+	const double	Ca_0		= 2.4E-4;			/* resting concentration 		*/
 
 	/* I_h activation parameters */
 	const double 	k1			= 2.5E7;
@@ -228,8 +228,8 @@ class Thalamic_Column {
 	const double 	N_tr		= 3;
     const double 	N_rt		= 5;
     const double 	N_rr		= 19;
-	const double 	N_et		= 5;
-	const double 	N_er		= 5;
+	const double 	N_et		= 2.5;
+	const double 	N_er		= 2.5;
 
 	/* Pointer to cortical column */
 	Cortical_Column* Cortex;