function dydt = gen_net_octave(y,t)


if y(1)>1
  r12 = -1; 
else
   r12=1;
end;

if y(2)>1
  r23 = -1;
else
  r23 = 1;  
end;

if y(3)>1
  r31 = -1;
else
  r31 = 1;  
end;



e1 = 0.1;
e2 = 0.1;
e3 = 0.1;

l1 = 0.00001;
l2 = 0.00001;
l3 = 0.00001;

d1 = r31+e1*randn()-l1*y(1);
d2 = r12+e2*randn()-l2*y(2);
d3 = r23+e3*randn()-l3*y(3);

if y(1)<0.1 && d1<0.0
  d1=0.0;
end;

if y(2)<0.1 && d2<0.0
  d2=0.0;
end;

if y(3)<0.1 && d3<0.0
  d3=0.0;
end;

if y(1)>1.9 && d1>0.0
  d1=0.0;
end;

if y(2)>1.9 && d2>0.0
  d2=0.0;
end;

if y(3)>1.9 && d3>0.0
  d3=0.0;
end;
				
dydt = [d1,d2,d3];

endfunction
