diff --git a/Figures.m b/Figures.m
new file mode 100644
index 0000000..669a734
--- /dev/null
+++ b/Figures.m
@@ -0,0 +1,59 @@
+% mex command is given by: 
+% mex CXXFLAGS="\$CXXFLAGS -std=c++0x" TC.cpp Cortical_Column.cpp Thalamic_Column.cpp
+
+function Figures(T)
+
+if nargin == 0
+    Input_N3    = [ 8.7;          % sigma_e
+                    2.6;        % alpha_Na
+                    3;          % tau_Na
+                    1.6;        % g_KNa
+                    60E-3];     % dphi
+                        
+                        
+    Input_N2    = [ 4.6;        % sigma_e
+                    2;          % alpha_Na
+                    1.2;        % tau_Na
+                    1.33;       % g_KNa
+                    60E-3];     % dphi
+
+    Connectivity= [2;           % N_et
+               	   2;           % N_er
+                   4;           % N_te
+                   4];          % N_ti    
+
+    % stimulation parameters
+    % first number is the mode of stimulation
+    % 0 == none
+    % 1 == periodic
+    % 2 == phase dependend up state
+    % 3 == phase dependend down state
+    
+    var_stim    = [ 1;          % mode of stimulation
+                    40;          % strength of the stimulus      in Hz (spikes per second)
+                    50;       	% duration of the stimulus      in ms
+                    5;          % time between stimuli          in s    
+                    0];         % time until stimuli after min 	in ms
+
+    T       	= 30;           % duration of the simulation
+end
+
+[Ve_N2, Vt_N2] = TC(T, Input_N2, Connectivity, var_stim);
+
+[Ve_N3, Vt_N3] = TC(T, Input_N3, Connectivity, var_stim);
+
+L        = max(size(Vt));
+timeaxis = linspace(0,T,L);
+
+figure(1)
+subplot(211), plot(timeaxis,Ve_N2)
+title('Pyramidal membrane voltage'), xlabel('time in s'), ylabel('Ve in mV')
+subplot(212), plot(timeaxis,Vt_N2)
+title('Thalamic relay membrane voltage'), xlabel('time in s'), ylabel('Vt in mV')
+
+figure(2)
+subplot(211), plot(timeaxis,Ve_N3)
+title('Pyramidal membrane voltage'), xlabel('time in s'), ylabel('Ve in mV')
+subplot(212), plot(timeaxis,Vt_N3)
+title('Thalamic relay membrane voltage'), xlabel('time in s'), ylabel('Vt in mV')
+%save('TC.mat','Ve','Vt')
diff --git a/Plots.m b/Plots.m
index b20624b..4af50d2 100644
--- a/Plots.m
+++ b/Plots.m
@@ -1,5 +1,5 @@
 % mex command is given by: 
-% mex CXXFLAGS="\$CXXFLAGS -std=gnu++0x -fpermissive" TC.cpp Cortical_Column.cpp Thalamic_Column.cpp
+% mex CXXFLAGS="\$CXXFLAGS -std=c++0x" TC.cpp Cortical_Column.cpp Thalamic_Column.cpp
 
 function Plots(T)
 
@@ -17,10 +17,10 @@ function Plots(T)
                     1.33;       % g_KNa
                     60E-3];     % dphi
 
-    Connectivity= [3;           % N_et
-               	   3;           % N_er
-                   15;           % N_te
-                   20];          % N_ti    
+    Connectivity= [2.5;           % N_et
+               	   2.5;           % N_er
+                   4;           % N_te
+                   4];          % N_ti    
 
     % stimulation parameters
     % first number is the mode of stimulation
@@ -30,8 +30,8 @@ function Plots(T)
     % 3 == phase dependend down state
     
     var_stim    = [ 0;          % mode of stimulation
-                    30;          % strength of the stimulus      in Hz (spikes per second)
-                    50;       	% duration of the stimulus      in ms
+                    40;          % strength of the stimulus      in Hz (spikes per second)
+                    150;       	% duration of the stimulus      in ms
                     5;          % time between stimuli          in s    
                     0];         % time until stimuli after min 	in ms
 
diff --git a/Thalamic_Column.h b/Thalamic_Column.h
index 7458a3e..10b5a32 100644
--- a/Thalamic_Column.h
+++ b/Thalamic_Column.h
@@ -181,7 +181,7 @@ class Thalamic_Column {
 	const double	g_T_r		= 2;
 
 	/* h current */
-	const double	g_h			= 0.07;
+	const double	g_h			= 0.06;
 
 	/* Reversal potentials in mV */
 	/* Synaptic */