function dydt = gen_net_octave(y,t)


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

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

e1 = 0.1;
e2=0.1;
l1 = 0.00001;
l2 = 0.00001;

d1 = r21+e1*randn()-l1*y(1);
d2 = r12+e2*randn()-l2*y(2);

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

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

dydt = [d1,d2];

endfunction
