If you want an uninitialized 61-by-61 array of double precision floating point numbers you can do Array(61, 61)
If you want the array initialized to zero you can call zeros(61, 61) A few for loops should accomplish what you're after here. Depending how expenaive your cost function is to evaluate this doesn't sound very big at all. How many state and input dimensions do you have?
