-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathZerlaut_adap.m
executable file
·428 lines (370 loc) · 46.2 KB
/
Zerlaut_adap.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
% Copyright 2021 Aix-Marseille Université
% "Licensed to the Apache Software Foundation (ASF) under one or more contributor license agreements; and to You under the Apache License, Version 2.0. "
function out = Zerlau_adap
out{1} = @init;
out{2} = @fun_eval;
out{3} = [];
out{4} = [];
out{5} = [];
out{6} = [];
out{7} = [];
out{8} = [];
out{9} = [];
out{10} = @Boundary_up_v_e;
out{11} = @Boundary_low_v_e;
out{12} = @Boundary_up_v_i;
out{13} = @Boundary_low_v_i;
out{14} = @Boundary_up_w_e;
out{15} = @Boundary_low_w_e;
out{16} = @Boundary_up_w_i;
out{17} = @Boundary_low_w_i;
out{18} = @Boundary_up_g_L;
out{19} = @Boundary_low_g_L;
out{20} = @Boundary_up_E_L_e;
out{21} = @Boundary_low_E_L_e;
out{22} = @Boundary_up_E_L_i;
out{23} = @Boundary_low_E_L_i;
out{24} = @Boundary_up_C_m;
out{25} = @Boundary_low_C_m;
out{26} = @Boundary_up_E_e;
out{27} = @Boundary_low_E_e;
out{28} = @Boundary_up_E_i;
out{29} = @Boundary_low_E_i;
out{30} = @Boundary_up_Q_e;
out{31} = @Boundary_low_Q_e;
out{32} = @Boundary_up_Q_i;
out{33} = @Boundary_low_Q_i;
out{34} = @Boundary_up_tau_e;
out{35} = @Boundary_low_tau_e;
out{36} = @Boundary_up_tau_i;
out{37} = @Boundary_low_tau_i;
out{38} = @Boundary_up_N_tot;
out{39} = @Boundary_low_N_tot;
out{40} = @Boundary_up_p_connect;
out{41} = @Boundary_low_p_connect;
out{42} = @Boundary_up_g;
out{43} = @Boundary_low_g;
out{44} = @Boundary_up_T;
out{45} = @Boundary_low_T;
out{46} = @Boundary_up_external_input_E_E;
out{47} = @Boundary_low_external_input_E_E;
out{48} = @Boundary_up_external_input_E_I;
out{49} = @Boundary_low_external_input_E_I;
out{50} = @Boundary_up_external_input_I_E;
out{51} = @Boundary_low_external_input_I_E;
out{52} = @Boundary_up_external_input_I_I;
out{53} = @Boundary_low_external_input_I_I;
out{54} = @Boundary_up_P_e_0;
out{55} = @Boundary_low_P_e_0;
out{56} = @Boundary_up_P_e_1;
out{57} = @Boundary_low_P_e_1;
out{58} = @Boundary_up_P_e_2;
out{59} = @Boundary_low_P_e_2;
out{60} = @Boundary_up_P_e_3;
out{61} = @Boundary_low_P_e_3;
out{62} = @Boundary_up_P_e_4;
out{63} = @Boundary_low_P_e_4;
out{64} = @Boundary_up_P_e_5;
out{65} = @Boundary_low_P_e_5;
out{66} = @Boundary_up_P_e_6;
out{67} = @Boundary_low_P_e_6;
out{68} = @Boundary_up_P_e_7;
out{69} = @Boundary_low_P_e_7;
out{70} = @Boundary_up_P_e_8;
out{71} = @Boundary_low_P_e_8;
out{72} = @Boundary_up_P_e_9;
out{73} = @Boundary_low_P_e_9;
out{74} = @Boundary_up_P_i_0;
out{75} = @Boundary_low_P_i_0;
out{76} = @Boundary_up_P_i_1;
out{77} = @Boundary_low_P_i_1;
out{78} = @Boundary_up_P_i_2;
out{79} = @Boundary_low_P_i_2;
out{80} = @Boundary_up_P_i_3;
out{81} = @Boundary_low_P_i_3;
out{82} = @Boundary_up_P_i_4;
out{83} = @Boundary_low_P_i_4;
out{84} = @Boundary_up_P_i_5;
out{85} = @Boundary_low_P_i_5;
out{86} = @Boundary_up_P_i_6;
out{87} = @Boundary_low_P_i_6;
out{88} = @Boundary_up_P_i_7;
out{89} = @Boundary_low_P_i_7;
out{90} = @Boundary_up_P_i_8;
out{91} = @Boundary_low_P_i_8;
out{92} = @Boundary_up_P_i_9;
out{93} = @Boundary_low_P_i_9;
% --------------------------------------------------------------------------
function dydt = fun_eval(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
% fprintf("%f, %f, %f, %f, %f\n", kmrgd(1),kmrgd(2),kmrgd(3),kmrgd(4),kmrgd(5));
if ~isreal(kmrgd)
warning('not real argument')
[v_e,v_i,c_ee,c_ei,c_ii,w_e]=Zerlaut_model_adap(abs(kmrgd(1)),abs(kmrgd(2)),abs(kmrgd(3)),abs(kmrgd(4)),abs(kmrgd(5)),abs(kmrgd(6)),par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9);
else
[v_e,v_i,c_ee,c_ei,c_ii,w_e]=Zerlaut_model_adap(kmrgd(1),kmrgd(2),kmrgd(3),kmrgd(4),kmrgd(5),kmrgd(6),par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9);
end
dydt=[v_e;v_i;c_ee;c_ei;c_ii;w_e];
% --------------------------------------------------------------------------
function [tspan,y0,options] = init
handles = feval(Zerlau_1);
y0=[0,0,0,0,0,0];
options = odeset('Jacobian',[],'JacobianP',[],'Hessians',[],'HessiansP',[]);
tspan = [0 10];
% --------------------------------------------------------------------------
function jac = jacobian(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
% --------------------------------------------------------------------------
function jacp = jacobianp(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
% --------------------------------------------------------------------------
function hess = hessians(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
% --------------------------------------------------------------------------
function hessp = hessiansp(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
%---------------------------------------------------------------------------
function tens3 = der3(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
%---------------------------------------------------------------------------
function tens4 = der4(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
%---------------------------------------------------------------------------
function tens5 = der5(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
%---------------------------------------------------------------------------
function userfun1=Boundary_up_v_e(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun1=kmrgd(1)-1.0;
function userfun2=Boundary_low_v_e(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun2=kmrgd(1)+1.e-7;
function userfun3=Boundary_up_v_i(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun3=kmrgd(2)-1.0;
function userfun4=Boundary_low_v_i(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun4=kmrgd(2)+1.0e-7;
%%
%%%%%%W_e
function userfun5=Boundary_up_w_e(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
% userfun5=kmrgd(6)+1.0e-7;
userfun5=kmrgd(6)+10000.0*1e-3;
function userfun6=Boundary_low_w_e(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun6=kmrgd(6)-10000.0*1e-3;
%%
%%%%%%G_L
function userfun7=Boundary_up_g_L(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun7=par_g_L-100.0;
function userfun8=Boundary_low_g_L(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun8=par_g_L-0.1;
%%
%%%%%%E_L_e
function userfun9=Boundary_up_E_L_e(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun9=par_E_L_e+50.0;
function userfun10=Boundary_low_E_L_e(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun10=par_E_L_e+90.0;
%%
%%%%%%E_L_i
function userfun11=Boundary_up_E_L_i(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun11=par_E_L_i+50.0;
function userfun12=Boundary_low_E_L_i(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun12=par_E_L_i+90.0;
%%
%%%%%%C_m
function userfun13=Boundary_up_C_m(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun13=par_C_m-500.0;
function userfun14=Boundary_low_C_m(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun14=par_C_m-10.0;
%%
%%%%%%E_e
function userfun15=Boundary_up_E_e(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun15=par_E_e-20.0;
function userfun16=Boundary_low_E_e(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun16=par_E_e+20.0;
%%
%%%%%%E_i
function userfun17=Boundary_up_E_i(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun17=par_E_i+60.0;
function userfun18=Boundary_low_E_i(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun18=par_E_i+100.0;
%%
%%%%%%Q_e
function userfun19=Boundary_up_Q_e(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun19=par_Q_e-5.0;
function userfun20=Boundary_low_Q_e(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun20=par_Q_e-0.0;
%%
%%%%%%Q_i
function userfun21=Boundary_up_Q_i(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun21=par_Q_i-10.0;
function userfun22=Boundary_low_Q_i(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun22=par_Q_i-0.0;
%%
%%%%%%tau_e
function userfun23=Boundary_up_tau_e(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun23=par_tau_e-10.0;
function userfun24=Boundary_low_tau_e(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun24=par_tau_e-1.0;
%%
%%%%%%tau_i
function userfun25=Boundary_up_tau_i(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun25=par_tau_i-10.0;
function userfun26=Boundary_low_tau_i(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun26=par_tau_i-0.5;
%%
%%%%%%b_e
function userfun27=Boundary_up_tau_w_e(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun27=par_tau_w_e-1.0;
function userfun28=Boundary_low_tau_w_e(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun28=par_tau_w_e-1000.0;
%%
%%%%%%b_e
function userfun29=Boundary_up_b_e(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun29=par_b_e-0.0;
function userfun30=Boundary_low_b_e(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun30=par_b_e-150.0;
%%
%%%%%%N_tot
function userfun31=Boundary_up_N_tot(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun31=par_N_tot-50000;
function userfun32=Boundary_low_N_tot(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun32=par_N_tot-1000;
%%
%%%%%%p_connect
function userfun33=Boundary_up_p_connect(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun33=par_p_connect-1.0;
function userfun34=Boundary_low_p_connect(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun34=par_p_connect-0.001;
%%
%%%%%%g
function userfun35=Boundary_up_g(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun35=par_g-0.6;
function userfun36=Boundary_low_g(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun36=par_g-0.001;
%%
%%%%%%T
function userfun37=Boundary_up_T(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun37=par_T-200.0;
function userfun38=Boundary_low_T(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun38=par_T-1.0;
%%
%%%%%% external input_E_E
function userfun39=Boundary_up_external_input_E_E(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun39=par_external_input_E_E-0.2;
function userfun40=Boundary_low_external_input_E_E(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun40=par_external_input_E_E+0.2;
%%
%%%%%% external input_E_I
function userfun41=Boundary_up_external_input_E_I(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun41=par_external_input_E_I-0.1;
function userfun42=Boundary_low_external_input_E_I(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun42=par_external_input_E_I+0.1;
%%
%%%%%% P_e_0
function userfun45=Boundary_up_P_e_0(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun45=par_P_e_0-5.0;
function userfun46=Boundary_low_P_e_0(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun46=par_P_e_0+5.0;
%%
%%%%%% P_e_1
function userfun47=Boundary_up_P_e_1(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun47=par_P_e_1-5.0;
function userfun48=Boundary_low_P_e_1(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun48=par_P_e_1+5.0;
%%
%%%%%% P_e_2
function userfun49=Boundary_up_P_e_2(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun49=par_P_e_2-5.0;
function userfun50=Boundary_low_P_e_2(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun50=par_P_e_2+5.0;
%%
%%%%%% P_e_3
function userfun51=Boundary_up_P_e_3(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun51=par_P_e_3-5.0;
function userfun52=Boundary_low_P_e_3(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun52=par_P_e_3+5.0;
%%
%%%%%% P_e_4
function userfun53=Boundary_up_P_e_4(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun53=par_P_e_4-5.0;
function userfun54=Boundary_low_P_e_4(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun54=par_P_e_4+5.0;
%%
%%%%%% P_e_5
function userfun55=Boundary_up_P_e_5(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun55=par_P_e_5-5.0;
function userfun56=Boundary_low_P_e_5(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun56=par_P_e_5+5.0;
%%
%%%%%% P_e_6
function userfun57=Boundary_up_P_e_6(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun57=par_P_e_6-5.0;
function userfun58=Boundary_low_P_e_6(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun58=par_P_e_6+5.0;
%%
%%%%%% P_e_7
function userfun59=Boundary_up_P_e_7(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun59=par_P_e_7-5.0;
function userfun60=Boundary_low_P_e_7(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun60=par_P_e_7+5.0;
%%
%%%%%% P_e_8
function userfun61=Boundary_up_P_e_8(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun61=par_P_e_8-5.0;
function userfun62=Boundary_low_P_e_8(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun62=par_P_e_8+5.0;
%%
%%%%%% P_e_9
function userfun63=Boundary_up_P_e_9(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun63=par_P_e_9-5.0;
function userfun64=Boundary_low_P_e_9(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun64=par_P_e_9+5.0;
%%
%%%%%% P_i_0
function userfun65=Boundary_up_P_i_0(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun65=par_P_i_0-5.0;
function userfun66=Boundary_low_P_i_0(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun66=par_P_i_0+5.0;
%%
%%%%%% P_i_1
function userfun67=Boundary_up_P_i_1(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun67=par_P_i_1-5.0;
function userfun68=Boundary_low_P_i_1(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun68=par_P_i_1+5.0;
%%
%%%%%% P_i_2
function userfun69=Boundary_up_P_i_2(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun69=par_P_i_2-5.0;
function userfun70=Boundary_low_P_i_2(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun70=par_P_i_2+5.0;
%%
%%%%%% P_i_3
function userfun71=Boundary_up_P_i_3(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun71=par_P_i_3-5.0;
function userfun72=Boundary_low_P_i_3(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun72=par_P_i_3+5.0;
%%
%%%%%% P_i_4
function userfun73=Boundary_up_P_i_4(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun73=par_P_i_4-5.0;
function userfun74=Boundary_low_P_i_4(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun74=par_P_i_4+5.0;
%%
%%%%%% P_i_5
function userfun75=Boundary_up_P_i_5(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun75=par_P_i_5-5.0;
function userfun76=Boundary_low_P_i_5(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun76=par_P_i_5+5.0;
%%
%%%%%% P_i_6
function userfun77=Boundary_up_P_i_6(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun77=par_P_i_6-5.0;
function userfun78=Boundary_low_P_i_6(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun78=par_P_i_6+5.0;
%%
%%%%%% P_i_7
function userfun79=Boundary_up_P_i_7(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun79=par_P_i_7-5.0;
function userfun80=Boundary_low_P_i_7(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun80=par_P_i_7+5.0;
%%
%%%%%% P_i_8
function userfun81=Boundary_up_P_i_8(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun81=par_P_i_8-5.0;
function userfun82=Boundary_low_P_i_8(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun82=par_P_i_8+5.0;
%%
%%%%%% P_i_9
function userfun83=Boundary_up_P_i_9(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun83=par_P_i_9-5.0;
function userfun84=Boundary_low_P_i_9(t,kmrgd,par_g_L,par_E_L_e,par_E_L_i,par_C_m,par_E_e,par_E_i,par_Q_e,par_Q_i,par_tau_e,par_tau_i,par_tau_w_e,par_b_e,par_N_tot,par_p_connect,par_g,par_T,par_external_input_E_E,par_external_input_E_I,par_P_e_0,par_P_e_1,par_P_e_2,par_P_e_3,par_P_e_4,par_P_e_5,par_P_e_6,par_P_e_7,par_P_e_8,par_P_e_9,par_P_i_0,par_P_i_1,par_P_i_2,par_P_i_3,par_P_i_4,par_P_i_5,par_P_i_6,par_P_i_7,par_P_i_8,par_P_i_9)
userfun84=par_P_i_9+5.0;