diff --git a/README.md b/README.md
index f488c9d..726a7ff 100644
--- a/README.md
+++ b/README.md
@@ -82,6 +82,16 @@ true     # Restrict the hole to the 1'st (i.e. half of natom) Family of states.
 mpiexec -n [nprocs] ./ex1 inpfile 
 ```
 
+_Sample Application_
+--------------------
+
+A 2D t-J model Hamiltonian description and setup for using DEHam to solve for few low lying states
+is provided in the notbooks folder. Please have a look about the details of using DEHam to study
+t-J Hamiltonians.
+
+![4x4 Lattice](https://raw.githubusercontent.com/v1j4y/DEHam/master/notebook/graph.png)
+
+
 _Publications using this code_
 -------------------------------
 
diff --git a/notebooks/graph.png b/notebooks/graph.png
new file mode 100644
index 0000000..d7e911a
Binary files /dev/null and b/notebooks/graph.png differ
diff --git a/notebooks/t_J_Model.ipynb b/notebooks/t_J_Model.ipynb
index a9259b6..d10d2f1 100644
--- a/notebooks/t_J_Model.ipynb
+++ b/notebooks/t_J_Model.ipynb
@@ -160,7 +160,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 138,
+   "execution_count": 150,
    "metadata": {},
    "outputs": [
     {
@@ -178,11 +178,14 @@
     "import matplotlib as mplt\n",
     "import networkx as nx\n",
     "import numpy as np\n",
+    "import matplotlib.pyplot as plt\n",
     "\n",
     "Lx = 4\n",
     "Ly = 4\n",
-    "\n",
-    "draw2DLattice(Lx,Ly)"
+    "#fig = plt.figure()\n",
+    "#draw2DLattice(Lx,Ly,fig)\n",
+    "draw2DLattice(Lx,Ly)\n",
+    "#fig.savefig(\"/tmp/graph.png\",dpi=1000)"
    ]
   },
   {
@@ -291,7 +294,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 17,
+   "execution_count": 149,
    "metadata": {
     "jupyter": {
      "source_hidden": true
@@ -364,6 +367,8 @@
     "        countcolor+=1\n",
     "\n",
     "\n",
+    "#    nx.draw(G, pos=pos, labels=labelscenter, node_color=color_map, font_color=\"white\", ax=fig.add_subplot(111));\n",
+    "#    nx.draw_networkx_edge_labels(G,pos,edge_labels=edgescenter,font_color=\"black\", ax=fig.add_subplot(111));\n",
     "    nx.draw(G, pos=pos, labels=labelscenter, node_color=color_map, font_color=\"white\");\n",
     "    nx.draw_networkx_edge_labels(G,pos,edge_labels=edgescenter,font_color=\"black\");\n"
    ]