function []=netcall(Einit,Iinit,PE)
%% Carol Lucas

[t,x]=ode45('netfun',[0 20],[Einit;Iinit],[],PE);

figure(1); plot(t,x);xlabel('time');
mess='Hit return to continue.'
pause;
a=16;b=12;c=16;d=5;PI=-4;tauE=1;tauI=1;

I=[.0001:.01:1];E=[.001:.01:1];

num=((1)./E)-1;
den=exp(-a*E-PE);
Ienull=log(num./den)/b;

num=((1)./I)-1;
den=exp(d*I-PI);
elnull=-log(num./den)./c;
last=length(t);

figure(2);plot(E,Ienull,elnull,I,x(:,1),x(:,2),x(last,1),x(last,2),'ro');
axis([0 1 -.2 1]);xlabel('E');ylabel('I');title('Phase Plane');
mess='Hit return to continue.'
pause;