Tupper's self-referencing function

Tupper.ipynb

{
  "cells": [
    {
      "metadata": {
        "trusted": true
      },
      "cell_type": "code",
      "source": "%matplotlib inline",
      "execution_count": 1,
      "outputs": []
    },
    {
      "metadata": {
        "trusted": true
      },
      "cell_type": "code",
      "source": "import matplotlib.pyplot as plt",
      "execution_count": 2,
      "outputs": []
    },
    {
      "metadata": {
        "trusted": true
      },
      "cell_type": "code",
      "source": "N = 4858450636189713423582095962494202044581400587983244549483093085061934704708809928450644769865524364849997247024915119110411605739177407856919754326571855442057210445735883681829823754139634338225199452191651284348332905131193199953502413758765239264874613394906870130562295813219481113685339535565290850023875092856892694555974281546386510730049106723058933586052544096664351265349363643957125565695936815184334857605266940161251266951421550539554519153785457525756590740540157929001765967965480064427829131488548259914721248506352686630476300\nH = 17\nW = 106",
      "execution_count": 3,
      "outputs": []
    },
    {
      "metadata": {
        "trusted": true
      },
      "cell_type": "code",
      "source": "def tupper(x,y):\n    return 0.5 < ((y//H) // (2**(H*x + y%H))) % 2",
      "execution_count": 4,
      "outputs": []
    },
    {
      "metadata": {
        "trusted": true
      },
      "cell_type": "code",
      "source": "plt.rc('patch', antialiased=True)\nfor x in range(W):\n    for yy in range(H):\n        y = N + yy\n        if tupper(x, y):\n            plt.bar(left=x, bottom=yy, height=1, width=1, linewidth=0, color='black')\n\nplt.axis('scaled')\nbuf = 10\nplt.xlim((-buf,W+buf))\nplt.ylim((-buf,H+buf))\nplt.xticks(range(0, W, 100))\nplt.yticks([0, H], ['k', f'k + {H}']);",
      "execution_count": 5,
      "outputs": [
        {
          "data": {
            "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAACFCAYAAABBqm9MAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAB+5JREFUeJzt3U+IJPUZxvHn0cVzYmYSgsmkIyyC5LCE1iQIkqMJwm4i\ni4oYswTmlEtu3uLRsxCEPUgMRCWQwwoJLLIXTwFnUXQXkSyi2cEQY0IWXA9L8M1ha9xm8k5XdVV1\n/evvB5qZrq6u39v/6unfr6qrHBECAOCw2/ouAAAwTAQEACBFQAAAUgQEACBFQAAAUgQEACBFQAAA\nUgQEACBFQAAAUsf6LqCJra2tmM1mfZcBAKNy8eLFTyJiu2y+UQfEbDbT3t5e32UAwKjY/rDKfAwx\nAQBSBAQAIEVAAABSBAQAIEVAAABSBAQAIEVAAABSBAQAIEVAAABSBAQAIEVAAABSBAQAIEVAAABS\nBAQAIEVAAABSBAQAIEVAAABSBAQAIEVAAABSBAQAIEVAAABSSwPC9sz2pTYbtH3a9mXbn9ueL0x/\nwvZbC5fPbZ9os20AQHWt9iBs/9z2MyWzXZL0U0mvL06MiN9HxImIOCHpSUkfRMRbbdYHAKiuckDY\nvtv2m7bva9JgRLwbEe+VzPa4pJebtAMAaOZYlZls3yPpFUlnOvpW/6ikkx20AwA4QpWA2JZ0TtIj\nEXH58I22vyLpQnH1Tkl32D5VXH8yIt5ZpSDb35P0WUSk2z5s70ralaSdnZ1VFj0otiVJEdFzJe05\neEzSeB5XVvPitEXZ7XUe5xifJ6n8PZvd3vR9Xuf+fbymZe2XzTvU90GVgLgm6aqkByT9X0BExL8k\nnZBuboOQNIuIZxrU9JiWDC9FxFlJZyVpPp8P81kFgAmoEhA3JJ2SdN72pxHx0rqKsX2bpNOSHlxX\nGwCAaiptpI6I65IelvQr2422Ddj+ie19ST+Q9Cfb5xduflDSfkS836QNDIftLy5DWtYyEfHFBdPA\na1qPx/yEzefz2Nvb67uMWoY+9lhH2dhvl2O7bS1zHY+JbRDttbnsPkfdb0jv077WA7YvRsS8bL5K\nezFtuuwb65g+2Os09JVdl6/dshVklXbLaq3acypbgZfdb5UeWp3Xf90r6DrqvmZVlzXEz0YVHGoD\nAJAiIAAAKYaYKqjT9UY1R3Xn2+qad7mNoOlwSdV5y4ZDytpfZdiqzhDXEPU9hDVW9CAAACkCAgCQ\nYogJaEnTvZi6UncvpbZMaS+fqaMHAQBI0YPAUl1uRG7zB2hVl1O2EXcVbS6ryfKParPORvB1WMeP\n0haXs+7XIWv/qLaG8juPuuhBAABSBAQAIMUQUwVZN3IMx1sZg3Uf12isXfuxqDPsOMTXoe4hUab+\nnqMHAQBIERAAgBRDTBV0dbiGVY7m2fTIoFXvs27rOpRGV8N6bR2BdepWOeXnKsvq+8ivB/quaV3o\nQQAAUvQgKqhzHP22vhmv8juEdZwoZV3fgJue26Bsmev8FreOA/BNQZPfH/T9nHZ5QqMxoQcBAEgR\nEACAFENMFbR1nH7c0tU5pYG2bdJ7jh4EACBFQAAAUgwxVdD0KKFd7elQtv9909NPtqlsL6wDdfdi\nykx9OADdm/p7jh4EACBFD6KCpr+cbLIveN1j+9f55tL3t506vZm+a8a0VT3fxFTRgwAApAgIAECK\nIaYVbNL+z2XGfipFoAtjPxQHPQgAQIqAAACkGGKqYOzdxLFguApDwPvwFnoQAIAUPYgVlJ0VC8D4\n8Xm+hR4EACBFQAAAUgwxVUCX82h1z4HR5PSUwJBNaSP3oHoQtme2L/VdBwBgYAEBABiOwQaE7btt\nv2n7vr5rWSfb6WWKVnmcm/KcYDOM9XM+yICwfY+kP0o6ExFv9F0PAGyiIW6k3pZ0TtIjEXH58I22\ndyXtStLOzk7HpbVn7BuvMm1ueJ7i84PNMKUdMIbYg7gm6aqkB7IbI+JsRMwjYr69vd1tZQCwQYbY\ng7gh6ZSk87Y/jYiX+i4IADbREANCEXHd9sOSXrN9PSLO9V0TAGyaQQVERHwg6TvF//+RNOk9mABg\nyIa4DQIAMAAEBAAgRUAAAFIEBAAgRUAAAFIEBAAgRUAAAFIEBAAgRUAAAFIEBAAgRUAAAFIEBAAg\nRUAAAFIEBAAgRUAAAFIEBAAgRUAAAFIEBAAgRUAAAFIEBAAgRUAAAFIEBAAg5Yjou4babP9T0ocd\nNbcl6ZOO2gIwTmNZT3wrIrbLZhp1QHTJ9l5EzPuuA8BwTW09wRATACBFQAAAUgREdWf7LgDA4E1q\nPcE2CABAih4EACBFQJSw/ZDt92xfsf103/UA6IftF2x/bPvSwrQ7bb9m+6/F3y8X0237uWK98bbt\n7/ZXeX0ExBK2b5f0G0k/knSvpMdt39tvVQB68ltJDx2a9rSkCxFxXNKF4rp0c51xvLjsSnq+oxpb\nRUAsd7+kKxHxfkTckPSKpJM91wSgBxHxuqR/H5p8UtKLxf8vSjq1MP13cdNfJH3J9te7qbQ9BMRy\nd0m6unB9v5gGAJL0tYj4uyQVf79aTJ/EuoOAWM7JNHb7AlBmEusOAmK5fUnfXLj+DUkf9VQLgOH5\nx8HQUfH342L6JNYdBMRyb0g6bvvbtu+Q9JikV3uuCcBwvCrpqeL/pySdW5j+s2Jvpu9LunYwFDUm\nx/ouYMgi4r+2fynpvKTbJb0QEZd7LgtAD2y/LOmHkrZs70v6taRnJf3B9i8k/U3S6WL2P0v6saQr\nkj6TdKbzglvAL6kBACmGmAAAKQICAJAiIAAAKQICAJAiIAAAKQICAJAiIAAAKQICAJD6H6+/ShMb\nlPNfAAAAAElFTkSuQmCC\n",
            "text/plain": "<matplotlib.figure.Figure at 0x7fef6a578b70>"
          },
          "metadata": {},
          "output_type": "display_data"
        }
      ]
    }
  ],
  "metadata": {
    "kernelspec": {
      "name": "conda-env-phd-py",
      "display_name": "Python [conda env:phd]",
      "language": "python"
    },
    "language_info": {
      "name": "python",
      "version": "3.6.2",
      "mimetype": "text/x-python",
      "codemirror_mode": {
        "name": "ipython",
        "version": 3
      },
      "pygments_lexer": "ipython3",
      "nbconvert_exporter": "python",
      "file_extension": ".py"
    },
    "gist": {
      "id": "",
      "data": {
        "description": "Tupper's self-referencing function",
        "public": true
      }
    }
  },
  "nbformat": 4,
  "nbformat_minor": 2
}