Hi! It's just interesting to know under what kind of conditions I get a convergence problem when computing a matrix exponent dimension 16x16? because sometimes it appears and I don't know a reason for that and sometimes Scilab computes the matrix exponent (I mean expm() function). thanks in advance. --------- best regards, Ivan