function udot=rnn_odeOctave(u,t)
% odefile for recurrent network
global dx w I_ext;
 tau_inv = 1.;      % inverse of membrane time constant
 r=1./(1+exp(-u)); 
 sum=w*r*dx;
 udot=tau_inv*(-u+sum+I_ext);
return
