function dy=traveling_wave_solve(t,y,i_electrode)


w1=0;%pi/4;
w2=.33;%pi/2;
w3=.67;%pi;
w4=1;%2*pi;

dy=zeros(4,1);

dy(1)=w1+sin(y(2)-y(1));

dy(2)=w2+sin(y(3)-y(2))+sin(y(1)-y(2));

dy(3)=w3+sin(y(2)-y(3))+sin(y(4)-y(3));

dy(4)=w4+sin(y(3)-y(4));