Wirte a program with c++, or c or java with the following detials.

The basic differential equation: dV Euler method of numerical integration (https://en.wikipedia.org/wiki/Euler method )gives: or m e where At is the integration time step, small compared to t. Define Vth as the firing threshold and initial condition V(O) -Vreset Use Euler's method to integrate the basic equation over a period of time T much greater than . If V(t) 2 Vh then register time t as a spike time and set V(t) to Vreset for the next integration time step. Firing rate curve calculation: Loop over electrode input level /e in small increments starting at zero. For each input level, integrate the above equation over a period of time T, computing the average inter-spike interval, tisi, by counting the number of spikes during the integration time T and then dividing by T. To achieve statistical significance, choose a T that gives a large spike count. Compute the firing rate as risi1/tisi Plot as a function of input level. Experiment with different At, Vth, EL, and Rmvalues. See Fig 5.6 in the text for initial parameter settings. The basic differential equation: dV Euler method of numerical integration (https://en.wikipedia.org/wiki/Euler method )gives: or m e where At is the integration time step, small compared to t. Define Vth as the firing threshold and initial condition V(O) -Vreset Use Euler's method to integrate the basic equation over a period of time T much greater than . If V(t) 2 Vh then register time t as a spike time and set V(t) to Vreset for the next integration time step. Firing rate curve calculation: Loop over electrode input level /e in small increments starting at zero. For each input level, integrate the above equation over a period of time T, computing the average inter-spike interval, tisi, by counting the number of spikes during the integration time T and then dividing by T. To achieve statistical significance, choose a T that gives a large spike count. Compute the firing rate as risi1/tisi Plot as a function of input level. Experiment with different At, Vth, EL, and Rmvalues. See Fig 5.6 in the text for initial parameter settings