clear all; close all;
global p;
global bp;
global bn;
 
  p1 = 0;
  p2 = 1;
  bp = 2;
  bn = 3;
  p = [p1;p2];
    
  [T,Y] = ode45('nndalay1',[0 0.1],[0;0]);
  Y = Y';
  T = T';
    plot(T,Y)  
xlabel('t'); ylabel('n_1^1 & n_2^1');
title('Response of Layer 1');
legend('n_1^1(t)','n_2^1(t)')