syms t x y z; a=sym('a(t)'); g=[-1,0,0,0;0,a^2,0
Sun Dec 08 14:33:17 CST 2019
Fri Aug 14 10:47:20 CST 2020