double spike_rate(double t, int i) {
  std::cout<<t<<std::endl;
  std::cout<<spikemon.t[0]<<std::endl;

  return ta2d(t, i);
}