diff --git a/doc/lqi-pdr/1000pkt-64len.csv b/doc/lqi-pdr/1000pkt-64len.csv
new file mode 100644
index 000000000..05520e6e2
--- /dev/null
+++ b/doc/lqi-pdr/1000pkt-64len.csv
@@ -0,0 +1,41 @@
+1,6
+5,4
+5,6
+5,6
+6,3
+30,5
+45,3
+51,9
+162,6
+170,7
+308,6
+317,6
+506,8
+511,7
+596,7
+602,7
+674,9
+743,8
+788,8
+842,9
+859,9
+951,9
+898,10
+976,12
+978,19
+981,13
+988,30
+991,17
+992,12
+993,16
+993,19
+994,15
+995,19
+995,19
+998,30
+999,38
+1000,34
+1000,40
+1000,46
+1000,47
+1000,54
diff --git a/doc/lqi-pdr/README b/doc/lqi-pdr/README
new file mode 100644
index 000000000..699c87b2d
--- /dev/null
+++ b/doc/lqi-pdr/README
@@ -0,0 +1,7 @@
+grep count 1000pkt-64len.txt | cut -d ' '  -f 2,5 | sed 's/ /,/g' |
+sort -n > 1000pkt-64len.csv
+
+then:
+
+ asy plot.asy
+ gv plot.eps
diff --git a/doc/lqi-pdr/plot.asy b/doc/lqi-pdr/plot.asy
new file mode 100644
index 000000000..7c2bef059
--- /dev/null
+++ b/doc/lqi-pdr/plot.asy
@@ -0,0 +1,13 @@
+import graph;
+size(350,250,IgnoreAspect);
+
+file fin=input("./1000pkt-64len.csv");
+real[][] A=dimension(csv(fin),0,2);
+real[][] pdr=transpose(A);
+
+draw(graph(pdr[1],pdr[0]));
+
+ylimits(0,1000);
+xlimits(0,100);
+xaxis("\rm Output Current (A)",BottomTop,LeftTicks("$%.1f$",10,begin=false,end=false,extend=true,pTick=dotted));
+yaxis("\rm Output Power (W)",LeftRight,RightTicks("$%#.1f$",8,begin=false,end=false,extend=true,pTick=dotted, ptick=dotted));