// Updated(18-7-07) // 2.2 F = [-1 0;1 0]; expm(F)