@tomleslie 

Dear Tom,

I need your help. I have a delay differential equation to solve and extract the value of the solution y(t) at a selected point of the independent variable t.  I am uploading a small sample code. 

Thanks.
 

Download delay-differential-equation.mw
 

Download delay-differential-equation.mw

I have 6 parameters, say p1 - p6, which enter my Monte Carlo Code and computes a 1-D Array of floats of size 1000 . I have an experimental array of floats of the same size. My code outputs an array which is the difference of these two arrays. The objective function is the sum of the sqares of the elements in this output array. I have initial guess for the parameters: p1i - p6i. I need to find the values of the 6 parameters which minimize the objective function.

Obviously I do not have an explicit mathematical form for the objective function and nor its Jacobian.

Could some one please help me in this task.  Thanks in advance. 

I have a spike train which is a series of time points, for example saved as a vector of floats. Could someone tell me how I can plot this spike train as a series of vertical bars at each of these time points ?

I would also like to know how to plot a series of spike trains one above the other[i.e. the raster plot] if I start with a 2-D matrix of floats where each row correspond to one spike train.

Thanks in advance.

The output of a solve command was:

solution := {p[1] = 2.788944999, p[2] = 4.940143518}, { p[1] = 15.29764736, p[2] = 4.946617373}

My question is: How to capture these 4 numbers in a 2 by 2 matrix ?

I tried assign, subs commands. Did not succeed. Could some one help, please?