function yprime = nndalay1(t,y)
%NNDALAY1 Calculates the derivative for layer 1 of the ART
%
%   NNDALAY1(t,y)
%     t - Current time
%     y - Current output
%   Returns dy.
 
 
global p;
global bp;
global bn;
 
w = [0  1; 1  1];
e = 0.02;
 
yprime = [-y(1) + (bp - y(1))*(p(1)+w(1,2)) - (y(1) + bn);
          -y(2) + (bp - y(2))*(p(2)+w(2,2)) - (y(2) + bn)] /e;