Updated the paramters to show better behavior.

Plots.m

@@ -4,9 +4,11 @@
 function Plots(T)
 
 if nargin == 0
-    Con     	= [10;		% N_tr
-               	   10;		% N_rt
-                   40];		% N_rr    
+    Con     	= [ 2;      % sigma
+                   5;		% N_tr
+               	   5;		% N_rt
+                   50];		% N_rr 
+               
 
     var_stim    = [ 0;          % strength of the stimulus in Hz (spikes per second)
                     0;          % time between   stimuli in s    

Thalamic_Column.h

@@ -54,7 +54,8 @@ class Thalamic_Column {
 
 	/* Constructor for simulation */
 	Thalamic_Column(double* Con)
-	: N_tr		(Con[0]),	N_rt		(Con[1]),	N_rr		(Con[2])
+	: 	sigma_t		(Con[0]),	sigma_r		(Con[0]),
+		N_tr		(Con[1]),	N_rt		(Con[2]),	N_rr		(Con[3])
 	{set_RNG();}
 
 	/* Set strength of input */
@@ -131,8 +132,8 @@ class Thalamic_Column {
 
 	/* Declaration and Initialization of parameters */
 	/* Membrane time in ms */
-	const double 	tau_t 		= 30;
-	const double 	tau_r 		= 30;
+	const double 	tau_t 		= 20;
+	const double 	tau_r 		= 20;
 
 	/* Maximum firing rate in ms^-1 */
 	const double 	Qt_max		= 400.E-3;
@@ -143,8 +144,8 @@ class Thalamic_Column {
 	const double 	theta_r		= -58.6;
 
 	/* Sigmoid gain in mV */
-	const double 	sigma_t		= 9;
-	const double 	sigma_r		= 9;
+	const double 	sigma_t		= 3;
+	const double 	sigma_r		= 3;
 
 	/* Scaling parameter for sigmoidal mapping (dimensionless) */
 	const double 	C1          = (3.14159265/sqrt(3));
@@ -160,14 +161,14 @@ class Thalamic_Column {
 
 	/* Potassium leak current */
 	const double 	g_LK_t 		= 0.02;
-	const double 	g_LK_r 		= 0.02;
+	const double 	g_LK_r 		= 0.012;
 
 	/* T current */
 	const double	g_T_t		= 3;
 	const double	g_T_r		= 2;
 
 	/* h current */
-	const double	g_h			= 0.04;
+	const double	g_h			= 0.07;
 
 	/* Reversal potentials in mV */
 	/* Synaptic */
@@ -206,9 +207,9 @@ class Thalamic_Column {
 	double			input		= 0.0;
 
 	/* Connectivities (dimensionless) */
-	const double 	N_tr		= 2;
-	const double 	N_rt		= 5.5;
-	const double 	N_rr		= 5;
+	const double 	N_tr		= 10;
+	const double 	N_rt		= 10;
+	const double 	N_rr		= 40;
 
 	friend class 	Stim;
 };

Thalamus.cpp

@@ -38,7 +38,7 @@
 /****************************************************************************************************/
 /*										Fixed simulation settings									*/
 /****************************************************************************************************/
-extern const int onset	= 5;								/* time until data is stored in  s		*/
+extern const int onset	= 15;								/* time until data is stored in  s		*/
 extern const int res 	= 1E4;								/* number of iteration steps per s		*/
 extern const int red 	= res/100;							/* number of iterations that is saved	*/
 extern const double dt 	= 1E3/res;							/* duration of a timestep in ms			*/