clear all; close all;
global p;
global bp;
global bn;
global A;
global i
  p1 = 1;
  p2 = 1;
  bp = 2;
  bn = 2;
  p = [p1;p2];
  A=10;  
  [T,Y] = ode45('nndalay2a',[0 0.2],[0;0]);
  Y = Y';
  T = T';
  plot(T,Y(1,:),'b',T,Y(2,:),'r:','LineWidth',2);hold on;
  plot([0 0.2],[i(1) i(1)],'c:',[0 0.2],[i(2) i(2)],'m-.');
xlabel('t'); ylabel('n_1^2 & n_2^2');
title('Response of Layer 2');
legend('n_1^2(t)','n_2^2(t)','(_1w^1^:^2)^Ta^1','(_2w^1^:^2)^Ta^1');hold off;



 