function V=nl_mat(t,A,x)
J=[[0, 1];[ -1, 0]];
B=J+[[0,0];[2*0.01*x(1),0]];
V=B*A;