009_2DFiniteElements_Asymptotic.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "from PIL import Image\n",
    "\n",
    "size = 10\n",
    "img = Image.new(\"RGB\", (size,size), \"white\") # create a new 15x15 image\n",
    "pixels = img.load() # create the pixel map\n",
    "\n",
    "black_2 = []\n",
    "for i in range(img.size[0]):\n",
    "    if i % 2 == 0:\n",
    "        black_2.append(i)\n",
    "\n",
    "black_1 = [i-1 for i in black_2 if i > 0]\n",
    "if img.size[0] % 2 == 0: # 'that' if statement\n",
    "    black_1.append(img.size[0]-1)\n",
    "\n",
    "\n",
    "for i in black_1:\n",
    "    for j in range(0, size, 2):\n",
    "        pixels[i,j] = (0,0,0)\n",
    "\n",
    "for k in black_2:\n",
    "    for l in range(1, size+1, 2):\n",
    "        pixels[k,l] = (0,0,0)\n",
    "\n",
    "img.save(\"trialSquare.png\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "img = plt.imread(\"trialSquare.png\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x21732217c88>"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAc4AAAHSCAYAAABl8itQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAPSElEQVR4nO3dwYtsB5XH8d+Z14omjiikNyZhEkGcCcIQuxA14MK40FF0M4sICuMmm1GjCBJn4z8gogsRQtSNQRcxCxFRB9TFbILdiaDxKYToJE8j9ixGxU0Mnlm8HsjE93x1tK5V1fP5QOB11c3lcHj0993b1VXV3QEA1vM32x4AAPaJcALAgHACwIBwAsCAcALAgHACwMDBEie94YYb+pZbblni1Bt3cnKy7RFGjo6Otj3C2ux2OXa7HLtdzr7ttrvrSo/XEr/HuVqt+vj4eOPnXULVFfeys/bp927tdjl2uxy7Xc4e7vaKA7tVCwADwgkAA8IJAAPCCQADwgkAA8IJAAPCCQADwgkAA8IJAAPCCQADwgkAA8IJAAPCCQADwgkAA8IJAAPCCQADa4Wzqt5aVT+pqser6t6lhwKAXXXNcFbVhSSfSfK2JLcleXdV3bb0YACwi9a54nxdkse7+4nufibJl5O8a9mxAGA3rRPOG5M89ZyvL509BgD/76wTzrrCY/1HB1XdXVXHVXV8enr6l08GADtonXBeSnLzc76+Kckvnn9Qd9/X3avuXh0eHm5qPgDYKeuE83tJXlVVt1bVC5PcleSry44FALvp4FoHdPezVfX+JN9MciHJ57v7scUnA4AddM1wJkl3fz3J1xeeBQB2nncOAoAB4QSAAeEEgAHhBIAB4QSAAeEEgAHhBIAB4QSAAeEEgAHhBIAB4QSAAeEEgAHhBIAB4QSAAeEEgAHhBICBtT7Ieurk5CRVtcSpN667tz3CyL7sNbHbJdntcux2Ofu029VqddXnXHECwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwMDBEic9OjrK8fHxEqfeuKra9ggj3b3tEdZmt8ux2+XY7XL2bbdX44oTAAaEEwAGhBMABoQTAAaEEwAGhBMABoQTAAaEEwAGhBMABoQTAAaEEwAGhBMABoQTAAaEEwAGhBMABoQTAAaEEwAGrhnOqrq5qr5TVRer6rGquuevMRgA7KKDNY55NslHuvuRqvrbJCdV9e/d/aOFZwOAnXPNK87ufrq7Hzn782+TXExy49KDAcAuGv2Ms6puSXJ7koeXGAYAdt3a4ayqlyT5SpIPdfdvrvD83VV1XFXHp6enm5wRAHbGWuGsqhfkcjQf6O6HrnRMd9/X3avuXh0eHm5yRgDYGeu8qraSfC7Jxe7+5PIjAcDuWueK844k703y5qr6/tl//7TwXACwk6756yjd/R9J6q8wCwDsPO8cBAADwgkAA8IJAAPCCQADwgkAA8IJAAPCCQADwgkAA8IJAAPCCQADwgkAA8IJAAPCCQADwgkAA8IJAAPCCQAD1/wg6z/HyclJqvbjs6+7e9sjjOzLXhO7XZLdLsdul7NPu12tVld9zhUnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwcLHHSo6OjHB8fL3HqjauqbY8w0t3bHmFtdrscu12O3S5n33Z7Na44AWBAOAFgQDgBYEA4AWBAOAFgQDgBYEA4AWBAOAFgQDgBYEA4AWBAOAFgQDgBYEA4AWBAOAFgQDgBYEA4AWBAOAFgYO1wVtWFqnq0qr625EAAsMsmV5z3JLm41CAAsA/WCmdV3ZTk7UnuX3YcANht615xfirJR5P8YcFZAGDnXTOcVfWOJL/q7pNrHHd3VR1X1fHp6enGBgSAXbLOFecdSd5ZVT9L8uUkb66qLz7/oO6+r7tX3b06PDzc8JgAsBuuGc7u/lh339TdtyS5K8m3u/s9i08GADvI73ECwMDB5ODu/m6S7y4yCQDsAVecADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAw+iDrdZ2cnKSqljj1xnX3tkcY2Ze9Jna7JLtdjt0uZ592u1qtrvqcK04AGBBOABgQTgAYEE4AGBBOABgQTgAYEE4AGBBOABgQTgAYEE4AGBBOABgQTgAYEE4AGBBOABgQTgAYEE4AGBBOABgQTgAYEE4AGBBOABgQTgAYEE4AGBBOABgQTgAYEE4AGBBOABgQTgAYEE4AGBBOABgQTgAYEE4AGBBOABgQTgAYEE4AGBBOABgQTgAYEE4AGBBOABgQTgAYEE4AGBBOABgQTgAYEE4AGDhY4qRHR0c5Pj5e4tQbV1XbHmGku7c9wtrsdjl2uxy7Xc6+7fZqXHECwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsDAWuGsqpdV1YNV9eOqulhVb1h6MADYRQdrHvfpJN/o7n+uqhcmuW7BmQBgZ10znFX10iRvSvIvSdLdzyR5ZtmxAGA3rXOr9pVJTpN8oaoerar7q+r6hecCgJ20TjgPkrw2yWe7+/Ykv0ty7/MPqqq7q+q4qo5PT083PCYA7IZ1wnkpyaXufvjs6wdzOaT/R3ff192r7l4dHh5uckYA2BnXDGd3/zLJU1X16rOH7kzyo0WnAoAdte6raj+Q5IGzV9Q+keR9y40EALtrrXB29/eTrBaeBQB2nncOAoAB4QSAAeEEgAHhBIAB4QSAAeEEgAHhBIAB4QSAAeEEgAHhBIAB4QSAAeEEgAHhBIAB4QSAAeEEgAHhBICBtT7Ieurk5CRVtcSpN667tz3CyL7sNbHbJdntcux2Ofu029VqddXnXHECwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwIBwAsCAcALAgHACwMDBEic9OjrK8fHxEqfeuKra9ggj3b3tEdZmt8ux2+XY7XL2bbdX44oTAAaEEwAGhBMABoQTAAaEEwAGhBMABoQTAAaEEwAGhBMABoQTAAaEEwAGhBMABoQTAAaEEwAGhBMABoQTAAaEEwAG1gpnVX24qh6rqh9W1Zeq6kVLDwYAu+ia4ayqG5N8MMmqu1+T5EKSu5YeDAB20bq3ag+SvLiqDpJcl+QXy40EALvrmuHs7p8n+USSJ5M8neTX3f2tpQcDgF20zq3alyd5V5Jbk7wiyfVV9Z4rHHd3VR1X1fHp6enmJwWAHbDOrdq3JPlpd5929++TPJTkjc8/qLvv6+5Vd68ODw83PScA7IR1wvlkktdX1XVVVUnuTHJx2bEAYDet8zPOh5M8mOSRJD84+3/uW3guANhJB+sc1N0fT/LxhWcBgJ3nnYMAYEA4AWBAOAFgQDgBYEA4AWBAOAFgQDgBYEA4AWBAOAFgQDgBYEA4AWBAOAFgQDgBYEA4AWBAOAFgQDgBYGCtD7KeOjk5SVUtceqN6+5tjzCyL3tN7HZJdrscu13OPu12tVpd9TlXnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwIJwAMCCcADAgnAAwcLDQef8ryX9u+Jw3nJ13o6pq06fcR3a7HLtdjt0ux26Tv7vaE9Xdf81B/mxVddzdq23PcR7Z7XLsdjl2uxy7/dPcqgWAAeEEgIF9Cud92x7gHLPb5djtcux2OXb7J+zNzzgBYBfs0xUnAGzdXoSzqt5aVT+pqser6t5tz3NeVNXNVfWdqrpYVY9V1T3bnum8qaoLVfVoVX1t27OcJ1X1sqp6sKp+fPb39w3bnum8qKoPn30/+GFVfamqXrTtmXbNzoezqi4k+UyStyW5Lcm7q+q27U51bjyb5CPd/Q9JXp/kX+124+5JcnHbQ5xDn07yje7++yT/GDveiKq6MckHk6y6+zVJLiS5a7tT7Z6dD2eS1yV5vLuf6O5nknw5ybu2PNO50N1Pd/cjZ3/+bS5/87lxu1OdH1V1U5K3J7l/27OcJ1X10iRvSvK5JOnuZ7r7v7c71blykOTFVXWQ5Lokv9jyPDtnH8J5Y5KnnvP1pfjmvnFVdUuS25M8vN1JzpVPJflokj9se5Bz5pVJTpN84ew2+P1Vdf22hzoPuvvnST6R5MkkTyf5dXd/a7tT7Z59COeV3qPJS4E3qKpekuQrST7U3b/Z9jznQVW9I8mvuvtk27OcQwdJXpvks919e5LfJfHahw2oqpfn8h29W5O8Isn1VfWe7U61e/YhnJeS3Pycr2+KWwcbU1UvyOVoPtDdD217nnPkjiTvrKqf5fKPF95cVV/c7kjnxqUkl7r7f++OPJjLIeUv95YkP+3u0+7+fZKHkrxxyzPtnH0I5/eSvKqqbq2qF+byD6q/uuWZzoW6/I7Ln0tysbs/ue15zpPu/lh339Tdt+Ty39lvd7d/uW9Ad/8yyVNV9eqzh+5M8qMtjnSePJnk9VV13dn3hzvjhVd/ZKlPR9mY7n62qt6f5Ju5/Aqvz3f3Y1se67y4I8l7k/ygqr5/9ti/dffXtzgTrOMDSR44+8f0E0net+V5zoXufriqHkzySC6/6v7ReBehP+KdgwBgYB9u1QLAzhBOABgQTgAYEE4AGBBOABgQTgAYEE4AGBBOABj4H2yBtky4OKujAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 576x576 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(8,8))\n",
    "plt.imshow(img)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}