Hi, this is on a Dec 12th snapshot on Windows, but I doubt that this problem has changed since then. Attention, the following example crashes (segfault/terminates) gretl: <hansl> matrices MM = defarray(I(2), I(3)) eval nelem(MM)     # 2, as expected MM = I(4)    # implicitly turn to a 1-elem array? actually legal? eval typestr(typeof(MM))    # array, OK eval nelem(MM)            # 1, semi-expected matrix m = MM[,1]    # illegal indexing -> crash </hansl>