Permalink
Browse files
Removed the matlab plot function from the pure code
- Loading branch information...
Showing
with
1 addition
and
26 deletions.
-
+0
−24
Plots.m
-
+1
−2
Thalamic_Column.h
|
|
@@ -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')
|
|
|
@@ -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 */
|
|
|
|
0 comments on commit
45fda04