syms x t
f1='sin(x)^2+cos(x)^3-3'
A = sym('A',[2 2])
det(A*t)
expm(A*t);

A=[1 2;3 5];B=[2;2];C=[1 0];
sys=ss(A,B,C,[]);
sys_dis=c2d(sys,1,'zoh');
sys_con=d2c(sys_dis);
G=sys_dis.a
A=sys_con.a

impulse(sys)