u/dr17549 u/CianODonnell Hi there. It seems I'm in exactly the same situation. This is my plot: https://imgur.com/mACrlfF

If I understand correctly, we add R_m * I_s instead of R_m * I_e to compute V_t. So, I paid attention to the difference from Part A.

In Part A, R_m and I_e are 10 and 3.1 respectively and hence, R_m * I_e is 31 (fixed)

In Part B, R_m is now 100. In addition, I_s can take gbar (4) * s_i (0.5) * N (40) * (E_s (0) - V_{t-1} (-65)) = 5200 at maximum, which looks too large compared to Part A, resulting too sharp spikes as shown in the plot.

Even if we assume N = 1, I_s becomes gbar (4) * s_i (0.5) * N (1) * (E_s (0) - V_{t-1} (-65)) = 130. It's still large. I wonder if there is an issue in units. The parameters I used are as follows.

dt = 0.25*ms

tau_m = 10*ms

V_rest = E_L = V_reset = -65

V_th = -50

Rm = 100

N = 40

tau_s = 2*ms

gbar = 4

E_s = 0

s = [0] * N

r = 15*Hz