The example script below crashes gretl every time (on linux with the cvs version from a couple of days ago). Note that the embedded function call in the string argument of fdjac is missing one required argument. Let me know if there is currently no time to address this, then I will file it in the bug tracker. thanks, sven <script> function getGGcoeffvec(const matrix coeffvec, const matrix beta_starmat,\ int numendo, int lagorder, int determcase, int numuexo, int numrexo) matrix result = vec(I(2)) return matrix result end function open denmark lagorder = 2 cirank = 1 determcase = 3 numendo = 4 numrexo = 0 numuexo = 0 matrix beta_starmat = ones(4,1) matrix coeffvec = vec(I(3)) # just a dummy for illustration matrix GGcoeffvec = getGGcoeffvec(coeffvec,beta_starmat,\ numendo,lagorder,determcase,numuexo,numrexo) matrix GGvecjacmat = fdjac(coeffvec, "getGGcoeffvec(coeffvec,beta_starmat,\ numendo,lagorder,determcase,numuexo)") # note the missing numrexo arg! </script>