diff --git a/ipynb/Examples/egfr-simple.ipynb b/ipynb/Examples/egfr-simple.ipynb
new file mode 100644
index 000000000..f227771fa
--- /dev/null
+++ b/ipynb/Examples/egfr-simple.ipynb
@@ -0,0 +1,371 @@
+{
+ "cells": [
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "# Simple EGFR model\n",
+    "\n",
+    "- http://bionetgen.org/index.php/Simple_EGFR_model\n",
+    "- M.L. Blinov, J.R. Faeder, B. Goldstein, W.S. Hlavacek, \"A network model of early events in epidermal growth factor receptor signaling that accounts for combinatorial complexity.\", Biosystems, 83(2-3), 136-151, 2006."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 1,
+   "metadata": {
+    "collapsed": true
+   },
+   "outputs": [],
+   "source": [
+    "%matplotlib inline\n",
+    "from ecell4 import *"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 2,
+   "metadata": {
+    "collapsed": true
+   },
+   "outputs": [],
+   "source": [
+    "NA = 6.02e23  # Avogadro's number (molecules/mol)\n",
+    "f = 1  # Fraction of the cell to simulate\n",
+    "Vo = f * 1.0e-10  # Extracellular volume=1/cell_density (L)\n",
+    "V = f * 3.0e-12  # Cytoplasmic volume (L)\n",
+    "\n",
+    "EGF_init = 20 * 1e-9 * NA * Vo  # Initial amount of ligand (20 nM) converted to copies per cell\n",
+    "\n",
+    "# Initial amounts of cellular components (copies per cell)\n",
+    "EGFR_init = f * 1.8e5\n",
+    "Grb2_init = f * 1.5e5\n",
+    "Sos1_init = f * 6.2e4\n",
+    "\n",
+    "# Rate constants\n",
+    "# Divide by NA*V to convert bimolecular rate constants\n",
+    "# from /M/sec to /(molecule/cell)/sec\n",
+    "kp1 = 9.0e7 / (NA * Vo)  # ligand-monomer binding\n",
+    "km1 = 0.06  # ligand-monomer dissociation\n",
+    "kp2 = 1.0e7 / (NA * V)  # aggregation of bound monomers\n",
+    "km2 = 0.1  # dissociation of bound monomers\n",
+    "kp3 = 0.5  # dimer transphosphorylation\n",
+    "km3 = 4.505  # dimer dephosphorylation\n",
+    "kp4 = 1.5e6 / (NA * V)  # binding of Grb2 to receptor\n",
+    "km4 = 0.05  # dissociation of Grb2 from receptor\n",
+    "kp5 = 1.0e7 / (NA * V)  # binding of Grb2 to Sos1\n",
+    "km5 = 0.06  # dissociation of Grb2 from Sos1\n",
+    "deg = 0.01  # degradation of receptor dimers"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 3,
+   "metadata": {
+    "collapsed": false
+   },
+   "outputs": [],
+   "source": [
+    "with reaction_rules():\n",
+    "  # R1: Ligand-receptor binding\n",
+    "  EGFR(L, CR1) + EGF(R) == EGFR(L^1, CR1).EGF(R^1) | (kp1, km1)\n",
+    "  \n",
+    "  # R2: Receptor-aggregation\n",
+    "  EGFR(L^_, CR1) + EGFR(L^_, CR1) == EGFR(L^_,CR1^1).EGFR(L^_,CR1^1) | (kp2, km2)\n",
+    "  \n",
+    "  # R3: Transphosphorylation of EGFR by RTK\n",
+    "  EGFR(CR1^_, Y1068=U) > EGFR(CR1^_, Y1068=P) | kp3\n",
+    "  \n",
+    "  # R4: Dephosphorylation\n",
+    "  EGFR(Y1068=P) > EGFR(Y1068=U) | km3\n",
+    "  \n",
+    "  # R5: Grb2 binding to pY1068\n",
+    "  EGFR(Y1068=P) + Grb2(SH2) == EGFR(Y1068=P^1).Grb2(SH2^1) | (kp4, km4)\n",
+    "  \n",
+    "  # R6: Grb2 binding to Sos1\n",
+    "  Grb2(SH3) + Sos1(PxxP) == Grb2(SH3^1).Sos1(PxxP^1) | (kp5, km5)\n",
+    "  \n",
+    "  # R7: Receptor dimer internalization/degradation\n",
+    "  (EGF(R^1).EGF(R^2).EGFR(L^1,CR1^3).EGFR(L^2,CR1^3) > ~EmptySet | deg\n",
+    "       | _policy(ReactionRule.STRICT | ReactionRule.DESTROY))\n",
+    "\n",
+    "m = get_model(is_netfree=True, effective=True)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 4,
+   "metadata": {
+    "collapsed": false
+   },
+   "outputs": [],
+   "source": [
+    "y0 = {\"EGF(R)\": EGF_init, \"EGFR(L, CR1, Y1068=U)\": EGFR_init, \"Grb2(SH2, SH3)\": Grb2_init, \"Sos1(PxxP)\": Sos1_init}"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 5,
+   "metadata": {
+    "collapsed": false
+   },
+   "outputs": [],
+   "source": [
+    "newm = m.expand([Species(serial) for serial in y0.keys()])"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 6,
+   "metadata": {
+    "collapsed": false
+   },
+   "outputs": [
+    {
+     "name": "stdout",
+     "output_type": "stream",
+     "text": [
+      "22 species and 86 reactions were generated.\n",
+      "1: EGF(R)\n",
+      "2: EGF(R^1).EGFR(CR1,L^1,Y1068=P)\n",
+      "3: EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)\n",
+      "4: EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)\n",
+      "5: EGF(R^1).EGFR(CR1,L^1,Y1068=U)\n",
+      "6: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3)\n",
+      "7: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3)\n",
+      "8: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5)\n",
+      "9: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)\n",
+      "10: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3)\n",
+      "11: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6)\n",
+      "12: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3^6).Sos1(PxxP^6).Grb2(SH2^3,SH3^7).Sos1(PxxP^7)\n",
+      "13: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3)\n",
+      "14: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5)\n",
+      "15: EGF(R^1).EGFR(CR1^2,L^1,Y1068=U).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)\n",
+      "16: EGFR(CR1,L,Y1068=P)\n",
+      "17: EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3)\n",
+      "18: EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3^2).Sos1(PxxP^2)\n",
+      "19: EGFR(CR1,L,Y1068=U)\n",
+      "20: Grb2(SH2,SH3)\n",
+      "21: Grb2(SH2,SH3^1).Sos1(PxxP^1)\n",
+      "22: Sos1(PxxP)\n",
+      "1: EGF(R)+EGFR(CR1,L,Y1068=U)>EGF(R^1).EGFR(CR1,L^1,Y1068=U)|1.49502e-06\n",
+      "2: Grb2(SH2,SH3)+Sos1(PxxP)>Grb2(SH2,SH3^1).Sos1(PxxP^1)|5.5371e-06\n",
+      "3: EGF(R^1).EGFR(CR1,L^1,Y1068=U)>EGF(R)+EGFR(CR1,L,Y1068=U)|0.06\n",
+      "4: EGF(R^1).EGFR(CR1,L^1,Y1068=U)+EGF(R^1).EGFR(CR1,L^1,Y1068=U)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=U).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)|2.76855e-06\n",
+      "5: Grb2(SH2,SH3^1).Sos1(PxxP^1)>Grb2(SH2,SH3)+Sos1(PxxP)|0.06\n",
+      "6: EGF(R^1).EGFR(CR1^2,L^1,Y1068=U).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)>EGF(R^1).EGFR(CR1,L^1,Y1068=U)+EGF(R^1).EGFR(CR1,L^1,Y1068=U)|0.1\n",
+      "7: EGF(R^1).EGFR(CR1^2,L^1,Y1068=U).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)|1\n",
+      "8: EGF(R^1).EGFR(CR1^2,L^1,Y1068=U).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)>|0.01\n",
+      "9: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)>EGF(R^1).EGFR(CR1,L^1,Y1068=P)+EGF(R^1).EGFR(CR1,L^1,Y1068=U)|0.1\n",
+      "10: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3)|0.5\n",
+      "11: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=U).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)|4.505\n",
+      "12: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)+Grb2(SH2,SH3^1).Sos1(PxxP^1)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5)|8.30565e-07\n",
+      "13: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)+Grb2(SH2,SH3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3)|8.30565e-07\n",
+      "14: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)>|0.01\n",
+      "15: EGF(R^1).EGFR(CR1,L^1,Y1068=P)>EGF(R)+EGFR(CR1,L,Y1068=P)|0.06\n",
+      "16: EGF(R^1).EGFR(CR1,L^1,Y1068=P)+EGF(R^1).EGFR(CR1,L^1,Y1068=P)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3)|2.76855e-06\n",
+      "17: EGF(R^1).EGFR(CR1,L^1,Y1068=P)+EGF(R^1).EGFR(CR1,L^1,Y1068=U)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)|5.5371e-06\n",
+      "18: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3)>EGF(R^1).EGFR(CR1,L^1,Y1068=P)+EGF(R^1).EGFR(CR1,L^1,Y1068=P)|0.1\n",
+      "19: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5)>EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)+EGF(R^1).EGFR(CR1,L^1,Y1068=U)|0.1\n",
+      "20: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3)>EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)+EGF(R^1).EGFR(CR1,L^1,Y1068=U)|0.1\n",
+      "21: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5)|0.5\n",
+      "22: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3)|0.5\n",
+      "23: EGF(R^1).EGFR(CR1,L^1,Y1068=P)>EGF(R^1).EGFR(CR1,L^1,Y1068=U)|4.505\n",
+      "24: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)|9.01\n",
+      "25: EGF(R^1).EGFR(CR1,L^1,Y1068=P)+Grb2(SH2,SH3^1).Sos1(PxxP^1)>EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)|8.30565e-07\n",
+      "26: EGF(R^1).EGFR(CR1,L^1,Y1068=P)+Grb2(SH2,SH3)>EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)|8.30565e-07\n",
+      "27: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3)+Grb2(SH2,SH3^1).Sos1(PxxP^1)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5)|1.66113e-06\n",
+      "28: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3)+Grb2(SH2,SH3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3)|1.66113e-06\n",
+      "29: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)+Grb2(SH2,SH3^1).Sos1(PxxP^1)|0.05\n",
+      "30: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=U).EGF(R^3)+Grb2(SH2,SH3)|0.05\n",
+      "31: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3)+Sos1(PxxP)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5)|5.5371e-06\n",
+      "32: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3)+Sos1(PxxP)|0.06\n",
+      "33: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3)>|0.01\n",
+      "34: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5)>|0.01\n",
+      "35: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3)>|0.01\n",
+      "36: EGF(R)+EGFR(CR1,L,Y1068=P)>EGF(R^1).EGFR(CR1,L^1,Y1068=P)|1.49502e-06\n",
+      "37: EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)>EGF(R)+EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3^2).Sos1(PxxP^2)|0.06\n",
+      "38: EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)>EGF(R)+EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3)|0.06\n",
+      "39: EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)+EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3^6).Sos1(PxxP^6).Grb2(SH2^3,SH3^7).Sos1(PxxP^7)|2.76855e-06\n",
+      "40: EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)+EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6)|5.5371e-06\n",
+      "41: EGF(R^1).EGFR(CR1,L^1,Y1068=P)+EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5)|5.5371e-06\n",
+      "42: EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)+EGF(R^1).EGFR(CR1,L^1,Y1068=U)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5)|5.5371e-06\n",
+      "43: EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)+EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3)|2.76855e-06\n",
+      "44: EGF(R^1).EGFR(CR1,L^1,Y1068=P)+EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3)|5.5371e-06\n",
+      "45: EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)+EGF(R^1).EGFR(CR1,L^1,Y1068=U)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3)|5.5371e-06\n",
+      "46: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5)>EGF(R^1).EGFR(CR1,L^1,Y1068=P)+EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)|0.1\n",
+      "47: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3)>EGF(R^1).EGFR(CR1,L^1,Y1068=P)+EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)|0.1\n",
+      "48: EGFR(CR1,L,Y1068=P)>EGFR(CR1,L,Y1068=U)|4.505\n",
+      "49: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3^5).Sos1(PxxP^5)|4.505\n",
+      "50: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=U).EGF(R^4).Grb2(SH2^3,SH3)|4.505\n",
+      "51: EGFR(CR1,L,Y1068=P)+Grb2(SH2,SH3^1).Sos1(PxxP^1)>EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3^2).Sos1(PxxP^2)|8.30565e-07\n",
+      "52: EGFR(CR1,L,Y1068=P)+Grb2(SH2,SH3)>EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3)|8.30565e-07\n",
+      "53: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5)+Grb2(SH2,SH3^1).Sos1(PxxP^1)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3^6).Sos1(PxxP^6).Grb2(SH2^3,SH3^7).Sos1(PxxP^7)|8.30565e-07\n",
+      "54: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5)+Grb2(SH2,SH3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6)|8.30565e-07\n",
+      "55: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3)+Grb2(SH2,SH3^1).Sos1(PxxP^1)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6)|8.30565e-07\n",
+      "56: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3)+Grb2(SH2,SH3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3)|8.30565e-07\n",
+      "57: EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)>EGF(R^1).EGFR(CR1,L^1,Y1068=P)+Grb2(SH2,SH3^1).Sos1(PxxP^1)|0.05\n",
+      "58: EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)>EGF(R^1).EGFR(CR1,L^1,Y1068=P)+Grb2(SH2,SH3)|0.05\n",
+      "59: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3)+Grb2(SH2,SH3^1).Sos1(PxxP^1)|0.05\n",
+      "60: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P).EGF(R^3)+Grb2(SH2,SH3)|0.05\n",
+      "61: EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)+Sos1(PxxP)>EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)|5.5371e-06\n",
+      "62: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3)+Sos1(PxxP)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5)|5.5371e-06\n",
+      "63: EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)>EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)+Sos1(PxxP)|0.06\n",
+      "64: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3)+Sos1(PxxP)|0.06\n",
+      "65: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5)>|0.01\n",
+      "66: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3)>|0.01\n",
+      "67: EGF(R)+EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3^2).Sos1(PxxP^2)>EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)|1.49502e-06\n",
+      "68: EGF(R)+EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3)>EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)|1.49502e-06\n",
+      "69: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3^6).Sos1(PxxP^6).Grb2(SH2^3,SH3^7).Sos1(PxxP^7)>EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)+EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)|0.1\n",
+      "70: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6)>EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)+EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3^3).Sos1(PxxP^3)|0.1\n",
+      "71: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3)>EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)+EGF(R^1).EGFR(CR1,L^1,Y1068=P^2).Grb2(SH2^2,SH3)|0.1\n",
+      "72: EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3^2).Sos1(PxxP^2)>EGFR(CR1,L,Y1068=P)+Grb2(SH2,SH3^1).Sos1(PxxP^1)|0.05\n",
+      "73: EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3)>EGFR(CR1,L,Y1068=P)+Grb2(SH2,SH3)|0.05\n",
+      "74: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3^6).Sos1(PxxP^6).Grb2(SH2^3,SH3^7).Sos1(PxxP^7)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5)+Grb2(SH2,SH3^1).Sos1(PxxP^1)|0.1\n",
+      "75: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3)+Grb2(SH2,SH3^1).Sos1(PxxP^1)|0.05\n",
+      "76: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3^5).Sos1(PxxP^5)+Grb2(SH2,SH3)|0.05\n",
+      "77: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P).EGFR(CR1^2,L^3,Y1068=P^4).EGF(R^3).Grb2(SH2^4,SH3)+Grb2(SH2,SH3)|0.1\n",
+      "78: EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3)+Sos1(PxxP)>EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3^2).Sos1(PxxP^2)|5.5371e-06\n",
+      "79: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6)+Sos1(PxxP)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3^6).Sos1(PxxP^6).Grb2(SH2^3,SH3^7).Sos1(PxxP^7)|5.5371e-06\n",
+      "80: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3)+Sos1(PxxP)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6)|1.10742e-05\n",
+      "81: EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3^2).Sos1(PxxP^2)>EGFR(CR1,L,Y1068=P^1).Grb2(SH2^1,SH3)+Sos1(PxxP)|0.06\n",
+      "82: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3^6).Sos1(PxxP^6).Grb2(SH2^3,SH3^7).Sos1(PxxP^7)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6)+Sos1(PxxP)|0.12\n",
+      "83: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6)>EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3)+Sos1(PxxP)|0.06\n",
+      "84: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3^6).Sos1(PxxP^6).Grb2(SH2^3,SH3^7).Sos1(PxxP^7)>|0.01\n",
+      "85: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3^6).Sos1(PxxP^6)>|0.01\n",
+      "86: EGF(R^1).EGFR(CR1^2,L^1,Y1068=P^3).EGFR(CR1^2,L^4,Y1068=P^5).EGF(R^4).Grb2(SH2^5,SH3).Grb2(SH2^3,SH3)>|0.01\n"
+     ]
+    }
+   ],
+   "source": [
+    "print(\"{} species and {} reactions were generated.\".format(len(newm.list_species()), len(newm.reaction_rules())))\n",
+    "\n",
+    "for i, sp in enumerate(newm.list_species()):\n",
+    "    print(\"{}: {}\".format(i + 1, sp.serial()))\n",
+    "\n",
+    "for i, rr in enumerate(newm.reaction_rules()):\n",
+    "    print(\"{}: {}\".format(i + 1, rr.as_string()))"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 7,
+   "metadata": {
+    "collapsed": false
+   },
+   "outputs": [
+    {
+     "data": {
+      "text/html": [
+       "<html>\n",
+       "  <head>\n",
+       "    <script type=\"text/javascript\">\n",
+       "     (function(){\n",
+       "       if(window[\"THREE\"] === undefined || window[\"d3\"] === undefined){\n",
+       "     require.config({\n",
+       "       paths:{\n",
+       "         d3: [\n",
+       "           'https://d3js.org/d3.v3.min',\n",
+       "           'custom/d3.min'\n",
+       "         ],\n",
+       "         THREE: [\n",
+       "           'https://cdnjs.cloudflare.com/ajax/libs/three.js/r66/three.min',\n",
+       "           'custom/three.min'\n",
+       "         ],\n",
+       "         Nyaplot: [\n",
+       "           'https://rawgit.com/domitry/Nyaplotjs/master/release/nyaplot',\n",
+       "           'custom/nyaplot'\n",
+       "         ],\n",
+       "         svg2png: 'https://cdn.rawgit.com/domitry/svg2png/850731291ef598b201687ed5b1e215df248d8aeb/src/main'\n",
+       "       },\n",
+       "       shim:{\n",
+       "         THREE: {\n",
+       "           exports: 'THREE'\n",
+       "         },\n",
+       "         Nyaplot: {\n",
+       "           deps: ['THREE'],\n",
+       "           exports: 'Nyaplot'\n",
+       "         }\n",
+       "       }\n",
+       "     });\n",
+       "\n",
+       "         require([\"d3\", \"THREE\"], function(d3){\n",
+       "           window[\"d3\"] = d3;\n",
+       "\n",
+       "           d3.select(document)\n",
+       "             .on(\"contextmenu\", function(){\n",
+       "               d3.selectAll(\"canvas\").style(\"z-index\", 10);\n",
+       "               d3.selectAll(\".download_menu\").style(\"z-index\", 100);\n",
+       "             });\n",
+       "         });\n",
+       "       }\n",
+       "\n",
+       "       require([\"svg2png\", \"Nyaplot\"], function(util){\n",
+       "         console.log(\"Begin rendering...\");\n",
+       "         var models = {\"data\": {\"data1\": [{\"y6\": 0.0, \"y4\": 0.0, \"y5\": 0.0, \"y2\": 180000.0, \"y1\": 1204000.0, \"x\": 0.0, \"y3\": 0.0}, {\"y6\": 34563.88146720779, \"y4\": 3895.399459291695, \"y5\": 40.33776119946268, \"y2\": 179778.63239455546, \"y1\": 1051939.9975628783, \"x\": 1.2, \"y3\": 48962.12627523841}, {\"y6\": 46164.41708806864, \"y4\": 9632.256561505948, \"y5\": 337.79736907140807, \"y2\": 178885.57215214949, \"y1\": 1031169.2811889593, \"x\": 2.4, \"y3\": 93744.63204801833}, {\"y6\": 50620.92423822618, \"y4\": 12939.590153575129, \"y5\": 823.2521490692546, \"y2\": 177634.48851418888, \"y1\": 1027335.054003944, \"x\": 3.6, \"y3\": 111975.97830265302}, {\"y6\": 52294.947031363256, \"y4\": 15048.913711080926, \"y5\": 1354.6781366006146, \"y2\": 176241.34100467025, \"y1\": 1026403.8254919862, \"x\": 4.8, \"y3\": 119092.0257303489}, {\"y6\": 52790.94693294999, \"y4\": 16597.94804461366, \"y5\": 1868.3622021353856, \"y2\": 174794.13014835457, \"y1\": 1026098.117250155, \"x\": 6.0, \"y3\": 121628.70018396988}, {\"y6\": 52777.33607410532, \"y4\": 17839.577224184155, \"y5\": 2343.0008126688263, \"y2\": 173330.09738863737, \"y1\": 1025970.7984679276, \"x\": 7.2, \"y3\": 122159.92501592467}, {\"y6\": 52549.140045520006, \"y4\": 18883.80308878195, \"y5\": 2774.0491536017566, \"y2\": 171865.8796898094, \"y1\": 1025908.8923385559, \"x\": 8.4, \"y3\": 121775.95509160103}, {\"y6\": 52239.66358550506, \"y4\": 19784.247558288123, \"y5\": 3162.870410049141, \"y2\": 170409.14388240955, \"y1\": 1025874.9665564192, \"x\": 9.6, \"y3\": 120966.5473226915}, {\"y6\": 51909.245548412946, \"y4\": 20570.16020388454, \"y5\": 3512.5923916087318, \"y2\": 168963.4568779353, \"y1\": 1025853.9399825004, \"x\": 10.8, \"y3\": 119959.77211231604}, {\"y6\": 51584.671334125254, \"y4\": 21259.44287285536, \"y5\": 3826.633586908932, \"y2\": 167530.45765435803, \"y1\": 1025839.0733024004, \"x\": 12.0, \"y3\": 118864.03971100145}, {\"y6\": 51277.19206170182, \"y4\": 21864.35753561827, \"y5\": 4108.2341183320505, \"y2\": 166110.8633618734, \"y1\": 1025827.2077324329, \"x\": 13.2, \"y3\": 117731.47736347724}, {\"y6\": 50990.88478857215, \"y4\": 22394.16724612147, \"y5\": 4360.342249537722, \"y2\": 164704.94673289137, \"y1\": 1025816.830011741, \"x\": 14.4, \"y3\": 116587.29931697063}, {\"y6\": 50726.55706449017, \"y4\": 22856.40067924043, \"y5\": 4585.614666019147, \"y2\": 163312.764840752, \"y1\": 1025807.2067495057, \"x\": 15.6, \"y3\": 115443.73494163548}, {\"y6\": 50483.56947475458, \"y4\": 23257.48055691496, \"y5\": 4786.446469558364, \"y2\": 161934.27049307313, \"y1\": 1025797.9801109537, \"x\": 16.8, \"y3\": 114306.70154225518}, {\"y6\": 50260.688799102754, \"y4\": 23603.047617878074, \"y5\": 4965.003397158172, \"y2\": 160569.36625820198, \"y1\": 1025788.9750304032, \"x\": 18.0, \"y3\": 113179.032622297}, {\"y6\": 50056.48316638921, \"y4\": 23898.149629919437, \"y5\": 5123.249531975854, \"y2\": 159217.9303740005, \"y1\": 1025780.1049366036, \"x\": 19.2, \"y3\": 112062.05902325027}, {\"y6\": 49869.49934515067, \"y4\": 24147.33282026996, \"y5\": 5262.970227810063, \"y2\": 157879.8308696269, \"y1\": 1025771.3275068069, \"x\": 20.4, \"y3\": 110956.3665849309}, {\"y6\": 49698.341693717186, \"y4\": 24354.72354968805, \"y5\": 5385.791263470774, \"y2\": 156554.931100133, \"y1\": 1025762.6214544062, \"x\": 21.6, \"y3\": 109862.17923076462}, {\"y6\": 49541.70325365259, \"y4\": 24524.079045465536, \"y5\": 5493.194905238117, \"y2\": 155243.0922312267, \"y1\": 1025753.9756783393, \"x\": 22.8, \"y3\": 108779.54510120797}, {\"y6\": 49398.37406330757, \"y4\": 24658.815542248132, \"y5\": 5586.533621536579, \"y2\": 153944.17562834496, \"y1\": 1025745.3850359592, \"x\": 24.0, \"y3\": 107708.41896666272}, {\"y6\": 49267.24131164062, \"y4\": 24762.059586635918, \"y5\": 5667.042122240137, \"y2\": 152658.0441487905, \"y1\": 1025736.8464935469, \"x\": 25.2, \"y3\": 106648.71625989168}, {\"y6\": 49147.28473302024, \"y4\": 24836.661016261005, \"y5\": 5735.8481082476665, \"y2\": 151384.5608582416, \"y1\": 1025728.3584689518, \"x\": 26.4, \"y3\": 105600.33016295906}, {\"y6\": 49037.56935907323, \"y4\": 24885.220216551847, \"y5\": 5793.9820554675525, \"y2\": 150123.59086003044, \"y1\": 1025719.9199959171, \"x\": 27.6, \"y3\": 104563.14490377449}, {\"y6\": 48937.239572849365, \"y4\": 24910.115127318815, \"y5\": 5842.385728490161, \"y2\": 148875.0007267741, \"y1\": 1025711.5306633129, \"x\": 28.8, \"y3\": 103537.03879328084}, {\"y6\": 48845.512245649516, \"y4\": 24913.520276281106, \"y5\": 5881.92005931755, \"y2\": 147638.65779427378, \"y1\": 1025703.1898158476, \"x\": 30.0, \"y3\": 102521.89111191881}, {\"y6\": 48761.67082741817, \"y4\": 24897.425765544052, \"y5\": 5913.372257154, \"y2\": 146414.43168995233, \"y1\": 1025694.8971487699, \"x\": 31.2, \"y3\": 101517.58038241926}, {\"y6\": 48685.05925193152, \"y4\": 24863.639897134093, \"y5\": 5937.462388820467, \"y2\": 145202.19289365315, \"y1\": 1025686.6522936783, \"x\": 32.4, \"y3\": 100523.98574838715}, {\"y6\": 48615.07685958917, \"y4\": 24813.82765346074, \"y5\": 5954.84901605262, \"y2\": 144001.8137060491, \"y1\": 1025678.4549801403, \"x\": 33.6, \"y3\": 99540.98732841594}, {\"y6\": 48551.17365536329, \"y4\": 24749.50551856357, \"y5\": 5966.134734840696, \"y2\": 142813.16732957636, \"y1\": 1025670.3048311098, \"x\": 34.8, \"y3\": 98568.46676897487}, {\"y6\": 48492.84603915497, \"y4\": 24672.06351418867, \"y5\": 5971.870796611711, \"y2\": 141636.1289504489, \"y1\": 1025662.2017393912, \"x\": 36.0, \"y3\": 97606.3056486952}, {\"y6\": 48439.63258504552, \"y4\": 24582.770427409563, \"y5\": 5972.561731203333, \"y2\": 140470.57492567264, \"y1\": 1025654.1452582716, \"x\": 37.2, \"y3\": 96654.3888525617}, {\"y6\": 48391.11061677976, \"y4\": 24482.788378361834, \"y5\": 5968.669305613184, \"y2\": 139316.38299695533, \"y1\": 1025646.1350929292, \"x\": 38.4, \"y3\": 95712.60198051616}, {\"y6\": 48346.89282019434, \"y4\": 24373.179467133014, \"y5\": 5960.616123816562, \"y2\": 138173.4323068438, \"y1\": 1025638.1710584968, \"x\": 39.6, \"y3\": 94780.83100733915}, {\"y6\": 48306.624043871365, \"y4\": 24254.91529483372, \"y5\": 5948.789111333126, \"y2\": 137041.60333474012, \"y1\": 1025630.252952958, \"x\": 40.8, \"y3\": 93858.9640455246}, {\"y6\": 48269.978841189724, \"y4\": 24128.885640039665, \"y5\": 5933.542420353773, \"y2\": 135920.7779085668, \"y1\": 1025622.3801860977, \"x\": 42.0, \"y3\": 92946.89002237584}, {\"y6\": 48236.65879940183, \"y4\": 23995.903863754793, \"y5\": 5915.200300237145, \"y2\": 134810.8392866944, \"y1\": 1025614.5527134696, \"x\": 43.2, \"y3\": 92044.50029595652}, {\"y6\": 48206.390248487485, \"y4\": 23856.712621311053, \"y5\": 5894.05962967369, \"y2\": 133711.67179101426, \"y1\": 1025606.7702865965, \"x\": 44.4, \"y3\": 91151.68561463087}, {\"y6\": 48178.92220254351, \"y4\": 23711.992918008433, \"y5\": 5870.392214708368, \"y2\": 132623.16099945578, \"y1\": 1025599.0325127972, \"x\": 45.6, \"y3\": 90268.33968618023}, {\"y6\": 48154.0244166268, \"y4\": 23562.366623612146, \"y5\": 5844.446913035946, \"y2\": 131545.19409720204, \"y1\": 1025591.3392545664, \"x\": 46.8, \"y3\": 89394.35587934082}, {\"y6\": 48131.485720460434, \"y4\": 23408.403071249726, \"y5\": 5816.451574199113, \"y2\": 130477.65947483394, \"y1\": 1025583.6901420479, \"x\": 48.0, \"y3\": 88529.6304338474}, {\"y6\": 48111.112252745166, \"y4\": 23250.62241567926, \"y5\": 5786.614829850334, \"y2\": 129420.446397881, \"y1\": 1025576.0849430693, \"x\": 49.2, \"y3\": 87674.05957923262}, {\"y6\": 48092.726144353, \"y4\": 23089.500056010136, \"y5\": 5755.127689436615, \"y2\": 128373.44575121025, \"y1\": 1025568.5233892206, \"x\": 50.4, \"y3\": 86827.54133598477}, {\"y6\": 48076.16414890803, \"y4\": 22925.471834301778, \"y5\": 5722.16500409589, \"y2\": 127336.54952767416, \"y1\": 1025561.0052422348, \"x\": 51.6, \"y3\": 85989.9743809501}, {\"y6\": 48061.27644409075, \"y4\": 22758.934175506263, \"y5\": 5687.886848715512, \"y2\": 126309.65089908589, \"y1\": 1025553.5302133014, \"x\": 52.8, \"y3\": 85161.25917427897}, {\"y6\": 48047.925409902455, \"y4\": 22590.25066613059, \"y5\": 5652.439773310293, \"y2\": 125292.6443354025, \"y1\": 1025546.097988835, \"x\": 54.0, \"y3\": 84341.29752103704}, {\"y6\": 48035.98473395826, \"y4\": 22419.752717192205, \"y5\": 5615.95788344553, \"y2\": 124285.4251791593, \"y1\": 1025538.7084082512, \"x\": 55.2, \"y3\": 83529.99090336318}, {\"y6\": 48025.338415006976, \"y4\": 22247.74352304974, \"y5\": 5578.563912810221, \"y2\": 123287.89035209756, \"y1\": 1025531.3611467123, \"x\": 56.4, \"y3\": 82727.24368728267}, {\"y6\": 48015.87991225786, \"y4\": 22074.49938282643, \"y5\": 5540.370167235847, \"y2\": 122299.93757007999, \"y1\": 1025524.0560201883, \"x\": 57.6, \"y3\": 81932.9597566808}, {\"y6\": 48007.51136994703, \"y4\": 21900.272394516345, \"y5\": 5501.4793917833495, \"y2\": 121321.46589128455, \"y1\": 1025516.7927291146, \"x\": 58.8, \"y3\": 81147.04530337936}, {\"y6\": 48000.14290555189, \"y4\": 21725.293886230367, \"y5\": 5461.985536083473, \"y2\": 120352.37541545193, \"y1\": 1025509.5709656149, \"x\": 60.0, \"y3\": 80369.40771015563}, {\"y6\": 47993.6919377598, \"y4\": 21549.772557338125, \"y5\": 5421.974542026925, \"y2\": 119392.56734169702, \"y1\": 1025502.3905461323, \"x\": 61.2, \"y3\": 79599.9542868602}, {\"y6\": 47988.08260252381, \"y4\": 21373.899859531062, \"y5\": 5381.524957465708, \"y2\": 118441.94406834751, \"y1\": 1025495.2512494243, \"x\": 62.4, \"y3\": 78838.59385690892}, {\"y6\": 47983.24515415825, \"y4\": 21197.85012372098, \"y5\": 5340.708586655985, \"y2\": 117500.40905962803, \"y1\": 1025488.1527754711, \"x\": 63.6, \"y3\": 78085.23683850546}, {\"y6\": 47979.11552963162, \"y4\": 21021.78165341416, \"y5\": 5299.591013079735, \"y2\": 116567.86677474319, \"y1\": 1025481.0949215452, \"x\": 64.8, \"y3\": 77339.79395369317}, {\"y6\": 47975.63482080916, \"y4\": 20845.838537515363, \"y5\": 5258.232158008264, \"y2\": 115644.22269360976, \"y1\": 1025474.0773602239, \"x\": 66.0, \"y3\": 76602.17792824}, {\"y6\": 47972.74892324408, \"y4\": 20670.152394445187, \"y5\": 5216.6866813329925, \"y2\": 114729.38353399943, \"y1\": 1025467.0999098099, \"x\": 67.2, \"y3\": 75872.30147009873}, {\"y6\": 47970.40805377507, \"y4\": 20494.841168255283, \"y5\": 5175.004506312871, \"y2\": 113823.25686349867, \"y1\": 1025460.1622849536, \"x\": 68.4, \"y3\": 75150.07905370682}, {\"y6\": 47968.56649413087, \"y4\": 20320.013502978443, \"y5\": 5133.231112129675, \"y2\": 112925.7513263502, \"y1\": 1025453.2643209542, \"x\": 69.6, \"y3\": 74435.42504734633}, {\"y6\": 47967.18217860369, \"y4\": 20145.766345108583, \"y5\": 5091.407983048831, \"y2\": 112036.77665519192, \"y1\": 1025446.4057941638, \"x\": 70.8, \"y3\": 73728.25562034325}, {\"y6\": 47966.21643060302, \"y4\": 19972.188377436963, \"y5\": 5049.572903593813, \"y2\": 111156.2435178479, \"y1\": 1025439.5863119647, \"x\": 72.0, \"y3\": 73028.48814419218}, {\"y6\": 47965.63369228028, \"y4\": 19799.359327747294, \"y5\": 5007.760257864902, \"y2\": 110284.0636491715, \"y1\": 1025432.8058786021, \"x\": 73.2, \"y3\": 72336.04032821291}, {\"y6\": 47965.40128039223, \"y4\": 19627.350972322165, \"y5\": 4966.00129776987, \"y2\": 109420.1496428947, \"y1\": 1025426.0639885112, \"x\": 74.4, \"y3\": 71650.83132132988}, {\"y6\": 47965.48911276381, \"y4\": 19456.228953036894, \"y5\": 4924.324456369468, \"y2\": 108564.41511553491, \"y1\": 1025419.3605613795, \"x\": 75.6, \"y3\": 70972.78089977225}, {\"y6\": 47965.86953519827, \"y4\": 19286.050739774775, \"y5\": 4882.755519296006, \"y2\": 107716.77462634257, \"y1\": 1025412.6953144062, \"x\": 76.8, \"y3\": 70301.81012196722}, {\"y6\": 47966.517112580266, \"y4\": 19116.8683564686, \"y5\": 4841.317871527521, \"y2\": 106877.14366739208, \"y1\": 1025406.0680715262, \"x\": 78.0, \"y3\": 69637.84005047343}, {\"y6\": 47967.40844568286, \"y4\": 18948.729476605266, \"y5\": 4800.032697236521, \"y2\": 106045.43872067519, \"y1\": 1025399.4785356295, \"x\": 79.2, \"y3\": 68980.79384617886}, {\"y6\": 47968.52201333792, \"y4\": 18781.674633656385, \"y5\": 4758.91915340107, \"y2\": 105221.57723201418, \"y1\": 1025392.9265945695, \"x\": 80.4, \"y3\": 68330.59407476058}, {\"y6\": 47969.838009889856, \"y4\": 18615.741333072667, \"y5\": 4717.994539961879, \"y2\": 104405.47743404549, \"y1\": 1025386.411937542, \"x\": 81.6, \"y3\": 67687.16598095116}, {\"y6\": 47971.338232452326, \"y4\": 18450.96216836304, \"y5\": 4677.274443126058, \"y2\": 103597.05847424174, \"y1\": 1025379.9342615751, \"x\": 82.8, \"y3\": 67050.43450711449}, {\"y6\": 47973.0059182337, \"y4\": 18287.36605468769, \"y5\": 4636.772898231111, \"y2\": 102796.2405131137, \"y1\": 1025373.493420967, \"x\": 84.0, \"y3\": 66420.32605964647}, {\"y6\": 47974.825657956746, \"y4\": 18124.97796407201, \"y5\": 4596.502496094191, \"y2\": 102002.94449583378, \"y1\": 1025367.0892642827, \"x\": 85.2, \"y3\": 65796.76662383205}, {\"y6\": 47976.78326895834, \"y4\": 17963.820028019505, \"y5\": 4556.474520119292, \"y2\": 101217.09225491666, \"y1\": 1025360.7214903343, \"x\": 86.4, \"y3\": 65179.68442759821}, {\"y6\": 47978.86570616174, \"y4\": 17803.911190518447, \"y5\": 4516.699039694679, \"y2\": 100438.60650053986, \"y1\": 1025354.3899067643, \"x\": 87.6, \"y3\": 64569.00795315626}, {\"y6\": 47981.06096459488, \"y4\": 17645.26773454019, \"y5\": 4477.185018349739, \"y2\": 99667.41077582118, \"y1\": 1025348.0942553151, \"x\": 88.8, \"y3\": 63964.666902524215}, {\"y6\": 47983.358006611306, \"y4\": 17487.903287476183, \"y5\": 4437.940404039772, \"y2\": 98903.42948261327, \"y1\": 1025341.8343631435, \"x\": 90.0, \"y3\": 63366.59107845313}, {\"y6\": 47985.74667197108, \"y4\": 17331.829319888733, \"y5\": 4398.972219633857, \"y2\": 98146.58785194805, \"y1\": 1025335.6099922063, \"x\": 91.2, \"y3\": 62774.711615795444}, {\"y6\": 47988.21761232437, \"y4\": 17177.05518761283, \"y5\": 4360.286630772571, \"y2\": 97396.81189261172, \"y1\": 1025329.4209192178, \"x\": 92.4, \"y3\": 62188.96037703229}, {\"y6\": 47990.76223880254, \"y4\": 17023.588185800363, \"y5\": 4321.889015385415, \"y2\": 96654.028566442, \"y1\": 1025323.2669810322, \"x\": 93.6, \"y3\": 61609.26952546516}, {\"y6\": 47993.37263575731, \"y4\": 16871.434104011165, \"y5\": 4283.784055229662, \"y2\": 95918.16545950208, \"y1\": 1025317.1479004006, \"x\": 94.8, \"y3\": 61035.57282689313}, {\"y6\": 47996.04152984132, \"y4\": 16720.596920291675, \"y5\": 4245.975757578923, \"y2\": 95189.1510770014, \"y1\": 1025311.0635106944, \"x\": 96.0, \"y3\": 60467.80412323293}, {\"y6\": 47998.76223770432, \"y4\": 16571.079272624116, \"y5\": 4208.467534526963, \"y2\": 94466.91471791417, \"y1\": 1025305.0136222941, \"x\": 97.2, \"y3\": 59905.898364448556}, {\"y6\": 48001.528602443286, \"y4\": 16422.882402349653, \"y5\": 4171.2622436633565, \"y2\": 93751.38633572504, \"y1\": 1025298.9979225373, \"x\": 98.4, \"y3\": 59349.791016756106}, {\"y6\": 48004.334966511386, \"y4\": 16276.006389679003, \"y5\": 4134.362250140845, \"y2\": 93042.49680168794, \"y1\": 1025293.0163018588, \"x\": 99.6, \"y3\": 58799.41853599947}, {\"y6\": 48007.17612648905, \"y4\": 16130.450275571822, \"y5\": 4097.769455729955, \"y2\": 92340.17763333325, \"y1\": 1025287.0684998889, \"x\": 100.8, \"y3\": 58254.71822842785}, {\"y6\": 48010.047301637656, \"y4\": 15986.211753061541, \"y5\": 4061.4853314832967, \"y2\": 91644.36115983225, \"y1\": 1025281.154384831, \"x\": 102.0, \"y3\": 57715.62707098667}, {\"y6\": 48012.94409176821, \"y4\": 15843.28809930378, \"y5\": 4025.510980233782, \"y2\": 90954.98042616021, \"y1\": 1025275.2736840596, \"x\": 103.2, \"y3\": 57182.08416979792}, {\"y6\": 48015.86246337032, \"y4\": 15701.675369228578, \"y5\": 3989.84713900076, \"y2\": 90271.96922703848, \"y1\": 1025269.4262399932, \"x\": 104.4, \"y3\": 56654.02833541349}, {\"y6\": 48018.79870146097, \"y4\": 15561.369036687729, \"y5\": 3954.4942310093893, \"y2\": 89595.26207991398, \"y1\": 1025263.6118413581, \"x\": 105.6, \"y3\": 56131.39939738839}, {\"y6\": 48021.749389059754, \"y4\": 15422.363776380551, \"y5\": 3919.4523944413486, \"y2\": 88924.79417449745, \"y1\": 1025257.8302556818, \"x\": 106.8, \"y3\": 55614.138345058294}, {\"y6\": 48024.71140846714, \"y4\": 15284.654209956247, \"y5\": 3884.721481728692, \"y2\": 88260.5015273902, \"y1\": 1025252.0813591847, \"x\": 108.0, \"y3\": 55102.1860313184}, {\"y6\": 48027.68187177452, \"y4\": 15148.233444605268, \"y5\": 3850.3011212483716, \"y2\": 87602.32074554755, \"y1\": 1025246.3648321411, \"x\": 109.2, \"y3\": 54595.48451005448}, {\"y6\": 48030.65813863728, \"y4\": 15013.094712884435, \"y5\": 3816.1907129720285, \"y2\": 86950.18916141632, \"y1\": 1025240.6805715781, \"x\": 110.4, \"y3\": 54093.976335116466}, {\"y6\": 48033.63779565599, \"y4\": 14879.230616536772, \"y5\": 3782.3894515937473, \"y2\": 86304.04480798432, \"y1\": 1025235.0283892527, \"x\": 111.6, \"y3\": 53597.6044585316}, {\"y6\": 48036.61861019281, \"y4\": 14746.633531958534, \"y5\": 3748.896364323512, \"y2\": 85663.82636970213, \"y1\": 1025229.4080335686, \"x\": 112.8, \"y3\": 53106.31320535179}, {\"y6\": 48039.59855691392, \"y4\": 14615.295219036936, \"y5\": 3715.7102915890387, \"y2\": 85029.47322385019, \"y1\": 1025223.8193534106, \"x\": 114.0, \"y3\": 52620.04676556372}, {\"y6\": 48042.57575907242, \"y4\": 14485.207212664653, \"y5\": 3682.8299485955763, \"y2\": 84400.92539413518, \"y1\": 1025218.2621338633, \"x\": 115.2, \"y3\": 52138.750414293856}, {\"y6\": 48045.54851525596, \"y4\": 14356.360951146504, \"y5\": 3650.2538922305002, \"y2\": 83778.12357889734, \"y1\": 1025212.7362275937, \"x\": 116.4, \"y3\": 51662.36948262714}, {\"y6\": 48048.51525990488, \"y4\": 14228.747262549043, \"y5\": 3617.980568183849, \"y2\": 83161.00911194028, \"y1\": 1025207.2414102762, \"x\": 117.6, \"y3\": 51190.850680757314}, {\"y6\": 48051.47456724052, \"y4\": 14102.356955762209, \"y5\": 3586.0083087312796, \"y2\": 82549.52396850112, \"y1\": 1025201.7775099704, \"x\": 118.8, \"y3\": 50724.14078552829}, {\"y6\": 48054.42511449156, \"y4\": 13977.180709497683, \"y5\": 3554.33535895599, \"y2\": 81943.61074874538, \"y1\": 1025196.3443207414, \"x\": 120.0, \"y3\": 50262.18762706345}]}, \"panes\": [{\"type\": \"rectangular\", \"options\": {\"width\": 600, \"height\": 400, \"zoom\": true, \"yrange\": [0.0, 1204000.0], \"legend\": true, \"xrange\": [0.0, 120.0]}, \"diagrams\": [{\"type\": \"line\", \"data\": \"data1\", \"options\": {\"stroke_width\": 2, \"title\": \"EGF(R)\", \"y\": \"y1\", \"x\": \"x\", \"color\": \"#a6cee3\"}}, {\"type\": \"line\", \"data\": \"data1\", \"options\": {\"stroke_width\": 2, \"title\": \"EGFR\", \"y\": \"y2\", \"x\": \"x\", \"color\": \"#1f78b4\"}}, {\"type\": \"line\", \"data\": \"data1\", \"options\": {\"stroke_width\": 2, \"title\": \"EGFR(CR1^_)\", \"y\": \"y3\", \"x\": \"x\", \"color\": \"#b2df8a\"}}, {\"type\": \"line\", \"data\": \"data1\", \"options\": {\"stroke_width\": 2, \"title\": \"EGFR(Y1068=P^_0)\", \"y\": \"y4\", \"x\": \"x\", \"color\": \"#33a02c\"}}, {\"type\": \"line\", \"data\": \"data1\", \"options\": {\"stroke_width\": 2, \"title\": \"EGFR(Y1068^1).Grb2(SH2^1,SH3^2).Sos1(PxxP^2)\", \"y\": \"y5\", \"x\": \"x\", \"color\": \"#e31a1c\"}}, {\"type\": \"line\", \"data\": \"data1\", \"options\": {\"stroke_width\": 2, \"title\": \"Grb2(SH2,SH3^1).Sos1(PxxP^1)\", \"y\": \"y6\", \"x\": \"x\", \"color\": \"#8dd3c7\"}}]}]};\n",
+       "         Nyaplot.core.parse(models, '#viz114766d6-d6c2-4ce5-80f1-f1eb6222259b');\n",
+       "\n",
+       "         var div = d3.select('#viz114766d6-d6c2-4ce5-80f1-f1eb6222259b');\n",
+       "         var svg = div.select('svg');\n",
+       "\n",
+       "         if(true){\n",
+       "           util.svg2uri(svg)\n",
+       "               .then(function(uri){\n",
+       "                 svg.remove();\n",
+       "                 div\n",
+       "                   .append(\"img\")\n",
+       "                   .attr(\"src\", uri);\n",
+       "               });\n",
+       "         }\n",
+       "     \n",
+       "       });\n",
+       "     })();\n",
+       "    </script>\n",
+       "  </head>\n",
+       "  <body>\n",
+       "    <div id=\"viz114766d6-d6c2-4ce5-80f1-f1eb6222259b\"></div>\n",
+       "  </body>\n",
+       "</html>"
+      ],
+      "text/plain": [
+       "<IPython.core.display.HTML object>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "species_list = [\"EGFR\",\n",
+    "                \"EGF(R)\",\n",
+    "                \"EGFR(CR1^_)\",\n",
+    "                \"EGFR(Y1068=P^_0)\",\n",
+    "                \"Grb2(SH2, SH3^1).Sos1(PxxP^1)\",\n",
+    "                \"EGFR(Y1068^1).Grb2(SH2^1, SH3^2).Sos1(PxxP^2)\"]\n",
+    "run_simulation(120, model=newm, y0=y0, species_list=species_list,\n",
+    "               opt_kwargs={'interactive': True, 'to_png': True})"
+   ]
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "Python 3",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.4.3"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 0
+}