function z=couple(t, x)
z=zeros(4,1);
z(1)=x(2);
z(2)=-2*x(1)+2*(x(3)-x(1));
z(3)=x(4);
z(4)=-2*(x(3)-x(1));