function dy=sin_osc_solve(t,y,i_electrode)

w=1;
a12=0;
a13=0;
a14=1;
a21=0;
a23=-1;
a24=0;
a31=0;
a32=1;
a34=0;
a41=0;
a42=0;
a43=-1;

% aa=y(1)
% bb=y(2)


dy = zeros(4,1); 

% a=y
% if find(y>2*pi)
% 	
% % 	evalin('base','y=mod(y,(2*pi))')
% end

dy(1)=w+a12*sin(y(2)-y(1))+a13*sin(y(3)-y(1)) +a14*sin(y(4)-y(1));

dy(2)=w+a21*sin(y(1)-y(2))+a23*sin(y(3)-y(2)) +a24*sin(y(4)-y(2));

dy(3)=w+a31*sin(y(1)-y(3))+a32*sin(y(2)-y(3)) +a34*sin(y(4)-y(3));

dy(4)=w+a41*sin(y(1)-y(4))+a42*sin(y(2)-y(4)) +a43*sin(y(3)-y(4));



% y(1)=v1
% y(2)=v2
% y(3)=s1
% y(4)=s2