Great. Thanks Nils. One more question. IF my y above was a 1D vector (different values generated at each iteration i) how would I concatenate this into the out matrix? Can I use the transpose command to transform y into a row 1D array (vector)? Something like?
out[i,:] = vcat(transpose(y) ,a[i],b[i])
Thanks Peter