clear, close all;

x = 0.05:1.1/24:1.15;
y = 1e-6:2e-6/24:3e-6;

z(1) = 1.113;
z(2) = 1.19;
z(3) = 1.268;
z(4) = 1.338;
z(5) = 1.411;
z(6) = 1.486;
z(7) = 1.563;
z(8) = 1.638;
z(9) = 1.714;
z(10) = 1.79;
z(11) = 1.866;
z(12) = 1.943;
z(13) = 2.018;
z(14) = 2.095;
z(15) = 2.171;
z(16) = 2.248;
z(17) = 2.324;
z(18) = 2.401;
z(19) = 2.477;
z(20) = 2.554;
z(21) = 2.63;
z(22) = 2.707;
z(23) = 2.784;
z(24) = 2.861;
z(25) = 2.944;

zr(1) = 1;
zr(2) = 1.117;
zr(3) = 1.201;
zr(4) = 1.284;
zr(5) = 1.367;
zr(6) = 1.451;
zr(7) = 1.534;
zr(8) = 1.618;
zr(9) = 1.702;
zr(10) = 1.785;
zr(11) = 1.868;
zr(12) = 1.952;
zr(13) = 2.036;
zr(14) = 2.119;
zr(15) = 2.202;
zr(16) = 2.287;
zr(17) = 2.371;
zr(18) = 2.457;
zr(19) = 2.546;
zr(20) = 2.641;
zr(21) = 2.745;
zr(22) = 2.879;
zr(23) = 3.09;
zr(24) = 3;
zr(25) = 3;

zr2(1) = 1.1;
zr2(2) = 1.178;
zr2(3) = 1.255;
zr2(4) = 1.332;
zr2(5) = 1.41;
zr2(6) = 1.487;
zr2(7) = 1.565;
zr2(8) = 1.642;
zr2(9) = 1.72;
zr2(10) = 1.798;
zr2(11) = 1.875;
zr2(12) = 1.953;
zr2(13) = 2.031;
zr2(14) = 2.108;
zr2(15) = 2.185;
zr2(16) = 2.261;
zr2(17) = 2.338;
zr2(18) = 2.413;
zr2(19) = 2.489;
zr2(20) = 2.565;
zr2(21) = 2.64;
zr2(22) = 2.716;
zr2(23) = 2.791;
zr2(24) = 2.864;
zr2(25) = 2.925;

zi(1) = 1.112;
zi(2) = 1.189;
zi(3) = 1.266;
zi(4) = 1.344;
zi(5) = 1.421;
zi(6) = 1.498;
zi(7) = 1.575;
zi(8) = 1.653;
zi(9) = 1.73;
zi(10) = 1.808;
zi(11) = 1.885;
zi(12) = 1.963;
zi(13) = 2.041;
zi(14) = 2.118;
zi(15) = 2.196;
zi(16) = 2.274;
zi(17) = 2.431;
zi(18) = 2.509;
zi(19) = 2.588;
zi(20) = 2.672;
zi(21) = 2.772;
zi(22) = 2.914;
zi(23) = 3.199;
zi(24) = 3;
zi(25) = 3;

zip(1) = 1.191;
zip(2) = 1.222;
zip(3) = 1.272;
zip(4) = 1.337;
zip(5) = 1.411;
zip(6) = 1.485;
zip(7) = 1.561;
zip(8) = 1.637;
zip(9) = 1.714;
zip(10) = 1.789;
zip(11) = 1.865;
zip(12) = 1.94;
zip(13) = 2.018;
zip(14) = 2.094;
zip(15) = 2.17;
zip(16) = 2.247;
zip(17) = 2.322;
zip(18) = 2.399;
zip(19) = 2.467;
zip(20) = 2.552;
zip(21) = 2.63;
zip(22) = 2.705;
zip(23) = 2.782;
zip(24) = 2.861;
zip(25) = 2.942;

resp = y-z*1e-6;
resp2 = y-zr*1e-6;
resp2b = y-zr2*1e-6;
resp3 = y-zi*1e-6;
resp4 = y-zip*1e-6;
hold on;

figure(1)
%plot(x,-resp,'r');
plot(x,-resp2,'b');
plot(x,-resp2b,'g');
%plot(x,-resp3,'g');
%plot(x,-resp4,'m');
xlabel('threshold mV');
ylabel('response time s');