clear all; close all;
a=0;
b=10;
 
input= [0 0];
 
hold on;
xlabel('time')
ylabel('neuron output')
 
[T,Y] = ode45('test_de_fun',[a b],input);
plot(T,Y(:,1),'r', T,Y(:,2),'b');
 
 
legend({'n1' 'n2'});

hold off;


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%