function dydt=coup_osc_func(t,y,A)

N=size(A,1);

A=[sparse(N,N) A;speye(N,N) sparse(N,N)];

dydt=A*y;

 
