def METHOD_NAME(self, i_mean=(400.0, 900.0, 50.0), i_std=(0.0, 600.0, 50.0)):
    #######################################################################
    # We loop through all possible combinations of `(I_mean, I_sigma)`
    # and measure the output rate of the neuron.
    self.i_range = numpy.arange(*i_mean)
    self.std_range = numpy.arange(*i_std)
    self.rate = numpy.zeros((self.i_range.size, self.std_range.size))
    nest.set_verbosity("M_WARNING")
    for n, i in enumerate(self.i_range):
        print("I  =  {0}".format(i))
        for m, std in enumerate(self.std_range):
            self.rate[n, m] = self.output_rate(i, std)