5-intermediate-custom-logic.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "d67e0524-e3e0-4010-a437-dbb2a7649412",
   "metadata": {},
   "outputs": [],
   "source": [
    "import dinkum\n",
    "from dinkum.vfg import Gene, Receptor\n",
    "from dinkum.vfn import Tissue"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "e61a38ef-a49c-4513-b358-3701cda4832d",
   "metadata": {},
   "source": [
    "# Custom activation logic\n",
    "\n",
    "dinkum supports custom cis-regulatory logic, in which you can write your arbitrary function of regulatory inputs.\n",
    "\n",
    "You have to write a function of the gene names & then returns a True or a False. See `custom_and`, `custom_or`, and `custom_not`, below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "1bc1fbf2-0f25-4b89-ad25-7562a6b54ecc",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "initializing: dinkum v0.3.2\n",
      "start=1 stop=3\n",
      "custom cis-regulatory logic: and! inputs: gene X=True, gene Y=False => gene active? False\n",
      "custom cis-regulatory logic: and! inputs: gene X=True, gene Y=True => gene active? True\n"
     ]
    },
    {
     "data": {
      "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wBDAAgGBgcGBQgHBwcJCQgKDBQNDAsLDBkSEw8UHRofHh0aHBwgJC4nICIsIxwcKDcpLDAxNDQ0Hyc5PTgyPC4zNDL/2wBDAQkJCQwLDBgNDRgyIRwhMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjL/wAARCACqAL4DASIAAhEBAxEB/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQAAAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWmp6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEAAwEBAQEBAQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAECAxEEBSExBhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRomJygpKjU2Nzg5OkNERUZHSElKU1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9oADAMBAAIRAxEAPwD3+iiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKAOa8ZeebKzXTNv9ufaAdOz0DgHeW/2Nm4H6jHOKv+GjZnw7Z/YRIIQhBEv+sEmT5m/wD29+7d/tZrWooAKKKKAOQ8Vaa82taYkMoji1h/7Nv1x/rIlSSYfjiORPpKfQV1wAUAAAAcADtS0UAFFFFAHE6e3hzUvH011CdKju7GSSCJEMa3E1wQRK5A+Y7RlB7+YfQ121FFABRRRQAUUUUAFFFFABRRRQAUUUUAFFQ3V1BY2c93cyLFbwRtLLI3RVUZJP0AqagAoqhpWtadrcdxJptytxHbzeTI6qQN+1W4JHzDaynIyDnrV+gAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAx/Femyax4Q1nTYYklnubKaKJHxgyFCF68D5sc9qgstEuItWhv0n+z2ioAthmYeUNm0LgTmLj2THHHPNb9FAGRo9jPaanr8sseyO6v0mgOQdyC2gTPHT5kYc+npiteiigDlPGmsX+k/YfsM/leb5m/5FbONuOoPqa5X/AITHXv8An/8A/IKf/E1ufEX/AJhv/bX/ANkrhq+ly+hSnhoylFN69F3Z8JnOMxFPHTjCpJLTRNrojd/4THXv+f8A/wDIKf8AxNavhrxLq+oeILW1urvzIX37l8tBnCEjkDPUVxtbvg7/AJGqy/4H/wCgNW+Jw1FUZtQWz6LscuBx2KliqcZVJNOS6vudLrmuajZ6xPBBcbIl24XYpxlQe496z/8AhJdX/wCfv/yGn+FHiX/kYLr/AIB/6AKya+TP0U1v+El1f/n7/wDIaf4V0HhjUrvUftX2qbzPL2bflAxndnoPYVxNdX4L/wCX7/tn/wCzUAZX/CS6v/z9/wDkNP8ACj/hJdX/AOfv/wAhp/hWTRQBrf8ACS6v/wA/f/kNP8K7DQ7ma80eCed98rbstgDOGI7fSvOa9B8Nf8i/a/8AA/8A0M0AatFFFABRRRQAUUUUAFFFFABRVDXNVj0PQb/VZo5JY7O3edkjBLMFBOBj6deg6niki13SprlbVL+3N4y7haFws/3d2DGfmBxzgjNAGhRWNo2tXOo3l1aXlgLOeCKGbaJvM+WTdhW+UbXGw5XkDI5Oa2aAOG+Iv/MN/wC2v/slcNXq3iLw7/b/ANm/0ryPI3f8s927dj3HpWF/wrr/AKiv/kv/APZV9BgcdQpUIwnKzV+j7nxmbZTjMRjJ1aULxduq7LzOGrd8Hf8AI1WX/A//AEBq3P8AhXX/AFFf/Jf/AOyq9o/gv+ydVhvv7Q83yt3yeTtzlSOu4+ta18ww06Uoxlq0+j7ehz4PJsdTxFOcoaKSb1XR+pleJf8AkYLr/gH/AKAKya7bUvDH9o6hLdfbPL8zHy+VnGAB1z7VU/4Qv/qIf+Qf/sq+aPuzlK6vwX/y/f8AbP8A9mo/4Qv/AKiH/kH/AOyrW0XRf7H8/wD0jzvN2/wbcYz7n1oA8+orq/8AhC/+oh/5B/8AsqP+EL/6iH/kH/7KgDlK9B8Nf8i/a/8AA/8A0M1k/wDCF/8AUQ/8g/8A2VdBptl/Z2nxWvmeZ5efm24zkk9PxoAt0UUUAFFFFABRRRQAUUUUAUda00azoOo6W0piW9tpLcyAZKb1K5x3xmmJpQW4F61xOb3ZyRcTCAttxnyfMKgd8frnmtGigDG8O6Re6Pbzpe31vezTSebJPHbNE8jn7zNl2zwFAAwAAAOMY2aKKAOG+Iv/ADDf+2v/ALJXDV3PxF/5hv8A21/9krhq+ry3/dY/P82fnme/8jCp8v8A0lBW74O/5Gqy/wCB/wDoDVhVu+Dv+Rqsv+B/+gNW+K/gT9H+Rx5f/vdL/FH80HjH/kar3/gH/oC1hVu+Mf8Akar3/gH/AKAtYVGF/gQ9F+QZh/vdX/FL82Fdz8Ov+Yl/2y/9nrhq7n4df8xL/tl/7PWGZf7rL5fmjsyL/kYU/n/6SzCooor5Q/Qwr0Hw1/yL9r/wP/0M159XoPhr/kX7X/gf/oZoA1aKKKACiiigAooooAKKKKACisjxVqF1pXhLWNRshGbm1spZ4/MPygqhOTwc9Onfpx1qqfFMA1NtJmENverGS7m7geONhHvJ2+YJCo9dg454HNAHQ0Vznh661L+2tW07ULmacWqQNG1wkSu+7eC6iMY8slPlB+bKtntXR0AcN8Rf+Yb/ANtf/ZK4avVvEX9g/wCjf23/ALXk/f8Abd938OtYX/FB/wCfPr6DA4r2dCMeST32Wm58Zm2Xe2xk6ntYRvbRys9l0scNW74O/wCRqsv+B/8AoDVuf8UH/nz6vaP/AMIl/asP9mf8fnzeX/rf7pz97jpmta+M5qUo+zkrp9PI58HlnJiKcvbU3aS0Utd+mm5yvjH/AJGq9/4B/wCgLWFXo+sf8Il/as39p/8AH58vmf63+6Mfd46Yqj/xQf8Anz6KGM5aUY+zk7JdPIMZlnPiKkvbU1eT0ctd+um5w1dz8Ov+Yl/2y/8AZ6P+KD/z59bvh3+wf9J/sT/Z877/AL7fvfj0rLHYr2lCUeSS23Wm50ZTl3scZCp7WErX0UrvZ9LHC0V0n2zwV/z0/Saj7Z4K/wCen6TV431Wv/I/uZ9R/aGE/wCfsf8AwJf5nN16D4a/5F+1/wCB/wDoZrC+2eCv+en6TV02lPZSabC2nHNod2w/N6nP3ueuaidGpBXnFr1RpSxVCq+WnNSfk0y5RRRWZuFFFFABRRRQAUUUUAQXlpBqFjcWV1GJbe4jaKVCSAyMMEceoJpfssJs/sjqZIDH5TLKxfcuMYYtktkdc9amooAoabo1jpPmfZI5A0gUO8szyuQv3Ruck4GTgZwMn1q/RRQBw3xF/wCYb/21/wDZK4au5+Iv/MN/7a/+yVw1fV5b/usfn+bPzzPf+RhU+X/pKCt3wd/yNVl/wP8A9AasKt3wd/yNVl/wP/0Bq3xX8Cfo/wAjjy//AHul/ij+aDxj/wAjVe/8A/8AQFrCrd8Y/wDI1Xv/AAD/ANAWsKjC/wACHovyDMP97q/4pfmwrufh1/zEv+2X/s9cNXc/Dr/mJf8AbL/2esMy/wB1l8vzR2ZF/wAjCn8//SWcNRRRXceQFereDv8AkVbL/gf/AKG1eU16t4O/5FWy/wCB/wDobV5WcfwF6/oz6Lhn/e5f4X+aNyiiivmz7gKKKKACiiigAooooAKKKKACiiigDhviL/zDf+2v/slcNXuVFerhsz9hSVPkvbz/AOAfO47IPreIlX9pa9tLX2SXc8Nrd8Hf8jVZf8D/APQGr1WitKub+0hKHJumt/8AgGWH4b9jWjV9rflafw9n6nlXjH/kar3/AIB/6AtYVe5UUUs39nCMOTZJb/8AADEcN+2rSq+1tzNv4e79Tw2u5+HX/MS/7Zf+z13NFZ4nM/b0nT5LX8/+Aa4HIPqmIjX9pe19LW3TXc8Nor3Kit/7a/ufj/wDk/1W/wCnv/kv/BPDa9W8Hf8AIq2X/A//AENq3KK5cZmH1mChy21vv/wD0csyb6jWdXn5rq21uq832CiiivNPcCiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooArf2jY+ZdR/bLffZgG5XzVzACu4b+fl4557c1LBPDdW8VxbypNBKgeOSNgyupGQQRwQR3rz280rUP+Es13yrK4NrrV5HaTyCM7TEsFsdxPQLt+0ru6biB1rrPB0E1r4I0C3uInimi023SSORSrIwjUEEHkEHtQBtUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFAH//2Q==",
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAL4AAACqCAIAAACcUWJHAAAI30lEQVR4Ae2dX2iOXxzAvcwWE4vccT2JEfIn1AqbzJ9Im2L+FLmxorlYKUphUVKuXCgkyp9l2i62xEhcUFyYf7XSliSU/242fl+/1esx55w953He1/Oc5/Ne6LznfM95zvl8P877vPKcN/Pjx48hvCBgT2CofRd6QOAnAdTBg4gEUCciOLqhDg5EJIA6EcHRDXVwICIB1IkIjm6ogwMRCaBORHB0Qx0ciEgAdSKCoxvq4EBEAqgTERzdUAcHIhJAnYjg6IY6Fg60tbWtWbOmv0NnZ+fs2bP7+vos+utDHz9+PHXq1O/fv/eHzJo169GjR/pwu5b6+vry/18LFy4sLCy062yIlv/qxSs8gaqqqps3b0r8kiVL7t69G77joJHr1q07f/68hLW0tFRXVw8aHyHg2LFjDQ0NEToquwxR1lKpI/D06dO5c+deunRp69atupho9bLxTJkyRbax+fPny5YWbRBDr5cvX5aVlX358sUQY9WUkWjDnkTTnwR27tzZ1NT04MGDcePG/dn6NzXr168fOXKkZPfcuXN/M46yb01NzcaNG2XXVLZGqCyI0CflXT5+/FhQUPD582fn6uzdu3fy5Mmy/TgnLHdpsp859ObnDK32KIJv374td8rt7e0rV67MBY0xY8Y4H/bbt2/Tp0/v6elxOzLfsCz+hvf29u7atevIkSNyjywbT3Nzs0Xnfxd64MCB2traCRMmuJ0C9zoWPI8ePfr27duDBw9Kn+7u7srKyvv37xcXF1sMMVhoSUnJ+/fvB4uyay8qKpJ/Rxg2bFh/N/kGN2rUKLshVNGoo6JCXQgCfGCFgESIigDqqKhQF4IA6oSARIiKAOqoqFAXggDqhIBEiIoA6qioUBeCAOqEgESIigDqqKhQF4IA6oSARIiKAOqoqFAXggDqhIBEiIoA6qioUBeCAOqEgESIigDqqKhQF4IA6oSARIiKAOqoqFAXggDqhIBEiIoA6qioUBeCAOqEgESIigDqqKhQF4IA6oSARIiKgJ/qNDY2qhb7s04eVdm8ebM8KacLoD4kgdSps2LFipkzZ2YymZCACNMR8FCdffv2yQPhFRUVyjVfvHixrq5O2USlFQE/H+Hrf4by2bNn27dvz+Lo6OjIlnPxkGV28JQUfFbHkMJo6mQ6DENaN/0o/9XlwoULv978dUlOdsqOkbuRfT4kxbDrZMlSiEzAT3XkUD55lZaWBj+kIjOio5KAh7fJsk45b1HOv1EumEpXBPzcdVpbW82AnB9EYr6cl61+7jpepipui0KduGUkMfNBncSkKm4TRZ24ZSQx80GdxKQqbhNFnbhlJDHzQZ3EpCpuE0WduGUkMfNBncSkKm4TRZ24ZSQx80GdxKQqbhNFnbhlJDHzQZ3EpCpuE0WduGUkMfNBncSkKm4TRZ24ZSQx80GdxKQqbhNFnbhlJDHzQZ3EpCpuE0WduGUkMfNBncSkKnYTdfsDxjEZ7dChQ7qZnDx5csGCBdOmTZOf/tbFUB+GQLoeHH7z5s3atWtv3Ljx/Pnz1atXP3nyJHZ/lZMzIQ8/sAwnXbx7927Hjh1Dhw6dOHGilJOTpjjO1Oddx/DM+enTp2/duiUfXnHMSULm5LM6uhR0dXWtWrVKPrbGjx+vi1HWZzqU1RErgydd5G5kTrqIkh7lriOnNtXU1Mh+Y+tNlBl43cfPZ851J13IF4dNmzbV19fPmTPH67TmY3Ee3iYLNt1JF6dOnZLv5CdOnCgvL1++fHk+APt7DT/vdXKUr9zdkeRu5Nzd6/i56+RIHYYNEkCdIA3KFgRQxwIWoUECqBOkQdmCAOpYwCI0SAB1gjQoWxBAHQtYhAYJoE6QBmULAqhjAYvQIAHUCdKgbEEAdSxgERokgDpBGpQtCKCOBSxCgwRQJ0iDsgUB1LGARWiQAOoEaVC2IIA6FrAIDRJAnSANyhYEUMcCFqFBAqgTpEHZggDqWMAiNEjAzyciGhsbGxoaguvsL3/9+nXLli2vX7/+9OnT/v37q6qq/oyhJiQBP9UpKSlR/jDs4cOH5ayC3bt3v3r1at68eS9evAiJibA/CXj49Gf2pIv29vYBC962bVtxcbFUdnZ2Dh8+fEArb60I+LzrKJ85Fzq1tbVNTU1Xr15dtGiRFSyCgwR8Vie4zgHl5ubmy5cvnzlzZkC9+W2mw9xu15qfky5yN2efv2HJriPPlmdfkti6urre3l4pyAPnra2tdqkm+ncCHt7ryAJ1J118+PDhypUrcibcnTt3SktLf0fBOzsCfn5gybfuTCbT0tIyAEZ3d7cckiJiFRYWHj9+fNKkSQMCzG9zt/kncmQ5csbMi9YsgUQmuCM7fQeF4P2Zz/c6DlAxhJ4A6ujZ0GIkgDpGPDTqCaCOng0tRgKoY8RDo54A6ujZ0GIkgDpGPDTqCaCOng0tRgKoY8RDo54A6ujZ0GIkgDpGPDTqCaCOng0tRgKoY8RDo54A6ujZ0GIkgDpGPDTqCaCOng0tRgKoY8RDo54A6ujZ0GIkgDpGPDTqCaCOng0tRgKoY8RDo56Anw/T6NdLizMC7DrOUKZtINRJW8adrRd1nKFM20Cok7aMO1sv6jhDmbaBUCdtGXe2XtRxhjJtA6FO2jLubL2o4wxl2gZCnbRl3Nl6UccZyrQNhDppy7iz9aKOM5RpGwh10pZxZ+tFHWco8znQiBEjqqurs1fcsGGD1GTf5qeAOvnh7PgqRUVFcqB4X1+fjCtn0HZ1dUmN42sMNhzqDEYoru0zZsy4d++ezO7hw4dlZWX5nybq5J+5mytWVla2tbXJWPKnlN0MajMK6tjQilNsRUXFtWvXZEbXr19fvHhx/qeGOvln7uaKY8eOlZ+F6+npkeFGjx7tZlCbUVDHhlbMYpcuXbpnz55/suUICdSJmQ4205FfUTl79uyyZctsOjmL5WEaZyjTNhC7Ttoy7my9qOMMZdoGQp20ZdzZelHHGcq0DYQ6acu4s/WijjOUaRsIddKWcWfrRR1nKNM2EOqkLePO1os6zlCmbaD/AEi/g1OPzmj7AAAAAElFTkSuQmCC",
      "text/plain": [
       "<PIL.Image.Image image mode=RGB size=190x170>"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dinkum.reset()\n",
    "\n",
    "def custom_and(X, Y):\n",
    "    is_active = (X and Y)\n",
    "    print(f\"custom cis-regulatory logic: and! inputs: gene X={X}, gene Y={Y} => gene active? {is_active}\")\n",
    "    return is_active\n",
    "\n",
    "x = Gene(name='X')\n",
    "y = Gene(name='Y')\n",
    "z = Gene(name='Z')\n",
    "m = Tissue(name='M')\n",
    "\n",
    "x.is_present(where=m, start=1)\n",
    "y.is_present(where=m, start=2)\n",
    "z.custom_activation(state_fn=custom_and, delay=1)\n",
    "\n",
    "dinkum.run_and_display(stop=3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "86ba9862-51c0-4941-b0fa-94f4ba9dee1c",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "initializing: dinkum v0.3.2\n",
      "start=1 stop=3\n",
      "custom cis-regulatory logic: or! inputs: gene X=True, gene Y=False => gene active? True\n",
      "custom cis-regulatory logic: or! inputs: gene X=True, gene Y=True => gene active? True\n"
     ]
    },
    {
     "data": {
      "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wBDAAgGBgcGBQgHBwcJCQgKDBQNDAsLDBkSEw8UHRofHh0aHBwgJC4nICIsIxwcKDcpLDAxNDQ0Hyc5PTgyPC4zNDL/2wBDAQkJCQwLDBgNDRgyIRwhMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjL/wAARCACqAL4DASIAAhEBAxEB/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQAAAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWmp6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEAAwEBAQEBAQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAECAxEEBSExBhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRomJygpKjU2Nzg5OkNERUZHSElKU1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9oADAMBAAIRAxEAPwD3+iiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKAOa8ZeebKzXTNv9ufaAdOz0DgHeW/2Nm4H6jHOKv+GjZnw7Z/YRIIQhBEv+sEmT5m/wD29+7d/tZrWooAKKKKAOQ8Vaa82taYkMoji1h/7Nv1x/rIlSSYfjiORPpKfQV1wAUAAAAcADtS0UAFFFFAHE6e3hzUvH011CdKju7GSSCJEMa3E1wQRK5A+Y7RlB7+YfQ121FFABRRRQAUUUUAFFFFABRRRQAUUUUAFFQ3V1BY2c93cyLFbwRtLLI3RVUZJP0AqagAoqhpWtadrcdxJptytxHbzeTI6qQN+1W4JHzDaynIyDnrV+gAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAx/Femyax4Q1nTYYklnubKaKJHxgyFCF68D5sc9qgstEuItWhv0n+z2ioAthmYeUNm0LgTmLj2THHHPNb9FAGRo9jPaanr8sseyO6v0mgOQdyC2gTPHT5kYc+npiteiigDlPGmsX+k/YfsM/leb5m/5FbONuOoPqa5X/AITHXv8An/8A/IKf/E1ufEX/AJhv/bX/ANkrhq+ly+hSnhoylFN69F3Z8JnOMxFPHTjCpJLTRNrojd/4THXv+f8A/wDIKf8AxNavhrxLq+oeILW1urvzIX37l8tBnCEjkDPUVxtbvg7/AJGqy/4H/wCgNW+Jw1FUZtQWz6LscuBx2KliqcZVJNOS6vudLrmuajZ6xPBBcbIl24XYpxlQe496z/8AhJdX/wCfv/yGn+FHiX/kYLr/AIB/6AKya+TP0U1v+El1f/n7/wDIaf4V0HhjUrvUftX2qbzPL2bflAxndnoPYVxNdX4L/wCX7/tn/wCzUAZX/CS6v/z9/wDkNP8ACj/hJdX/AOfv/wAhp/hWTRQBrf8ACS6v/wA/f/kNP8K7DQ7ma80eCed98rbstgDOGI7fSvOa9B8Nf8i/a/8AA/8A0M0AatFFFABRRRQAUUUUAFFFFABRVDXNVj0PQb/VZo5JY7O3edkjBLMFBOBj6deg6niki13SprlbVL+3N4y7haFws/3d2DGfmBxzgjNAGhRWNo2tXOo3l1aXlgLOeCKGbaJvM+WTdhW+UbXGw5XkDI5Oa2aAOG+Iv/MN/wC2v/slcNXq3iLw7/b/ANm/0ryPI3f8s927dj3HpWF/wrr/AKiv/kv/APZV9BgcdQpUIwnKzV+j7nxmbZTjMRjJ1aULxduq7LzOGrd8Hf8AI1WX/A//AEBq3P8AhXX/AFFf/Jf/AOyq9o/gv+ydVhvv7Q83yt3yeTtzlSOu4+ta18ww06Uoxlq0+j7ehz4PJsdTxFOcoaKSb1XR+pleJf8AkYLr/gH/AKAKya7bUvDH9o6hLdfbPL8zHy+VnGAB1z7VU/4Qv/qIf+Qf/sq+aPuzlK6vwX/y/f8AbP8A9mo/4Qv/AKiH/kH/AOyrW0XRf7H8/wD0jzvN2/wbcYz7n1oA8+orq/8AhC/+oh/5B/8AsqP+EL/6iH/kH/7KgDlK9B8Nf8i/a/8AA/8A0M1k/wDCF/8AUQ/8g/8A2VdBptl/Z2nxWvmeZ5efm24zkk9PxoAt0UUUAFFFFABRRRQAUUUUAUda00azoOo6W0piW9tpLcyAZKb1K5x3xmmJpQW4F61xOb3ZyRcTCAttxnyfMKgd8frnmtGigDG8O6Re6Pbzpe31vezTSebJPHbNE8jn7zNl2zwFAAwAAAOMY2aKKAOG+Iv/ADDf+2v/ALJXDV3PxF/5hv8A21/9krhq+ry3/dY/P82fnme/8jCp8v8A0lBW74O/5Gqy/wCB/wDoDVhVu+Dv+Rqsv+B/+gNW+K/gT9H+Rx5f/vdL/FH80HjH/kar3/gH/oC1hVu+Mf8Akar3/gH/AKAtYVGF/gQ9F+QZh/vdX/FL82Fdz8Ov+Yl/2y/9nrhq7n4df8xL/tl/7PWGZf7rL5fmjsyL/kYU/n/6SzhqKKK7jyAr1bwd/wAirZf8D/8AQ2rymvVvB3/Iq2X/AAP/ANDavKzj+AvX9GfRcM/73L/C/wA0blFFFfNn3AUUUUAFFFFABRRRQAUVkeKtQutK8JaxqNkIzc2tlLPH5h+UFUJyeDnp079OOtVT4pgGptpMwht71YyXc3cDxxsI95O3zBIVHrsHHPA5oA6Giuc8PXWpf21q2nahczTi1SBo2uEiV33bwXURjHlkp8oPzZVs9q6OgDhviL/zDf8Atr/7JXDV6t4i/sH/AEb+2/8Aa8n7/tu+7+HWsL/ig/8APn19BgcV7OhGPJJ77LTc+MzbLvbYydT2sI3to5Wey6WOGrd8Hf8AI1WX/A//AEBq3P8Aig/8+fV7R/8AhEv7Vh/sz/j8+by/9b/dOfvcdM1rXxnNSlH2cldPp5HPg8s5MRTl7am7SWilrv003OV8Y/8AI1Xv/AP/AEBawq9H1j/hEv7Vm/tP/j8+XzP9b/dGPu8dMVR/4oP/AD59FDGctKMfZydkunkGMyznxFSXtqavJ6OWu/XTc4au5+HX/MS/7Zf+z0f8UH/nz63fDv8AYP8ApP8AYn+z533/AH2/e/HpWWOxXtKEo8kltutNzoynLvY4yFT2sJWvopXez6WPKaK7n/ig/wDPn0f8UH/nz66/r3/Tuf3f8E8/+yf+n9P/AMC/4Bw1ereDv+RVsv8Agf8A6G1YX/FB/wCfPrq9H+wf2VD/AGZ/x5/N5f3v7xz97nrmvPzLE+1pKPJJa9Vboz2cjwP1fESl7SEvd2i7vdfgXqKKK8Q+qCiiigAooooAKKKKAILy0g1CxuLK6jEtvcRtFKhJAZGGCOPUE0v2WE2f2R1MkBj8pllYvuXGMMWyWyOuetTUUAUNN0ax0nzPskcgaQKHeWZ5XIX7o3OScDJwM4GT61foooA4b4i/8w3/ALa/+yVw1dz8Rf8AmG/9tf8A2SuGr6vLf91j8/zZ+eZ7/wAjCp8v/SUFbvg7/karL/gf/oDVhVu+Dv8AkarL/gf/AKA1b4r+BP0f5HHl/wDvdL/FH80HjH/kar3/AIB/6AtYVbvjH/kar3/gH/oC1hUYX+BD0X5BmH+91f8AFL82Fdz8Ov8AmJf9sv8A2euGrufh1/zEv+2X/s9YZl/usvl+aOzIv+RhT+f/AKSzhqKKK7jyAr1bwd/yKtl/wP8A9DavKa9W8Hf8irZf8D/9DavKzj+AvX9GfRcM/wC9y/wv80blFFFfNn3AUUUUAFFFFABRRRQAUUUUAFFFFAHDfEX/AJhv/bX/ANkrhq9yor1cNmfsKSp8l7ef/APncdkH1vESr+0te2lr7JLueG1u+Dv+Rqsv+B/+gNXqtFaVc39pCUOTdNb/APAMsPw37GtGr7W/K0/h7P1PKvGP/I1Xv/AP/QFrCr3Kiilm/s4RhybJLf8A4AYjhv21aVX2tuZt/D3fqeG13Pw6/wCYl/2y/wDZ67mis8Tmft6Tp8lr+f8AwDXA5B9UxEa/tL2vpa26a7nhtFe5UVv/AG1/c/H/AIByf6rf9Pf/ACX/AIJ4bXq3g7/kVbL/AIH/AOhtW5RXLjMw+swUOW2t9/8AgHo5Zk31Gs6vPzXVtrdV5vsFFFFeae4FFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQBW/tGx8y6j+2W++zANyvmrmAFdw38/Lxzz25qWCeG6t4ri3lSaCVA8ckbBldSMggjggjvXnt5pWof8ACWa75VlcG11q8jtJ5BGdpiWC2O4noF2/aV3dNxA611ng6Ca18EaBb3ETxTRabbpJHIpVkYRqCCDyCD2oA2qKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiigAAAAAAKAP/9k=",
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAL4AAACqCAIAAACcUWJHAAAI0UlEQVR4Ae2dW4hNXxjAnTFmYsRE3ngeiSHkEmoKMzIukWYU41LkxRSNhylFKUyUlCcPColymYxmHmYSQ+KB4sG41ZRmkoRy9zLD//OfOrY96+xZa9v72Hut33nQ2mt/6ztr/b6fdfaRvU/m58+fw3hBwJxAgfkQRkDgFwHUwYOQBFAnJDiGoQ4OhCSAOiHBMQx1cCAkAdQJCY5hqIMDIQmgTkhwDEMdHAhJAHVCgmMY6uBASAKoExIcw1AHB0ISQJ2Q4BiGOgYOtLe3r127dmBAV1fXnDlz+vv7DcbnDn3y5Mm0adN+/PgxEDJ79uzHjx/nDjc709DQUPH/a9GiRUVFRWaDA6Llv3rx0idQXV1969YtiV+6dOm9e/f0Bw4ZuX79+gsXLkhYa2trTU3NkPEhAo4fP97Y2BhioHLIMGUvnbkIPHv2bN68eZcvX962bVuumHD9svFMnTpVtrEFCxbIlhYuScCoV69elZeXf/36NSDG6FRGogP2JE4NJrBr167m5uaHDx+OHz9+8Nm/6dmwYcOoUaOkuufPn/+bPMqxtbW1mzZtkl1TeTZEZ2GIMY4P+fTpU2Fh4ZcvXyJXZ9++fVOmTJHtJ3LCcpUm+1mE3vyaodEeRfCdO3fkSrmjo2PVqlVx0Bg7dmzkab9//z5jxoze3t5oM/MNy+BveF9f3+7du48ePSrXyLLxtLS0GAz+d6EHDx6sq6ubOHFitFPgWseA57Fjx969e3fo0CEZ09PTU1VV9eDBg5KSEoMUQ4WWlpZ++PBhqCiz88XFxfLvCMOHDx8YJt/gRo8ebZZCFY06Kir0aRDgA0sDEiEqAqijokKfBgHU0YBEiIoA6qio0KdBAHU0IBGiIoA6Kir0aRBAHQ1IhKgIoI6KCn0aBFBHAxIhKgKoo6JCnwYB1NGARIiKAOqoqNCnQQB1NCARoiKAOioq9GkQQB0NSISoCKCOigp9GgRQRwMSISoCqKOiQp8GAdTRgESIigDqqKjQp0EAdTQgEaIiYKc6TU1NqsX+6pNbVbZs2SJ3yuUKoF+TgHPqrFy5ctasWZlMRhMQYbkIWKjO/v375YbwyspK5ZovXbpUX1+vPEWnEQE7b+EbuIfy+fPnO3bsyOLo7OzMtuO4yTKb3JGGzeoElDCcOpnOgJTGp35W/B5y8eLF3wd/3ZInO2VzxJfZ5oekBOw6WbI0QhOwUx15KJ+8ysrKvB9SoRkxUEnAwstkWac8b1Gef6NcMJ1REbBz12lrawsGFPmDSILfzsqzdu46VpYqaYtCnaRVJDXzQZ3UlCppE0WdpFUkNfNBndSUKmkTRZ2kVSQ180Gd1JQqaRNFnaRVJDXzQZ3UlCppE0WdpFUkNfNBndSUKmkTRZ2kVSQ180Gd1JQqaRNFnaRVJDXzQZ3UlCppE0WdpFUkNfNBndSUKmkTRZ2kVSQ180Gd1JQqaRNFnaRVJDXzQZ3UlCpxE432B4wTku3w4cO5ZnLq1KmFCxdOnz5dfvo7Vwz9OgTcunH47du369atu3nz5osXL9asWfP06dPE/VVOz4Qs/MAKeNLF+/fvd+7cWVBQMGnSJGmnp0xJnKnNu07APednzpy5ffu2fHglsSYpmZPN6uQqQXd39+rVq+Vja8KECblilP2ZTmV3yE7vky7SmNnOG4cHiqncdeSpTbW1tbLfmHoTUhB7h9mpTq4nXcgXh82bNzc0NMydO9femuZpZRZeJgu5XE+6OH36tHwnP3nyZEVFxYoVK/LE2NK3sfNaJ6ZipfGKJL4527nrxKQOab0EUMdLg7YBAdQxgEWolwDqeGnQNiCAOgawCPUSQB0vDdoGBFDHABahXgKo46VB24AA6hjAItRLAHW8NGgbEEAdA1iEegmgjpcGbQMCqGMAi1AvAdTx0qBtQAB1DGAR6iWAOl4atA0IoI4BLEK9BFDHS4O2AQHUMYBFqJcA6nhp0DYggDoGsAj1ErDzjoimpqbGxkbvOgfa375927p165s3bz5//nzgwIHq6urBMfRoErBTndLSUuUPwx45ckSeVbBnz57Xr1/Pnz//5cuXmpgIG0zAwrs/s0+66Ojo8C14+/btJSUl0tnV1TVixAjfWQ6NCNi86yjvORc6dXV1zc3N165dW7x4sREsgr0EbFbHu05fu6Wl5cqVK2fPnvX1Bx9mOoPPm51N+5MubP6GJbuO3FuefUlh6+vr+/r6pCE3nLe1tZmVmug/CVh4rSMLzPWki48fP169elWeCXf37t2ysrI/UXBkRsDODyz51p3JZFpbW30wenp65CEpIlZRUdGJEycmT57sCwg+5APrjw9ZeeRMMC/OZgmgjlcdm691siWnEQcB1ImDqhM5UceJMsexSNSJg6oTOVHHiTLHsUjUiYOqEzlRx4kyx7FI1ImDqhM5UceJMsexSNSJg6oTOVHHiTLHsUjUiYOqEzlRx4kyx7FI1ImDqhM5UceJMsexSNSJg6oTOVHHiTLHsUjUiYOqEzlRx4kyx7FI1ImDqhM5UceJMsexSDtvpomDFDl9BNh1fEA41CWAOrqkiPMRQB0fEA51CaCOLinifARQxweEQ10CqKNLijgfAdTxAeFQlwDq6JIizkcAdXxAONQlgDq6pIjzEUAdHxAOdQmgji4p4nwEUMcHhENdAqijS4o4HwHU8QFJx+HIkSNramqyc924caP0ZA/z00Cd/HCO+F2Ki4vlgeL9/f2SV55B293dLT0Rv8dQ6VBnKEJJPT9z5sz79+/L7B49elReXp7/aaJO/plH845VVVXt7e2SS/6UdjRJTbKgjgmtJMVWVlZev35dZnTjxo0lS5bkf2qok3/m0bzjuHHj5Gfhent7Jd2YMWOiSWqSBXVMaCUsdtmyZXv37v0nW46QQJ2E6WAyHfkVlXPnzi1fvtxkUGSx3EwTGUrXErHruFbxyNaLOpGhdC0R6rhW8cjWizqRoXQtEeq4VvHI1os6kaF0LRHquFbxyNaLOpGhdC0R6rhW8cjWizqRoXQt0X9CyIFTJaErswAAAABJRU5ErkJggg==",
      "text/plain": [
       "<PIL.Image.Image image mode=RGB size=190x170>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dinkum.reset()\n",
    "\n",
    "def custom_or(X, Y):\n",
    "    is_active = (X or Y)\n",
    "    print(f\"custom cis-regulatory logic: or! inputs: gene X={X}, gene Y={Y} => gene active? {is_active}\")\n",
    "    return is_active\n",
    "\n",
    "x = Gene(name='X')\n",
    "y = Gene(name='Y')\n",
    "z = Gene(name='Z')\n",
    "m = Tissue(name='M')\n",
    "\n",
    "x.is_present(where=m, start=1)\n",
    "y.is_present(where=m, start=2)\n",
    "z.custom_activation(state_fn=custom_or, delay=1)\n",
    "\n",
    "dinkum.run_and_display(stop=3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "664964a5-baf1-4d5d-95e2-2776aa02eca4",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "initializing: dinkum v0.3.2\n",
      "start=1 stop=3\n",
      "custom cis-regulatory logic: not! inputs: gene X=False => gene active? True\n",
      "custom cis-regulatory logic: not! inputs: gene X=True => gene active? False\n"
     ]
    },
    {
     "data": {
      "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wBDAAgGBgcGBQgHBwcJCQgKDBQNDAsLDBkSEw8UHRofHh0aHBwgJC4nICIsIxwcKDcpLDAxNDQ0Hyc5PTgyPC4zNDL/2wBDAQkJCQwLDBgNDRgyIRwhMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjL/wAARCACqAKADASIAAhEBAxEB/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQAAAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWmp6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEAAwEBAQEBAQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAECAxEEBSExBhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRomJygpKjU2Nzg5OkNERUZHSElKU1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9oADAMBAAIRAxEAPwD3+iiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigDmvGXnmys10zb/bn2gHTs9A4B3lv9jZuB+oxzir/ho2Z8O2f2ESCEIQRL/rBJk+Zv8A9vfu3f7Wa1qKACiiigDkPFWmvNrWmJDKI4tYf+zb9cf6yJUkmH44jkT6Sn0FdcAFAAAAHAA7UtFABRRRQAUUUUAFFFFABRRRQAUUUUAFFQ3V1BY2c93cyLFbwRtLLI3RVUZJP0AqagAoqhpWtadrcdxJptytxHbzeTI6qQN+1W4JHzDaynIyDnrV+gAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKAMfxXpsmseENZ02GJJZ7mymiiR8YMhQhevA+bHPaoLLRLiLVob9J/s9oqALYZmHlDZtC4E5i49kxxxzzW/RQBkaPYz2mp6/LLHsjur9JoDkHcgtoEzx0+ZGHPp6YrXoooA5/xPqV3p32X7LN5fmb93yg5xtx1Hua5/8A4SXV/wDn7/8AIaf4Vq+NP+XH/tp/7LXKUAa3/CS6v/z9/wDkNP8ACtDQ9c1G81iCCe43xNuyuxRnCk9h7VzNa3hr/kYLX/gf/oBoA0Nc1zUbPWJ4ILjZEu3C7FOMqD3HvWf/AMJLq/8Az9/+Q0/wo8S/8jBdf8A/9AFZNAGt/wAJLq//AD9/+Q0/wroPDGpXeo/avtU3meXs2/KBjO7PQewria6vwX/y/f8AbP8A9moA6uiiigAooooAKKKKACiiigAoqhrmqx6HoN/qs0cksdnbvOyRglmCgnAx9OvQdTxSRa7pU1ytql/bm8ZdwtC4Wf7u7BjPzA45wRmgDQorG0bWrnUby6tLywFnPBFDNtE3mfLJuwrfKNrjYcryBkcnNbNAHKeNP+XH/tp/7LXKV6DrWi/2x5H+keT5W7+DdnOPcelZP/CF/wDUQ/8AIP8A9lQByla3hr/kYLX/AIH/AOgGtX/hC/8AqIf+Qf8A7Krem+GP7O1CK6+2eZ5efl8rGcgjrn3oA5/xL/yMF1/wD/0AVk122peGP7R1CW6+2eX5mPl8rOMADrn2qp/whf8A1EP/ACD/APZUAcpXV+C/+X7/ALZ/+zUf8IX/ANRD/wAg/wD2Va2i6L/Y/n/6R53m7f4NuMZ9z60AatFFFABRRRQAUUUUAFFFFAFHWtNGs6DqOltKYlvbaS3MgGSm9Sucd8ZpiaUFuBetcTm92ckXEwgLbcZ8nzCoHfH655rRooAxvDukXuj286Xt9b3s00nmyTx2zRPI5+8zZds8BQAMAAADjGNmiigDhviL/wAw3/tr/wCyVw1dz8Rf+Yb/ANtf/ZK4avq8t/3WPz/Nn55nv/IwqfL/ANJQVu+Dv+Rqsv8Agf8A6A1YVbvg7/karL/gf/oDVviv4E/R/kceX/73S/xR/NB4x/5Gq9/4B/6AtYVbvjH/AJGq9/4B/wCgLWFRhf4EPRfkGYf73V/xS/NhXc/Dr/mJf9sv/Z64au5+HX/MS/7Zf+z1hmX+6y+X5o7Mi/5GFP5/+ks7miiivlD9DCiiigAooooAKKKKACisjxVqF1pXhLWNRshGbm1spZ4/MPygqhOTwc9Onfpx1qqfFMA1NtJmENverGS7m7geONhHvJ2+YJCo9dg454HNAHQ0Vznh661L+2tW07ULmacWqQNG1wkSu+7eC6iMY8slPlB+bKtntXR0AcN8Rf8AmG/9tf8A2SuGr1bxF/YP+jf23/teT9/23fd/DrWF/wAUH/nz6+gwOK9nQjHkk99lpufGZtl3tsZOp7WEb20crPZdLHDVu+Dv+Rqsv+B/+gNW5/xQf+fPq9o//CJf2rD/AGZ/x+fN5f8Arf7pz97jpmta+M5qUo+zkrp9PI58HlnJiKcvbU3aS0Utd+mm5yvjH/kar3/gH/oC1hV6bqtt4Xk1OZtRTN2du85l9Bj7vHTFUvsfgr/nn+s1Y0c1owpxg09El0/zOrFcPYqrXnUjKNpNvd9X6Hn9dz8Ov+Yl/wBsv/Z6n+x+Cv8Ann+s1bOgw6LF9o/sdcZ2+by/vj7349KzxeZUq1GVOKd36d/U3y3I8RhcVGtOUbK+1+qa7GzRRRXiH1QUUUUAFFFFABRRRQBBeWkGoWNxZXUYlt7iNopUJIDIwwRx6gml+ywmz+yOpkgMflMsrF9y4xhi2S2R1z1qaigChpujWOk+Z9kjkDSBQ7yzPK5C/dG5yTgZOBnAyfWr9FFAHDfEX/mG/wDbX/2SuGrufiL/AMw3/tr/AOyVw1fV5b/usfn+bPzzPf8AkYVPl/6Sgrd8Hf8AI1WX/A//AEBqwq3fB3/I1WX/AAP/ANAat8V/An6P8jjy/wD3ul/ij+aNXxL/AMjBdf8AAP8A0AVk1reJf+Rguv8AgH/oArJr40/Twrq/Bf8Ay/f9s/8A2auUrq/Bf/L9/wBs/wD2agDq6KKKACiiigAooooAKKKKACiiigAooooA4b4i/wDMN/7a/wDslcNXuVFerhsz9hSVPkvbz/4B87jsg+t4iVf2lr20tfZJdzw2t3wd/wAjVZf8D/8AQGr1WitKub+0hKHJumt/+AZYfhv2NaNX2t+Vp/D2fqef+Jf+Rguv+Af+gCsmvVqK8Y+nPKa6vwX/AMv3/bP/ANmrq6KACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAg+2WpvjYi5h+1iMTGDePMCE4Dbeu3IIz0yKW0vLW/tlubO5huIGJCywuHUkEg4I44II+orhvElrqtt4zutZ0yzuJJhp1vZI6RFgfMa6Hb+7J9nZvReTW94J006R4aFh5MkKQXt4saSKQfL+0ylDz2KkEHuCDQB0NFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAEf/9k=",
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAKAAAACqCAIAAAClNtMcAAAHhklEQVR4Ae2dTUgVfRTGu2ZJShaGS924SKQSMvzA/IBMQ9NAohYVuSjaqCS6C3ITJQkudNVCMZEWUaKhC0NMBN0Y1CKtFkLkQoiEFN2IH+95ESa5Odc5di79z/G5C5k788y55zy/e2ZGuGcmtLm5eQAvuw7E2C0Nlf3vAAAb/x4AMAAbd8B4eehgADbugPHy0MEAbNwB4+WhgwHYuAPGy0MHA7BxB4yXhw4GYOMOGC8PHQzAxh1glDc8PFxdXb21w/T0dHZ29vr6OmN/f+nMzMzp06c3Nja2JOfOnfv06ZO/nLEFHcwwq6ysbHV1dXx8nPZpaGhob28/ePAgY39/aUZGxqlTp16+fEmSoaGhtLQ0eusv52yhX3TgFdyBL1++5Obmvnr16s6dO8H3CqKkJiaodEjIz8+nw0OQXYJoQiTifB+gPXD//v2+vr4PHz6cOHFC1o4bN27Ex8evrKy8ePFCKnKsVKD9E2dpaSk2NnZ5eVkc8MOHD+lYTa0saCbOwTwzJyYmFhcXnz17Vl9fz9szgPrkyZNHjx6lvwG0QSUAHNQp0q2trdG1VWtr68WLF6mJBwYGGDv/IynOwQzj29rafv78+fjxY9rn+/fvdFH9/v37hIQERojdpMePH//169duKsZ2AGaYpVGKQ7RGaoycAZhhlkYpAGukxsgZgBlmaZQCsEZqjJwBmGGWRikAa6TGyBmAGWZplAKwRmqMnAGYYZZGKQBrpMbIGYAZZmmUArBGaoycAZhhlkYpAGukxsgZgBlmaZQCsEZqjJwBmGGWRikAa6TGyBmAGWZplAKwRmqMnG0Cbmlp8fOAfpRaU1Nz7NgxP4Gx9fsOcGVlZVZWVigUMgbSt5wgE2q6NDThQ1OdNHywY9rz8/O0njp4x632Vtr84fvWfMDXr1/v3bvnfbXHxsa8ZfEBAi+yawuWAUfwem+Atwa0I4Rlbbp27Zqnj15ky+OjETrYc9b8gk3AdLMLetEc5vbDsnmWOxZo8yq6oKCgqqpqx4L320qbHUz3MYkMUnZEM/Jn/dutNjv433rq1KcDsFM45JMBYHlPnYoIwE7hkE8GgOU9dSoiADuFQz4ZAJb31KmIAOwUDvlkAFjeU6ciArBTOOSTAWB5T52KCMBO4ZBPBoDlPXUqIgA7hUM+GQCW99SpiADsFA75ZABY3lOnIgKwUzjkkwFgeU+digjATuGQTwaA5T11K6K9aRyq6MmTJ351dXZ2nj9/PjMzkx5E6KextP6ApWK8Wvxmy378+FFYWEiPj/v8+XN6erqnN7xg8BDd3NxMjyUrLS3981C5sLBQW1sbExOTkpJCy38K7K2xPHwWYTbp+fPn9BBROlzbIxpWkWXAYaV6b2dnZ69cufLu3bvk5GRvZZCF0FgQVVDNZvFvZfQi2xxd2XJuxw6mo/f169epd7l0f9NQtWQTsN90IV1M3b59u7GxMScnRxWmvSdr8CKLzPCbLuzu7qb/jujZocXFxZcvX967bXr2tHkOjpL/0TtTRi+yzQ6OEmCNYQFYIzVGzgDMMEujFIA1UmPkDMAMszRKAVgjNUbOAMwwS6MUgDVSY+QMwAyzNEoBWCM1Rs4AzDBLoxSANVJj5AzADLM0SgFYIzVGzgDMMEujFIA1UmPkDMAMszRKAVgjNUbOAMwwS6MUgDVSY+QMwAyzVEpNzl35TReurKzQ04qKiorOnj07ODhosvawomz+bNbvuVdPnz6lybOmpiZ6wF1eXt63b99UNiUnaYOTDd504du3b8OsuHv3bkJCAq2cnp4+dOhQ2FaTby138I6zSUTx1q1bfX19b968uXDhgkmo24uyDHh7nWHLAwMDr1+/7unpCVsf+W1oLPJ23tbt04XRe3ah5ato6mCaQfJeZH9dXd3a2hot0GDSrg/P4uFyVW3wHExW+00XLi4u9vf3X716dXJykp5s6CoUybxsAt6aLqR/hMKsevToEY2PdnR0HD58uKurK2yrybc2AfsdflNTU2mw3yRIv6Isn4P9at5X6wHYOG4ABmDjDhgvDx0MwMYdMF4eOhiAjTtgvDx0MAAbd8B4eehgADbugPHy0MEAbNwB4+WhgwHYuAPGy0MHA7BxB4yXhw4GYOMOGC8PHWwcsM3JBuPQOOWhgzluKdQCsEJonJQBmOOWQi0AK4TGSRmAOW4p1AKwQmiclAGY45ZCLQArhMZJGYA5binUArBCaJyUAZjjlkItACuExkkZgDluKdQCsEJonJQBmONWNLVHjhyhG6V6n3Dz5k1a473d8wIA79k64R3j4uLozm3r6+sUl+4YOzs7S2v+/jMA+O89FItANzmempqicB8/fjxz5oxIXAAWsVEmSFlZ2fDwMMWiv7QsEhSARWyUCVJaWjoyMkKxRkdHS0pKRIICsIiNMkGSkpLofuVzc3MULjExUSQoAIvYKBbk0qVLDx48kGpfSguAxdiIBKqoqOjt7S0vLxeJRkHws1kpJx2Ngw52FIxUWgAs5aSjcQDYUTBSaQGwlJOOxgFgR8FIpQXAUk46GgeAHQUjlRYASznpaBwAdhSMVFoALOWko3H+A6stvXkaOOnSAAAAAElFTkSuQmCC",
      "text/plain": [
       "<PIL.Image.Image image mode=RGB size=160x170>"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dinkum.reset()\n",
    "\n",
    "def custom_not(X):\n",
    "    is_active = not X\n",
    "    print(f\"custom cis-regulatory logic: not! inputs: gene X={X} => gene active? {is_active}\")\n",
    "    return is_active\n",
    "\n",
    "x = Gene(name='X')\n",
    "y = Gene(name='Y')\n",
    "m = Tissue(name='M')\n",
    "\n",
    "x.is_present(where=m, start=2)\n",
    "y.custom_activation(state_fn=custom_not, delay=1)\n",
    "\n",
    "dinkum.run_and_display(stop=3)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "6afbbfa3-699d-4d1e-9dce-54210b978769",
   "metadata": {},
   "source": [
    "## Building more complex input functions\n",
    "\n",
    "You can build input functions that take arbitrary numbers of genes (not just two) and apply more complex logic operations to them.\n",
    "\n",
    "Here's one that activates the downstream gene if one of A and B are active, and one of C and D are active: (A or B) and (C or D)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "490c804a-6c18-4a1c-8b2b-767d3d1f396a",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "initializing: dinkum v0.3.2\n",
      "start=1 stop=7\n",
      "custom cis-regulatory logic! inputs: gene A=True, gene B=False, gene C=False, gene D=False => gene active? False\n",
      "custom cis-regulatory logic! inputs: gene A=True, gene B=True, gene C=False, gene D=False => gene active? False\n",
      "custom cis-regulatory logic! inputs: gene A=False, gene B=True, gene C=True, gene D=False => gene active? True\n",
      "custom cis-regulatory logic! inputs: gene A=False, gene B=True, gene C=True, gene D=True => gene active? True\n",
      "custom cis-regulatory logic! inputs: gene A=False, gene B=False, gene C=True, gene D=True => gene active? False\n",
      "custom cis-regulatory logic! inputs: gene A=False, gene B=False, gene C=False, gene D=False => gene active? False\n"
     ]
    },
    {
     "data": {
      "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wBDAAgGBgcGBQgHBwcJCQgKDBQNDAsLDBkSEw8UHRofHh0aHBwgJC4nICIsIxwcKDcpLDAxNDQ0Hyc5PTgyPC4zNDL/2wBDAQkJCQwLDBgNDRgyIRwhMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjL/wAARCAEiAPoDASIAAhEBAxEB/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQAAAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWmp6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEAAwEBAQEBAQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAECAxEEBSExBhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRomJygpKjU2Nzg5OkNERUZHSElKU1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9oADAMBAAIRAxEAPwD3+iiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigDm/GUUmoWFtokMMU76hOFkillMatEnzvlgrEA7QvQ/fq74Yup7vw7aG6IN3CDbXJBzmWJjG5+hZSR7GteigAooooA8u8dvbr4l1hWs9NubptGtfs/n3IiuVbzbn/j2Gxi0nTGCOQvXNeo0UUAFFFFAHCeFjozPpBf8A5Gb5/wC0PKx5wl2N5vn9/L3fdzxny8cV3dFFABWD44/5EDxJ/wBgu5/9FNW9RQByWijSv+Erz4b+y/2d9if7Z9h2+R529PK+78u/b5ucc4257V1tFFABXPeMr3TrLQh/aEVhKZplhtkvwpi85s4LbuAFAZj3wpxzXQ0UAZPhq002w8PWdrpE8M9lGpCywspWRtxLt8vGSxYnHGSa1qKKACiiigAooooAKKKKACiiigAooooAKKhurqCxs57u5kWK3gjaWWRuiqoySfoBU1ABRVDSta07W47iTTblbiO3m8mR1Ugb9qtwSPmG1lORkHPWr9ABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAY/ivTZNY8IazpsMSSz3NlNFEj4wZChC9eB82Oe1QWWiXEWrQ36T/Z7RUAWwzMPKGzaFwJzFx7Jjjjnmt+igDI0exntNT1+WWPZHdX6TQHIO5BbQJnjp8yMOfT0xWvRRQBynjTWL/SfsP2GfyvN8zf8itnG3HUH1Ncr/wAJjr3/AD//APkFP/ia3PiL/wAw3/tr/wCyVw1fS5fQpTw0ZSim9ei7s+EznGYinjpxhUklpom10Ru/8Jjr3/P/AP8AkFP/AImtXw14l1fUPEFra3V35kL79y+WgzhCRyBnqK42t3wd/wAjVZf8D/8AQGrfE4aiqM2oLZ9F2OXA47FSxVOMqkmnJdX3Ol1zXNRs9YngguNkS7cLsU4yoPce9Z//AAkur/8AP3/5DT/CjxL/AMjBdf8AAP8A0AVk18mfoprf8JLq/wDz9/8AkNP8K6DwxqV3qP2r7VN5nl7NvygYzuz0HsK4mur8F/8AL9/2z/8AZqAMr/hJdX/5+/8AyGn+FH/CS6v/AM/f/kNP8KyaKANb/hJdX/5+/wDyGn+FdB/aV3/wh/2/zv8ASf7+0f8APTHTGOlcTXV/80//AM/89aAMr/hJdX/5+/8AyGn+FH/CS6v/AM/f/kNP8KyaKANb/hJdX/5+/wDyGn+FdB4n1K7077L9lm8vzN+75Qc4246j3NcTXV+NP+XH/tp/7LQBlf8ACS6v/wA/f/kNP8KP+El1f/n7/wDIaf4Vk0UAdNoeuajeaxBBPcb4m3ZXYozhSew9q7CvP/DX/IwWv/A//QDXoFABRRRQAUUUUAFFFFABRRRQAUVQ1zVY9D0G/wBVmjkljs7d52SMEswUE4GPp16DqeKSLXdKmuVtUv7c3jLuFoXCz/d3YMZ+YHHOCM0AaFFY2ja1c6jeXVpeWAs54IoZtom8z5ZN2Fb5RtcbDleQMjk5rZoA4b4i/wDMN/7a/wDslcNXq3iLw7/b/wBm/wBK8jyN3/LPdu3Y9x6Vhf8ACuv+or/5L/8A2VfQYHHUKVCMJys1fo+58Zm2U4zEYydWlC8Xbquy8zhq3fB3/I1WX/A//QGrc/4V1/1Ff/Jf/wCyq9o/gv8AsnVYb7+0PN8rd8nk7c5UjruPrWtfMMNOlKMZatPo+3oc+DybHU8RTnKGikm9V0fqZXiX/kYLr/gH/oArJrttS8Mf2jqEt19s8vzMfL5WcYAHXPtVT/hC/wDqIf8AkH/7Kvmj7s5Sur8F/wDL9/2z/wDZqP8AhC/+oh/5B/8Asq1tF0X+x/P/ANI87zdv8G3GM+59aAPPqK6v/hC/+oh/5B/+yo/4Qv8A6iH/AJB/+yoA5Sur/wCaf/5/560f8IX/ANRD/wAg/wD2Va39i/8AFP8A9lfaP+2uz/b3dM/h1oA8+orq/wDhC/8AqIf+Qf8A7Kj/AIQv/qIf+Qf/ALKgDlK6vxp/y4/9tP8A2Wj/AIQv/qIf+Qf/ALKtbWtF/tjyP9I8nyt38G7Oce49KAPPqK6v/hC/+oh/5B/+yo/4Qv8A6iH/AJB/+yoAyvDX/IwWv/A//QDXoFc/pvhj+ztQiuvtnmeXn5fKxnII65966CgAooooAKKKKACiiigAooooAo61po1nQdR0tpTEt7bSW5kAyU3qVzjvjNMTSgtwL1ric3uzki4mEBbbjPk+YVA74/XPNaNFAGN4d0i90e3nS9vre9mmk82SeO2aJ5HP3mbLtngKABgAAAcYxs0UUAcN8Rf+Yb/21/8AZK4au5+Iv/MN/wC2v/slcNX1eW/7rH5/mz88z3/kYVPl/wCkoK3fB3/I1WX/AAP/ANAasKt3wd/yNVl/wP8A9Aat8V/An6P8jjy//e6X+KP5oPGP/I1Xv/AP/QFrCrd8Y/8AI1Xv/AP/AEBawqML/Ah6L8gzD/e6v+KX5sK7n4df8xL/ALZf+z1w1dz8Ov8AmJf9sv8A2esMy/3WXy/NHZkX/Iwp/P8A9JZhUUUV8ofoYV1f/NP/APP/AD1rlK6v/mn/APn/AJ60AcpRRRQAV1fjT/lx/wC2n/stcpXV+NP+XH/tp/7LQBylFFFAGt4a/wCRgtf+B/8AoBr0CvP/AA1/yMFr/wAD/wDQDXoFABRRRQAUUUUAFFFFABRRRQAUVkeKtQutK8JaxqNkIzc2tlLPH5h+UFUJyeDnp079OOtVT4pgGptpMwht71YyXc3cDxxsI95O3zBIVHrsHHPA5oA6Giuc8PXWpf21q2nahczTi1SBo2uEiV33bwXURjHlkp8oPzZVs9q6OgDhviL/AMw3/tr/AOyVw1es69Dosv2f+2Fzjd5XL+2fu/h1rG+x+Cv+ef6zV7eEzKlRoqnJO69O/qfK5lkeIxWKlWhKNnbe/RJdjz+t3wd/yNVl/wAD/wDQGro/sfgr/nn+s1XdKtvC8epwtpyYuxu2HMvoc/e46ZrStmtGdOUEnqmun+ZhheHsVSrwqSlG0Wnu+j9Dj/GP/I1Xv/AP/QFrCr0fWP8AhEv7Vm/tP/j8+XzP9b/dGPu8dMVR/wCKD/z59bUMZy0ox9nJ2S6eRy4zLOfEVJe2pq8no5a79dNzhq7n4df8xL/tl/7PR/xQf+fPrd8O/wBg/wCk/wBif7Pnff8Afb978elZY7Fe0oSjySW2603OjKcu9jjIVPawla+ild7PpY4Wiuk+2eCv+en6TUfbPBX/AD0/SavG+q1/5H9zPqP7Qwn/AD9j/wCBL/M5uur/AOaf/wCf+etQfbPBX/PT9Jq20fRpPDoZTnSz3+f+/wD99feqJ0akFecWvVGlLFUKr5ac1J+TTOAorq/+KS/z5tH/ABSX+fNrM3OUrq/Gn/Lj/wBtP/ZaP+KS/wA+bWtrX9kfuP7V/wBry/v+2fu/hQB59RXV/wDFJf582j/ikv8APm0AZXhr/kYLX/gf/oBr0Cuf03/hHf7Qi+wf8fPOz/Wehz146ZroKACiiigAooooAKKKKACiiigCC8tINQsbiyuoxLb3EbRSoSQGRhgjj1BNL9lhNn9kdTJAY/KZZWL7lxjDFslsjrnrU1FAFDTdGsdJ8z7JHIGkCh3lmeVyF+6NzknAycDOBk+tX6KKAOU8af8ALj/20/8AZa5Sur8af8uP/bT/ANlrlKACtbw1/wAjBa/8D/8AQDWTWt4a/wCRgtf+B/8AoBoAyvGP/I1Xv/AP/QFrCrd8Y/8AI1Xv/AP/AEBawq+ywv8AAh6L8j8wzD/e6v8Ail+bCu5+HX/MS/7Zf+z1w1dz8Ov+Yl/2y/8AZ6wzL/dZfL80dmRf8jCn8/8A0lnDUUUV3HkBXoFn/wAk1T8f/R1ef16BZ/8AJNU/H/0dXlZx/AXr+jPouGf97l/hf5o5uiiivmz7gK6vxp/y4/8AbT/2WuUrq/Gn/Lj/ANtP/ZaAOUooooA1vDX/ACMFr/wP/wBANegV5/4a/wCRgtf+B/8AoBr0CgAooooAKKKKACiiigAooooAKKwvGwm/4QTxAYJZIpRp1wyNF97IjY8cd+nr6YNMtbzWTqcSohu9K8sMl0IYszDy8ht4nHU+kQHPYc0AdBRXC+EDqJ8a63Jqdjc295c2FnNMJJI2RD5lyFRdrn5QuADxkoxIBIz3VAHKeNP+XH/tp/7LXKV6DrWtf2P5H+j+d5u7+PbjGPY+tZP/AAmn/UP/API3/wBjQByla3hr/kYLX/gf/oBrV/4TT/qH/wDkb/7Grem+J/7R1CK1+x+X5mfm83OMAnpj2oA4rxj/AMjVe/8AAP8A0Bawq9H1jxp/ZOqzWP8AZ/m+Vt+fztucqD02n1qj/wALF/6hX/kx/wDY19LQr4lUoqNK6svtLsfCYzC4GWIqOWIs+Z3XI3bXY4au5+HX/MS/7Zf+z0f8LF/6hX/kx/8AY1u+HfEX9v8A2n/RfI8jb/y03bt2fYelZY6tXlQkp07LTW6fU6Mpw2DhjISpV+aWunK10fU8poruf+Fi/wDUK/8AJj/7Gj/hYv8A1Cv/ACY/+xrr+sYr/nz/AOTI8/6nl/8A0E/+SSOGr0Cz/wCSap+P/o6oP+Fi/wDUK/8AJj/7Gt3/AISL/ilf7b+y/wDbHzP9vZ97H49K4MfKvVhGEqdrtdU9dT1snp4PD1p1adbmtF391qyurv8A4BwtFbv/AAsX/qFf+TH/ANjR/wALF/6hX/kx/wDY1w/2biv5fxX+Z6/9u5f/AM/Pwl/kYVdX40/5cf8Atp/7LVH/AIWL/wBQr/yY/wDsa3fEXiL+wPs3+i+f5+7/AJabdu3HsfWolga8ZKDjq9tV0NYZtg505VYz92Nr6Prt0OFord/4WL/1Cv8AyY/+xo/4WL/1Cv8AyY/+xq/7NxX8v4r/ADMv7dy//n5+Ev8AIg8Nf8jBa/8AA/8A0A16BXKaP40/tbVYbH+z/K83d8/nbsYUnptHpXV1zVqFSjLlqKzO7DYujioOdF3V7dV+YUUUVkdAUUUUAFFFFABRRRQAUUUUAFFFFAHKeNP+XH/tp/7LXKV1fjT/AJcf+2n/ALLXKUAFa3hr/kYLX/gf/oBrJrW8Nf8AIwWv/A//AEA0AZXjH/kar3/gH/oC1hVu+Mf+Rqvf+Af+gLWFX2WF/gQ9F+R+YZh/vdX/ABS/NhXc/Dr/AJiX/bL/ANnrhq7n4df8xL/tl/7PWGZf7rL5fmjsyL/kYU/n/wCks4aiiiu48gK7n/mlf+f+e9cNXc/80r/z/wA964cd/wAu/wDGv1PXyn/l/wD9e5focNRRRXceQFdz8Rf+Yb/21/8AZK4au5+Iv/MN/wC2v/slcOI/3qj/ANvfkevg/wDkX4n/ALc/9KOGoooruPIN3wd/yNVl/wAD/wDQGr1WvKvB3/I1WX/A/wD0Bq9Vr5vOP469P1Z9xwz/ALpL/E/yQUUUV5R9EFFFFABRRRQAUUUUAQ3V1BY2c93cyLFbwRtLLI3RVUZJP0Ap8siQxPLI21EUsx9AOtZXivTZNY8IazpsMSSz3NlNFEj4wZChC9eB82Oe1Z03hy7maWaGbyLaS1aNNP3zJ5eYigTiZolwT2jIGOOeaANzTdVtdWhM1qLny8Agz2ssO4EZBG9Rn8Ku1zvhez1C0kuvPhvbeyMcKw297cieRZAG8xgwZgFIKADPVTwM10VAHP8AifTbvUfsv2WHzPL37vmAxnbjqfY1z/8AwjWr/wDPp/5ET/Gug8T6ld6d9l+yzeX5m/d8oOcbcdR7muf/AOEl1f8A5+//ACGn+FAB/wAI1q//AD6f+RE/xrQ0PQ9Rs9Ygnnt9kS7stvU4ypHY+9Z//CS6v/z9/wDkNP8ACtDQ9c1G81iCCe43xNuyuxRnCk9h7UAZ/iXw1q+oeILq6tbTzIX2bW8xBnCAHgnPUVlf8Idr3/Ph/wCRk/8Aiq6XXNc1Gz1ieCC42RLtwuxTjKg9x71n/wDCS6v/AM/f/kNP8K9OnmtaEFBJaadf8zwq3D2Fq1JVJSleTb3XX5GV/wAIdr3/AD4f+Rk/+KrqvBej3+k/bvt0Hleb5ez51bON2ehPqKyv+El1f/n7/wDIaf4V0HhjUrvUftX2qbzPL2bflAxndnoPYVFfMqtam6ckrP1/zNMJkeHwtZVoSldd7dVbscV/wh2vf8+H/kZP/iqP+EO17/nw/wDIyf8AxVH/AAmOvf8AP/8A+QU/+Jo/4THXv+f/AP8AIKf/ABNe3/t39z8T5X/hJ/6ef+Sh/wAIdr3/AD4f+Rk/+Krqv7Hv/wDhAP7M8j/TP+ee9f8Anru65x0965X/AITHXv8An/8A/IKf/E11X9sX/wDwgH9p+f8A6Z/z02L/AM9dvTGOntXJivrXuc/L8Sta+/n5HoZd/Z3732XP8Er35dtL2t17dDlf+EO17/nw/wDIyf8AxVH/AAh2vf8APh/5GT/4qj/hMde/5/8A/wAgp/8AE0f8Jjr3/P8A/wDkFP8A4muv/bv7n4nn/wDCT/08/wDJQ/4Q7Xv+fD/yMn/xVdV400e/1b7D9hg83yvM3/Oq4ztx1I9DXK/8Jjr3/P8A/wDkFP8A4muq8aaxf6T9h+wz+V5vmb/kVs4246g+prkrfWvb0+blvrbe22tz0MN/Z31Ovyc/L7t78t99LdN97nK/8Idr3/Ph/wCRk/8AiqP+EO17/nw/8jJ/8VR/wmOvf8//AP5BT/4mj/hMde/5/wD/AMgp/wDE11/7d/c/E8//AISf+nn/AJKavhrw1q+n+ILW6urTy4U37m8xDjKEDgHPU16BXn/hrxLq+oeILW1urvzIX37l8tBnCEjkDPUV6BXiZl7X2q9ra9ul+77n1WR/V/q8vq3Ny832rXvZdugUUUV557IUUUUAFFFFABRRRQAUUUUAFFFFAHKeNP8Alx/7af8AstcpXV+NP+XH/tp/7LXKUAFa3hr/AJGC1/4H/wCgGsmtbw1/yMFr/wAD/wDQDQAeJf8AkYLr/gH/AKAKya1vEv8AyMF1/wAA/wDQBWTQAV1fgv8A5fv+2f8A7NXKV1fgv/l+/wC2f/s1AHnFFFFfcn5KFdz/AM0r/wA/8964au5/5pX/AJ/571w47/l3/jX6nr5T/wAv/wDr3L9DhqKKK7jyArufiL/zDf8Atr/7JXDV3PxF/wCYb/21/wDZK4cR/vVH/t78j18H/wAi/E/9uf8ApRw1FFFdx5Bu+Dv+Rqsv+B/+gNXqteVeDv8AkarL/gf/AKA1eq183nH8den6s+44Z/3SX+J/kgoooryj6IKKKKACiiigAooooAoa5qseh6Df6rNHJLHZ27zskYJZgoJwMfTr0HU8VWufFOiWkE7y6jbrLDA1w9sXCzhFQucxnDA7RnBANWta00azoOo6W0piW9tpLcyAZKb1K5x3xmql54bttStrhL2e6eS6t3gn8u6mWI70KMViLlBwTgYOOvXmgBPD+uSa0LjzLeCMwlfmt7oXCHIOVLADDrjlccZHJzW1WPpGjTWF5dXt1dRT3NxFFCfIt/JRUj3bfl3Nz87ZOemBgYrYoAyta0X+2PI/0jyfK3fwbs5x7j0rJ/4Qv/qIf+Qf/sqPGn/Lj/20/wDZa5SgDq/+EL/6iH/kH/7Krem+GP7O1CK6+2eZ5efl8rGcgjrn3ria1vDX/IwWv/A//QDQB0GpeGP7R1CW6+2eX5mPl8rOMADrn2qp/wAIX/1EP/IP/wBlWV4l/wCRguv+Af8AoArJoA6v/hC/+oh/5B/+yrW0XRf7H8//AEjzvN2/wbcYz7n1rz6ur8F/8v3/AGz/APZqAKP/AArr/qK/+S//ANlR/wAK6/6iv/kv/wDZVhUV3f2liv5vwX+R5P8AYWX/APPv8Zf5m7/wrr/qK/8Akv8A/ZVu/wDCO/8AFK/2J9q/7beX/t7/ALufw61wtdX/AM0//wA/89aieOr1Lc0tnfZbmtHKcHR5vZwtzJp6vZ79Sj/wrr/qK/8Akv8A/ZUf8K6/6iv/AJL/AP2VYVFX/aWK/m/Bf5GX9hZf/wA+/wAZf5m7/wAK6/6iv/kv/wDZVu+IvDv9v/Zv9K8jyN3/ACz3bt2PcelcLXV+NP8Alx/7af8AstRLHV5SU3LVbaLqawynBwpypRh7srX1fTbqUf8AhXX/AFFf/Jf/AOyo/wCFdf8AUV/8l/8A7KsKir/tLFfzfgv8jL+wsv8A+ff4y/zOr0fwX/ZOqw339oeb5W75PJ25ypHXcfWurrz/AMNf8jBa/wDA/wD0A16BXNWr1K0uao7s7sNhKOFg4UVZXv1f5hRRRWR0BRRRQAUUUUAFFFFABRRRQAUUUUAcp40/5cf+2n/stcpXV+NP+XH/ALaf+y1ylABWt4a/5GC1/wCB/wDoBrJrW8Nf8jBa/wDA/wD0A0AHiX/kYLr/AIB/6AKya1vEv/IwXX/AP/QBWTQAV1fgv/l+/wC2f/s1cpXV+C/+X7/tn/7NQBylFFFABXV/80//AM/89a5Sur/5p/8A5/560AcpRRRQAV1fjT/lx/7af+y1yldX40/5cf8Atp/7LQBylFFFAGt4a/5GC1/4H/6Aa9Arz/w1/wAjBa/8D/8AQDXoFABRRRQAUUUUAFFFFABRRRQBkeKtQutK8JaxqNkIzc2tlLPH5h+UFUJyeDnp079OOtV18RxT6v8A2E48i/ZCrPDd27+W2zcdqlvMOO2Y/cjFbF5aQahY3FldRiW3uI2ilQkgMjDBHHqCakjiWKBIVLlUUKCzlmwBjlick+5OaAMPw097PLqU82pXF7ZifyLUzpED+7ysjZRF4L7l5zwgI61v1BZ2cGn2UFnaxiO3gQRxpknCgYHJ5P1NT0AZWtf2R+4/tX/a8v7/ALZ+7+FZP/FJf582jxp/y4/9tP8A2WuUoA6v/ikv8+bVvTf+Ed/tCL7B/wAfPOz/AFnoc9eOma4mtbw1/wAjBa/8D/8AQDQB0Gpf8I7/AGhL9v8A+Pnjf/rPQY6cdMVU/wCKS/z5tZXiX/kYLr/gH/oArJoA6v8A4pL/AD5ta2i/2R+//sr/AGfM+/74+9+NefV1fgv/AJfv+2f/ALNQAf8AFJf582j/AIpL/Pm1ylFAHV/8Ul/nza1v+JR/wj//AFDf+B/3/wDvr71efV1f/NP/APP/AD1oAP8Aikv8+bR/xSX+fNrlKKAOr/4pL/Pm1ra1/ZH7j+1f9ry/v+2fu/hXn1dX40/5cf8Atp/7LQAf8Ul/nzaP+KS/z5tcpRQB22m/8I7/AGhF9g/4+edn+s9Dnrx0zXQV5/4a/wCRgtf+B/8AoBr0CgAooooAKKKKACiiigAooooAKKKKACiiigDlPGn/AC4/9tP/AGWuUrq/Gn/Lj/20/wDZa5SgArW8Nf8AIwWv/A//AEA1k1reGv8AkYLX/gf/AKAaADxL/wAjBdf8A/8AQBWTWt4l/wCRguv+Af8AoArJoAK6vwX/AMv3/bP/ANmrlK6vwX/y/f8AbP8A9moA5SiiigArq/8Amn/+f+etcpXV/wDNP/8AP/PWgDlKKKKACur8af8ALj/20/8AZa5Sur8af8uP/bT/ANloA5SiiigDW8Nf8jBa/wDA/wD0A16BXn/hr/kYLX/gf/oBr0CgAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiqSaxpklpb3aajZtbXLrHBMJ1KSuxwFVs4Yk8ACrtABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFAHj+k6Hqseg+HtKk027WDTpdOv1DQsNrvLbbx06qftZYfwggmvYKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigKKKKAP/2Q==",
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPoAAAEiCAIAAABflFicAAAQZUlEQVR4Ae2dfYhV5RaHm+GO9v1hjNKHU2gx9EdaBqVUzqjhVFPZF1qRlsRghIoxQ00UBGYRE32ACI2aTWUQhZbjJI2ZTVKZ2B8SlWUJkaRJRVZ+pGnnru6B07kzZ+9z9t7vlr3e9cwfcs7e+31nref33N0+A3edqlwudww/ELBBoNpGm3QJgX8IoDseGCKA7obCplV0xwFDBNDdUNi0iu44YIgAuhsKm1bRHQcMEUB3Q2HTKrrjgCEC6G4obFpFdxwwRADdDYVNq+iOA4YIoLuhsGkV3SM7sHTp0sGDB+/evTvyytAFxx9/fGNjY0NDw8UXX7x69erQa6OdXLJkyZgxY2Tn5ubmHTt2RFscenW+5gkTJlx++eVdXV2h10Y72draKjTk58orrxw0aFC0xSFXy/+9g59IBG644Ya2trZly5ZFWlX24lNOOSV/zZYtW+rq6speX+EFa9euFR33798v169Zs2bixIkVLqzkskLNe/fulZ1ff/31SlZFuua5555rb2+PtCTk4mNCznFqIIF9+/ZNmjTpq6++uvnmmweeTXKkoM7ff/89YsSIJFsVr508efLGjRsLR1paWg4dOlR4m/BFoWbZ59NPP5V7fMIN+y3/4YcfRo0aJcz7HY/9Ft2joVu5cuVTTz0la+Tx4ODBg9EWh15dUGfdunXd3d2h10Y4edZZZ/35558RFkS5tFCzLJL/FQ0bNizK6vLXTp06taenp/x1FV/Bs3vIg16JU6tWrVq+fPnYsWN37tz5wQcflLgi7iHRRR5Vx40b19TUtHDhwrjb9F935MiR/ofSeX/48OGamhqHe/f29krx8nnD4Z7c3Su+M+RykqjomF/wzjvvzJ07N8LicpcW7pSfffZZ4XW5ReXPjx8//pNPPslfJ49JM2bMKL+m4iuK6+zr67vmmmsqXlrmwgMHDlx00UXywbrMdRFPc3ePcO/46KOPRo8enV8gfzGQT4ERFld86emnnz5y5MiKLy9z4X333ffII4/Ic5dc99prr+VflFkT/fSvv/764IMPPvDAA9GXll7x+OOPT58+/eyzzy59Ou7R/8RdaHGdPMnI3x/yncvf4IYOHbp169YLLrjACYv8w0x19T83oM7OTid7yibTpk375ptvLrnkktraWil40aJFrnaWffI1V1VV/fXXX6K7PIy52ryjo+PSSy+VzzD5DeUJ/sQTT0y+eZX81yD5LuwAARUEeJhRERNFuiGA7m44sosKAuiuIiaKdEMA3d1wZBcVBNBdRUwU6YYAurvhyC4qCKC7ipgo0g0BdHfDkV1UEEB3FTFRpBsC6O6GI7uoIIDuKmKiSDcE0N0NR3ZRQQDdVcREkW4IoLsbjuyiggC6q4iJIt0QQHc3HNlFBQF0VxETRbohgO5uOLKLCgLoriIminRDAN3dcGQXFQTQXUVMFOmGgJ+6P/nkk0F49uzZc/fdd8s8oKALOO4xAXO6X3/99TJ0RUajeBwqrQUR8FD3Rx99VOYvy+Tbkj2/8cYbc+bMKXmKg94T8HOs0qmnnioPLV9//fWsWbMKEcoQw8Lr/AWFt7wwQsBn3UMijKd7VV/IlpFP5Rr/XSJfBPDvm8SvZE50YQ92Lqbh84zIkLt7wQZemCLgp+4y2Vl+6uvrix9gTOVKsyUJePhRVfqUadTyDUolG+agZQJ+3t3ffvvt8FDlg2z4BZz1koCfd3cvo6Kp5ATQPTlDdlBDAN3VREWhyQmge3KG7KCGALqriYpCkxNA9+QM2UENAXRXExWFJieA7skZsoMaAuiuJioKTU4A3ZMzZAc1BNBdTVQUmpwAuidnyA5qCKC7mqgoNDkBdE/OkB3UEEB3NVFRaHIC6J6cITuoIYDuaqKi0OQE0D05Q3ZQQwDd1URFockJoHtyhuyghgC6q4mKQpMT8HOKmEwAbm9vL0ln2bJlL7744h9//NHR0RE0R7LkQg56QMBP3YNm4v3000+33nrr+++/v23btptuumnr1q0eREgLlRPw8GEmZALwL7/8Mnv27Orq6uHDh8vryjFxpR8EfL67h8yIfOmllzZs2PDCCy/4kSJdVEjAZ92DEGzfvn3KlCnySFNbWxt0TcnjVX0lD8c8WDwBOL2dmQBsegKwfNPBtGnT5L4e1fWYUrMsSwT8nBEZNAE4l8vdddddra2tl112WZZSoJajRMDDj6pCLmgCcFdXV29vb2dnZ2Nj43XXXXeUGPNrMkPAz7t70ATgmf/7yQx8CjnaBPy8ux9tivw+JQTQXUlQlOmCALq7oMgeSgigu5KgKNMFAXR3QZE9lBBAdyVBUaYLAujugiJ7KCGA7kqCokwXBNDdBUX2UEIA3ZUERZkuCKC7C4rsoYQAuisJijJdEEB3FxTZQwkBdFcSFGW6IIDuLiiyhxIC6K4kKMp0QQDdXVBkDyUE0F1JUJTpggC6u6DIHkoIoLuSoCjTBQF0d0GRPZQQ8HOKWNAE4P3798ssgt27d8sE4Pnz5zc3NyuJiTLdEPBT96AJwDLkWuahtrW17dq1a9y4cd99950biuyihICHc2YKE4DXrl3bL4WWlpYTTjhBDn7xxRc1NTX9zvLWewI+392DJgBPnz595cqV3d3dkyZN8j5gGiwm4LPuxX32e71q1aoVK1a8/PLL/Y6Hv01vmm5VX/hvjnaW2cLFvIonAPv8lxm5u8ssyMKPIJgzZ87hw4flhQyIDBqsV0yK154R8PDZXRIKmgD822+/vfXWW/J9NR9//HF9fb1nWdJOWQJ+6p6fANzT09Ov/wULFsjA64ULFw4aNEi+k6zfWd56T8BP3YMeVOrq6uRLO7wPlQaDCPj87B7UM8fNEkB3s9FbbBzdLaZutmd0Nxu9xcbR3WLqZntGd7PRW2wc3S2mbrZndDcbvcXG0d1i6mZ7Rnez0VtsHN0tpm62Z3Q3G73FxtHdYupme0Z3s9FbbBzdLaZutmd0Nxu9xcbR3WLqZntGd7PRW2wc3S2mbrZndDcbvcXG0d1i6mZ79nOKWNAE4ELMa9asueWWWw4cOFA4wgsLBPzUPWgCcD5RmXbd1NT0+eef//777xYypscCAQ8fZgoTgAtN9nvR3t4+b948mXzd7zhvvSfg89295ATgDz/88Omnn37zzTfD/wvgffA2G/RZ94GJHjx4sKGhQaZdn3nmmTF0ZwLw0ZktXNU3MLr4R4pr9nNoXp7NwLu7DLmWB/c77rhDLti7d++MGTOizryOT52VGSDgp+5BE4BF9LzrQl7u7rieAQOPagl+flzLTwA+qiD5ZRoI+Hl3D5oAXJzInj17it/y2gIBP+/uFpKjxxgE0D0GNJZoJYDuWpOj7hgE0D0GNJZoJYDuWpOj7hgE0D0GNJZoJYDuWpOj7hgE0D0GNJZoJYDuWpOj7hgE0D0GNJZoJYDuWpOj7hgE0D0GNJZoJYDuWpOj7hgE0D0GNJZoJYDuWpOj7hgE0D0GNJZoJYDuWpOj7hgE0D0GNJZoJYDuWpOj7hgE0D0GNJZoJYDuWpOj7hgE/JwiFjQB+N133501a1ZdXZ2QuuKKKxYsWBADGUv0EvBT96CBeK+88sq+ffvuvfdevYFReRICHj7MhEwA3rVr1xlnnJGEF2tVE/D57j5wRuT9998vd/cvv/xyyJAhzz777MiRI1WHR/FRCfg5RSxPob6+vq+vr5hIVVXV6NGjFy9eLEOAW1pa1q9fX3y27Ov0JgCnt3PV/wEo22KZC4qn6Wrc2WfdB97d586dO3z4cIl0ypQp8pm1TLac9o6An7oHTQCW7+24/fbbxfVNmzZdeOGF3qVJQ2UI+Kl7fgJwT09Pv+4fe+yxmTNnylP7scceu2TJkn5nees9AT91D5oAfP7558uX1XgfKg0GEfDwD5FBrXIcAuiOA4YIoLuhsGkV3XHAEAF0NxQ2raI7DhgigO6GwqZVdMcBQwTQ3VDYtIruOGCIALobCptW0R0HDBFAd0Nh0yq644AhAuhuKGxaRXccMEQA3Q2FTavojgOGCKC7obBpFd1xwBABdDcUNq36OTQvaAKwzJ+ZN2/e5s2ba2pqurq6RowYgQGmCPh5dxfdS6b4/PPPn3TSSRs3bpRhkW1tbSWv4aDHBDzUPWQC8KuvvipjlSTO5ubmsWPHepwrrZUk4OfDTH6++8AZkbW1tQ899FB3d/dpp50ms8TOPffcklA46CsBP6eI5dMaOAH40KFD55xzjowFXrFixT333PPee+9FyjW9Ob3sPHXq1EIW6dHwWfeBd/dhw4bdeOONglX+5Ts8CnrZeeGn7kETgCdOnLhhw4YJEybIvzLo3U7MdJon4OFHVWksPwF4YMbz589/5plnRPcnnniis7Nz4AUc8ZuAn3f3oAnAQ4cOXb16td+J0l0IAT/v7iENc8oyAXS3nL653tHdXOSWG0Z3y+mb6x3dzUVuuWF0t5y+ud7R3VzklhtGd8vpm+sd3c1FbrlhdLecvrne0d1c5JYbRnfL6ZvrHd3NRW65YXS3nL653tHdXOSWG0Z3y+mb6x3dzUVuuWF0t5y+ud7R3VzklhtGd8vpm+sd3c1FbrlhP4fmBU0Avu2223788UfJW8aJbdu27eeff7acvcHe/dQ9PyMyJM6lS5d+//33MnYm5BpO+UfAw4eZkAnA+fxyudyiRYtmz57tX5x0FE7A57v7wBmReRYyAbinp2fx4sXhaDjrHwGfdQ9Kq6GhQVyX+cBBFwQdT28yLTszATjIukqPl7y7b9q0SZ7sY7he6W/lugwT8HNGZNAEYAmio6OjtbU1w4lQWooEPPyoKrSCJgB/++23O3fuHD9+fIpE2TrDBPy8uwdNAD7vvPPke8gyHAelpUvAz7t7uszYXS0BdFcbHYVHJ4Du0ZmxQi0BdFcbHYVHJ4Du0ZmxQi0BdFcbHYVHJ4Du0ZmxQi0BdFcbHYVHJ4Du0ZmxQi0BdFcbHYVHJ4Du0ZmxQi0BdFcbHYVHJ4Du0ZmxQi0BdFcbHYVHJ4Du0ZmxQi0BdFcbHYVHJ4Du0ZmxQi0BdFcbHYVHJ4Du0ZmxQi0BdFcbHYVHJ4Du0ZmxQi0BP6eIqY2DwtMlwN09Xb7snikC6J6pOCgmXQLoni5fds8UAXTPVBwUky4BdE+XL7tnigC6ZyoOikmXALqny5fdM0UA3TMVB8WkSwDd0+XL7pkigO6ZioNi0iWA7unyZfdMEUD3TMVBMekSQPd0+bJ7pgige6bioJh0CaB7unzV7X7ccccVf6PvnXfeKUfUdRFUMLoHkTF6fPDgwfL1y0eOHJH+c7nc9u3b5Yg3LNDdmyidNTJmzJjNmzfLdlu2bBk1apSzfTOwEbpnIISMldDU1NTb2ytFyb/yOmPVJSoH3RPh83Lx5MmT161bJ62tX7/+qquu8qlHdPcpTTe9DBkypLq6eseOHbLdySef7GbTbOyC7tnIIWNVXH311Q8//LBnt3ZhjO4ZEy0b5TQ3Ny9fvvzaa6/NRjnOqmDwhjOUbJR9Atzds58RFTojgO7OULJR9gmge/YzokJnBNDdGUo2yj4BdM9+RlTojAC6O0PJRtkngO7Zz4gKnRFAd2co2Sj7BNA9+xlRoTMC6O4MJRtln8B/AVo5GsUAjQQjAAAAAElFTkSuQmCC",
      "text/plain": [
       "<PIL.Image.Image image mode=RGB size=250x290>"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dinkum.reset()\n",
    "\n",
    "def custom_complex(A, B, C, D):\n",
    "    is_active = (A or B) and (C or D)\n",
    "    print(f\"custom cis-regulatory logic! inputs: gene A={A}, gene B={B}, gene C={C}, gene D={D} => gene active? {is_active}\")\n",
    "    return is_active\n",
    "\n",
    "a = Gene(name='A')\n",
    "b = Gene(name='B')\n",
    "c = Gene(name='C')\n",
    "d = Gene(name='D')\n",
    "z = Gene(name='Z')\n",
    "m = Tissue(name='M')\n",
    "\n",
    "a.is_present(where=m, start=1, duration=2)\n",
    "b.is_present(where=m, start=2, duration=3)\n",
    "c.is_present(where=m, start=3, duration=3)\n",
    "d.is_present(where=m, start=4, duration=2)\n",
    "z.custom_activation(state_fn=custom_complex, delay=1)\n",
    "\n",
    "dinkum.run_and_display(stop=7)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "4647d206-4982-4d75-8a98-75783dd0607e",
   "metadata": {},
   "source": [
    "## Building input functions with stochasticity\n",
    "\n",
    "You can use the `random` module to introduce stochastic gene activation (either unbiased - 50/50 - or biased).\n",
    "\n",
    "Here's a function that activates the gene based on an even coin flip (`random.choice([True, False])`)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "2c18bdc3-9339-46df-9e30-1ae838178bdb",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "initializing: dinkum v0.3.2\n",
      "start=1 stop=5\n",
      "custom cis-regulatory logic: random! gene active? False\n",
      "custom cis-regulatory logic: random! gene active? True\n",
      "custom cis-regulatory logic: random! gene active? True\n",
      "custom cis-regulatory logic: random! gene active? False\n",
      "custom cis-regulatory logic: random! gene active? False\n",
      "custom cis-regulatory logic: random! gene active? True\n",
      "custom cis-regulatory logic: random! gene active? True\n",
      "custom cis-regulatory logic: random! gene active? True\n",
      "custom cis-regulatory logic: random! gene active? True\n",
      "custom cis-regulatory logic: random! gene active? False\n",
      "custom cis-regulatory logic: random! gene active? False\n",
      "custom cis-regulatory logic: random! gene active? True\n"
     ]
    },
    {
     "data": {
      "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wBDAAgGBgcGBQgHBwcJCQgKDBQNDAsLDBkSEw8UHRofHh0aHBwgJC4nICIsIxwcKDcpLDAxNDQ0Hyc5PTgyPC4zNDL/2wBDAQkJCQwLDBgNDRgyIRwhMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjL/wAARCADmAL4DASIAAhEBAxEB/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQAAAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWmp6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEAAwEBAQEBAQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAECAxEEBSExBhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRomJygpKjU2Nzg5OkNERUZHSElKU1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9oADAMBAAIRAxEAPwD3+iiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKAOa8ZeebKzXTNv9ufaAdOz0DgHeW/2Nm4H6jHOKv+GjZnw7Z/YRIIQhBEv+sEmT5m/wD29+7d/tZrWooAKKKKAOQ8Vaa82taYkMoji1h/7Nv1x/rIlSSYfjiORPpKfQV1wAUAAAAcADtS0UAFFFFAHE6e3hzUvH011CdKju7GSSCJEMa3E1wQRK5A+Y7RlB7+YfQ121FFABRRRQAUUUUAFFFFABRRRQAUUUUAFFQ3V1BY2c93cyLFbwRtLLI3RVUZJP0AqagAoqhpWtadrcdxJptytxHbzeTI6qQN+1W4JHzDaynIyDnrV+gAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAx/Femyax4Q1nTYYklnubKaKJHxgyFCF68D5sc9qgstEuItWhv0n+z2ioAthmYeUNm0LgTmLj2THHHPNb9FAGRo9jPaanr8sseyO6v0mgOQdyC2gTPHT5kYc+npiteiigDn/ABPqV3p32X7LN5fmb93yg5xtx1Hua5//AISXV/8An7/8hp/hWr40/wCXH/tp/wCy1ylAGt/wkur/APP3/wCQ0/wrQ0PXNRvNYggnuN8TbsrsUZwpPYe1czWt4a/5GC1/4H/6AaANDXNc1Gz1ieCC42RLtwuxTjKg9x71n/8ACS6v/wA/f/kNP8KPEv8AyMF1/wAA/wDQBWTQBrf8JLq//P3/AOQ0/wAK6DwxqV3qP2r7VN5nl7NvygYzuz0HsK4mur8F/wDL9/2z/wDZqAMr/hJdX/5+/wDyGn+FH/CS6v8A8/f/AJDT/CsmigDW/wCEl1f/AJ+//Iaf4V2Gh3M15o8E8775W3ZbAGcMR2+lec16D4a/5F+1/wCB/wDoZoA1aKKKACiiigAooooAKKKKACiqGuarHoeg3+qzRySx2du87JGCWYKCcDH069B1PFJFrulTXK2qX9ubxl3C0LhZ/u7sGM/MDjnBGaANCisbRtaudRvLq0vLAWc8EUM20TeZ8sm7Ct8o2uNhyvIGRyc1s0Acp40/5cf+2n/stcpXoOtaL/bHkf6R5Plbv4N2c49x6Vk/8IX/ANRD/wAg/wD2VAHKVreGv+Rgtf8Agf8A6Aa1f+EL/wCoh/5B/wDsqt6b4Y/s7UIrr7Z5nl5+XysZyCOufegDn/Ev/IwXX/AP/QBWTXbal4Y/tHUJbr7Z5fmY+Xys4wAOufaqn/CF/wDUQ/8AIP8A9lQByldX4L/5fv8Atn/7NR/whf8A1EP/ACD/APZVraLov9j+f/pHnebt/g24xn3PrQB59RXV/wDCF/8AUQ/8g/8A2VH/AAhf/UQ/8g//AGVAHKV6D4a/5F+1/wCB/wDoZrJ/4Qv/AKiH/kH/AOyroNNsv7O0+K18zzPLz823Gcknp+NAFuiiigAooooAKKKKACiiigCjrWmjWdB1HS2lMS3ttJbmQDJTepXOO+M0xNKC3AvWuJze7OSLiYQFtuM+T5hUDvj9c81o0UAY3h3SL3R7edL2+t72aaTzZJ47Zonkc/eZsu2eAoAGAAABxjGzRRQBynjT/lx/7af+y1yldX40/wCXH/tp/wCy1ylABWt4a/5GC1/4H/6Aaya1vDX/ACMFr/wP/wBANAGV4x/5Gq9/4B/6AtYVbvjH/kar3/gH/oC1hV9lhf4EPRfkfmGYf73V/wAUvzYV3Pw6/wCYl/2y/wDZ64au5+HX/MS/7Zf+z1hmX+6y+X5o7Mi/5GFP5/8ApLOGoooruPICvVvB3/Iq2X/A/wD0Nq8pr1bwd/yKtl/wP/0Nq8rOP4C9f0Z9Fwz/AL3L/C/zRuUUUV82fcBRRRQAUUUUAFFFFABRWR4q1C60rwlrGo2QjNza2Us8fmH5QVQnJ4OenTv0461VPimAam2kzCG3vVjJdzdwPHGwj3k7fMEhUeuwcc8DmgDoaK5zw9dal/bWradqFzNOLVIGja4SJXfdvBdRGMeWSnyg/NlWz2ro6AOU8af8uP8A20/9lrlK9B1r+yP3H9q/7Xl/f9s/d/Csn/ikv8+bQByla3hr/kYLX/gf/oBrV/4pL/Pm1b03/hHf7Qi+wf8AHzzs/wBZ6HPXjpmgDivGP/I1Xv8AwD/0Bawq9N1W28LyanM2opm7O3ecy+gx93jpiqX2PwV/zz/Wavfo5rRhTjBp6JLp/mfHYrh7FVa86kZRtJt7vq/Q8/rufh1/zEv+2X/s9T/Y/BX/ADz/AFmrZ0GHRYvtH9jrjO3zeX98fe/HpWeLzKlWoypxTu/Tv6m+W5HiMLio1pyjZX2v1TXY8moruf8Aig/8+fR/xQf+fPr0Pr3/AE7n93/BPG/sn/p/T/8AAv8AgHDV6t4O/wCRVsv+B/8AobVhf8UH/nz66vR/sH9lQ/2Z/wAefzeX97+8c/e565rz8yxPtaSjySWvVW6M9nI8D9XxEpe0hL3dou73X4F6iiivEPqgooooAKKKKACiiigCC8tINQsbiyuoxLb3EbRSoSQGRhgjj1BNL9lhNn9kdTJAY/KZZWL7lxjDFslsjrnrU1FAFDTdGsdJ8z7JHIGkCh3lmeVyF+6NzknAycDOBk+tX6KKAOU8af8ALj/20/8AZa5Sur8af8uP/bT/ANlrlKACtbw1/wAjBa/8D/8AQDWTWt4a/wCRgtf+B/8AoBoAPEv/ACMF1/wD/wBAFZNa3iX/AJGC6/4B/wCgCsmgArq/Bf8Ay/f9s/8A2auUrq/Bf/L9/wBs/wD2agDziiiivuT8lCvVvB3/ACKtl/wP/wBDavKa9W8Hf8irZf8AA/8A0Nq8rOP4C9f0Z9Fwz/vcv8L/ADRuUUUV82fcBRRRQAUUUUAFFFFABRWF42E3/CCeIDBLJFKNOuGRovvZEbHjjv09fTBplreaydTiVEN3pXlhkuhDFmYeXkNvE46n0iA57DmgDoKK4XwgdRPjXW5NTsbm3vLmws5phJJGyIfMuQqLtc/KFwAeMlGJAJGe6oA5Txp/y4/9tP8A2WuUruvEXiL+wPs3+i+f5+7/AJabdu3HsfWsL/hYv/UK/wDJj/7Guulga9WKnCN0/NHnV82weHqOlVnaS8n/AJGFWt4a/wCRgtf+B/8AoBqf/hYv/UK/8mP/ALGr2j+NP7W1WGx/s/yvN3fP527GFJ6bR6U5ZfiYRcpR0Xmv8yKec4GpNQjPVuy0fX5GV4l/5GC6/wCAf+gCsmur1jxp/ZOqzWP9n+b5W35/O25yoPTafWqP/Cxf+oV/5Mf/AGNEcvxM4qUY6PzX+YVM5wNObhKeqdno+nyMKur8F/8AL9/2z/8AZqo/8LF/6hX/AJMf/Y1u+HfEX9v/AGn/AEXyPI2/8tN27dn2HpSq4GvSi5zjZLzRdDNsHiKipUp3k/J/5HlNFdz/AMLF/wCoV/5Mf/Y0f8LF/wCoV/5Mf/Y19B9YxX/Pn/yZHxv1PL/+gn/ySRw1ereDv+RVsv8Agf8A6G1YX/Cxf+oV/wCTH/2NdXo+o/2tpUN95Xlebu+TduxhiOuB6V5+ZVa06SVSnyq/dPoz2cjoYWniJOjW53y7crXVa6l6iiivEPqgooooAKKKKACiiigAooooAKKKKAOG+Iv/ADDf+2v/ALJXDV3PxF/5hv8A21/9krhq+ry3/dY/P82fnme/8jCp8v8A0lBW74O/5Gqy/wCB/wDoDVhVu+Dv+Rqsv+B/+gNW+K/gT9H+Rx5f/vdL/FH80HjH/kar3/gH/oC1hVu+Mf8Akar3/gH/AKAtYVGF/gQ9F+QZh/vdX/FL82Fdz8Ov+Yl/2y/9nrhq7n4df8xL/tl/7PWGZf7rL5fmjsyL/kYU/n/6SzhqKKK7jyAr1bwd/wAirZf8D/8AQ2rymvVvB3/Iq2X/AAP/ANDavKzj+AvX9GfRcM/73L/C/wA0blFFFfNn3AUUUUAFFFFABRRRQBDdXUFjZz3dzIsVvBG0ssjdFVRkk/QCnyyJDE8sjbURSzH0A61leK9Nk1jwhrOmwxJLPc2U0USPjBkKEL14HzY57VnTeHLuZpZoZvItpLVo00/fMnl5iKBOJmiXBPaMgY455oA3NN1W11aEzWoufLwCDPayw7gRkEb1Gfwq7XO+F7PULSS68+G9t7IxwrDb3tyJ5FkAbzGDBmAUgoAM9VPAzXRUAcp400e/1b7D9hg83yvM3/Oq4ztx1I9DXK/8Idr3/Ph/5GT/AOKrtfE+pXenfZfss3l+Zv3fKDnG3HUe5rn/APhJdX/5+/8AyGn+FehQzKrRpqnFKy9f8zxsXkeHxVZ1pyld9rdFbsZX/CHa9/z4f+Rk/wDiq1fDXhrV9P8AEFrdXVp5cKb9zeYhxlCBwDnqaP8AhJdX/wCfv/yGn+FaGh65qN5rEEE9xvibdldijOFJ7D2q6ma1pwcGlrp1/wAzOjw9haVSNSMpXi0910+Rn+JfDWr6h4gurq1tPMhfZtbzEGcIAeCc9RWV/wAIdr3/AD4f+Rk/+Krpdc1zUbPWJ4ILjZEu3C7FOMqD3HvWf/wkur/8/f8A5DT/AAop5rWhBQSWmnX/ADCtw9hatSVSUpXk2911+Rlf8Idr3/Ph/wCRk/8Aiq6rwXo9/pP277dB5Xm+Xs+dWzjdnoT6isr/AISXV/8An7/8hp/hXQeGNSu9R+1fapvM8vZt+UDGd2eg9hUV8yq1qbpySs/X/M0wmR4fC1lWhKV13t1VuxxX/CHa9/z4f+Rk/wDiqP8AhDte/wCfD/yMn/xVH/CY69/z/wD/AJBT/wCJo/4THXv+f/8A8gp/8TXt/wC3f3PxPlf+En/p5/5KH/CHa9/z4f8AkZP/AIqvQPDVncaf4ftbW6j8uZN+5cg4y5I5HHQ15/8A8Jjr3/P/AP8AkFP/AImvQPDV5cah4ftbq6k8yZ9+5sAZw5A4HHQV5+ZfWfZL2vLa/S/Z9z2cj+o/WJfVufm5ftWta67dTVooorxD6oKKKKACiiigAooooAKKKKACiiigDlPGn/Lj/wBtP/Za5Sur8af8uP8A20/9lrlKACtbw1/yMFr/AMD/APQDWTWt4a/5GC1/4H/6AaADxL/yMF1/wD/0AVk1reJf+Rguv+Af+gCsmgArq/Bf/L9/2z/9mrlK6vwX/wAv3/bP/wBmoA84ooor7k/JQr1bwd/yKtl/wP8A9DavKa9W8Hf8irZf8D/9DavKzj+AvX9GfRcM/wC9y/wv80blFFFfNn3AUUUUAFFFFABRRRQAUUUUAFFFFAGVrWi/2x5H+keT5W7+DdnOPcelZP8Awhf/AFEP/IP/ANlXV0UAcp/whf8A1EP/ACD/APZVb03wx/Z2oRXX2zzPLz8vlYzkEdc+9dBRQBz+peGP7R1CW6+2eX5mPl8rOMADrn2qp/whf/UQ/wDIP/2VdXRQByn/AAhf/UQ/8g//AGVa2i6L/Y/n/wCked5u3+DbjGfc+tatFAHDf8K6/wCor/5L/wD2VH/Cuv8AqK/+S/8A9lXc0V3f2liv5vwX+R5P9hZf/wA+/wAZf5nDf8K6/wCor/5L/wD2VdXo+nf2TpUNj5vm+Vu+fbtzliemT61eorKti61aPLUldfI6MNluFws3OjGztbdv82FFFFcx3BRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAcbL4uvodZ8T2DQW/+giNNN4bM0jRxEq/PPzzxDjHDVv+HNRm1fwxpOp3Cos95ZQ3EixghQzoGIAJJxk+pqrN4T0+fV/7SeS4843gvCoYbWYRxoFIx939zE+M53KDnHFaWladDpGj2WmW7O0Nnbx28bSEFiqKFBOABnA9KALdFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFEUUAf/2Q==",
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAL4AAADmCAIAAADky6hQAAALyElEQVR4Ae2dW0hVSxjHz7bSyjApeiiylwghUqOiCxkFlUZWUoR20xISe8hQ7EEICropRRfoqaLoRoGVZeiDFrULKaKgHtQuIEQSEhWUaRFonq822DrumTUzq9Gz1sx/PxyWM998zvy+n7PXjrNmh3p7e//BCwTUCcSoD8EIEPhFAOrAA48EoI5HcBgGdeCARwJQxyM4DIM6cMAjAajjERyGQR044JEA1PEIDsOgDhzwSADqeASHYVAHDngkAHU8gsMwqAMHPBKAOh7BYRjUUXCgvr5+zZo1kQHNzc2zZ8/u6elRGM8PbWlpSUlJ+fnzZyRk1qxZTU1N/HC1nrKyskW/XwsWLIiNjVUb7BJN/6sXXvIEsrKy7t+/T/FLly599OiR/EBh5Lp1665cuUJhtbW1OTk5wngPAcePHy8vL/cwkDnkH2YrGnkEXr58OXfu3GvXrm3dupUX462dNp5p06bRNjZ//nza0rwlcRn17t271NTUrq4ulxilrhBFu+xJ6IomUFJSUl1d/ezZs7Fjx0b3/k3Lxo0bR44cSdW9fPny3+Rhjs3Nzc3Pz6ddk9nroXGohzGWD+no6Bg6dGhnZ6d2dXbv3j116lTafrQTprs02s80evNrhkp7FIIbGxvpTrmhoWHVqlUDQWP06NHa037//n369OltbW16M+MTlsJfeHd3d2lp6eHDh+kemTaempoahcH/X+iBAwfy8vImTpyodwq411HgefTo0Y8fPx48eJDGvH37NjMz8+nTp/Hx8QopRKGJiYmfP38WRan1x8XF0b8jDBkyJDKMPsGNGjVKLQUrGuqwqKBNggDesCQgIYRFAOqwqKBNggDUkYCEEBYBqMOigjYJAlBHAhJCWASgDosK2iQIQB0JSAhhEYA6LCpokyAAdSQgIYRFAOqwqKBNggDUkYCEEBYBqMOigjYJAlBHAhJCWASgDosK2iQIQB0JSAhhEYA6LCpokyAAdSQgIYRFAOqwqKBNggDUkYCEEBYBqMOigjYJAlBHAhJCWATMVKeyspK12F9t9KjKli1b6Ek5XgDaJQlYp87KlStnzpwZCoUkASGMR8BAdfbs2UMPhGdkZDDXfPXq1eLiYmYXGpUImPkIX+QZylevXhUVFfXhCIfDfdcD8ZBlX3JLLkxWx6WE3tSpqqpyyanaRecv9Q0JYmaTD0lx2XX6aoYLzwTMVIcO5aNXcnKy803KMyMMZBIw8DaZ1knnLdL5N8wFo1EXATN3nbq6OndA2g8icf91RvaauesYWSq/LQrq+K0igZkP1AlMqfw2Uajjt4oEZj5QJzCl8ttEoY7fKhKY+UCdwJTKbxOFOn6rSGDmA3UCUyq/TRTq+K0igZkP1AlMqfw2Uajjt4oEZj5QJzCl8ttEoY7fKhKY+UCdwJTKbxOFOn6rSGDmA3UCUyq/TRTq+K0igZkP1AlMqfw2Uajjt4oEZj5QJzCl8t1E9X6BsU+yVVRU8GZy5syZ9PT0tLQ0+upvXgzaZQjY9eDwhw8f1q5de+/evdevX69evfrFixe++1MOzoQMfMNyOeni06dP27dvj4mJSUpKouvglMmPMzV513F55vz8+fMPHjygNy8/1iQgczJZHV4JWltbs7Oz6W1r3LhxvBhm+8CdRxEKM3+hx8beRX8GDlxmMx8cjpBj7jp0alNubi7tN6re/KkGrn4TMFMd3kkX9MFh8+bNZWVlc+bMgQB/ScDA22Qiwjvp4ty5c/SZ/OTJk4sWLVqxYsVfsrN8uJn3OgNUVNzrOO+izNx1BkgdpHUSgDpOGrhWIAB1FGAh1EkA6jhp4FqBANRRgIVQJwGo46SBawUCUEcBFkKdBKCOkwauFQhAHQVYCHUSgDpOGrhWIAB1FGAh1EkA6jhp4FqBANRRgIVQJwGo46SBawUCUEcBFkKdBKCOkwauFQhAHQVYCHUSgDpOGrhWIAB1FGAh1EkA6jhp4FqBANRRgIVQJwEzn4iorKwsLy93rjNy/e3bt4KCgvfv33/9+nXv3r1ZWVnRMWiRJGCmOomJicwvhj106BCdVbBz58729vZ58+a9efNGEhPCogkY+PRn30kXDQ0N/RZcWFgYHx9Pjc3NzcOGDevXix+VCJi86zCfOSc6eXl51dXVt27dWrx4sRIsBDsJmKyOc539rmtqaq5fv37hwoV+7e4/DtzTnwOXORR2X5Nary1Pf9KuQ8+W970IUnFxcXd3N13QA+d1dXVq2BD9XwIG3uvQAnknXXz58uXmzZt0JtzDhw+Tk5P/iwI/qREwU53ISRe1tbX9YOzfv58OSTlx4kRsbOzZs2f79eJHJQJmqsN7M5o0aRId5qUECME8AvjXZB4ZtAsIQB0BIHTzCEAdHhm0CwhAHQEgdPMIQB0eGbQLCEAdASB08whAHR4ZtAsIQB0BIHTzCEAdHhm0CwhAHQEgdPMIQB0eGbQLCEAdASB08whAHR4ZtAsIQB0BIHTzCEAdHhm0CwhAHQEgdPMIQB0eGbQLCEAdASB08whAHR4ZtAsIQB0BIHTzCEAdHhm0iwjQl/Ca96qoqHBfFD0yMXz4cPcY9LoTsPHBYTohJTMzs6mpqaOjQ/SXhX4uAQPfsPpOuuAtmo7eKSkpodNSeAFolyFg8q7DPOmisbHxyJEjN27c4J3BI0MNMUTAZHWiC/zjx4+FCxfSCSkTJkzwoE4oHJ3Se4vz1IggZjbzweFIPaN3HToYhW50NmzYQAGdnZ35+fmq56R4N8W4kWbuOgkJCXQgnPvdDHYdDzI7d0ozbxUjJ114QIMh8gTMfMPinXTh5MI8p9IZgGt3AmbuOu5rRq8WAlBHC0Ybk0AdG6uuZc1QRwtGG5NAHRurrmXNUEcLRhuTQB0bq65lzVBHC0Ybk0AdG6uuZc1QRwtGG5NAHRurrmXNUEcLRhuTQB0bq65lzVBHC0Ybk0AdG6uuZc1QRwtGG5NAHRurrmXNUEcLRhuTQB0bq65lzVBHC0Ybk0AdG6uuZc1mPodVWVlJD5ZHA7p9+3ZRURF9Ayh1paen07fIRsegRZKAmerwHs+7ePFiV1fXtm3bJOkgzIWAgW9YLiddtLe3jx8/3gUHuuQJmLzrRD9zXlpaSrtOS0vLmDFjjh07NnnyZHlSiOxHwMynPyOLTE5ODofDzgWHQqG0tLRTp07RYReFhYV379519gqvq6qqhDHyATk5OX3BA5c5FO77JRounM+cm6xO9K6zY8eOpKQkQpidnU33yxpYWpzCTHV+/n5F7zr0sWv9+vXkzePHj1NSUiyuu4alm6lO5KSL2trafoT27dtXUFBAdzl0BuXp06f79eJHJQJmqsM76WLKlCl0IJwSIATzCBj44Zy3VLTrJQB19PK0KBvUsajYepcKdfTytCgb1LGo2HqXCnX08rQoG9SxqNh6lwp19PK0KBvUsajYepcKdfTytCgb1LGo2HqXCnX08rQoG9SxqNh6lwp19PK0KBvUsajYepcKdfTytCgb1LGo2HqXCnX08rQoG9SxqNh6lwp19PK0KBvUsajYepdq5oPDehkhG5MAdh0mFjSKCUAdMSNEMAlAHSYWNIoJQB0xI0QwCUAdJhY0iglAHTEjRDAJQB0mFjSKCUAdMSNEMAlAHSYWNIoJQB0xI0QwCUAdJhY0iglAHTEjRDAJQB0mFjSKCUAdMSNEMAlAHSYWvzeOGDHCeWL3pk2bqGWQJw11Bhm4nl8XFxdHB4r39PRQut7e3tbWVmrRk1o6C9SRRuWzwBkzZjx58oQm9fz589TU1MGfHdQZfOZ6fmNmZmZ9fT3lov/StZ6kKlmgjgotP8VmZGTcuXOHZkRfkrJkyZLBnxrUGXzmen4jfS9TTExMW1sbpUtISNCTVCUL1FGh5bPYZcuW7dq163/ZcogE1PGZDirTycrKunTp0vLly1UGaYvFwzTaUNqWCLuObRXXtl6oow2lbYmgjm0V17ZeqKMNpW2JoI5tFde2XqijDaVtiaCObRXXtl6oow2lbYmgjm0V17ZeqKMNpW2J/gV3CIHL9vw1HwAAAABJRU5ErkJggg==",
      "text/plain": [
       "<PIL.Image.Image image mode=RGB size=190x230>"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dinkum.reset()\n",
    "\n",
    "def custom_random(X, Y):\n",
    "    import random\n",
    "    \n",
    "    is_active = random.choice([True, False])\n",
    "    print(f\"custom cis-regulatory logic: random! gene active? {is_active}\")\n",
    "    return is_active\n",
    "\n",
    "x = Gene(name='X')\n",
    "y = Gene(name='Y')\n",
    "z = Gene(name='Z')\n",
    "m = Tissue(name='M')\n",
    "\n",
    "x.custom_activation(state_fn=custom_random, delay=1)\n",
    "y.custom_activation(state_fn=custom_random, delay=1)\n",
    "z.custom_activation(state_fn=custom_random, delay=1)\n",
    "\n",
    "dinkum.run_and_display(stop=5)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "0ed084cd-352e-41a6-9ebf-21c95d889cd0",
   "metadata": {},
   "source": [
    "### Biased stochastic activation\n",
    "\n",
    "You can also use other functions from the Python random library, e.g. `randrange`, which lets you select randomly from a range of numbers.\n",
    "\n",
    "Using this you can implement a biased stochastic activation function that (for example) activates a gene 80% of the time.\n",
    "\n",
    "See [the docs for the `random` module for more information](https://docs.python.org/3/library/random.html)."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1a6ad700-94b8-4543-92dc-e608bbe03c02",
   "metadata": {},
   "source": [
    "### Stochastic but reproducible\n",
    "\n",
    "When debugging, it can be annoying to get new random behavior every time.\n",
    "\n",
    "You can use `random.seed(NUMBER)` to force the function to return random numbers reproducibly. Use `random.seed(1)` (or any other number) to initialize the [pseudo-random number generator](https://en.wikipedia.org/wiki/Pseudorandom_number_generator)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "26605325-27ac-42e0-873f-94f72d193662",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "initializing: dinkum v0.3.2\n",
      "start=1 stop=5\n",
      "custom cis-regulatory logic: random 80%! gene active? True\n",
      "custom cis-regulatory logic: random 80%! gene active? True\n",
      "custom cis-regulatory logic: random 80%! gene active? False\n",
      "custom cis-regulatory logic: random 80%! gene active? True\n",
      "custom cis-regulatory logic: random 80%! gene active? True\n",
      "custom cis-regulatory logic: random 80%! gene active? True\n",
      "custom cis-regulatory logic: random 80%! gene active? True\n",
      "custom cis-regulatory logic: random 80%! gene active? False\n",
      "custom cis-regulatory logic: random 80%! gene active? True\n",
      "custom cis-regulatory logic: random 80%! gene active? True\n",
      "custom cis-regulatory logic: random 80%! gene active? False\n",
      "custom cis-regulatory logic: random 80%! gene active? True\n"
     ]
    },
    {
     "data": {
      "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wBDAAgGBgcGBQgHBwcJCQgKDBQNDAsLDBkSEw8UHRofHh0aHBwgJC4nICIsIxwcKDcpLDAxNDQ0Hyc5PTgyPC4zNDL/2wBDAQkJCQwLDBgNDRgyIRwhMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjL/wAARCADmAL4DASIAAhEBAxEB/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQAAAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWmp6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEAAwEBAQEBAQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAECAxEEBSExBhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRomJygpKjU2Nzg5OkNERUZHSElKU1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9oADAMBAAIRAxEAPwD3+iiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKAOa8ZeebKzXTNv9ufaAdOz0DgHeW/2Nm4H6jHOKv+GjZnw7Z/YRIIQhBEv+sEmT5m/wD29+7d/tZrWooAKKKKAOQ8Vaa82taYkMoji1h/7Nv1x/rIlSSYfjiORPpKfQV1wAUAAAAcADtS0UAFFFFAHE6e3hzUvH011CdKju7GSSCJEMa3E1wQRK5A+Y7RlB7+YfQ121FFABRRRQAUUUUAFFFFABRRRQAUUUUAFFQ3V1BY2c93cyLFbwRtLLI3RVUZJP0AqagAoqhpWtadrcdxJptytxHbzeTI6qQN+1W4JHzDaynIyDnrV+gAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAx/Femyax4Q1nTYYklnubKaKJHxgyFCF68D5sc9qgstEuItWhv0n+z2ioAthmYeUNm0LgTmLj2THHHPNb9FAGRo9jPaanr8sseyO6v0mgOQdyC2gTPHT5kYc+npiteiigDn/ABPqV3p32X7LN5fmb93yg5xtx1Hua5//AISXV/8An7/8hp/hWr40/wCXH/tp/wCy1ylAGt/wkur/APP3/wCQ0/wrQ0PXNRvNYggnuN8TbsrsUZwpPYe1czWt4a/5GC1/4H/6AaANDXNc1Gz1ieCC42RLtwuxTjKg9x71n/8ACS6v/wA/f/kNP8KPEv8AyMF1/wAA/wDQBWTQBrf8JLq//P3/AOQ0/wAK6DwxqV3qP2r7VN5nl7NvygYzuz0HsK4mur8F/wDL9/2z/wDZqAMr/hJdX/5+/wDyGn+FH/CS6v8A8/f/AJDT/CsmigDW/wCEl1f/AJ+//Iaf4V2Gh3M15o8E8775W3ZbAGcMR2+lec16D4a/5F+1/wCB/wDoZoA1aKKKACiiigAooooAKKKKACiqGuarHoeg3+qzRySx2du87JGCWYKCcDH069B1PFJFrulTXK2qX9ubxl3C0LhZ/u7sGM/MDjnBGaANCisbRtaudRvLq0vLAWc8EUM20TeZ8sm7Ct8o2uNhyvIGRyc1s0Acp40/5cf+2n/stcpXoOtaL/bHkf6R5Plbv4N2c49x6Vk/8IX/ANRD/wAg/wD2VAHKVreGv+Rgtf8Agf8A6Aa1f+EL/wCoh/5B/wDsqt6b4Y/s7UIrr7Z5nl5+XysZyCOufegDn/Ev/IwXX/AP/QBWTXbal4Y/tHUJbr7Z5fmY+Xys4wAOufaqn/CF/wDUQ/8AIP8A9lQByldX4L/5fv8Atn/7NR/whf8A1EP/ACD/APZVraLov9j+f/pHnebt/g24xn3PrQB59RXV/wDCF/8AUQ/8g/8A2VH/AAhf/UQ/8g//AGVAHKV6D4a/5F+1/wCB/wDoZrJ/4Qv/AKiH/kH/AOyroNNsv7O0+K18zzPLz823Gcknp+NAFuiiigAooooAKKKKACiiigCjrWmjWdB1HS2lMS3ttJbmQDJTepXOO+M0xNKC3AvWuJze7OSLiYQFtuM+T5hUDvj9c81o0UAY3h3SL3R7edL2+t72aaTzZJ47Zonkc/eZsu2eAoAGAAABxjGzRRQBw3xF/wCYb/21/wDZK4au5+Iv/MN/7a/+yVw1fV5b/usfn+bPzzPf+RhU+X/pKCt3wd/yNVl/wP8A9AasKt3wd/yNVl/wP/0Bq3xX8Cfo/wAjjy//AHul/ij+aDxj/wAjVe/8A/8AQFrCrd8Y/wDI1Xv/AAD/ANAWsKjC/wACHovyDMP97q/4pfmwrufh1/zEv+2X/s9cNXc/Dr/mJf8AbL/2esMy/wB1l8vzR2ZF/wAjCn8//SWYVFFFfKH6GFeg+Gv+Rftf+B/+hmvPq9B8Nf8AIv2v/A//AEM0AatFFFABRRRQAUUUUAFFFFABRWR4q1C60rwlrGo2QjNza2Us8fmH5QVQnJ4OenTv0461VPimAam2kzCG3vVjJdzdwPHGwj3k7fMEhUeuwcc8DmgDoaK5zw9dal/bWradqFzNOLVIGja4SJXfdvBdRGMeWSnyg/NlWz2ro6AOG+Iv/MN/7a/+yVw1ereIv7B/0b+2/wDa8n7/ALbvu/h1rC/4oP8Az59fQYHFezoRjySe+y03PjM2y722MnU9rCN7aOVnsuljhq3fB3/I1WX/AAP/ANAatz/ig/8APn1e0f8A4RL+1Yf7M/4/Pm8v/W/3Tn73HTNa18ZzUpR9nJXT6eRz4PLOTEU5e2pu0lopa79NNzlfGP8AyNV7/wAA/wDQFrCr0fWP+ES/tWb+0/8Aj8+XzP8AW/3Rj7vHTFUf+KD/AM+fRQxnLSjH2cnZLp5BjMs58RUl7amryejlrv103OGrufh1/wAxL/tl/wCz0f8AFB/58+t3w7/YP+k/2J/s+d9/32/e/HpWWOxXtKEo8kltutNzoynLvY4yFT2sJWvopXez6WOForpPtngr/np+k1H2zwV/z0/SavG+q1/5H9zPqP7Qwn/P2P8A4Ev8zm69B8Nf8i/a/wDA/wD0M1hfbPBX/PT9Jq6bSnspNNhbTjm0O7Yfm9Tn73PXNROjUgrzi16o0pYqhVfLTmpPyaZcooorM3CiiigAooooAKKKKAILy0g1CxuLK6jEtvcRtFKhJAZGGCOPUE0v2WE2f2R1MkBj8pllYvuXGMMWyWyOuetTUUAUNN0ax0nzPskcgaQKHeWZ5XIX7o3OScDJwM4GT61foooA4b4i/wDMN/7a/wDslcNXc/EX/mG/9tf/AGSuGr6vLf8AdY/P82fnme/8jCp8v/SUFbvg7/karL/gf/oDVhVu+Dv+Rqsv+B/+gNW+K/gT9H+Rx5f/AL3S/wAUfzQeMf8Akar3/gH/AKAtYVbvjH/kar3/AIB/6AtYVGF/gQ9F+QZh/vdX/FL82Fdz8Ov+Yl/2y/8AZ64au5+HX/MS/wC2X/s9YZl/usvl+aOzIv8AkYU/n/6SzhqKKK7jyAr1bwd/yKtl/wAD/wDQ2rymvVvB3/Iq2X/A/wD0Nq8rOP4C9f0Z9Fwz/vcv8L/NG5RRRXzZ9wFFFFABRRRQAUUUUAFFYXjYTf8ACCeIDBLJFKNOuGRovvZEbHjjv09fTBplreaydTiVEN3pXlhkuhDFmYeXkNvE46n0iA57DmgDoKK4XwgdRPjXW5NTsbm3vLmws5phJJGyIfMuQqLtc/KFwAeMlGJAJGe6oA4b4i/8w3/tr/7JXDV6t4i8Rf2B9m/0Xz/P3f8ALTbt249j61hf8LF/6hX/AJMf/Y19Bga1eNCKhTutdbpdT4zNsNg54ycqtfllppyt9F1OGrd8Hf8AI1WX/A//AEBq3P8AhYv/AFCv/Jj/AOxq9o/jT+1tVhsf7P8AK83d8/nbsYUnptHpWteviXSkpUrKz+0uxz4PC4GOIpuOIu+ZWXI1fXY5Xxj/AMjVe/8AAP8A0Bawq9S1LxP/AGdqEtr9j8zy8fN5uM5APTHvVT/hNP8AqH/+Rv8A7GuOlm/s4RhybJLf/gHqYjhv21aVX2tuZt/D3fqecV3Pw6/5iX/bL/2er3/Caf8AUP8A/I3/ANjWtoutf2x5/wDo/k+Vt/j3Zzn2HpWeJzP29J0+S1/P/gGuByD6piI1/aXtfS1t013PIqK7n/hYv/UK/wDJj/7Gj/hYv/UK/wDJj/7GvW+sYr/nz/5Mj576nl//AEE/+SSOGr1bwd/yKtl/wP8A9DasL/hYv/UK/wDJj/7Gur0fUf7W0qG+8ryvN3fJu3YwxHXA9K8/MqtadJKpT5Vfun0Z7OR0MLTxEnRrc75duVrqtdS9RRRXiH1QUUUUAFFFFABRRRQAUUUUAFFFFAHDfEX/AJhv/bX/ANkrhq7n4i/8w3/tr/7JXDV9Xlv+6x+f5s/PM9/5GFT5f+koK3fB3/I1WX/A/wD0Bqwq3fB3/I1WX/A//QGrfFfwJ+j/ACOPL/8Ae6X+KP5o1fEv/IwXX/AP/QBWTWt4l/5GC6/4B/6AKya+NP08K6vwX/y/f9s//Zq5Sur8F/8AL9/2z/8AZqAPOKKKK+5PyUK9W8Hf8irZf8D/APQ2rymvVvB3/Iq2X/A//Q2rys4/gL1/Rn0XDP8Avcv8L/NG5RRRXzZ9wFFFFABRRRQAUUUUAQ3V1BY2c93cyLFbwRtLLI3RVUZJP0Ap8siQxPLI21EUsx9AOtZXivTZNY8IazpsMSSz3NlNFEj4wZChC9eB82Oe1Z03hy7maWaGbyLaS1aNNP3zJ5eYigTiZolwT2jIGOOeaANzTdVtdWhM1qLny8Agz2ssO4EZBG9Rn8Ku1zvhez1C0kuvPhvbeyMcKw297cieRZAG8xgwZgFIKADPVTwM10VAHKeNNHv9W+w/YYPN8rzN/wA6rjO3HUj0Ncr/AMIdr3/Ph/5GT/4quq8aaxf6T9h+wz+V5vmb/kVs4246g+prlf8AhMde/wCf/wD8gp/8TX0GB+tewj7Pltrve+58Zm39nfXJ+35+bS9uW2y767B/wh2vf8+H/kZP/iq1fDXhrV9P8QWt1dWnlwpv3N5iHGUIHAOeprK/4THXv+f/AP8AIKf/ABNavhrxLq+oeILW1urvzIX37l8tBnCEjkDPUVrX+ueylzctrPv2OfB/2Z9Yp+z9pzcytflte+l/I0Nc0PUbzWJ54LffE23Db1GcKB3PtWf/AMI1q/8Az6f+RE/xrQ1zXNRs9YngguNkS7cLsU4yoPce9Z//AAkur/8AP3/5DT/Cvmj7sP8AhGtX/wCfT/yIn+NdB4Y027077V9qh8vzNm35gc43Z6H3Fc//AMJLq/8Az9/+Q0/wroPDGpXeo/avtU3meXs2/KBjO7PQewoA4r/hDte/58P/ACMn/wAVR/wh2vf8+H/kZP8A4qj/AITHXv8An/8A/IKf/E0f8Jjr3/P/AP8AkFP/AImvq/8Abv7n4n55/wAJP/Tz/wAlD/hDte/58P8AyMn/AMVXoHhqzuNP8P2trdR+XMm/cuQcZckcjjoa8/8A+Ex17/n/AP8AyCn/AMTXoHhq8uNQ8P2t1dSeZM+/c2AM4cgcDjoK8/MvrPsl7Xltfpfs+57OR/UfrEvq3Pzcv2rWtddupq0UUV4h9UFFFFABRRRQAUUUUAFFFFABRRRQBw3xF/5hv/bX/wBkrhq7n4i/8w3/ALa/+yVw1fV5b/usfn+bPzzPf+RhU+X/AKSgrd8Hf8jVZf8AA/8A0Bqwq3fB3/I1WX/A/wD0Bq3xX8Cfo/yOPL/97pf4o/mjV8S/8jBdf8A/9AFZNa3iX/kYLr/gH/oArJr40/Twrq/Bf/L9/wBs/wD2auUrq/Bf/L9/2z/9moA84ooor7k/JQr1bwd/yKtl/wAD/wDQ2rymvVvB3/Iq2X/A/wD0Nq8rOP4C9f0Z9Fwz/vcv8L/NG5RRRXzZ9wFFFFABRRRQAUUUUAFFFFABRRRQBh+IvDv9v/Zv9K8jyN3/ACz3bt2PcelYX/Cuv+or/wCS/wD9lXc0V10sdXpRUISsl5I86vlODxFR1asLyfm/8zhv+Fdf9RX/AMl//sqvaP4L/snVYb7+0PN8rd8nk7c5UjruPrXV0U5ZhiZxcZS0fkv8iKeTYGnNTjDVO61fT5nP6l4Y/tHUJbr7Z5fmY+Xys4wAOufaqn/CF/8AUQ/8g/8A2VdXRXGeocp/whf/AFEP/IP/ANlWtoui/wBj+f8A6R53m7f4NuMZ9z61q0UAcN/wrr/qK/8Akv8A/ZUf8K6/6iv/AJL/AP2VdzRXd/aWK/m/Bf5Hk/2Fl/8Az7/GX+Zw3/Cuv+or/wCS/wD9lXV6Pp39k6VDY+b5vlbvn27c5Ynpk+tXqKyrYutWjy1JXXyOjDZbhcLNzoxs7W3b/NhRRRXMdwUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFAHGy+Lr6HWfE9g0Fv/oIjTTeGzNI0cRKvzz888Q4xw1b/AIc1GbV/DGk6ncKiz3llDcSLGCFDOgYgAknGT6mqs3hPT59X/tJ5LjzjeC8KhhtZhHGgUjH3f3MT4zncoOccVpaVp0OkaPZaZbs7Q2dvHbxtIQWKooUE4AGcD0oAt0UUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFAAAUUUAf/9k=",
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAL4AAADmCAIAAADky6hQAAAL40lEQVR4Ae2dW0hUTxzHWyutDJOihyJ7iTAiNSq6kJJQaWQpRWg3LSGxhwzFHoSgoJtSdIGeKoJuFHSxDH3QoraQIgrqQe0CQighUUGZFoHm/1cL22l3Zs6Z0+z5nzPnuw9/zs78Zs7M5/dx9mz8ZzYwNDQ0DC8QkCcQJ98ELUDgFwGoAw9sEoA6NsGhGdSBAzYJQB2b4NAM6sABmwSgjk1waAZ14IBNAlDHJjg0gzpwwCYBqGMTHJpBHThgkwDUsQkOzaAOHLBJAOrYBIdmUEfCgebm5rVr14YatLe3z58/f3BwUKI9P7SjoyMtLe3nz5+hkHnz5rW1tfHD5Wqqq6uzf7+ysrLi4+PlGgui6X/1wss6gby8vAcPHlD88uXLHz9+bL2haeT69euvXLlCYY2NjYWFhabxNgJOnDhRU1NjoyGzyTBmKQp5BF69erVw4cLr169v27aNF2OvnBaeWbNm0TK2ePFiWtLsdSJo9e7du/T09P7+fkGMVFWAogVrEqqiCVRWVtbX1z9//nzChAnRtf9SsmnTpjFjxlB2L1++/C/9MNsWFRWVlJTQqsmstVE4wkYbnzfp7e0dMWJEX1+fcnX27Nkzc+ZMWn6UE6anNFrPFHrza4RSaxSCW1tb6Um5paUlPz8/FjTGjRunvNvv37/Pnj27u7tbbc/4hiXxFz4wMFBVVXXkyBF6RqaFp6GhQaLx/xd68ODB4uLiKVOmqB0CnnUkeB47duzjx4+HDh2iNl1dXbm5uc+ePUtMTJTowiw0OTn58+fPZlFy9QkJCfTvCMOHDw81o29wY8eOleuCFQ11WFRQZoEAPrAsQEIIiwDUYVFBmQUCUMcCJISwCEAdFhWUWSAAdSxAQgiLANRhUUGZBQJQxwIkhLAIQB0WFZRZIAB1LEBCCIsA1GFRQZkFAlDHAiSEsAhAHRYVlFkgAHUsQEIIiwDUYVFBmQUCUMcCJISwCEAdFhWUWSAAdSxAQgiLANRhUUGZBQJQxwIkhLAIQB0WFZRZIAB1LEBCCIuAnurU1dWxJvurjLaqbN26lXbK8QJQbpGA79RZvXr13LlzA4GARUAI4xHQUJ29e/fShvCcnBzmnK9du1ZRUcGsQqEUAT238IX2UL5+/bq8vDyMIxgMhq9jscky3LlPLnRWR5BCe+pcvXpV0KdsFZ2/FG7ixZ51PiRFsOqEc4YL2wT0VIcO5aNXamqq8UPKNiM0ZBLQ8DGZ5knnLdL5N8wJo1AVAT1XnaamJjEg5QeRiG+nZa2eq46WqXLbpKCO2zLimfFAHc+kym0DhTpuy4hnxgN1PJMqtw0U6rgtI54ZD9TxTKrcNlCo47aMeGY8UMczqXLbQKGO2zLimfFAHc+kym0DhTpuy4hnxgN1PJMqtw0U6rgtI54ZD9TxTKrcNlCo47aMeGY8UMczqXLbQKGO2zLimfFAHc+kym0DhTpuy4hnxgN1PJMq1w1U7Q8Yu6S32tpa3kjOnj2bmZmZkZFBP/3Ni0G5FQL+2jj84cOHdevW3b9//82bN2vWrHn58qXr/pS9MyANP7AEJ118+vRpx44dcXFxKSkpdO2dNLlxpDqvOoI95+fPn3/48CF9eLkxJx4Zk87q8FLQ2dlZUFBAH1sTJ07kxTDLA0Fmsc3Coew/DWPXc+zO0NBz43AoJ8xVh05tKioqovVG1ps/ecbVbwJ6qsM76YK+OGzZsqW6unrBggUQ4B8JaPiYTER4J12cO3eOvpOfOnUqOzt71apV/8jO5831fNaJUVJj90QSu55j96yj56oTI3XQrZEA1DHSwLUEAagjAQuhRgJQx0gD1xIEoI4ELIQaCUAdIw1cSxCAOhKwEGokAHWMNHAtQQDqSMBCqJEA1DHSwLUEAagjAQuhRgJQx0gD1xIEoI4ELIQaCUAdIw1cSxCAOhKwEGokAHWMNHAtQQDqSMBCqJEA1DHSwLUEAagjAQuhRgJQx0gD1xIEoI4ELIQaCei5I6Kurq6mpsY4z9D1t2/fSktL379///Xr13379uXl5UXHoMQiAT3VSU5OZv4w7OHDh+msgl27dvX09CxatOjt27cWMSEsmoCGuz/DJ120tLRETLisrCwxMZEK29vbR44cGVGLt1IEdF51mHvOiU5xcXF9ff3t27eXLl0qBQvBRgI6q2OcZ8R1Q0PDjRs3Lly4EFEufhsIiuvlap056SJ2Y9b5GxatOrS3PPyixFZUVAwMDNAFbThvamqSSzWi/yag4bMOTZB30sWXL19u3bpFZ8I9evQoNTX1bxR4J0dAzw8s+tYdCAQaGxsjYHR1ddEhKSRWfHz8yZMnZ8yYEREgfhu7xd+TPdORM2JeqA0T8GSCg+HhK7gwPp/p/KyjABW64BOAOnw2qBESgDpCPKjkE4A6fDaoERKAOkI8qOQTgDp8NqgREoA6Qjyo5BOAOnw2qBESgDpCPKjkE4A6fDaoERKAOkI8qOQTgDp8NqgREoA6Qjyo5BOAOnw2qBESgDpCPKjkE4A6fDaoERKAOkI8qOQTgDp8NqgREoA6Qjyo5BOAOnw2qBESgDpCPKgUEKAdEfq9amtrxZOi/XujRo0Sx6BWTEDPfVi8ky5Cf0J0Qkpubm5bW1tvb6/gjwpVYgIafmCFT7rgzZyO3qmsrKTTUngBKLdCQOdVh3nSRWtr69GjR2/evClemayw83mMzupEp/bHjx9LliyhE1ImT55sQ51AMLpL+yXGnZSx+zXy2I1Zz+MKQvmMXnXoYBR60Nm4cSMF9PX1lZSUyJ6TYt8U7VrqueokJSXRgXDipxmsOjZkNq6Uej4qZmVl5efn20CDJtYJ6PmBZeXYJeY5ldbBIVLPVQd5dYAA1HEAsp63gDp65tWBWUEdByDreQuoo2deHZgV1HEAsp63gDp65tWBWUEdByDreQuoo2deHZgV1HEAsp63gDp65tWBWUEdByDreQuoo2deHZgV1HEAsp63gDp65tWBWUEdByDreQuoo2deHZgV1HEAsp63gDp65tWBWUEdByDreQuoo2deHZiVnvuw6urqaGN5NL47d+6Ul5dPnTqVqjIzMw8cOBAdgxKLBPRUh7c97+LFi/39/du3b7dIB2ECAhp+YAlOuujp6Zk0aZIAB6qsE9B51Ynec15VVUWrTkdHx/jx448fPz5t2jTrpBAZQUDP3Z+hSaampgaDQeOEA4FARkbG6dOn6bCLsrKye/fuGWtNrwN/dWYabhJg3L+Nky5MYDlcHb3q7Ny5MyUlhYZRUFBAz8sOj0ez2+m56vz8/Ypedehr14YNG8ibJ0+epKWlaZZLh6ejpzqhky4aGxsjaO7fv7+0tJSecugMyjNnzkTU4q0UAT3V4Z10MX36dDoQTgoQgnkENPxyzpsqytUSgDpqefqoN6jjo2SrnSrUUcvTR71BHR8lW+1UoY5anj7qDer4KNlqpwp11PL0UW9Qx0fJVjtVqKOWp496gzo+SrbaqUIdtTx91BvU8VGy1U4V6qjl6aPeoI6Pkq12qlBHLU8f9QZ1fJRstVOFOmp5+qg3qOOjZKudKtRRy9NHvUEdHyVb7VT13DislhF6YxLAqsPEgkJzAlDHnBEimASgDhMLCs0JQB1zRohgEoA6TCwoNCcAdcwZIYJJAOowsaDQnADUMWeECCYBqMPEgkJzAlDHnBEimASgDhMLCs0JQB1zRohgEoA6TCwoNCcAdcwZIYJJAOowsbi9cPTo0YWFheFRbt68mUrCb525gDrOcFZ8l4SEBDpQfHBwkPodGhrq7OykEsX3MOsO6pgRcmv9nDlznj59SqN78eJFenq688OEOs4zV3PH3Nzc5uZm6ov+S9dqOpXpBerI0HJTbE5Ozt27d2lE9CMpy5Ytc35oUMd55mruSL/LFBcX193dTd0lJSWp6VSmF6gjQ8tlsStWrNi9e/f/suQQCajjMh1khpOXl3fp0qWVK1fKNFIWi800ylD6rSOsOn7LuLL5Qh1lKP3WEdTxW8aVzRfqKEPpt46gjt8yrmy+UEcZSr91BHX8lnFl84U6ylD6rSOo47eMK5sv1FGG0m8d/QcZsY/+L2zNJAAAAABJRU5ErkJggg==",
      "text/plain": [
       "<PIL.Image.Image image mode=RGB size=190x230>"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dinkum.reset()\n",
    "\n",
    "import random\n",
    "random.seed(1)\n",
    "\n",
    "def custom_random_biased(X, Y):\n",
    "    # activate 20% of the time\n",
    "    is_active = random.randrange(0, 100) < 80\n",
    "    print(f\"custom cis-regulatory logic: random 80%! gene active? {is_active}\")\n",
    "    return is_active\n",
    "\n",
    "x = Gene(name='X')\n",
    "y = Gene(name='Y')\n",
    "z = Gene(name='Z')\n",
    "m = Tissue(name='M')\n",
    "\n",
    "x.custom_activation(state_fn=custom_random_biased, delay=1)\n",
    "y.custom_activation(state_fn=custom_random_biased, delay=1)\n",
    "z.custom_activation(state_fn=custom_random_biased, delay=1)\n",
    "\n",
    "dinkum.run_and_display(stop=5)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c6191902-2ecf-470a-9628-771fadfab0f5",
   "metadata": {},
   "source": [
    "## Using lambda functions\n",
    "\n",
    "Python supports defining functions using `lambda` instead of `def`, and you can use this syntax too; see below.\n",
    "\n",
    "This can make things shorter, at the expense of readability. It's your call!\n",
    "\n",
    "Here's a simple example of the `AND` logic used above:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "e7fc53be-78e8-4938-bda5-b37d3030720d",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "initializing: dinkum v0.3.2\n",
      "start=1 stop=5\n"
     ]
    },
    {
     "data": {
      "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wBDAAgGBgcGBQgHBwcJCQgKDBQNDAsLDBkSEw8UHRofHh0aHBwgJC4nICIsIxwcKDcpLDAxNDQ0Hyc5PTgyPC4zNDL/2wBDAQkJCQwLDBgNDRgyIRwhMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjIyMjL/wAARCADmAL4DASIAAhEBAxEB/8QAHwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUFBAQAAAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkKFhcYGRolJicoKSo0NTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWmp6ipqrKztLW2t7i5usLDxMXGx8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEAAwEBAQEBAQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAECAxEEBSExBhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRomJygpKjU2Nzg5OkNERUZHSElKU1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOEhYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9oADAMBAAIRAxEAPwD3+iiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKAOa8ZeebKzXTNv9ufaAdOz0DgHeW/2Nm4H6jHOKv+GjZnw7Z/YRIIQhBEv+sEmT5m/wD29+7d/tZrWooAKKKKAOQ8Vaa82taYkMoji1h/7Nv1x/rIlSSYfjiORPpKfQV1wAUAAAAcADtS0UAFFFFAHE6e3hzUvH011CdKju7GSSCJEMa3E1wQRK5A+Y7RlB7+YfQ121FFABRRRQAUUUUAFFFFABRRRQAUUUUAFFQ3V1BY2c93cyLFbwRtLLI3RVUZJP0AqagAoqhpWtadrcdxJptytxHbzeTI6qQN+1W4JHzDaynIyDnrV+gAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAx/Femyax4Q1nTYYklnubKaKJHxgyFCF68D5sc9qgstEuItWhv0n+z2ioAthmYeUNm0LgTmLj2THHHPNb9FAGRo9jPaanr8sseyO6v0mgOQdyC2gTPHT5kYc+npiteiigDlPGmsX+k/YfsM/leb5m/5FbONuOoPqa5X/AITHXv8An/8A/IKf/E1ufEX/AJhv/bX/ANkrhq+ly+hSnhoylFN69F3Z8JnOMxFPHTjCpJLTRNrojd/4THXv+f8A/wDIKf8AxNavhrxLq+oeILW1urvzIX37l8tBnCEjkDPUVxtbvg7/AJGqy/4H/wCgNW+Jw1FUZtQWz6LscuBx2KliqcZVJNOS6vudLrmuajZ6xPBBcbIl24XYpxlQe496z/8AhJdX/wCfv/yGn+FHiX/kYLr/AIB/6AKya+TP0U1v+El1f/n7/wDIaf4V0HhjUrvUftX2qbzPL2bflAxndnoPYVxNdX4L/wCX7/tn/wCzUAZX/CS6v/z9/wDkNP8ACj/hJdX/AOfv/wAhp/hWTRQBrf8ACS6v/wA/f/kNP8K7DQ7ma80eCed98rbstgDOGI7fSvOa9B8Nf8i/a/8AA/8A0M0AatFFFABRRRQAUUUUAFFFFABRVDXNVj0PQb/VZo5JY7O3edkjBLMFBOBj6deg6niki13SprlbVL+3N4y7haFws/3d2DGfmBxzgjNAGhRWNo2tXOo3l1aXlgLOeCKGbaJvM+WTdhW+UbXGw5XkDI5Oa2aAOG+Iv/MN/wC2v/slcNXq3iLw7/b/ANm/0ryPI3f8s927dj3HpWF/wrr/AKiv/kv/APZV9BgcdQpUIwnKzV+j7nxmbZTjMRjJ1aULxduq7LzOGrd8Hf8AI1WX/A//AEBq3P8AhXX/AFFf/Jf/AOyq9o/gv+ydVhvv7Q83yt3yeTtzlSOu4+ta18ww06Uoxlq0+j7ehz4PJsdTxFOcoaKSb1XR+pleJf8AkYLr/gH/AKAKya7bUvDH9o6hLdfbPL8zHy+VnGAB1z7VU/4Qv/qIf+Qf/sq+aPuzlK6vwX/y/f8AbP8A9mo/4Qv/AKiH/kH/AOyrW0XRf7H8/wD0jzvN2/wbcYz7n1oA8+orq/8AhC/+oh/5B/8AsqP+EL/6iH/kH/7KgDlK9B8Nf8i/a/8AA/8A0M1k/wDCF/8AUQ/8g/8A2VdBptl/Z2nxWvmeZ5efm24zkk9PxoAt0UUUAFFFFABRRRQAUUUUAUda00azoOo6W0piW9tpLcyAZKb1K5x3xmmJpQW4F61xOb3ZyRcTCAttxnyfMKgd8frnmtGigDG8O6Re6Pbzpe31vezTSebJPHbNE8jn7zNl2zwFAAwAAAOMY2aKKAOG+Iv/ADDf+2v/ALJXDV3PxF/5hv8A21/9krhq+ry3/dY/P82fnme/8jCp8v8A0lBW74O/5Gqy/wCB/wDoDVhVu+Dv+Rqsv+B/+gNW+K/gT9H+Rx5f/vdL/FH80HjH/kar3/gH/oC1hVu+Mf8Akar3/gH/AKAtYVGF/gQ9F+QZh/vdX/FL82Fdz8Ov+Yl/2y/9nrhq7n4df8xL/tl/7PWGZf7rL5fmjsyL/kYU/n/6SzCooor5Q/Qwr0Hw1/yL9r/wP/0M159XoPhr/kX7X/gf/oZoA1aKKKACiiigAooooAKKKKACisjxVqF1pXhLWNRshGbm1spZ4/MPygqhOTwc9Onfpx1qqfFMA1NtJmENverGS7m7geONhHvJ2+YJCo9dg454HNAHQ0Vznh661L+2tW07ULmacWqQNG1wkSu+7eC6iMY8slPlB+bKtntXR0AcN8Rf+Yb/ANtf/ZK4avVvEX9g/wCjf23/ALXk/f8Abd938OtYX/FB/wCfPr6DA4r2dCMeST32Wm58Zm2Xe2xk6ntYRvbRys9l0scNW74O/wCRqsv+B/8AoDVuf8UH/nz6vaP/AMIl/asP9mf8fnzeX/rf7pz97jpmta+M5qUo+zkrp9PI58HlnJiKcvbU3aS0Utd+mm5yvjH/AJGq9/4B/wCgLWFXo+sf8Il/as39p/8AH58vmf63+6Mfd46Yqj/xQf8Anz6KGM5aUY+zk7JdPIMZlnPiKkvbU1eT0ctd+um5w1dz8Ov+Yl/2y/8AZ6P+KD/z59bvh3+wf9J/sT/Z877/AL7fvfj0rLHYr2lCUeSS23Wm50ZTl3scZCp7WErX0UrvZ9LHC0V0n2zwV/z0/Saj7Z4K/wCen6TV431Wv/I/uZ9R/aGE/wCfsf8AwJf5nN16D4a/5F+1/wCB/wDoZrC+2eCv+en6TV02lPZSabC2nHNod2w/N6nP3ueuaidGpBXnFr1RpSxVCq+WnNSfk0y5RRRWZuFFFFABRRRQAUUUUAQXlpBqFjcWV1GJbe4jaKVCSAyMMEceoJpfssJs/sjqZIDH5TLKxfcuMYYtktkdc9amooAoabo1jpPmfZI5A0gUO8szyuQv3Ruck4GTgZwMn1q/RRQBw3xF/wCYb/21/wDZK4au5+Iv/MN/7a/+yVw1fV5b/usfn+bPzzPf+RhU+X/pKCt3wd/yNVl/wP8A9AasKt3wd/yNVl/wP/0Bq3xX8Cfo/wAjjy//AHul/ij+aDxj/wAjVe/8A/8AQFrCrd8Y/wDI1Xv/AAD/ANAWsKjC/wACHovyDMP97q/4pfmwrufh1/zEv+2X/s9cNXc/Dr/mJf8AbL/2esMy/wB1l8vzR2ZF/wAjCn8//SWcNRRRXceQFereDv8AkVbL/gf/AKG1eU16t4O/5FWy/wCB/wDobV5WcfwF6/oz6Lhn/e5f4X+aNyiiivmz7gKKKKACiiigAooooAKKwvGwm/4QTxAYJZIpRp1wyNF97IjY8cd+nr6YNMtbzWTqcSohu9K8sMl0IYszDy8ht4nHU+kQHPYc0AdBRXC+EDqJ8a63Jqdjc295c2FnNMJJI2RD5lyFRdrn5QuADxkoxIBIz3VAHDfEX/mG/wDbX/2SuGr1bxF4i/sD7N/ovn+fu/5abdu3HsfWsL/hYv8A1Cv/ACY/+xr6DA1q8aEVCnda63S6nxmbYbBzxk5Va/LLTTlb6LqcNW74O/5Gqy/4H/6A1bn/AAsX/qFf+TH/ANjV7R/Gn9rarDY/2f5Xm7vn87djCk9No9K1r18S6UlKlZWf2l2OfB4XAxxFNxxF3zKy5Gr67HK+Mf8Akar3/gH/AKAtYVej6x40/snVZrH+z/N8rb8/nbc5UHptPrVH/hYv/UK/8mP/ALGihXxKpRUaV1ZfaXYMZhcDLEVHLEWfM7rkbtrscNXc/Dr/AJiX/bL/ANno/wCFi/8AUK/8mP8A7Gt3w74i/t/7T/ovkeRt/wCWm7duz7D0rLHVq8qElOnZaa3T6nRlOGwcMZCVKvzS105Wuj6nlNFdz/wsX/qFf+TH/wBjR/wsX/qFf+TH/wBjXX9YxX/Pn/yZHn/U8v8A+gn/AMkkcNXq3g7/AJFWy/4H/wChtWF/wsX/AKhX/kx/9jXV6PqP9raVDfeV5Xm7vk3bsYYjrgelefmVWtOklUp8qv3T6M9nI6GFp4iTo1ud8u3K11WupeooorxD6oKKKKACiiigAooooAKKKKACiiigDhviL/zDf+2v/slcNXc/EX/mG/8AbX/2SuGr6vLf91j8/wA2fnme/wDIwqfL/wBJQVu+Dv8AkarL/gf/AKA1YVbvg7/karL/AIH/AOgNW+K/gT9H+Rx5f/vdL/FH80HjH/kar3/gH/oC1hVu+Mf+Rqvf+Af+gLWFRhf4EPRfkGYf73V/xS/NhXc/Dr/mJf8AbL/2euGrufh1/wAxL/tl/wCz1hmX+6y+X5o7Mi/5GFP5/wDpLOGoooruPICvVvB3/Iq2X/A//Q2rymvVvB3/ACKtl/wP/wBDavKzj+AvX9GfRcM/73L/AAv80blFFFfNn3AUUUUAFFFFABRRRQBDdXUFjZz3dzIsVvBG0ssjdFVRkk/QCnyyJDE8sjbURSzH0A61leK9Nk1jwhrOmwxJLPc2U0USPjBkKEL14HzY57VnTeHLuZpZoZvItpLVo00/fMnl5iKBOJmiXBPaMgY455oA3NN1W11aEzWoufLwCDPayw7gRkEb1Gfwq7XO+F7PULSS68+G9t7IxwrDb3tyJ5FkAbzGDBmAUgoAM9VPAzXRUAcp400e/wBW+w/YYPN8rzN/zquM7cdSPQ1yv/CHa9/z4f8AkZP/AIquq8aaxf6T9h+wz+V5vmb/AJFbONuOoPqa5X/hMde/5/8A/wAgp/8AE19BgfrXsI+z5ba73vufGZt/Z31yft+fm0vbltsu+uwf8Idr3/Ph/wCRk/8Aiq1fDXhrV9P8QWt1dWnlwpv3N5iHGUIHAOeprK/4THXv+f8A/wDIKf8AxNavhrxLq+oeILW1urvzIX37l8tBnCEjkDPUVrX+ueylzctrPv2OfB/2Z9Yp+z9pzcytflte+l/IPEvhrV9Q8QXV1a2nmQvs2t5iDOEAPBOeorK/4Q7Xv+fD/wAjJ/8AFVq+JfEur6f4gurW1u/LhTZtXy0OMoCeSM9TWV/wmOvf8/8A/wCQU/8AiaKH1z2UeXltZd+wYz+zPrFT2ntObmd7ctr31t5B/wAIdr3/AD4f+Rk/+KrqvBej3+k/bvt0Hleb5ez51bON2ehPqK5X/hMde/5//wDyCn/xNdV4L1i/1b7d9un83yvL2fIq4zuz0A9BWWO+tewl7Tltpte+50ZT/Z31yHsOfm1tflts+2uxyv8Awh2vf8+H/kZP/iqP+EO17/nw/wDIyf8AxVH/AAmOvf8AP/8A+QU/+Jo/4THXv+f/AP8AIKf/ABNdf+3f3PxPP/4Sf+nn/kof8Idr3/Ph/wCRk/8Aiq9A8NWdxp/h+1tbqPy5k37lyDjLkjkcdDXn/wDwmOvf8/8A/wCQU/8Aia9A8NXlxqHh+1urqTzJn37mwBnDkDgcdBXn5l9Z9kva8tr9L9n3PZyP6j9Yl9W5+bl+1a1rrt1NWiiivEPqgooooAKKKKACiiigAooooAKKKKAOG+Iv/MN/7a/+yVw1dz8Rf+Yb/wBtf/ZK4avq8t/3WPz/ADZ+eZ7/AMjCp8v/AElBW74O/wCRqsv+B/8AoDVhVu+Dv+Rqsv8Agf8A6A1b4r+BP0f5HHl/+90v8UfzQeMf+Rqvf+Af+gLWFW74x/5Gq9/4B/6AtYVGF/gQ9F+QZh/vdX/FL82Fdz8Ov+Yl/wBsv/Z64au5+HX/ADEv+2X/ALPWGZf7rL5fmjsyL/kYU/n/AOks4aiiiu48gK9W8Hf8irZf8D/9DavKa9W8Hf8AIq2X/A//AENq8rOP4C9f0Z9Fwz/vcv8AC/zRuUUUV82fcBRRRQAUUUUAFFFFABRRRQAUUUUAYfiLw7/b/wBm/wBK8jyN3/LPdu3Y9x6Vhf8ACuv+or/5L/8A2VdzRXXSx1elFQhKyXkjzq+U4PEVHVqwvJ+b/wAzhv8AhXX/AFFf/Jf/AOyq9o/gv+ydVhvv7Q83yt3yeTtzlSOu4+tdXRTlmGJnFxlLR+S/yIp5Ngac1OMNU7rV9PmcprHgv+1tVmvv7Q8rzdvyeTuxhQOu4elUf+Fdf9RX/wAl/wD7Ku5oojmGJhFRjLReS/yCpk2BqTc5Q1bu9X1+Zw3/AArr/qK/+S//ANlW74d8O/2B9p/0rz/P2/8ALPbt259z61uUUquOr1YuE5XT8kXQynB4eoqtKFpLzf8AmcN/wrr/AKiv/kv/APZUf8K6/wCor/5L/wD2VdzRV/2liv5vwX+Rl/YWX/8APv8AGX+Zw3/Cuv8AqK/+S/8A9lXV6Pp39k6VDY+b5vlbvn27c5Ynpk+tXqKyrYutWjy1JXXyOjDZbhcLNzoxs7W3b/NhRRRXMdwUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFAHGy+Lr6HWfE9g0Fv/oIjTTeGzNI0cRKvzz888Q4xw1b/hzUZtX8MaTqdwqLPeWUNxIsYIUM6BiACScZPqaqzeE9Pn1f+0nkuPON4LwqGG1mEcaBSMfd/cxPjOdyg5xxWlpWnQ6Ro9lpluztDZ28dvG0hBYqihQTgAZwPSgC3RRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQEUUUAf/2Q==",
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAL4AAADmCAIAAADky6hQAAAL7ElEQVR4Ae2dW0hVSxjHz1bTyjApgorspQcjUqOiCykFlUaWUoR20xISe0hQ7EEICqJSii4QPRQE3SjoYhn6oEXtQoooqAe1G4IoIWFBmRaB1vlK2GflnjV7ZjV7nbVm/vvhsPbMN7Nmft/P2WvHmdmBnz9//oMXCMgTiJFvghYg8IsA1IEHDglAHYfg0AzqwAGHBKCOQ3BoBnXggEMCUMchODSDOnDAIQGo4xAcmkEdOOCQANRxCA7NoA4ccEgA6jgEh2ZQBw44JAB1HIJDM6gj4UBTU9P69euHG7S1tS1YsGBoaEiivX1oe3t7Wlrajx8/hkPmz5/f2tpqHy5XU1VVtez3KysrKz4+Xq4xJ5r+Vy+8xAnk5uY+ePCA4leuXPn48WPxhhEjN27ceOXKFQpraGgoKCiIGO8g4MSJE9XV1Q4aMpv8wyxFoR2BV69eLVq06Pr16zt27LCLcVZOC8/s2bNpGVuyZAktac464bR69+5denr6wMAAJ0aqKkDRnDUJVeEEKioq6urqnj9/PnHixPDavynZsmXL2LFjKbuXL1/+m36YbQsLC4uLi2nVZNY6KIxz0MbwJn19fXFxcf39/crV2bt376xZs2j5UU6YntJoPVPoza8RSq1RCG5paaEn5ebm5ry8vGjQGD9+vPJuv337NmfOnO7ubrU94xuWxF/44OBgZWXlkSNH6BmZFp76+nqJxv9f6MGDB4uKiqZNm6Z2CHjWkeB57NixDx8+HDp0iNp0dXXl5OQ8e/YsMTFRootIocnJyZ8+fYoUJVefkJBA/44QGxs73Iy+wY0bN06uC1Y01GFRQZkAAXxgCUBCCIsA1GFRQZkAAagjAAkhLAJQh0UFZQIEoI4AJISwCEAdFhWUCRCAOgKQEMIiAHVYVFAmQADqCEBCCIsA1GFRQZkAAagjAAkhLAJQh0UFZQIEoI4AJISwCEAdFhWUCRCAOgKQEMIiAHVYVFAmQADqCEBCCIsA1GFRQZkAAagjAAkhLAJQh0UFZQIEoI4AJISwCOipTm1tLWuyv8poq8r27dtpp5xdAMoFCRinztq1a+fNmxcIBAQBIcyOgIbq7Nu3jzaEZ2dnM+d87dq18vJyZhUKpQjouYVveA/l69evy8rKQjiCwWDoOhqbLEOdG3KhszqcFDpTJxDkdCld9XPZf02uXr3635u/vqKTnUJ9RK9nnQ9J4aw6IbK4cExAT3XoUD56paamWj+kHDNCQyYBDR+TaZ503iKdf8OcMApVEdBz1WlsbOQDUn4QCf92WtbquepomSqvTQrqeC0jvhkP1PFNqrw2UKjjtYz4ZjxQxzep8tpAoY7XMuKb8UAd36TKawOFOl7LiG/GA3V8kyqvDRTqeC0jvhkP1PFNqrw2UKjjtYz4ZjxQxzep8tpAoY7XMuKb8UAd36TKawOFOl7LiG/GA3V8kyqvDRTqeC0jvhkP1PFNqrw2UKjjtYz4ZjxQxzep8txA1f6AsUd6q6mpsRvJ2bNnMzMzMzIy6Ke/7WJQLkLArI3Dvb29GzZsuH///ps3b9atW/fy5UvP/Sn7Z0AafmBxTrr4+PHjrl27YmJiUlJS6No/afLiSHVedTh7zs+fP//w4UP68PJiTnwyJp3VsUtBR0dHfn4+fWxNmjTJLoZZHggyix0WWk+6iF7POOnCSXqYqw6d2lRYWEjrjaw3TkagdRs995zbnXRBXxy2bdtWVVW1cOFCrdPqxuQ0fEwmbHYnXZw7d46+k58+fXrZsmVr1qxxA7C+99DzWSdK+YreE0n0eo7es46eq06U1EG3VgJQx0oD1xIEoI4ELIRaCUAdKw1cSxCAOhKwEGolAHWsNHAtQQDqSMBCqJUA1LHSwLUEAagjAQuhVgJQx0oD1xIEoI4ELIRaCUAdKw1cSxCAOhKwEGolAHWsNHAtQQDqSMBCqJUA1LHSwLUEAagjAQuhVgJQx0oD1xIEoI4ELIRaCUAdKw1cSxCAOhKwEGoloOeOiNra2urqaus8h6+/fv1aUlLy/v37L1++7N+/Pzc3NzwGJYIE9FQnOTmZ+cOwhw8fprMKdu/e3dPTs3jx4s7OTkFMCAsnoOHuz9BJF83NzSMmXFpampiYSIVtbW2jRo0aUYu3UgR0XnWYe86JTlFRUV1d3e3bt5cvXy4FC8FWAjqrY53niOv6+vobN25cuHBhRDn/bSDIr5erdeeki+iNWedvWLTq0N7y0IsSW15ePjg4SBe04byxsVEu1Yj+k4CGzzo0QbuTLj5//nzr1i06E+7Ro0epqal/osA7OQJ6fmDRt+5AINDQ0DACRldXFx2SQmLFx8efPHly5syZIwL4b6O3+PuyZzpyhs8LtSECvkxwMDR8BRfW5zOdn3UUoEIX9gSgjj0b1HAJQB0uHlTaE4A69mxQwyUAdbh4UGlPAOrYs0ENlwDU4eJBpT0BqGPPBjVcAlCHiweV9gSgjj0b1HAJQB0uHlTaE4A69mxQwyUAdbh4UGlPAOrYs0ENlwDU4eJBpT0BqGPPBjVcAlCHiweV9gSgjj0b1HAJQB0uHlTaE4A69mxQwyUAdbh4UMkhQDsi9HvV1NTwJ0X790aPHs2PQS2fgJ77sOxOuhj+E6ITUnJyclpbW/v6+jh/VKjiE9DwAyt00oXdzOnonYqKCjotxS4A5SIEdF51mCddtLS0HD169ObNm/yVSYSd4TE6qxOe2u/fvy9dupROSJk6daoDdQLB8C6dl1h3UvqxZz2PKxjOZ/iqQwej0IPO5s2bKaC/v7+4uFj2nBTnpmjXUs9VJykpiQ6E4z/NYNVxILN1pdTzUTErKysvL88BGjQRJ6DnB5bIsUvMcyrFwSFSz1UHeXWBANRxAbKet4A6eubVhVlBHRcg63kLqKNnXl2YFdRxAbKet4A6eubVhVlBHRcg63kLqKNnXl2YFdRxAbKet4A6eubVhVlBHRcg63kLqKNnXl2YFdRxAbKet4A6eubVhVlBHRcg63kLqKNnXl2YFdRxAbKet4A6eubVhVlBHRcg63kLqKNnXl2YlZ77sGpra2ljeTi+O3fulJWVTZ8+naoyMzMPHDgQHoMSQQJ6qmO3Pe/ixYsDAwM7d+4UpIMwDgENP7A4J1309PRMmTKFgwNV4gR0XnXC95xXVlbSqtPe3j5hwoTjx4/PmDFDnBQiRxDQc/fn8CRTU1ODwaB1woFAICMj48yZM3TYRWlp6b1796y1Ea8Df3QWMTxCgHX/ti97pkO/IkzRh9XDzzrhq05nZ2dKSkpsbOzQ0NDkyZN7e3ulJufLBAelphgh2Kq7nqvOj9+v8FWHvnZt2rQpPz//yZMnaWlpETihmktAz2ed3Nxc+mxqaGgYMfe3b9+WlJTExcXRGZSnTp2SfdbBqmNddfRUZ4Qxqt5CHas6Maqwoh/TCEAd0zKubL5QRxlK0zqCOqZlXNl8oY4ylKZ1BHVMy7iy+UIdZShN6wjqmJZxZfOFOspQmtYR1DEt48rmC3WUoTStI6hjWsaVzRfqKENpWkdQx7SMK5sv1FGG0rSOoI5pGVc2X6ijDKVpHUEd0zKubL5QRxlK0zqCOqZlXNl8oY4ylKZ1hB0RpmVc2Xyx6ihDaVpHUMe0jCubL9RRhtK0jqCOaRlXNl+oowylaR1BHdMyrmy+UEcZStM6gjqmZVzZfKGOMpSmdQR1TMu4svlCHWUoTesI6piWcWXzhTrKUJrWEdQxLePK5gt1lKF0s6MxY8YUFBSE7rh161YqCb115wLquMNZ8V0SEhLoQHE6OJz6pUPTOzo6qETxPSJ1B3UiEfJq/dy5c58+fUqje/HiRXp6uvvDhDruM1dzx5ycnKamJuqL/kvXajqV6QXqyNDyUmx2dvbdu3dpRPQjKStWrHB/aFDHfeZq7ki/yxQTE9Pd3U3dJSUlqelUpheoI0PLY7GrVq3as2fP/7LkEAmo4zEdZIZDv6Jy6dKl1atXyzRSFovNNMpQmtYRVh3TMq5svlBHGUrTOoI6pmVc2XyhjjKUpnUEdUzLuLL5Qh1lKE3rCOqYlnFl84U6ylCa1hHUMS3jyuYLdZShNK2jfwHK16lQ6dKLiwAAAABJRU5ErkJggg==",
      "text/plain": [
       "<PIL.Image.Image image mode=RGB size=190x230>"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dinkum.reset()\n",
    "\n",
    "custom_lambda_and = lambda X, Y: X and Y\n",
    "\n",
    "x = Gene(name='X')\n",
    "y = Gene(name='Y')\n",
    "z = Gene(name='Z')\n",
    "m = Tissue(name='M')\n",
    "\n",
    "x.is_present(where=m, start=1)\n",
    "y.is_present(where=m, start=2)\n",
    "\n",
    "z.custom_activation(state_fn=custom_lambda_and, delay=1)\n",
    "\n",
    "dinkum.run_and_display(stop=5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "5fb402c6-849a-4a29-932a-fb5d3dd112b4",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.11.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}