Warning

This version of the documentation is NOT an official release. You are looking at ‘latest’, which is in active and ongoing development. You can change versions on the bottom left of the screen.

aeif_cond_beta_multisynapse

import nest
import numpy as np
import matplotlib.pyplot as plt

neuron = nest.Create('aeif_cond_beta_multisynapse')
nest.SetStatus(neuron, {"V_peak": 0.0, "a": 4.0, "b": 80.5})
nest.SetStatus(neuron, {'E_rev': [0.0, 0.0, 0.0, -85.0],
                        'tau_decay': [50.0, 20.0, 20.0, 20.0],
                        'tau_rise': [10.0, 10.0, 1.0, 1.0]})

spike = nest.Create('spike_generator', params={'spike_times':
                                               np.array([10.0])})

voltmeter = nest.Create('voltmeter')

delays = [1.0, 300.0, 500.0, 700.0]
w = [1.0, 1.0, 1.0, 1.0]
for syn in range(4):
    nest.Connect(spike, neuron, syn_spec={'synapse_model': 'static_synapse',
                                          'receptor_type': 1 + syn,
                                          'weight': w[syn],
                                          'delay': delays[syn]})

nest.Connect(voltmeter, neuron)

nest.Simulate(1000.0)

Vms = voltmeter.get("events", "V_m")
ts = voltmeter.get("events", "times")

plt.plot(ts, Vms)
plt.show()