function V=lin_rw_mat(t,A,x)
J=[[0, 1];[ -1, 0]];
B=[[0,1];[25*cosh(5*x(1)),0]];
V=B*A;