0-Intro.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "0af4e1b0-f590-44cf-8798-f9bcdb3243f5",
   "metadata": {},
   "source": [
    "This is our first lecture for the introduction of finite element mesh in 1d and 2d.\n",
    "\n",
    "- [x] mesh1d\n",
    "- [x] mesh2d\n",
    "- [x] order\n",
    "\n",
    "Author: Yang Bai\n",
    "Date  : 2021.12.19"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "177b7b73-43ac-4e29-b7c3-3b0fb3cce025",
   "metadata": {},
   "outputs": [],
   "source": [
    "from FEToy.mesh.lagrange1dmesh import mesh1d # before we create the mesh, you must import mesh1d"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "263c7b64-034c-48c5-a4be-ebbc35f1d736",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "*** node coordinates of the mesh (total nodes=6, nodes per element=2)\n",
      "     1-th node: x=  0.000000e+00\n",
      "     2-th node: x=  2.000000e-01\n",
      "     3-th node: x=  4.000000e-01\n",
      "     4-th node: x=  6.000000e-01\n",
      "     5-th node: x=  8.000000e-01\n",
      "     6-th node: x=  1.000000e+00\n",
      "*** element connectivity information(bulk elements=5)\n",
      "     1-th element    1     2 \n",
      "     2-th element    2     3 \n",
      "     3-th element    3     4 \n",
      "     4-th element    4     5 \n",
      "     5-th element    5     6 \n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAUQ0lEQVR4nO3df5DcdX3H8ec7OWOARElMgjEBLtq0JgGB5AypUvyBwSQyYkfaAVEQZDKkjYPTOhJxpjBjZ4oyVKBFaQYZwTplHKs1raEUQRtbQTlCCMQUiRHNwQkRq8FQmh+8+8cudln3cpvcZ39c7vmY2bn9fj+f73ff79mbfe33+93bi8xEkqSRGtfpAiRJhwcDRZJUhIEiSSrCQJEkFWGgSJKK6Ol0AZ00bdq07O3t7XQZkjSqPPDAAz/PzOn168d0oPT29tLf39/pMiRpVImInzRa7ykvSVIRBookqQgDRZJUhIEiSSrCQJEkFWGgSJKKMFAkSUUYKJKkIgwUSVIRBookqQgDRZJUhIEiSSrCQJEkFWGgSJKKMFAkSUUYKJKkIgwUSVIRBookqQgDRZJUhIEiSSrCQJEkFWGgSJKKMFAkSUUYKJKkIgwUSVIRXRUoEbEsIh6NiG0RsabBeETEDdXxzRGxsG58fEQ8GBH/0r6qJUnQRYESEeOBG4HlwHzgvIiYXzdtOTC3elsJfK5u/DJga4tLlSQ10DWBAiwGtmXm9szcA9wOnF0352zgtqy4Dzg6ImYCRMRs4F3Aze0sWpJU0U2BMgvYUbM8UF3X7JzrgI8BLxzoQSJiZUT0R0T/zp07R1SwJOn/dVOgRIN12cyciDgLeDozHxjuQTJzbWb2ZWbf9OnTD6VOSVID3RQoA8CxNcuzgSebnPNm4N0R8TiVU2Vvj4i/b12pkqR63RQo9wNzI2JOREwAzgXW1c1ZB1xQ/bTXEuBXmTmYmR/PzNmZ2Vvd7p7MfH9bq5ekMa6n0wW8KDP3RcRq4E5gPHBLZm6JiEur4zcB64EVwDbgOeCiTtUrSXqpyKy/TDF29PX1ZX9/f6fLkKRRJSIeyMy++vXddMpLkjSKGSiSpCIMFElSEQaKJKkIA0WSVISBIkkqwkCRJBVhoEiSijBQJElFGCiSpCIMFElSEQaKJKkIA0WSVISBIkkqwkCRJBVhoEiSijBQJElFGCiSpCIMFElSEQaKJKkIA0WSVISBIkkqwkCRJBVhoEiSijBQJElFGCiSpCIMFElSEQaKJKkIA0WSVISBIkkqoqsCJSKWRcSjEbEtItY0GI+IuKE6vjkiFlbXHxsR34qIrRGxJSIua3/1kjS2dU2gRMR44EZgOTAfOC8i5tdNWw7Mrd5WAp+rrt8H/HlmzgOWAH/aYFtJUgt1TaAAi4Ftmbk9M/cAtwNn1805G7gtK+4Djo6ImZk5mJkbATLzWWArMKudxUvSWNdNgTIL2FGzPMBvh8KwcyKiFzgF+F75EiVJQ+mmQIkG6/Jg5kTEJOAfgY9k5q6GDxKxMiL6I6J/586dh1ysJOmluilQBoBja5ZnA082OyciXkYlTL6UmV8d6kEyc21m9mVm3/Tp04sULknqrkC5H5gbEXMiYgJwLrCubs464ILqp72WAL/KzMGICODzwNbM/Ov2li1JAujpdAEvysx9EbEauBMYD9ySmVsi4tLq+E3AemAFsA14DriouvmbgQ8AD0fEpuq6KzJzfRtbkKQxLTLrL1OMHX19fdnf39/pMiRpVImIBzKzr359N53ykiSNYgaKJKkIA0WSVISBIkkqwkCRJBVhoEiSijBQJElFGCiSpCIMFElSEQaKJKkIA0WSVISBIkkqwkCRJBVhoEiSijBQJElFGCiSpCIMFElSEQaKJKkIA0WSVISBIkkqwkCRJBVhoEiSijBQJElFGCiSpCIMFElSEQaKJKkIA0WSVISBMgIXX3wxM2bM4IQTTuh0KW2zY8cO3va2tzFv3jwWLFjA9ddf3+mSWu75559n8eLFnHTSSSxYsIArr7yy0yW1xf79+znllFM466yzOl1KW/T29nLiiSdy8skn09fX1+ly2uKXv/wl55xzDq9//euZN28e995774j211OorjHpgx/8IKtXr+aCCy7odClt09PTw7XXXsvChQt59tlnWbRoEUuXLmX+/PmdLq1lXv7yl3PPPfcwadIk9u7dy2mnncby5ctZsmRJp0trqeuvv5558+axa9euTpfSNt/61reYNm1ap8tom8suu4xly5bxla98hT179vDcc8+NaH8eoYzA6aefztSpUztdRlvNnDmThQsXAjB58mTmzZvHE0880eGqWisimDRpEgB79+5l7969RESHq2qtgYEBvvGNb3DJJZd0uhS1yK5du9iwYQMf+tCHAJgwYQJHH330iPZpoOiQPf744zz44IOceuqpnS6l5fbv38/JJ5/MjBkzWLp06WHf80c+8hE+/elPM27c2HmJiAjOPPNMFi1axNq1aztdTstt376d6dOnc9FFF3HKKadwySWXsHv37hHt84C/LRHR1mP6iFgWEY9GxLaIWNNgPCLihur45ohY2Oy2RQ0OwlveAj/7WUsfpqvU9fzrX/+a9773vVx33XW84hWv6HBxLVLT8/jx49m0aRMDAwN8//vf55FHHul0da0xOMi/nHACMyZNYtGiRZ2upj2qz/N/fu1rbNy4kTvuuIMbb7yRDRs2dLqy1hkcZN/FF7Nx40ZWrVrFgw8+yFFHHcXVV189sv1m5pA3YA/wSaDnQPNK3IDxwI+A1wITgIeA+XVzVgB3AAEsAb7X7LaNbosWLcpDsmpV5rhxmatW5Y9//ONcsGDBoe1nNKnpec+ePXnmmWfmtdde2+mqWqum51pXXXVVXnPNNR0qqsVWrco1kLOOOiqPP/74POaYY/KII47I888/v9OVtU6D5/nKK688fJ/jzMxVq3IwIo+fPPk3qzZs2JArVqxoanOgPxu8pkZlrLGIOAP4PPAL4P2Z+YORxdfQIuL3gasy853V5Y8DZOZf1cz5O+DbmfkP1eVHgbcCvcNt20hfX1/29/c3X+QRR8Dzz79k1ePAu4Dpb3lL8/sZRe78znd4+Qsv/GY5gQuBVwIPj5GedwIvA44YN453vOlNbN68meOOO45XvepVHauxtPqeX/RvEZw3ZQonnnhiB6pqrdqedwMvAJOBX0QwZ/Jkjj/++MPuGmn98/wHwM3A702cyFWXX87u3bu55pprht1PRDyQmb/1UbgDnvLKzLuBE4EHgf6I+LODrP9gzAJ21CwPVNc1M6eZbQGIiJUR0R8R/Tt37jy4Crdvh/e9jz09lQ/H/TGwaNw4tkZw7733Mjg4eHD7GwXOW7yYu2bM4PnqufR7IvgicNuRR9Lf309/fz/PPPNMZ4ssrL7nxyPo6+lhxsSJbNy4kSlTphxWYQK/3fPz48Zx14wZfHLevA5X1jq1PT8FvBl4XU8PvRMnMnXq1MMuTOC3n+e/mTiR86dO5Q1z5rBp0yauuOKKkT1Ao8OWRjfgvcA+4NfArtpbs/sYZv9/BNxcs/wB4G/q5nwDOK1m+W5gUTPbNrod0imvSy+tHB5PnNjwdMhhyZ7t+XBlz4fUM0Oc8mrqIxwR0Qf8JfAYcBnw4bpbCQPAsTXLs4Enm5zTzLZlPPUUXHop3Hdf5edYuDBvz/Z8uLLnoj0Pdw2lB7gSuBz4LLAmM58fcoORFFJ5rB8CZwBPAPcD78vMLTVz3gWspnJx/lTghsxc3My2jRz0NRRJ0pDXUIb7S/n7ganA8qxcT2mZzNwXEauBO6l8auuWzNwSEZdWx28C1lMJk23Ac8BFB9q2lfVKkl5quCOULwKrM/NX7SupfTxCkaSDd0hHKJn5gdaVJEk6nIyd71WQJLWUgSJJKsJAkSQVYaBIkoowUCRJRRgokqQiDBRJUhEGiiSpCANFklSEgSJJKsJAkSQVYaBIkoowUCRJRRgokqQiDBRJUhEGiiSpCANFklSEgSJJKsJAkSQVYaBIkoowUCRJRRgokqQiDBRJUhEGiiSpCANFklSEgSJJKsJAkSQVYaBIkoowUCRJRXRFoETE1Ii4KyIeq/6cMsS8ZRHxaERsi4g1NeuviYj/iojNEfG1iDi6bcVLkoAuCRRgDXB3Zs4F7q4uv0REjAduBJYD84HzImJ+dfgu4ITMfAPwQ+DjbalakvQb3RIoZwO3Vu/fCrynwZzFwLbM3J6Ze4Dbq9uRmf+Wmfuq8+4DZre2XElSvW4JlGMycxCg+nNGgzmzgB01ywPVdfUuBu4oXqEk6YB62vVAEfFN4NUNhj7R7C4arMu6x/gEsA/40gHqWAmsBDjuuOOafGhJ0nDaFiiZ+Y6hxiLiqYiYmZmDETETeLrBtAHg2Jrl2cCTNfu4EDgLOCMzkyFk5lpgLUBfX9+Q8yRJB6dbTnmtAy6s3r8Q+HqDOfcDcyNiTkRMAM6tbkdELAMuB96dmc+1oV5JUp1uCZSrgaUR8RiwtLpMRLwmItYDVC+6rwbuBLYCX87MLdXt/xaYDNwVEZsi4qZ2NyBJY13bTnkdSGY+A5zRYP2TwIqa5fXA+gbzfqelBUqShtUtRyiSpFHOQJEkFWGgSJKKMFAkSUUYKJKkIgwUSVIRBookqQgDRZJUhIEiSSrCQJEkFWGgSJKKMFAkSUUYKJKkIgwUSVIRBookqQgDRZJUhIEiSSrCQJEkFWGgSJKKMFAkSUUYKJKkIgwUSVIRBookqQgDRZJUhIEiSSrCQJEkFWGgSJKKMFAkSUUYKJKkIgwUSVIRBookqYiuCJSImBoRd0XEY9WfU4aYtywiHo2IbRGxpsH4RyMiI2Ja66uWJNXqikAB1gB3Z+Zc4O7q8ktExHjgRmA5MB84LyLm14wfCywFftqWiiVJL9EtgXI2cGv1/q3AexrMWQxsy8ztmbkHuL263Ys+A3wMyBbWKUkaQrcEyjGZOQhQ/TmjwZxZwI6a5YHqOiLi3cATmfnQcA8UESsjoj8i+nfu3DnyyiVJAPS064Ei4pvAqxsMfaLZXTRYlxFxZHUfZzazk8xcC6wF6Ovr82hGkgppW6Bk5juGGouIpyJiZmYORsRM4OkG0waAY2uWZwNPAq8D5gAPRcSL6zdGxOLM/FmxBiRJB9Qtp7zWARdW718IfL3BnPuBuRExJyImAOcC6zLz4cyckZm9mdlLJXgWGiaS1F7dEihXA0sj4jEqn9S6GiAiXhMR6wEycx+wGrgT2Ap8OTO3dKheSVKdtp3yOpDMfAY4o8H6J4EVNcvrgfXD7Ku3dH2SpOF1yxGKJGmUM1AkSUUYKJKkIgwUSVIRBookqQgDRZJUhIEiSSrCQJEkFWGgSJKKMFAkSUUYKJKkIgwUSVIRBookqQgDRZJUhIEiSSrCQJEkFWGgSJKKMFAkSUUYKJKkIgwUSVIRBookqQgDRZJUhIEiSSrCQJEkFRGZ2ekaOiYidgI/OcTNpwE/L1jOaGDPY4M9jw0j6fn4zJxev3JMB8pIRER/ZvZ1uo52suexwZ7Hhlb07CkvSVIRBookqQgD5dCt7XQBHWDPY4M9jw3Fe/YaiiSpCI9QJElFGCiSpCIMlGFExLKIeDQitkXEmgbjERE3VMc3R8TCTtRZUhM9n1/tdXNEfDciTupEnSUN13PNvDdGxP6IOKed9ZXWTL8R8daI2BQRWyLi39tdY2lN/F6/MiL+OSIeqvZ8USfqLCkibomIpyPikSHGy75+Zaa3IW7AeOBHwGuBCcBDwPy6OSuAO4AAlgDf63Tdbej5TcCU6v3lY6Hnmnn3AOuBczpdd4uf46OBHwDHVZdndLruNvR8BfCp6v3pwC+ACZ2ufYR9nw4sBB4ZYrzo65dHKAe2GNiWmdszcw9wO3B23Zyzgduy4j7g6IiY2e5CCxq258z8bmb+d3XxPmB2m2ssrZnnGeDDwD8CT7ezuBZopt/3AV/NzJ8CZOZY6DmByRERwCQqgbKvvWWWlZkbqPQxlKKvXwbKgc0CdtQsD1TXHeyc0eRg+/kQlXc4o9mwPUfELOAPgZvaWFerNPMc/y4wJSK+HREPRMQFbauuNZrp+W+BecCTwMPAZZn5QnvK65iir189Iy7n8BYN1tV/zrqZOaNJ0/1ExNuoBMppLa2o9Zrp+Trg8szcX3kDO6o1028PsAg4AzgCuDci7svMH7a6uBZppud3ApuAtwOvA+6KiO9k5q4W19ZJRV+/DJQDGwCOrVmeTeXdy8HOGU2a6ici3gDcDCzPzGfaVFurNNNzH3B7NUymASsiYl9m/lNbKiyr2d/rn2fmbmB3RGwATgJGa6A00/NFwNVZubiwLSJ+DLwe+H57SuyIoq9fnvI6sPuBuRExJyImAOcC6+rmrAMuqH5aYgnwq8wcbHehBQ3bc0QcB3wV+MAofsdaa9ieM3NOZvZmZi/wFeBPRmmYQHO/118H/iAieiLiSOBUYGub6yypmZ5/SuWIjIg4Bvg9YHtbq2y/oq9fHqEcQGbui4jVwJ1UPiVyS2ZuiYhLq+M3UfnEzwpgG/AclXc5o1aTPf8F8Crgs9V37PtyFH9Ta5M9Hzaa6Tczt0bEvwKbgReAmzOz4UdPR4Mmn+NPAl+IiIepnAq6PDNH9VfaR8Q/AG8FpkXEAHAl8DJozeuXX70iSSrCU16SpCIMFElSEQaKJKkIA0WSVISBIkkqwkCRJBVhoEhdICLGRcSGiKj/I9Ijq1+5/rlO1SY1y0CRukD1Swg/CLw9Ii6uGfoUlT9A/mgn6pIOhn/YKHWR6l9ufxo4EfgdKn/Z/dbM/I+OFiY1wUCRukxE3EnlG357gdsz82OdrUhqjoEidZmImEPlvwv+CDghM/+3wyVJTfEaitR9Lgb+h8pXib+2w7VITfMIReoiEfFG4LvAu4FVwDHAmzJzf0cLk5rgEYrUJSJiInAb8IXMvANYSeXCvNdQNCp4hCJ1iYj4DPAe4A2Z+Wx13bnArcCi0fz/SDQ2GChSF4iI04F7gHdk5rfrxr5M5VrKkszc14HypKYYKJKkIryGIkkqwkCRJBVhoEiSijBQJElFGCiSpCIMFElSEQaKJKkIA0WSVMT/AYh/D+RdfAXbAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "mymesh=mesh1d()\n",
    "mymesh.createmesh()\n",
    "mymesh.plotmesh(withnode=True,withnodeid=True)\n",
    "mymesh.printnodes() # print out all node's coordinates\n",
    "mymesh.printelements() # print out all elements' connectivity"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "id": "a627c757-250f-47bc-8849-5554aec512e0",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "*** node coordinates of the mesh (total nodes=21, nodes per element=3)\n",
      "     1-th node: x= -1.000000e+00\n",
      "     2-th node: x= -7.000000e-01\n",
      "     3-th node: x= -4.000000e-01\n",
      "     4-th node: x= -1.000000e-01\n",
      "     5-th node: x=  2.000000e-01\n",
      "     6-th node: x=  5.000000e-01\n",
      "     7-th node: x=  8.000000e-01\n",
      "     8-th node: x=  1.100000e+00\n",
      "     9-th node: x=  1.400000e+00\n",
      "    10-th node: x=  1.700000e+00\n",
      "    11-th node: x=  2.000000e+00\n",
      "    12-th node: x=  2.300000e+00\n",
      "    13-th node: x=  2.600000e+00\n",
      "    14-th node: x=  2.900000e+00\n",
      "    15-th node: x=  3.200000e+00\n",
      "    16-th node: x=  3.500000e+00\n",
      "    17-th node: x=  3.800000e+00\n",
      "    18-th node: x=  4.100000e+00\n",
      "    19-th node: x=  4.400000e+00\n",
      "    20-th node: x=  4.700000e+00\n",
      "    21-th node: x=  5.000000e+00\n",
      "*** element connectivity information(bulk elements=10)\n",
      "     1-th element    1     2     3 \n",
      "     2-th element    3     4     5 \n",
      "     3-th element    5     6     7 \n",
      "     4-th element    7     8     9 \n",
      "     5-th element    9    10    11 \n",
      "     6-th element   11    12    13 \n",
      "     7-th element   13    14    15 \n",
      "     8-th element   15    16    17 \n",
      "     9-th element   17    18    19 \n",
      "    10-th element   19    20    21 \n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAblklEQVR4nO3dfXRV9b3n8fc3JOFBwoNCMAZIkBsQwkMw4UGdii2Gglptr9XKWGWkLirr6tJWRBzbpe2MxfGhYCvXDq1Ue3XpwuVVWRWHWtSiHakERRQxCogQg4DypERIQr7zRw6ZkJ6EfZIfOSfyea11VrL3/n3377th53xy9j5JzN0RERFpq7RkNyAiIl8PChQREQlCgSIiIkEoUEREJAgFioiIBJGe7AaSqU+fPp6fn5/sNkREOpQ1a9Z85u59m64/oQMlPz+fsrKyZLchItKhmNnH8dbrkpeIiAShQBERkSAUKCIiEoQCRUREglCgiIhIEAoUEREJQoEiIiJBKFBERCQIBYqIiAShQBERkSAUKCIiEoQCRUREglCgiIhIEAoUEREJQoEiIiJBKFBERCQIBYqIiAShQBERkSAUKCIiEoQCRUREglCgiIhIEAoUEREJQoEiIiJBKFBERCQIBYqIiASRUoFiZlPMrNzMNprZ3Djbzcx+E9u+zszObLK9k5m9ZWZ/br+uRUQEUihQzKwTsBCYCgwHppnZ8CbDpgIFscdM4KEm228ENhznVkVEJI6UCRRgHLDR3Te7ezXwJHBJkzGXAH/yequAXmaWA2Bm/YELgT+0Z9MiIlIvlQIlF9jWaLkiti7qmAXAHKCupUnMbKaZlZlZ2a5du9rUsIiI/H+pFCgWZ51HGWNmFwE73X3NsSZx90XuXuLuJX379m1NnyIiEkcqBUoFMKDRcn+gMuKYc4CLzWwL9ZfKvmVmjx2/VkVEpKlUCpTVQIGZDTKzTOAKYGmTMUuBq2Pv9poA7HP37e5+m7v3d/f8WN1L7v7Ddu1eROQEl57sBo5w91ozux5YDnQCFrv7ejO7Lrb9d8Ay4AJgI1AFXJOsfkVE5Gjm3vQ2xYmjpKTEy8rKkt2GiEiHYmZr3L2k6fpUuuQlIiIdmAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJIiUChQzm2Jm5Wa20czmxtluZvab2PZ1ZnZmbP0AM3vZzDaY2Xozu7H9uxcRObGlTKCYWSdgITAVGA5MM7PhTYZNBQpij5nAQ7H1tcDN7j4MmAD8W5xaERE5jlImUIBxwEZ33+zu1cCTwCVNxlwC/MnrrQJ6mVmOu2939zcB3P0LYAOQ257Ni4ic6FIpUHKBbY2WK/jnUDjmGDPLB8YA/wjfooiINCeVAsXirPNExphZd+Bp4CZ33x93ErOZZlZmZmW7du1qdbMiInK0VAqUCmBAo+X+QGXUMWaWQX2YPO7u/9ncJO6+yN1L3L2kb9++QRoXEZHUCpTVQIGZDTKzTOAKYGmTMUuBq2Pv9poA7HP37WZmwMPABnf/dfu2LSIiAOnJbuAId681s+uB5UAnYLG7rzez62LbfwcsAy4ANgJVwDWx8nOAq4B3zGxtbN1/d/dl7XgIIiInNHNvepvixFFSUuJlZWXJbkNEpEMxszXuXtJ0fSpd8hIRkQ5MgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFShvMmDGD7OxsRowYkXDttm3b+OY3v8mwYcMoLCzkgQceiFx78OBBxo0bx+jRoyksLOSOO+5IeP7Dhw8zZswYLrroooRr8/PzGTlyJEVFRZSUlCRUu3fvXr7//e9zxhlnMGzYMF5//fVIdeXl5RQVFTU8evTowYIFCyLPO3/+fAoLCxkxYgTTpk3j4MGDCfX9wAMPMGLECAoLCxOaN945snv3bkpLSykoKKC0tJQ9e/ZErn3qqacoLCwkLS2NsrKyhOe+5ZZbOOOMMxg1ahTf+9732Lt3b+Tan//854waNYqioiImT55MZWVl5Noj7rvvPsyMzz77LHLtnXfeSW5ubsP//bJlyxI6ZoDf/va3DB06lMLCQubMmRO59gc/+EHDvPn5+RQVFUWuXbt2LRMmTGj4OnnjjTci17799tucddZZjBw5ku985zvs378/bm1zzyNRzrHmahM5x/6Ju5+wj+LiYm+Lv/3tb75mzRovLCxMuLaystLXrFnj7u779+/3goICX79+faTauro6/+KLL9zdvbq62seNG+evv/56QvPff//9Pm3aNL/wwgsTa9zd8/LyfNeuXQnXubtfffXV/vvf/97d3Q8dOuR79uxJeB+1tbXer18/37JlS6TxFRUVnp+f71VVVe7uftlll/kf//jHyPO98847XlhY6AcOHPCamhqfNGmSf/DBB5Fq450jt9xyi8+bN8/d3efNm+dz5syJXPvee+/5+++/7xMnTvTVq1cnPPfy5cu9pqbG3d3nzJmT0Nz79u1r+PyBBx7wH//4x5Fr3d23bt3qkydP9oEDBzZ7/sSrveOOO/zee+9t8Vhbqn/ppZd80qRJfvDgQXd337FjR0J9H/HTn/7Uf/GLX0SuLS0t9WXLlrm7+/PPP+8TJ06MXFtSUuKvvPKKu7s//PDD/rOf/SxubXPPI1HOseZqo5xjQJnHeU7VK5Q2OPfcczn55JNbVZuTk8OZZ54JQFZWFsOGDeOTTz6JVGtmdO/eHYCamhpqamows8hzV1RU8Pzzz3Pttdcm3ngb7N+/n5UrV/KjH/0IgMzMTHr16pXwflasWMHgwYPJy8uLXFNbW8tXX31FbW0tVVVVnHbaaZFrN2zYwIQJE+jWrRvp6elMnDiRZ555JlJtvHPkueeeY/r06QBMnz6dZ599NnLtsGHDGDp0aKvnnjx5Munp6QBMmDCBioqKyLU9evRo+PzAgQPNnnPNfV385Cc/4Z577mnxXG3L11Rz9Q899BBz586lc+fOAGRnZyc8t7uzZMkSpk2bFrnWzBpeWezbt6/Zcy5ebXl5Oeeeey4ApaWlPP3003Frm3seiXKONVebyDnWlAIlBWzZsoW33nqL8ePHR645fPgwRUVFZGdnU1pamlDtTTfdxD333ENaWuv++82MyZMnU1xczKJFiyLXbd68mb59+3LNNdcwZswYrr32Wg4cOJDw/E8++WSzX9jx5ObmMnv2bAYOHEhOTg49e/Zk8uTJketHjBjBypUr+fzzz6mqqmLZsmVs27Yt4b6P2LFjBzk5OUD9F/XOnTtbva+2WLx4MVOnTk2o5vbbb2fAgAE8/vjj/PKXv4xct3TpUnJzcxk9enSibQLw4IMPMmrUKGbMmNHsJcLmfPDBB7z66quMHz+eiRMnsnr16oTnf/XVV+nXrx8FBQWRaxYsWMAtt9zCgAEDmD17NvPmzYtcO2LECJYuXQrUX4KKcr41fh5J9BxrzXNQPC0+o5jZhDbtPUFmNsXMys1so5nNjbPdzOw3se3rzOzMqLVBbd8OEyfCp5+2ufbLL7/k0ksvZcGCBUd9B3is+k67drF27VoqKip44403ePfddyPV/nnECLK7d6e4uLjVff/973/nzTff5IUXXmDhwoWsXLkyUm3tp5/y5ptvMmvWLN566y1OOukk7r777sjzAlRXV7N06VIuu+yyyH3vOftsnnvqKT766CMqKys5cOAAjz32WORjHta7N7feeiulpaVMmTKF0aNHN3yXH7XvhLSlNmL9XXfdRXp6OldeeWVCtXfddRfbtm3jyiuv5MEHH4xUW1VVxV133dVyALUw76xZs9i0aRNr164lJyeHm2++OaGea2tr2bNnD6tWreLee+/l8ssvp/6qTbR6gCeeeCL+NzEt1D700EPMnz+fbdu2MX/+/IZX5lFqFy9ezMKFCykuLuaLL74gMzOzxdrWPo/w6aeJ17Yk3nWwIw+gGvgfQHpL40I8gE7AJuB0IBN4GxjeZMwFwAuAAROAf0Stjfdo9T2UWbPc09LcZ83yjz76KLF7KI1qq6urffLkyX7//fe3qv6IO++8M9o15lmzfC547kkneV5envfr18+7du3qV155ZavmdY94fTtWu336dM/Ly2tYvXLlSr/gggsSmvfZZ5/10tLSY/fbqH6Jmc8444yGVY8++qjPanIcUeY+4rbbbvOFCxdGrm16jgwZMsQrKyvdvf469pAhQyLXHtHiPZRj1D/yyCM+YcIEP3DgQMK1R2zZsuWftzVTu27dOu/bt6/n5eV5Xl6ed+rUyQcMGODbt29PeN6425r8PzUd8+1vf9tffvnlhuXTTz/dd+7cGXnumpoaz87O9m3btiX079WjRw+vq6tz9/r7nllZWZFrGysvL/exY8c2WxvveSTqOVY9c2azz0GtuYdi3jipmzCzScDDwG7gh+7+Xtviq3lmdhZwp7t/O7Z8G4C7z2s05n8Dr7j7E7HlcuA8IP9YtfGUlJR4Qu9i6NoVmrw7aAtwIdB34sQWS5e/+iqd6+oalh2YDvQE3jlGbdP6XUAG0AvYa0ZeVhYDBw7klFNOiTT3EX8xY1rv3owcOTJy3weAOiAzLY3zzz6bdevWkZeXF/fac7x5vwH8uxk3nHsuW7Zs4fDhwwwePDhyz5cD64YO5dRTT22256b1/wBmAKuBNDPys7PJysoiNzc30jED7AR6pqUxcexY1q1bx5gxY8jIyIhUu4Wjz5FNmzaRkZHBwIED2bp1KzU1NQwePDhS7RFr165l8ODBZGVlJTT37t272bRpE6NHjz7qu94otVVVVXTr1g2ATz75hL1791JYWJhQ3wCrVq2iuLiYjIyMSLWHDh1quP9RUVHB/v37GT58eLPnSLkZZ3btytixYwGorKzk0KFDDBo0iKqqKtatW8f48eP5y2uvRep79+7dbN269ah3eEXpe/Xq1RQUFNCrVy/27NnD5s2bKS4ujlRbXV1NZmYm7k55eTk9e/YkJycn8vNIlHPsSO3JwIIuXeCrr47q6bzzzuO+++6L+05OM1vj7v+0ocVAiRVmAQuAacDP3P3XLRa0kpl9H5ji7tfGlq8Cxrv79Y3G/Bm4291fiy2vAG6lPlBarG20j5nATICBAwcWf/zxx9Gb3L4dZs+meskSMmtruRxYkZbGHncyMjLIz89vuG7Z1MmHDjFr82a+8dlndKmrY4UZ57vTo1s3DsduUg4aNKjZUGhc/0FdHVcD+zt14rPMTE7OziY/P7/ZtpvOfTAtjVf79OGXffrw3qefthgoTWs3mHFRbN5a6m9wNndzPN68j/TsyW3V1dQAXbp0YejQoXGfmOPVvnjKKVy6Zw/jxo8/5iWnpvW3A3/q1Il9mZl0zcpi6NChzd5Dijf32E6d+Dg9nbq0NAYPHkzv3r0j1cY7R/r06cN7773X8EQ5fPhwMjIyItVmZGTw4YcfUlNTQ3p6Ot27d2fUqFGR5966dSvu3vDv16NHD4YMGRKpdvfu3VRVVWFmdO7cmSFDhtC5c+dItY2/LhoHSpTavXv3Ntxr69KlCwUFBXHnPZiWxtSMDF6vq6O6trahvl+/fpSXl/Pll1+SlpbG6aefTu/evSP3/f7779OjR4+jbqpHqe3WrRsbN26sf/dTWhoFBQVkZWVFqj18+HDD27L79OnDoEGDMLPIzyM9evQ45jlWVlfHN4CRvXqRlpsL6en86le/4tChQ9xwww3s2rWLXr16UVRUxPLly486z5sLlEQuSV0K1AJfAvsbP6Lu4xj7vwz4Q6Plq4DfNhnzPPBfGi2vAIqj1MZ7tOqS13XX1b/U7NIl7uWQ41abzLk7Ym1H7VvHrH+v41kbot7b+LZhMysB/ifwIXAjcEOTRwgVwIBGy/2Bpj851dyYKLVh7NgB110Hq1bVf0zkxmlbapM5d0es7ah965j173U8a0PUt+BY91DSgTuov6z078Bcd0/sR4yjNlI/1wfAJOAT6i95/1d3X99ozIXA9dTfnB8P/Mbdx0WpjSfheygiItLsJa9jvPeR1dTfs5nq7iuOS2cx7l5rZtcDy6l/19Zid19vZtfFtv8OWEZ9mGwEqoBrWqo9nv2KiMjRjvUK5T+A6919X/u11H70CkVEJHGteoXi7lcdv5ZEROTrRL96RUREglCgiIhIEAoUEREJQoEiIiJBKFBERCQIBYqIiAShQBERkSAUKCIiEoQCRUREglCgiIhIEAoUEREJQoEiIiJBKFBERCQIBYqIiAShQBERkSAUKCIiEoQCRUREglCgiIhIEAoUEREJQoEiIiJBKFBERCQIBYqIiAShQBERkSAUKCIiEoQCRUREglCgiIhIEAoUEREJQoEiIiJBKFBERCSIlAgUMzvZzF40sw9jH3s3M26KmZWb2UYzm9to/b1m9r6ZrTOzZ8ysV7s1LyIiQIoECjAXWOHuBcCK2PJRzKwTsBCYCgwHppnZ8NjmF4ER7j4K+AC4rV26FhGRBqkSKJcAj8Y+fxT4bpwx44CN7r7Z3auBJ2N1uPtf3L02Nm4V0P/4tisiIk2lSqD0c/ftALGP2XHG5ALbGi1XxNY1NQN4IXiHIiLSovT2msjM/gqcGmfT7VF3EWedN5njdqAWeLyFPmYCMwEGDhwYcWoRETmWdgsUdz+/uW1mtsPMctx9u5nlADvjDKsABjRa7g9UNtrHdOAiYJK7O81w90XAIoCSkpJmx4mISGJS5ZLXUmB67PPpwHNxxqwGCsxskJllAlfE6jCzKcCtwMXuXtUO/YqISBOpEih3A6Vm9iFQGlvGzE4zs2UAsZvu1wPLgQ3AEndfH6t/EMgCXjSztWb2u/Y+ABGRE127XfJqibt/DkyKs74SuKDR8jJgWZxx/3JcGxQRkWNKlVcoIiLSwSlQREQkCAWKiIgEoUAREZEgFCgiIhKEAkVERIJQoIiISBAKFBERCUKBIiIiQShQREQkCAWKiIgEoUAREZEgFCgiIhKEAkVERIJQoIiISBAKFBERCUKBIiIiQShQREQkCAWKiIgEoUAREZEgFCgiIhKEAkVERIJQoIiISBAKFBERCUKBIiIiQShQREQkCAWKiIgEoUAREZEgFCgiIhKEAkVERIJQoIiISBApEShmdrKZvWhmH8Y+9m5m3BQzKzezjWY2N8722WbmZtbn+HctIiKNpUSgAHOBFe5eAKyILR/FzDoBC4GpwHBgmpkNb7R9AFAKbG2XjkVE5CipEiiXAI/GPn8U+G6cMeOAje6+2d2rgSdjdUfMB+YAfhz7FBGRZqRKoPRz9+0AsY/ZccbkAtsaLVfE1mFmFwOfuPvbx5rIzGaaWZmZle3atavtnYuICADp7TWRmf0VODXOptuj7iLOOjezbrF9TI6yE3dfBCwCKCkp0asZEZFA2i1Q3P385raZ2Q4zy3H37WaWA+yMM6wCGNBouT9QCQwGBgFvm9mR9W+a2Th3/zTYAYiISItS5ZLXUmB67PPpwHNxxqwGCsxskJllAlcAS939HXfPdvd8d8+nPnjOVJiIiLSvVAmUu4FSM/uQ+ndq3Q1gZqeZ2TIAd68FrgeWAxuAJe6+Pkn9iohIE+12yasl7v45MCnO+krggkbLy4Blx9hXfuj+RETk2FLlFYqIiHRwChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCMHdPdg9JY2a7gI9bWd4H+CxgO8mkY0k9X5fjAB1LqmrLseS5e9+mK0/oQGkLMytz95Jk9xGCjiX1fF2OA3Qsqep4HIsueYmISBAKFBERCUKB0nqLkt1AQDqW1PN1OQ7QsaSq4MeieygiIhKEXqGIiEgQChQREQlCgdIGZnaZma03szoz63BvJTSzKWZWbmYbzWxusvtpCzNbbGY7zezdZPfSFmY2wMxeNrMNsXPrxmT31Fpm1sXM3jCzt2PH8otk99QWZtbJzN4ysz8nu5e2MLMtZvaOma01s7KQ+1agtM27wL8CK5PdSKLMrBOwEJgKDAemmdnw5HbVJo8AU5LdRAC1wM3uPgyYAPxbB/5/OQR8y91HA0XAFDObkNyW2uRGYEOymwjkm+5epJ9DSSHuvsHdy5PdRyuNAza6+2Z3rwaeBC5Jck+t5u4rgd3J7qOt3H27u78Z+/wL6p/AcpPbVet4vS9jixmxR4d8F5CZ9QcuBP6Q7F5SmQLlxJULbGu0XEEHfeL6ujKzfGAM8I8kt9JqsctEa4GdwIvu3lGPZQEwB6hLch8hOPAXM1tjZjND7jg95M6+jszsr8CpcTbd7u7PtXc/AVmcdR3yu8evIzPrDjwN3OTu+5PdT2u5+2GgyMx6Ac+Y2Qh371D3uczsImCnu68xs/OS3E4I57h7pZllAy+a2fuxV/htpkA5Bnc/P9k9HCcVwIBGy/2ByiT1Io2YWQb1YfK4u/9nsvsJwd33mtkr1N/n6lCBApwDXGxmFwBdgB5m9pi7/zDJfbWKu1fGPu40s2eov/wdJFB0yevEtRooMLNBZpYJXAEsTXJPJzwzM+BhYIO7/zrZ/bSFmfWNvTLBzLoC5wPvJ7WpVnD329y9v7vnU/918lJHDRMzO8nMso58DkwmYMArUNrAzL5nZhXAWcDzZrY82T1F5e61wPXAcupv/C5x9/XJ7ar1zOwJ4HVgqJlVmNmPkt1TK50DXAV8K/a2zrWx74w7ohzgZTNbR/03MC+6e4d+y+3XQD/gNTN7G3gDeN7d/0+onetXr4iISBB6hSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIpAAzSzOzlWa2tMn6brE/MfBQsnoTiUqBIpIC3L0O+G/U/0DjjEab/hf1vyJpdjL6EkmEfrBRJIWY2XXAPcBI4F+o/00G57n7a0ltTCQCBYpIion9Cp+uQD7wpLvPSW5HItEoUERSjJkNAjbFHiPc/VCSWxKJRPdQRFLPDOAr6v+kwOlJ7kUkMr1CEUkhZjYW+L/AxcAs6n877NmxP1QlktL0CkUkRZhZF+BPwCPu/gIwk/ob87qHIh2CXqGIpAgzmw98Fxjl7l/E1l0BPAoUd7Q/nSsnHgWKSAows3OBl4Dz3f2VJtuWUH8vZULsD6OJpCQFioiIBKF7KCIiEoQCRUREglCgiIhIEAoUEREJQoEiIiJBKFBERCQIBYqIiAShQBERkSD+H4Olq7y0kohcAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "mymesh.setnx(10)\n",
    "mymesh.setmeshtype('edge3') # set to 2nd order\n",
    "mymesh.update()\n",
    "mymesh.plotmesh(withnode=True,withnodeid=True)\n",
    "mymesh.printnodes() # print out all node's coordinates\n",
    "mymesh.printelements() # print out all elements' connectivity"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "fac40411-c7e2-42b3-8a86-d32240ffa529",
   "metadata": {},
   "source": [
    "The order of node index depends your implement, here, we use: 1->2, 1->2->3, 1->2->3->4\n",
    "In other code, it may look like: 1->2 (linear), 1->3->2 (2nd), 1->3->4->2 (3rd)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "9eaa986c-ead2-4ad7-bfcb-5193decbea1a",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "*** node coordinates of the mesh (total nodes=31, nodes per element=4)\n",
      "     1-th node: x= -1.000000e+00\n",
      "     2-th node: x= -8.000000e-01\n",
      "     3-th node: x= -6.000000e-01\n",
      "     4-th node: x= -4.000000e-01\n",
      "     5-th node: x= -2.000000e-01\n",
      "     6-th node: x=  0.000000e+00\n",
      "     7-th node: x=  2.000000e-01\n",
      "     8-th node: x=  4.000000e-01\n",
      "     9-th node: x=  6.000000e-01\n",
      "    10-th node: x=  8.000000e-01\n",
      "    11-th node: x=  1.000000e+00\n",
      "    12-th node: x=  1.200000e+00\n",
      "    13-th node: x=  1.400000e+00\n",
      "    14-th node: x=  1.600000e+00\n",
      "    15-th node: x=  1.800000e+00\n",
      "    16-th node: x=  2.000000e+00\n",
      "    17-th node: x=  2.200000e+00\n",
      "    18-th node: x=  2.400000e+00\n",
      "    19-th node: x=  2.600000e+00\n",
      "    20-th node: x=  2.800000e+00\n",
      "    21-th node: x=  3.000000e+00\n",
      "    22-th node: x=  3.200000e+00\n",
      "    23-th node: x=  3.400000e+00\n",
      "    24-th node: x=  3.600000e+00\n",
      "    25-th node: x=  3.800000e+00\n",
      "    26-th node: x=  4.000000e+00\n",
      "    27-th node: x=  4.200000e+00\n",
      "    28-th node: x=  4.400000e+00\n",
      "    29-th node: x=  4.600000e+00\n",
      "    30-th node: x=  4.800000e+00\n",
      "    31-th node: x=  5.000000e+00\n",
      "*** element connectivity information(bulk elements=10)\n",
      "     1-th element    1     2     3     4 \n",
      "     2-th element    4     5     6     7 \n",
      "     3-th element    7     8     9    10 \n",
      "     4-th element   10    11    12    13 \n",
      "     5-th element   13    14    15    16 \n",
      "     6-th element   16    17    18    19 \n",
      "     7-th element   19    20    21    22 \n",
      "     8-th element   22    23    24    25 \n",
      "     9-th element   25    26    27    28 \n",
      "    10-th element   28    29    30    31 \n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAejUlEQVR4nO3dfXRU9b3v8fc3PFPkQSAYeUjQgwgBDBgxthejxURQW+2pbaVVqdSL0CPWKiAePWo9XUdrtWpbj16rVtvjqsfWqrTipYoKdl0fAq2AiAgiaggYFBAKQh7me/+YIWeY7EkC+ZGZyOe11iyy9/7Ob3/3/GbPJ7NnAHN3REREWisn0w2IiMjngwJFRESCUKCIiEgQChQREQlCgSIiIkF0zHQDmdSvXz8vKCjIdBsiIu3KsmXLPnb3/qnrD+tAKSgoYOnSpZluQ0SkXTGz96PW65KXiIgEoUAREZEgFCgiIhKEAkVERIJQoIiISBAKFBERCUKBIiIiQShQREQkCAWKiIgEoUAREZEgFCgiIhKEAkVERIJQoIiISBAKFBERCUKBIiIiQShQREQkCAWKiIgEoUAREZEgFCgiIhKEAkVERIJQoIiISBAKFBERCUKBIiIiQShQREQkCAWKiIgEkVWBYmaTzGyNma0zs3kR283Mfp7YvsLMxqVs72BmfzezP7dd1yIiAlkUKGbWAbgHmAyMBKaY2ciUssnAsMRtOnBvyvYfAKsPcasiIhIhawIFGA+sc/f17l4DPAacm1JzLvAbj3sV6G1meQBmNgg4G3igLZsWEZG4bAqUgcCHScuViXUtrbkLmAvEmtqJmU03s6VmtnTLli2talhERP5HNgWKRazzltSY2TlAtbsva24n7n6/uxe7e3H//v0Ppk8REYmQTYFSCQxOWh4EVLWw5kvAV81sA/FLZV82s/86dK2KiEiqbAqUCmCYmQ01s87ABcD8lJr5wMWJb3uVAJ+6+yZ3v9bdB7l7QeJ+L7j7hW3avYjIYa5jphvYx93rzOxyYCHQAXjI3VeZ2YzE9vuABcBZwDpgN3BJpvoVEZH9mXvqxxSHj+LiYl+6dGmm2xARaVfMbJm7F6euz6ZLXiIi0o4pUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCyKpAMbNJZrbGzNaZ2byI7WZmP09sX2Fm4xLrB5vZi2a22sxWmdkP2r57EZHDW9YEipl1AO4BJgMjgSlmNjKlbDIwLHGbDtybWF8HXO3uI4AS4F8i7isiIodQ1gQKMB5Y5+7r3b0GeAw4N6XmXOA3Hvcq0NvM8tx9k7v/DcDddwKrgYFt2byIyOEumwJlIPBh0nIljUOh2RozKwDGAq+Fb1FERNLJpkCxiHV+IDVm1gN4ArjS3XdE7sRsupktNbOlW7ZsOehmRURkf9kUKJXA4KTlQUBVS2vMrBPxMHnU3f+Ybifufr+7F7t7cf/+/YM0LiIi2RUoFcAwMxtqZp2BC4D5KTXzgYsT3/YqAT51901mZsCDwGp3/1nbti0iIgAdM93APu5eZ2aXAwuBDsBD7r7KzGYktt8HLADOAtYBu4FLEnf/EnARsNLM3kis+1d3X9CGhyAiclgz99SPKQ4fxcXFvnTp0ky3ISLSrpjZMncvTl2fTZe8RESkHVOgiIhIEAoUEREJQoEiIiJBKFBERCQIBYqIiAShQBERkSAUKCIiEoQCRUREglCgiIhIEAoUEREJQoEiIiJBKFBERCQIBYqIiAShQBERkSAUKCIiEoQCRUREglCgiIhIEAoUEREJQoEiIiJBKFBERCQIBYqIiAShQBERkSAUKCIiEoQCRUREglCgiIhIEAoUEREJQoHSCtOmTSM3N5dRo0Y1Wffhhx9y+umnM2LECAoLC7n77rsj6/bs2cP48eM54YQTKCws5MYbb2xy3Pr6esaOHcs555zTZF1BQQGjR4+mqKiI4uLitHXbt2/n/PPP5/jjj2fEiBG88sorjWrWrFlDUVFRw61nz57cddddkePdeeedFBYWMmrUKKZMmcKePXsi6+6++25GjRpFYWFh2rEg+vHeunUrZWVlDBs2jLKyMrZt2xZZ9/vf/57CwkJycnJYunRp2vHmzJnD8ccfz5gxY/ja177G9u3bI+v+7d/+jTFjxlBUVER5eTlVVVVNPh9uv/12zIyPP/44su6mm25i4MCBDY/rggUL0o73i1/8guHDh1NYWMjcuXPTHsu3vvWthvEKCgooKiqKrHvjjTcoKSlpeH68/vrrkXXLly/nlFNOYfTo0XzlK19hx44daZ/bqfOycuXKyLrUeUk3Xuq8vPnmm5F1qfNSUVHR5Lm3b16WL18eWRc1L02dz8lzM2PGjMi61HkZOXJkZF3qvPzpT3+KrEudl+rq6sjXkdQ52bRpU2Rd1LnSYu5+2N5OPPFEb43Fixf7smXLvLCwsMm6qqoqX7Zsmbu779ixw4cNG+arVq1qVBeLxXznzp3u7l5TU+Pjx4/3V155Je24d9xxh0+ZMsXPPvvsJvefn5/vW7Zsae5w/OKLL/Zf/epX7u6+d+9e37ZtW5P1dXV1PmDAAN+wYUOjbZWVlV5QUOC7d+92d/dvfOMb/utf/7pR3cqVK72wsNB37drltbW1PnHiRH/nnXci9xf1eM+ZM8dvueUWd3e/5ZZbfO7cuZF1b731lr/99tteWlrqFRUVacdbuHCh19bWurv73Llz04736aefNvx89913+2WXXZb2+fDBBx94eXm5DxkyxLds2RJZd+ONN/pPf/rTZo/3hRde8IkTJ/qePXvc3f2jjz5KW5vsqquu8h/96EeRdWVlZb5gwQJ3d3/mmWe8tLQ0sq64uNhfeukld3d/8MEH/frrr0/73E6dl+9///uRdanzkm681HlJN17qvFx44YVpz73keVm5cmVkXdS8pOsxdW5WrFjR7Hl/1VVX+ezZsyPrUuelpKQksi51Xq677rrI15HUOZkzZ05kXdS5kgpY6hGvqXqH0gqnnnoqRx55ZLN1eXl5jBs3DoAjjjiCESNGsHHjxkZ1ZkaPHj0AqK2tpba2FjOLHLOyspJnnnmGSy+9tBVH8D927NjBkiVL+N73vgdA586d6d27d5P3WbRoEcceeyz5+fmR2+vq6vjss8+oq6tj9+7dHH300Y1qVq9eTUlJCd27d6djx46Ulpby5JNPRo4X9Xg//fTTTJ06FYCpU6fy1FNPRdaNGDGC4cOHNzteeXk5HTt2BKCkpITKysrIup49ezb8vGvXLsws7fPhhz/8IbfddlvDXLb0eRNVd++99zJv3jy6dOkCQG5ubrNjujuPP/44U6ZMiawzM3bs2AHAp59+ytFHHx1Zt2bNGk499VQAysrKeOKJJ9I+t1Pn5fnnn4+sS52XdOOlzsv27dsj61LnpUePHmnPveR5Oeqoo1p0jjbVY+rcjB49uskx983L9OnTI+tS52Xo0KGRdanz8sc//jHydSR1Tp5++unIuqhzpaUUKG1sw4YN/P3vf+fkk0+O3F5fX09RURG5ubmUlZWlrbvyyiu57bbbyMlpfgrNjPLyck488UTuv//+yJr169fTv39/LrnkEsaOHcull17Krl27mhz3scceY8qUKZHbBg4cyOzZsxkyZAh5eXn06tWL8vLyRnWjRo1iyZIlfPLJJ+zevbvhkkJLffTRR+Tl5QHxE726urrF923OQw89xOTJk9Nuv+666xg8eDCPPvooN998c2TN/PnzGThwICeccEKz+/vlL3/JmDFjmDZtGtu2bYuseeedd3j55Zc5+eSTKS0tpaKiotlxX375ZQYMGMCwYcMit991113MmTOHwYMHM3v2bG655ZbIulGjRjF//nwgflkkdZ6Sn9tNzUtz50BzdanzklqXbl6S65qal9TxmpqX5Nqm5ibqWKLmJbmuqXlJroual6jXkag5aenrTYtFvW3ZdwNKmtoe+gZMAtYA64B5EdsN+Hli+wpgXEvvG3U76EteVVXup57qvmmTv/fee+kveSXVubvv3LnTx40b50888USTde7u27Zt89NOO81XrlzZqPZPhYU+c+pUd3d/8cUXoy95JY25ceNGd49fHhkzZowvXry4UV3Fs896hw4d/NVXX3V39yuuuMKvv/76tD3u3bvX+/bt65s3b47c79bVq/3000/36upqr6mp8XPPPdd/+9vfRo73wAMP+NixY33ChAl+2WWX+ZVXXpl2v6mPd69evfar7d2xY5PzUlpa6hXPPtvs/P34xz/28yZN8tiECc3O839ce63fMGRIo7pdu3b5+PHjffv27e7unj9okG855ZTI8TZv3ux1dXVeX1/v/3rFFX7JgAGRdYWFhT5r1iyPxWL+2jPPeEGXLh6rqop8bPaZMWOG337DDWmPedasWf6HP/zB3d3/+777fGLv3pF1q1ev9rKyMh83bpzfdPXVfmTisXZv/NxumJfE/PVOLKc7B0pPOcUrxo5t9lz58TXX+Hl9+zYcc9pzKmVekuuampfU8dLNS9S+G+Zm40Z/rajICwYPbriMHdXjjIsv9tuPOSbteA3zUlXl/z1ihE+cMCGyrmFeRo/2m/Lz/cg+fRr2kfw6ku5cSa1rmJODuOTV3At8DfDvQMem6kLcgA7Au8AxQGdgOTAypeYs4NlEsJQAr7X0vlG3gw6UmTPdc3LcZ85sOlCS6mpqary8vNzvuOOOJuuS3XTTTY2u3/rMmT4PfOAXvuD5+fk+YMAA79atm3/nO99p0ZiNrgkn6jZNner5+fkNq5csWeJnnXVW2vGeeuopLysrS3ssj5eX+7Rp0xpWP/LIIz4zuZc0/V177bV+zz33pK1LfbyPO+44r0q8uFRdfLEfB03OS2lpqVd8/etNzt/DDz/sJSUlvuvSS1s0zxsuvNALI/a7YsUK79+/v+fn53t+fr53MPPB4JumTm1yvPe+853I8dzdzzzzTH/xxRcbHptjwKu/+93Ix8bdvba21nNzc/3Diy5Keyw9e/b0WCzm7u6xGTP8iGYeQ3f3NVOm+EmJuqjndsO8zJzpVWZ+XK9eTZ4DpUcf7RVmTZ4rDz/8sJcMGOC7mqnbZ9+81Eyfvl9dunn54KKLmhwveV6i9t0wN4nn7DE9e/rGjRsjx6ytrfXcbt38wyaOpWFeZs70mJkf0alTs68ja8z8pNzc/Vbvex1Jd66k1jXMyUEEisW3RTOzicCDwFbgQnd/qzXvhppiZqcAN7n7mYnlawHc/Zakmv8DvOTuv0ssrwFOAwqau2+U4uJiP6BvMXTrBinfVNoAnA30Ly1tWLfw5ZfpEos1LDswFegFrGyibgvQCeiWk8MZX/wiK1asYMiQIfTt27dR7T5/MWNKnz6MHj06csxdQAzonDRmfn4+FatWNRpvAvCfZsw69VQ2bNhAfX0966qqIvf7TWDF8OEcddRRkft9DZgG/NWM8yZMYM2aNRxxxBG8uX59o/GqgV45OZSedBIrVqxg7NixvPDqq5H7XWPGuG7dOOmkkwB49913uWzjRq5351biT9TbSD8vZ8Zi3A7s+65bat3WrVupXbmSJUD/pP2m1u3evZuXly2jSyzGL4DFwB+a2G+XWIwCYCnQL6Ju7969vPT663SJxbgz8fg9FlFXVVXF+WvX8h/AO8BE4APiv2GlPjYAVy9Zwu3uLG7iWCoqKnjys88od2cRMBdYFlFXU1PDi6+9RqdYjO8SP/EuIfq5PW7xYnKBecCtwCfARxF1+x6b04DbgRPTjHf1kiVckziO/qQ/p1Ln5SXgC8CRwE9ycjhzwoTIeakAro6oi5qX36XZ9/lLllDtzs2Jufly4hZ1zC/GYtxC/HmT7lg2L17MfYnHeREwBxgVUffbJUsY7E4M+G7SY5j6OrJ9+/b9zpVK4MdAb+CzLl0oP+kkrrnmmoZvjZ522mncfvvtkd8MNbNl7t5oQ5OBkrjjEcBdwBTgenf/WZN3OEhmdj4wyd0vTSxfBJzs7pcn1fwZuNXd/5pYXgRcQzxQmrxv0hjTgekAQ4YMOfH9999veZObNsHs2dQ8/jid6+r4JrAoJ4dt7nTq1ImCggLy8vI4cu9eZq5fz4SPP6ZrLMYiM85wp2f37tQnPpgdOnQow3r02K+uwowpHTqwpXNn6oD+/ftTUFAA0GjMPTk5vNyvHzf368dbmzc3BEpq3WozzunQgY8TY+bm5pKfnx853sO9enFtTQ21QNeuXRk+fDgDYrFGdc/17cvXt21j/MknN3xQGjXe/+7aladiMWI5OfTo0YPhw4fTr7a2Ud1JHTrwfseOxHJyOPbYY+nTp0/keJM7deKVWIyaurqGx7tfv36sffNNOu3axXH19TwBzIiYl06dOvHu2rXU19TQBygifmKm1n3wwQdYLEZeLEaPujpOIR5SqXVbt25l765d9KupYUR9PfcDP2zm+XBpdTXLgO9H1G3fvp09O3fSt6aGkfX1/Ar4QUTdgAED2LB6NTnbttGnvp47gC+meWzy8vJ4b9Uqztu7l5/s2kXXWCzyOdu9e3c2rF3LkXv3cmRdHfcBP4moq6+vZ/PGjfSureWC+np+CryQ5rk9pEsXqlasYGdtLfnALOBr0KjuCzU1vLtuHTtiMXoDQ4iHbmrde2vX0q22lqNjMQwYDPw5om7z5s37zculwD8DQzt0YGuXLtSZMXToUPr27bvfvNwDfCWirrq6utG8vJXmmI/t3p2dy5ezee9eugDfJh4CUef9P954g6/s3css97SvD31iMTa98w7d6uroRjwsvh8xnu3cydbKSrrV1/P1xPF+u2PHRq8jtbW1+50rNwNXmFHfqxexAQP45re/zQ033MCTTz7JrFmz2LJlC71796aoqIiFCxfu91KYLlAO5JLU14E64B/AjuRbS8doZvxvAA8kLV8E/CKl5hngfyUtLyIeyM3eN+p2UJe8ZsyIXzro2jXyks0hq8vkvrO9rj30qMdGx5wNdQdamwat+dqwmRUTf3e0lvgvTrNSbiFUEv/FY59BQFULa1py3zA++ghmzIBXX43/uXlz29Rlct/ZXtceetRjo2POhroDrT1AzX2G0hG4kfhlpf8k/u2pPcH23nhf+y4LbyR+WfPb7r4qqeZs4HLiH86fDPzc3ce35L5RDvgzFBERSXvJq2Mz96sg/jnVZHdfdEg6S3D3OjO7HFhI/FtbD7n7KjObkdh+H7CAeJisA3YT/0ww7X0PZb8iIrK/5t6h/Ba43N0/bbuW2o7eoYiIHLiDeofi7hcdupZEROTzRP/0ioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQWREoZnakmT1nZmsTf/ZJUzfJzNaY2Tozm5e0/qdm9raZrTCzJ82sd5s1LyIiQJYECjAPWOTuw4BFieX9mFkH4B5gMjASmGJmIxObnwNGufsY4B3g2jbpWkREGmRLoJwLPJL4+RHgvIia8cA6d1/v7jXAY4n74e5/cfe6RN2rwKBD266IiKTKlkAZ4O6bABJ/5kbUDAQ+TFquTKxLNQ14NniHIiLSpI5ttSMzex44KmLTdS0dImKdp+zjOqAOeLSJPqYD0wGGDBnSwl2LiEhz2ixQ3P2MdNvM7CMzy3P3TWaWB1RHlFUCg5OWBwFVSWNMBc4BJrq7k4a73w/cD1BcXJy2TkREDky2XPKaD0xN/DwVeDqipgIYZmZDzawzcEHifpjZJOAa4KvuvrsN+hURkRTZEii3AmVmthYoSyxjZkeb2QKAxIfulwMLgdXA4+6+KnH/XwJHAM+Z2Rtmdl9bH4CIyOGuzS55NcXdPwEmRqyvAs5KWl4ALIio+6dD2qCIiDQrW96hiIhIO6dAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkFkRaCY2ZFm9pyZrU382SdN3SQzW2Nm68xsXsT22WbmZtbv0HctIiLJsiJQgHnAIncfBixKLO/HzDoA9wCTgZHAFDMbmbR9MFAGfNAmHYuIyH6yJVDOBR5J/PwIcF5EzXhgnbuvd/ca4LHE/fa5E5gL+CHsU0RE0siWQBng7psAEn/mRtQMBD5MWq5MrMPMvgpsdPflze3IzKab2VIzW7ply5bWdy4iIgB0bKsdmdnzwFERm65r6RAR69zMuifGKG/JIO5+P3A/QHFxsd7NiIgE0maB4u5npNtmZh+ZWZ67bzKzPKA6oqwSGJy0PAioAo4FhgLLzWzf+r+Z2Xh33xzsAEREpEnZcslrPjA18fNU4OmImgpgmJkNNbPOwAXAfHdf6e657l7g7gXEg2ecwkREpG1lS6DcCpSZ2Vri39S6FcDMjjazBQDuXgdcDiwEVgOPu/uqDPUrIiIp2uySV1Pc/RNgYsT6KuCspOUFwIJmxioI3Z+IiDQvW96hiIhIO6dAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkEoUEREJAhz90z3kDFmtgV4/yDv3g/4OGA7maRjyT6fl+MAHUu2as2x5Lt7/9SVh3WgtIaZLXX34kz3EYKOJft8Xo4DdCzZ6lAciy55iYhIEAoUEREJQoFy8O7PdAMB6Viyz+flOEDHkq2CH4s+QxERkSD0DkVERIJQoIiISBAKlFYws2+Y2Sozi5lZu/sqoZlNMrM1ZrbOzOZlup/WMLOHzKzazN7MdC+tYWaDzexFM1udeG79INM9HSwz62pmr5vZ8sSx/CjTPbWGmXUws7+b2Z8z3UtrmNkGM1tpZm+Y2dKQYytQWudN4J+BJZlu5ECZWQfgHmAyMBKYYmYjM9tVqzwMTMp0EwHUAVe7+wigBPiXdjwve4Evu/sJQBEwycxKMttSq/wAWJ3pJgI53d2L9PdQsoi7r3b3NZnu4yCNB9a5+3p3rwEeA87NcE8Hzd2XAFsz3Udrufsmd/9b4uedxF/ABma2q4Pjcf9ILHZK3Nrlt4DMbBBwNvBApnvJZgqUw9dA4MOk5Ura6QvX55WZFQBjgdcy3MpBS1wmegOoBp5z9/Z6LHcBc4FYhvsIwYG/mNkyM5secuCOIQf7PDKz54GjIjZd5+5Pt3U/AVnEunb52+PnkZn1AJ4ArnT3HZnu52C5ez1QZGa9gSfNbJS7t6vPuczsHKDa3ZeZ2WkZbieEL7l7lZnlAs+Z2duJd/itpkBphrufkekeDpFKYHDS8iCgKkO9SBIz60Q8TB519z9mup8Q3H27mb1E/HOudhUowJeAr5rZWUBXoKeZ/Ze7X5jhvg6Ku1cl/qw2syeJX/4OEii65HX4qgCGmdlQM+sMXADMz3BPhz0zM+BBYLW7/yzT/bSGmfVPvDPBzLoBZwBvZ7Spg+Du17r7IHcvIH6evNBew8TMvmBmR+z7GSgnYMArUFrBzL5mZpXAKcAzZrYw0z21lLvXAZcDC4l/8Pu4u6/KbFcHz8x+B7wCDDezSjP7XqZ7OkhfAi4Cvpz4Wucbid+M26M84EUzW0H8F5jn3L1df+X2c2AA8FczWw68Djzj7v831OD6p1dERCQIvUMREZEgFCgiIhKEAkVERIJQoIiISBAKFBERCUKBIiIiQShQRLKAmeWY2RIzm5+yvnvivxi4N1O9ibSUAkUkC7h7DPgu8b/QOC1p00+I/xNJszPRl8iB0F9sFMkiZjYDuA0YDfwT8X/J4DR3/2tGGxNpAQWKSJZJ/BM+3YAC4DF3n5vZjkRaRoEikmXMbCjwbuI2yt33ZrglkRbRZygi2Wca8Bnx/1LgmAz3ItJieocikkXM7CTg/wFfBWYS/9dhv5j4j6pEspreoYhkCTPrCvwGeNjdnwWmE/9gXp+hSLugdygiWcLM7gTOA8a4+87EuguAR4AT29t/nSuHHwWKSBYws1OBF4Az3P2llG2PE/8spSTxH6OJZCUFioiIBKHPUEREJAgFioiIBKFAERGRIBQoIiIShAJFRESCUKCIiEgQChQREQlCgSIiIkH8f0/ZWS82AoEoAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "mymesh.setmeshtype('edge4') # default option is edge2, you can choose edge2~edge4\n",
    "mymesh.setdomainsize(xmin=-1.0,xmax=5.0)\n",
    "mymesh.update()\n",
    "mymesh.plotmesh(withnode=True,withnodeid=True)\n",
    "mymesh.printnodes() # print out all node's coordinates\n",
    "mymesh.printelements() # print out all elements' connectivity"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d84242a3-c88e-4a04-b482-c777aec3504a",
   "metadata": {},
   "source": [
    "Now, lets try the 2d mesh case.\n",
    "Before we run the code, you must import mesh2d"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "c37fc2fb-ccc8-4eea-90d3-f2473ba81ef9",
   "metadata": {},
   "outputs": [],
   "source": [
    "from FEToy.mesh.lagrange2dmesh import mesh2d"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "7420d43a-bfe5-4bc6-9726-d319811991c4",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "*** node coordinates of the mesh (total nodes=9, nodes per element=4)\n",
      "     1-th node: x=  0.000000e+00,y=  0.000000e+00\n",
      "     2-th node: x=  1.000000e+00,y=  0.000000e+00\n",
      "     3-th node: x=  2.000000e+00,y=  0.000000e+00\n",
      "     4-th node: x=  0.000000e+00,y=  1.000000e+00\n",
      "     5-th node: x=  1.000000e+00,y=  1.000000e+00\n",
      "     6-th node: x=  2.000000e+00,y=  1.000000e+00\n",
      "     7-th node: x=  0.000000e+00,y=  2.000000e+00\n",
      "     8-th node: x=  1.000000e+00,y=  2.000000e+00\n",
      "     9-th node: x=  2.000000e+00,y=  2.000000e+00\n",
      "*** element connectivity information(bulk elements=4)\n",
      "     1-th element :    0     1     4     3 \n",
      "     2-th element :    1     2     5     4 \n",
      "     3-th element :    3     4     7     6 \n",
      "     4-th element :    4     5     8     7 \n"
     ]
    }
   ],
   "source": [
    "mynewmesh=mesh2d(nx=2,ny=2,meshtype='quad4') # order=1\n",
    "mynewmesh.createmesh()\n",
    "mynewmesh.printnodes()\n",
    "mynewmesh.printelements()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "id": "6821fba6-73df-4b9f-bedf-501b83a95de6",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEKCAYAAAAB0GKPAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAcZUlEQVR4nO3df5TVdb3v8edrGIZfgkgMiIAOt8syhABxAstSuF0IUq9HpRZoWSqLZMVaeu+tq/WHuuyu1bl5W6mFERXH0+kEt1WRVIPozbz4i2TAUUGjOISHEZUJUwIEHHjfP/ZGd+Oemc84s797mHk91tpr9vfz47vf7vVxXny/3z37q4jAzMysPRXlLsDMzE4ODgwzM0viwDAzsyQODDMzS+LAMDOzJJXlLqCUhg8fHjU1NeUuw8zspLF58+a/RER1sb4eHRg1NTXU19eXuwwzs5OGpBdb6/MpKTMzS+LAaMP27duZOnXq248hQ4Zw1113lbss60G+9a1vMXHiRCZNmsTChQs5fPhwuUuyHuTuu+9m0qRJTJw4sUt+dzkw2nD22WfT0NBAQ0MDmzdvZuDAgVx++eXlLst6iJdeeol77rmH+vp6tm7dyrFjx1i9enW5y7IeYuvWrXz/+9/nqaee4plnnuHXv/41f/rTnzq1TwdGot/+9re8//3v56yzzip3KdaDNDc38+abb9Lc3MyhQ4c444wzyl2S9RAvvPAC559/PgMHDqSyspKLLrqINWvWdGqfDoxEq1evZuHCheUuw3qQ0aNH86UvfYkzzzyTUaNGceqppzJnzpxyl2U9xKRJk9iwYQP79u3j0KFD1NXVsXv37k7tM7PAkDRW0u8kvSBpm6Qbi4yRpHsk7ZD0rKRpBX1zJW3P991S0mJffhkuugheeQWAo0ePsnbtWj71qU+V9GWtl8ivr7/+4Q/cf//9/PnPf2bPnj0cPHiQH//4x+Wuzk52+fU14bTTuPnmm5k9ezZz585lypQpVFZ27oOxWR5hNAP/PSImAOcDX5R0Tosx84Dx+cdi4LsAkvoAy/L95wALi8ztOl/7Gjz2GNxxBwDr1q1j2rRpjBw5smQvab1Ifn393xtvZNy4cVRXV9O3b1+uuOIKnnjiiXJXZye7gt9f119/PVu2bGHDhg0MGzaM8ePHd2rXKtfXm0u6H/hORDxU0PY94JGIWJXf3g7MBGqA2yPiE/n2rwBExNfbeo3a2tro0N9hDBgART6l8mng2bPP5vTTT0/fl1kL6x99lH7Hj7+9/XvgOuAxiX/42MfYvn07gwcPZvTo0WWr0U5eLdcXwF5gRP/+/Pv27cyZM4cnn3yS0047rc39SNocEbXF+spyDUNSDXAuuf9nCo0GCk+yNebbWmsvtu/Fkuol1Tc1NXWssJ074aqrOFyRe1sOV1Twq+pqfllZyfDhwzu2L7MWFk6fzkMjRnAwvz2looJpAwdyZr9+1NfXExGMGjWqrDXayevE+jp0omHgQK6sruacsWO59NJLWbZsWbth0Z7M/9Jb0inAz4GbImJ/y+4iU6KN9nc3RqwAVkDuCKNDxY0aBUOGUHX8OEcqKugPXDp/PkfvvbdDuzFr1ZIlHFu+nDeBAcC/fO5z4PVlXWXJEo4vX86Rigr6HT7Mo128vjINDEl9yYXFv0bEL4oMaQTGFmyPAfYAVa20d71XX2XtGWfwq1Gj+OGMGbkLSGZd5dVX+aeqKu6rquKxa67x+rKuVeLfX5ldw5Ak4J+B1yLiplbGXAwsBT4JzADuiYjpkiqBPwIfB14CNgFXRcS2tl6zw9cw8mbOnAnAI4880uG5Zu0ZOnQoAK+//npZ67CeqbO/v9q6hpHlEcYFwGeB5yQ15Nu+CpwJEBHLgTpyYbEDOARcm+9rlrQUWA/0AVa2FxZmZta1MguMiHiM4tciCscE8MVW+urIBYqZmZWB/9LbzMySODDMzCyJA8PMzJI4MMzMLIkDw8zMkjgwzMwsiQPDzMySODDMzCyJA8PMzJI4MMzMLIkDw8zMkjgwzMwsiQPDzMySODDMzCyJA8PMzJI4MMzMLElmN1CStBK4BNgbEZOK9H8ZuLqgrglAdUS8JmkX8DfgGNDc2u0DzcysdLI8wrgPmNtaZ0TcGRFTI2Iq8BXg/0XEawVDZuX7HRZmZmWQWWBExAbgtXYH5iwEVpWwHDMz66Budw1D0kByRyI/L2gO4EFJmyUtbmf+Ykn1kuqbmppKWaqZWa/S7QIDuBR4vMXpqAsiYhowD/iipAtbmxwRKyKiNiJqq6urS12rmVmv0R0DYwEtTkdFxJ78z73AGmB6GeoyM+vVulVgSDoVuAi4v6BtkKTBJ54Dc4Ct5anQzKz3yvJjtauAmcBwSY3AbUBfgIhYnh92OfBgRBwsmDoSWCPpRL0/iYgHsqrbzMxyMguMiFiYMOY+ch+/LWzbCUwpTVVmZpaqW52SMjOz7suBYWZmSRwYZmaWxIFhZmZJHBhmZpbEgWFmZkkcGGZmlsSBYWZmSRwYZmaWxIFhZmZJHBhmZpbEgWFmZkkcGGZmlsSBYWZmSRwYZmaWxIFhZmZJMgsMSSsl7ZVU9PaqkmZKekNSQ/5xa0HfXEnbJe2QdEtWNZuZ2TuyPMK4D5jbzphHI2Jq/nEHgKQ+wDJgHnAOsFDSOSWt1MzM3iWzwIiIDcBr72HqdGBHROyMiKPAauCyLi3OzMza1d2uYXxY0jOS1kmamG8bDewuGNOYbytK0mJJ9ZLqm5qaSlmrmVmv0p0CYwtwVkRMAb4N/DLfriJjo7WdRMSKiKiNiNrq6uqur9LMrJfqNoEREfsj4kD+eR3QV9JwckcUYwuGjgH2lKFEM7NerdsEhqTTJSn/fDq52vYBm4DxksZJqgIWAGvLV6mZWe9UmdULSVoFzASGS2oEbgP6AkTEcmA+sERSM/AmsCAiAmiWtBRYD/QBVkbEtqzqNjOznMwCIyIWttP/HeA7rfTVAXWlqMvMzNJ0m1NSZmbWvTkwzMwsiQPDzMySODDMzCyJA8PMzJI4MMzMLIkDw8zMkjgwzMwsiQPDzMySODDMzCyJA8PMzJI4MMzMLIkDw8zMkjgwzMwsiQPDzMySZBYYklZK2itpayv9V0t6Nv94QtKUgr5dkp6T1CCpPquazczsHVkeYdwHzG2j/8/ARRExGfgasKJF/6yImBoRtSWqz8zM2pDlHfc2SKppo/+Jgs2NwJiSF2VmZsm66zWM64F1BdsBPChps6TFbU2UtFhSvaT6pqamkhZpZtabZHaEkUrSLHKB8dGC5gsiYo+kEcBDkv4QERuKzY+IFeRPZ9XW1kbJCzYz6yW61RGGpMnAD4DLImLfifaI2JP/uRdYA0wvT4VmZr1XtwkMSWcCvwA+GxF/LGgfJGnwiefAHKDoJ63MzKx0MjslJWkVMBMYLqkRuA3oCxARy4FbgfcB90oCaM5/ImoksCbfVgn8JCIeyKpuMzPLyfJTUgvb6V8ELCrSvhOY8u4ZZmaWpW5zSsrMzLo3B4aZmSVxYJiZWRIHhpmZJXFgmJlZEgeGmZklcWCYmVkSB4aZmSVxYJiZWRIHhpmZJXFgmJlZEgeGmZklcWC049ixY5x77rlccskl5S7FeqCamho++MEPMnXqVGprfbt661qvv/468+fP5wMf+AATJkzgySef7NT+ut0d97qbu+++mwkTJrB///5yl2I91O9+9zuGDx9e7jKsB7rxxhuZO3cuP/vZzzh69CiHDh3q1P58hNGGxsZGfvOb37Bo0bu+dd3MrFvbv38/GzZs4PrrrwegqqqKoUOHdmqfDow23HTTTXzjG9+gosJvk5WGJObMmcN5553HihUryl2O9SA7d+6kurqaa6+9lnPPPZdFixZx8ODBTu2zzd+Eks7v1N7/fl8rJe2VVPT2qsq5R9IOSc9KmlbQN1fS9nzfLV1VU2uGHTnCFZs2MeKUUzjvvPNK/XLWy4w8fpzfHDgAr7zC448/zpYtW1i3bh3Lli1jw4YN5S7PTnLDjhzhroYGml95hS1btrBkyRKefvppBg0axD/+4z92at/t/dN5g6SvSeqKax33AXPb6J8HjM8/FgPfBZDUB1iW7z8HWCjpnC6op1XXvPgiew4dYu3PfkZNTQ0LFizg4Ycf5jOf+UwpX9Z6iS8fOcL5x47BHXdwxhlnADBixAguv/xynnrqqTJXZye7a158kclvvMGY1asZM2YMM2bMAGD+/Pls2bKlU/tWRLTeKX0c+CHwGvCZiHi+Uy8m1QC/johJRfq+BzwSEavy29vJ3QO8Brg9Ij6Rb/8KQER8vb3Xq62tjfr6+vQCBwyAw4ff1fwgcHFlJYMGDUrfl1kLr7zxBv0Ltg8Cx8nd2H7kkCEcPHiQfv360bdv3/IUaCe1lusL4GPAD6qqOPvIEW6//XYOHjzInXfe2eZ+JG2OiKIf2WvzCCMifgt8EHgaqJf039LL77DRwO6C7cZ8W2vtRUlaLKleUn1TU1PHKti5E666ihNn+Q4BP+3bl+sGDOjYfsyKmDJ4MD/t2/ft9bULmAScWlHBgQMHqKysdFjYe9ZyfTFwIN+eN4+rP/ABJk+eTENDA1/96lc79RrtnmqKiL8B10uqA/6PpDvI/cOocMyQTlWRo2Iv30Z7URGxAlgBuSOMDlUwahQMGcIA4EhFBQOBTy9axKfvvbdDuzFr1ZIlHFu+nDeBiRUVvPiFL4DXl3WVJUs4vnw5Ryoq6Hf4MFNraqivq+uy3Sddm5BUC/xP4E/A/waau6yCdzQCYwu2xwB7gKpW2kvj1VdZe8YZ/GrUKH44Ywa8/HLJXsp6oVdf5Z+qqrivqorHrrnG68u6Vol/f7V3DaMSuA24GbgXuCUi3n2SP/XF2r6GcTGwFPgkMAO4JyKm52v4I/Bx4CVgE3BVRGxr7/U6fA0jb+bMmQA88sgjHZ5r1p4Tn4V//fXXy1qH9Uyd/f3V1jWM9o4wNgHDgHn56xnvmaRV5C5iD5fUSC6I+gJExHKgjlxY7CB3+eDafF+zpKXAeqAPsDIlLMzMrGu1FxhbgaUR8UZnXygiFrbTH8AXW+mrIxcoZmZWJm0GRkR8NqtCzMyse/N3XpiZWRIHhpmZJXFgmJlZEgeGmZklcWCYmVkSB4aZmSVxYJiZWRIHhpmZJXFgmJlZEgeGmZklcWCYmVkSB4aZmSVxYJiZWRIHhpmZJXFgmJlZkkwDQ9JcSdsl7ZB0S5H+L0tqyD+2SjomaVi+b5ek5/J9Hb/vqpmZdUp7d9zrMpL6AMuA2UAjsEnS2oh4/sSYiLgTuDM//lLgv0bEawW7mRURf8mqZjMze0eWRxjTgR0RsTMijgKrgcvaGL8QWJVJZWZm1q4sA2M0sLtguzHf9i6SBgJzgZ8XNAfwoKTNkha39iKSFkuql1Tf1NTUBWWbmRlkGxgq0hatjL0UeLzF6agLImIaMA/4oqQLi02MiBURURsRtdXV1Z2r2MzM3pZlYDQCYwu2xwB7Whm7gBanoyJiT/7nXmANuVNcZmaWkSwDYxMwXtI4SVXkQmFty0GSTgUuAu4vaBskafCJ58AcYGsmVZuZGZDhp6QiolnSUmA90AdYGRHbJN2Q71+eH3o58GBEHCyYPhJYI+lEzT+JiAeyqt3MzDIMDICIqAPqWrQtb7F9H3Bfi7adwJQSl2dmZm3wX3qbmVkSB4aZmSVxYJiZWRIHhpmZJXFgmJlZEgeGmZklcWCYmVkSB4aZmSVxYJiZWRIHhpmZJXFgmJlZEgeGmZklcWCYmVkSB4aZmSVxYJiZWZJMA0PSXEnbJe2QdEuR/pmS3pDUkH/cmjrXzMxKK7MbKEnqAywDZpO7v/cmSWsj4vkWQx+NiEve41wzMyuRLI8wpgM7ImJnRBwFVgOXZTDXzMy6QJaBMRrYXbDdmG9r6cOSnpG0TtLEDs5F0mJJ9ZLqm5qauqJuMzMj28BQkbZosb0FOCsipgDfBn7Zgbm5xogVEVEbEbXV1dXvtVYzM2shy8BoBMYWbI8B9hQOiIj9EXEg/7wO6CtpeMpcMzMrrSwDYxMwXtI4SVXAAmBt4QBJp0tS/vn0fH37UuaamVlpZfYpqYholrQUWA/0AVZGxDZJN+T7lwPzgSWSmoE3gQUREUDRuVnVbmZmGQYGvH2aqa5F2/KC598BvpM618zMsuO/9DYzsyQODDMzS+LAMDOzJA4MMzNL4sAwM7MkDgwzM0viwDAzsyQODDMzS+LAMDOzJA4MMzNL4sAwM7MkDgwzM0viwDAzsyQODDMzS+LAMDOzJA4MMzNLkmlgSJorabukHZJuKdJ/taRn848nJE0p6Nsl6TlJDZLqs6zbzMwyvOOepD7AMmA20AhskrQ2Ip4vGPZn4KKI+KukecAKYEZB/6yI+EtWNZuZ2TuyPMKYDuyIiJ0RcRRYDVxWOCAinoiIv+Y3NwJjMqzPzMzakGVgjAZ2F2w35ttacz2wrmA7gAclbZa0uLVJkhZLqpdU39TU1KmCzczsHZmdkgJUpC2KDpRmkQuMjxY0XxAReySNAB6S9IeI2PCuHUasIHcqi9ra2qL7NzOzjsvyCKMRGFuwPQbY03KQpMnAD4DLImLfifaI2JP/uRdYQ+4Ul5mZZSTLwNgEjJc0TlIVsABYWzhA0pnAL4DPRsQfC9oHSRp84jkwB9iaWeVmZpbdKamIaJa0FFgP9AFWRsQ2STfk+5cDtwLvA+6VBNAcEbXASGBNvq0S+ElEPJBV7WZmlu01DCKiDqhr0ba84PkiYFGReTuBKS3bzcwsO/5LbzMzS+LAMDOzJA4MMzNL4sAwM7MkDgwzM0viwDAzsyQODDMzS+LAMDOzJA4MMzNL4sAwM7MkDgwzM0viwDAzsyQOjDZcd911jBgxgkmTJpW7FOuBdu/ezaxZs5gwYQITJ07k7rvvLndJ1oMcPnyY6dOnM2XKFCZOnMhtt93W6X06MNrw+c9/ngce8LeoW2lUVlbyzW9+kxdeeIGNGzeybNkynn/++XKXZT1Ev379ePjhh3nmmWdoaGjggQceYOPGjZ3apwOjDRdeeCHDhg0rdxnWQ40aNYpp06YBMHjwYCZMmMBLL71U5qqsp5DEKaecAsBbb73FW2+9Rf6eQu+ZA8OsG9i1axdPP/00M2bMKHcp1oMcO3aMqVOnMmLECGbPnt3p9ZVpYEiaK2m7pB2SbinSL0n35PuflTQtdW5XGnbkCHc1NMArr5TyZayXGnn8OL85cODt9XXgwAGuvPJK7rrrLoYMGVLm6uxkV/j7q0+fPjQ0NNDY2MhTTz3F1q2du7N1ZoEhqQ+wDJgHnAMslHROi2HzgPH5x2Lgux2Y22WuefFFJr/xBtxxR6lewnqxLx85wvnHjsEdd/DWW29x5ZVXcvXVV3PFFVeUuzTrAYr9/ho6dCgzZ87s9DVZRURn60t7IenDwO0R8Yn89lcAIuLrBWO+BzwSEavy29uBmUBNe3OLqa2tjfr6+vQiBwyAw4f/rmkXcDFQfdFF6fsxK2L9o4/S7/jxt7cD+BxwKvCc15d1Usv11QT0BYb278+br73GnDlzuPnmm7nkkkva3I+kzRFRW6wvy1NSo4HdBduN+baUMSlzAZC0WFK9pPqmpqaOVbhzJ1x1FQfzm58Gzquo4AWJJ598kpdffrlj+zMrsHD6dB4aMeLt9fWwxL8APxo4kPr6eurr69m3b185S7ST2In1dbgi92v95f79mXXaaUweN44PfehDzJ49u92waE9lVxSaqNjl+ZaHN62NSZmba4xYAayA3BFGRwpk1CgYMoRBFRVQVcVPjx6FL3wB7r23Q7sxa9WSJbBiBVRV8fGjR4kbbvD6sq5zYn3178/ko0d5+tpru3R9ZXmE0QiMLdgeA+xJHJMyt2u8+irccANs3Jj76Qvf1pW8vqyUSry+sryGUQn8Efg48BKwCbgqIrYVjLkYWAp8EpgB3BMR01PmFtPhaxhmZr1cW9cwMjslFRHNkpYC64E+wMqI2Cbphnz/cqCOXFjsAA4B17Y1N6vazcwswyOMcvARhplZx3SXT0mZmdlJzIFhZmZJHBhmZpbEgWFmZkl69EVvSU3Ai+9x+nDgL11YTldxXR3jujrGdXVMT6zrrIioLtbRowOjMyTVt/ZJgXJyXR3jujrGdXVMb6vLp6TMzCyJA8PMzJI4MFq3otwFtMJ1dYzr6hjX1TG9qi5fwzAzsyQ+wjAzsyQODDMzS9LrAkPSXEnbJe2QdEuRfkm6J9//rKRpqXNLXNfV+XqelfSEpCkFfbskPSepQVKXfttiQl0zJb2Rf+0GSbemzi1xXV8uqGmrpGOShuX7Svl+rZS0V9LWVvrLtb7aq6tc66u9usq1vtqrq1zra6yk30l6QdI2STcWGVO6NRYRveZB7qvR/w34D0AV8AxwTosxnwTWkbvL3/nA71PnlriujwCn5Z/PO1FXfnsXMLxM79dM4NfvZW4p62ox/lLg4VK/X/l9XwhMA7a20p/5+kqsK/P1lVhX5usrpa4yrq9RwLT888Hk7hOU2e+w3naEMR3YERE7I+IosBq4rMWYy4AfRc5GYKikUYlzS1ZXRDwREX/Nb24kd9fBUuvMf3NZ368WFgKruui12xQRG4DX2hhSjvXVbl1lWl8p71dryvp+tZDl+no5Irbkn/8NeAEY3WJYydZYbwuM0cDugu1G3v1mtzYmZW4p6yp0Pbl/QZwQwIOSNkta3EU1daSuD0t6RtI6SRM7OLeUdSFpIDAX+HlBc6nerxTlWF8dldX6SpX1+kpWzvUlqQY4F/h9i66SrbHM7rjXTahIW8vPFbc2JmXue5W8b0mzyP0P/dGC5gsiYo+kEcBDkv6Q/xdSFnVtIffdMwckfRL4JTA+cW4p6zrhUuDxiCj812Kp3q8U5VhfyTJeXynKsb46oizrS9Ip5ELqpojY37K7yJQuWWO97QijERhbsD0G2JM4JmVuKetC0mTgB8BlEbHvRHtE7Mn/3AusIXfomUldEbE/Ig7kn9cBfSUNT5lbyroKLKDF6YISvl8pyrG+kpRhfbWrTOurIzJfX5L6kguLf42IXxQZUro1VooLM931Qe6Iaicwjncu+kxsMeZi/v6C0VOpc0tc15nk7nX+kRbtg4DBBc+fAOZmWNfpvPMHoNOBf8+/d2V9v/LjTiV3HnpQFu9XwWvU0PpF3MzXV2Jdma+vxLoyX18pdZVrfeX/238E3NXGmJKtsV51SioimiUtBdaT+8TAyojYJumGfP9yoI7cpwx2AIeAa9uam2FdtwLvA+6VBNAcuW+jHAmsybdVAj+JiAcyrGs+sERSM/AmsCByq7Pc7xfA5cCDEXGwYHrJ3i8ASavIfbJnuKRG4Dagb0Fdma+vxLoyX1+JdWW+vhLrgjKsL+AC4LPAc5Ia8m1fJRf4JV9j/moQMzNL0tuuYZiZ2XvkwDAzsyQODDMzS+LAMDOzJA4MMzNL4sAwM7MkDgyzDEiqkLRB0toW7QPzXzf93XLVZpbKgWGWgYg4Dnwe+E+Srivo+l/k/sDrS+Woy6wj/Id7ZhnK/zX6N4APAv+R3F/dzoyIx8pamFkCB4ZZxiStBwaQ+66i1RHxP8pbkVkaB4ZZxiSNI3fns38DJkXEkTKXZJbE1zDMsncduS/SG0PudplmJwUfYZhlSNKHyH3l9X8BlpD7dtOPRMSxshZmlsBHGGYZkdSf3L0M7ouIdcBiche+fQ3DTgo+wjDLiKRvAf8ATI6Iv+XbFgD/DJwXEVvLWJ5ZuxwYZhmQdCHwMPCfI+KRFn0/JXct4/yIaC5DeWZJHBhmZpbE1zDMzCyJA8PMzJI4MMzMLIkDw8zMkjgwzMwsiQPDzMySODDMzCyJA8PMzJL8f+sLS1u8EIIwAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "mynewmesh.plotmesh(withnode=True,withnodeid=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "id": "dce8a39f-41dd-423e-8fcc-d93f2b072eb2",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "*** node coordinates of the mesh (total nodes=25, nodes per element=9)\n",
      "     1-th node: x=  0.000000e+00,y=  0.000000e+00\n",
      "     2-th node: x=  2.500000e-01,y=  0.000000e+00\n",
      "     3-th node: x=  5.000000e-01,y=  0.000000e+00\n",
      "     4-th node: x=  7.500000e-01,y=  0.000000e+00\n",
      "     5-th node: x=  1.000000e+00,y=  0.000000e+00\n",
      "     6-th node: x=  0.000000e+00,y=  2.500000e-01\n",
      "     7-th node: x=  2.500000e-01,y=  2.500000e-01\n",
      "     8-th node: x=  5.000000e-01,y=  2.500000e-01\n",
      "     9-th node: x=  7.500000e-01,y=  2.500000e-01\n",
      "    10-th node: x=  1.000000e+00,y=  2.500000e-01\n",
      "    11-th node: x=  0.000000e+00,y=  5.000000e-01\n",
      "    12-th node: x=  2.500000e-01,y=  5.000000e-01\n",
      "    13-th node: x=  5.000000e-01,y=  5.000000e-01\n",
      "    14-th node: x=  7.500000e-01,y=  5.000000e-01\n",
      "    15-th node: x=  1.000000e+00,y=  5.000000e-01\n",
      "    16-th node: x=  0.000000e+00,y=  7.500000e-01\n",
      "    17-th node: x=  2.500000e-01,y=  7.500000e-01\n",
      "    18-th node: x=  5.000000e-01,y=  7.500000e-01\n",
      "    19-th node: x=  7.500000e-01,y=  7.500000e-01\n",
      "    20-th node: x=  1.000000e+00,y=  7.500000e-01\n",
      "    21-th node: x=  0.000000e+00,y=  1.000000e+00\n",
      "    22-th node: x=  2.500000e-01,y=  1.000000e+00\n",
      "    23-th node: x=  5.000000e-01,y=  1.000000e+00\n",
      "    24-th node: x=  7.500000e-01,y=  1.000000e+00\n",
      "    25-th node: x=  1.000000e+00,y=  1.000000e+00\n",
      "*** element connectivity information(bulk elements=4)\n",
      "     1-th element :    0     2    12    10     1     7    11     5     6 \n",
      "     2-th element :    2     4    14    12     3     9    13     7     8 \n",
      "     3-th element :   10    12    22    20    11    17    21    15    16 \n",
      "     4-th element :   12    14    24    22    13    19    23    17    18 \n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEKCAYAAAAW8vJGAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAlYklEQVR4nO3dfXRUhZ3/8fc3iSEgIKKEDY9BGyQmJCGJQPdXAXWDYB+s1e4KVFF0aeh2T7ddQX7tqq3WarUewEqxVmndXz1y6mqVVdR1UQrWIgYNFFQUQU1ICKA8BnlI8v39MQPNDQlMJJk7IZ/XOXPMfZg73/vlej9z7507Y+6OiIjIEUlhFyAiIolFwSAiIgEKBhERCVAwiIhIgIJBREQCUsIu4GSdffbZnpmZGXYZIiIdyurVq3e4e5/mpnX4YMjMzKSsrCzsMkREOhQz+6ilaTqVJCIiAZ0+GCoqKrjooovIzs4mJyeHefPmAfDEE0+Qk5NDUlKSjkgaaalfM2fOZNiwYeTl5XHFFVewa9eucAtNEC3165ZbbiEvL4+CggLGjx9PVVVVyJUmhpb6dcQvfvELzIwdO3aEVGHiaalnP/7xj+nfvz8FBQUUFBSwZMmS2Bfq7h36UVRU5CejqqrKV69e7e7ue/bs8aysLF+/fr2//fbb/u677/rYsWP9jTfeOKnXOJW01K8XX3zRDx8+7O7us2bN8lmzZoVZZsJoqV+7d+8+Os+8efP829/+dlglJpSW+uXu/vHHH/v48eN90KBBvn379jDLTCgt9ey2227ze++9t8XnAWXewn610x8xZGRkUFhYCECPHj3Izs5my5YtZGdnc95554VcXeJpqV/jx48nJSVyyWr06NFUVlaGWWbCaKlfPXv2PDpPbW0tZhZWiQmlpX4BfP/73+eee+5Rr5o4Xs8+r04fDI19+OGHvPXWW4waNSrsUjqElvq1cOFCJk6cGFJViatpv370ox8xcOBAHnvsMW6//faQq0s8jfu1ePFi+vfvT35+fthlJbSm29gDDzxAXl4e06ZNY+fOnbEvqKVDibZ+AAuBbcC6FqYbcD+wEVgLFMay3M99Kqmqyn3MGPfqand337t3rxcWFvqTTz4ZmE2nkhpp1LOW+vXTn/7Uv/71r3tDQ0NIRSaQGPrl7v6zn/3Mb7311hAKTDAt9Ku2ttZHjhzpu3btcnf3wYMH61SS+wn3YVu3bvW6ujqvr6/3H/7wh3799dcHns5xTiXFMxjGAIXHCYbLgOejATEaeD2W5X7uYJgxwz0pyX3GDD906JCPHz/e77vvvmNmUzA0Eu3ZoenTm+3X7373Ox89erTX1taGVGCCOUG/jvjwww89JycnzsUloBb6tXbtWu/Tp48PHjzYBw8e7MnJyT5w4ECvju4QO60Y92Hu7ps3bz5mGzteMFhkenyYWSbwrLvnNjPt18Ayd388OrwBGOfu1cdbZnFxsbfqU0Ndu8KBA0cHHZgKnAH8dezYY2YvLy/n3HPPpUePHrG/xinmxRUr6NLQAPytX72BnyclcemFFwLw6aef8sEHH5Cfn09qampotSaCWPq1f/9+unXrBsCWLVvYtWsXOTk54RQcslj61djKlSspKiritNNOi2udiaJxv6BRz5KTmVtXd3R8dXU1GRkZAMyZM4fXX3+dRYsWHZ1uZqvdvbi510ikYHgWuNvdX40OLwVudvdj9vpmNh2YDjBo0KCijz5q8T6NY1VXw003cWDRItIaGlhqxj+407NbN+qjF7WGDBmCu/P+++9z+PBhUlJS6N69O3l5ea1f6VNA74MHmbFpExfu2EFZQwMXAkOSk/m0SxfqzBgyZAgbN27E3Y9egO7ZsydDhw4Nt/CQHOnX32/bxlvQbL+2bt3K/v37MTO6dOnC0KFD6dKlS9ilhyKW7euss846On9nD4Yj/fo/27bRDXi1SxcuPHiQ4dnZJEXflP3sZz/j8ccfp7y8HDMjMzOTX//610eDAo4fDHH9aCmQScunkp4DvtRoeClQdKJlfq5TSaWlXg9+ICnp6KGYnEBpaaRXaWnqWSxKS70OfD+oX7HQ9tU6bbAPo4N8XLUSGNhoeADQPnf91NSwuF8/vjNiBJSWwtat7fIyp5SamkivVq5Uz2JRU8NvU1Mp6d5d/YqFtq/Waed9WCKdSvoy8F0iF6FHAfe7+8gTLbPV1xiixo0bB8CyZcta/VyRWPTq1QtAd4FLuzjZfdjxTiXF7Uv0zOxxYBxwtplVArcBpwG4+4PAEiKhsBHYD1wfr9pERORv4hYM7j7pBNMd+Jc4lSMiIi1IpGsMIiKSABQMIiISoGAQEZEABYOIiAQoGEREJEDBICIiAQoGEREJUDCIiEiAgkFERAIUDCIiEqBgEBGRAAWDiIgEKBhERCRAwSAiIgEKBhERCVAwiIhIgIJBREQCFAwiIhKgYBARkQAFg4iIBCgYREQkQMEgIiIBCgYREQlQMIiISICCQUREAhQMIiISoGAQEZEABYOIiAQoGEREJEDBICIiAQoGEREJiGswmNkEM9tgZhvNbHYz088ws/82szVmtt7Mro9nfSIiEsdgMLNkYD4wETgfmGRm5zeZ7V+At909HxgH3GdmqfGqUURE4nvEMBLY6O6b3P0QsAi4vMk8DvQwMwO6A58CdXGsUUSk04tnMPQHKhoNV0bHNfYAkA1UAX8FvufuDU0XZGbTzazMzMq2b9/eXvWKiHRK8QwGa2acNxm+FCgH+gEFwANm1vOYJ7k/5O7F7l7cp0+ftq5TRKRTi2cwVAIDGw0PIHJk0Nj1wFMesRHYDAyLU30iIkJ8g+ENIMvMhkQvKF8NLG4yz8fAJQBm1hc4D9jU3oVNmzaN9PR0cnNzA+N/+ctfct5555GTk8OsWbPau4wOo7l+/dM//RMFBQUUFBSQmZlJQUFBeAUmmOb6VV5ezujRoykoKKC4uJhVq1aFWGFiaa5fa9as4Ytf/CLDhw/nq1/9Knv27AmxwsRSUVHBRRddRHZ2Njk5OcybNw+ATz/9lJKSErKysigpKWHnzp2xL9Td4/YALgPeAz4AfhQdVwqURv/uB/wPkesL64BvnWiZRUVF/nmMHTvWx44d6+7uf/rTn3z16tWek5NzdPrLL7/sl1xyiR84cMDd3Wtqaj7X65yKmutXYz/4wQ/8Jz/5SZyrSjxnnHGGn3HGGc32q6SkxJcsWeLu7s8999zRbVGa376Ki4t92bJl7u7+yCOP+H/8x3+EVV7COLIPq6qq8tWrV7u7+549ezwrK8vXr1/vM2fO9Lvuusvd3e+66y6fNWtW4PlAmbewX43rfQzuvsTdh7r7ue5+Z3Tcg+7+YPTvKncf7+7D3T3X3X8fj7rGjBlD7969A+MWLFjA7Nmz6dKlCwDp6enxKKVDaK5fR7g7f/jDH5g0aVKcq0pczfXLzI6+6929ezf9+vULo7SE1Fy/NmzYwJgxYwAoKSnhySefDKO0hJSRkUFhYSEAPXr0IDs7my1btvDMM88wdepUAKZOncrTTz8d8zJ153ML3nvvPVasWMGoUaMYO3Ysb7zxRtgldQgrVqygb9++ZGVlhV1KQps7dy4zZ85k4MCB3HTTTdx1111hl5TQcnNzWbw4cub5iSeeoKKi4gTP6Jw+/PBD3nrrLUaNGkVNTQ0ZGRlAJDy2bdsW83I6bTD0PniQueXlsHVrs9Pr6urYuXMnK1eu5N577+Uf//Efj5wO67yqq2Hs2BZ7BvD444/raCGqb0MDz+3b12y/FixYwJw5c6ioqGDOnDnccMMNIVSYYI6zfS1cuJD58+dTVFTE3r17SU3Vfa9N92H79u3jyiuvZO7cufTsecyHOVunpXNMHeXxea8x/DEjw+vBfcYMd3ffvHlz4JzmpZde6q+88srR4XPOOce3bdv2uV7rlDFjhntSkvuMGcf0y9398OHDnp6e7hUVFSEVmFh+k5rqddFtrGm/evbs6Q0NDe7u3tDQ4D169AirzMRxgu3riA0bNvgFF1wQ5+IST+N92KFDh3z8+PF+3333HZ0+dOhQr6qqcnf3qqoqHzp0aOD5HOcaQ+g79pN9tDoY0tIiq93ksTk1NbAhLliwwG+55RZ3j2yIAwYMOPo/cqfTTM82g+eYBWZ7/vnnfcyYMSEVmUBi6NewYcOOvvH43//9Xy8sLAyp2AQQQ7+OfPijvr7er7nmGn/kkUfCqjZ8TfrVAH4N+PeSkwOz3XTTTYGLzzNnzgxMVzA0VlXlPnmyf5aUFFn9bt386sGD/e/S0z0lJcX79+/vDz/8sB88eNCnTJniOTk5PmLECF+6dGnrXudUEu2Zd+vmDn51crL/XVpaoF/u7lOnTvUFCxaEXGwCiParNvo/bnP9WrFihRcWFnpeXp6PHDnSy8rKwq46PDFsX3PnzvWsrCzPysrym2++ufO+SXM/Zh+2oksXB3x4drbn5+d7fn6+P/fcc75jxw6/+OKL/Qtf+IJffPHF/sknnwQWo2BoqrTU68EPJCUdPXSVEygtjfQqLU09i0VpqdeB7wf1KxbavlqnDfZhxwuGznnxuaaGxf368Z0RI6C09LgXUyWqpibSq5Ur1bNY1NTw29RUSrp3V79ioe2rddp5H2aR4Oi4iouLvaysrNXPGzduHADLli1r24JEonr16gXArl27Qq1DTk0nuw8zs9XuXtzctM55xCAiIi1SMIiISICCQUREAhQMIiISoGAQEZEABYOIiAQoGEREJEDBICIiAQoGEREJUDCIiEiAgkFERAIUDCIiEqBgEBGRAAWDiIgEKBhERCRAwSAiIgEKBhERCVAwiIhIgIJBREQCFAwiIhKgYBARkQAFg4iIBMQ1GMxsgpltMLONZja7hXnGmVm5ma03sz/Fsz4REYGUeL2QmSUD84ESoBJ4w8wWu/vbjebpBfwKmODuH5tZerzqExGRiHgeMYwENrr7Jnc/BCwCLm8yz2TgKXf/GMDdt8WxPhERIb7B0B+oaDRcGR3X2FDgTDNbZmarzeza5hZkZtPNrMzMyrZv395O5YqIdE7xDAZrZpw3GU4BioAvA5cCt5jZ0GOe5P6Quxe7e3GfPn3avlIRkU4sbtcYiBwhDGw0PACoamaeHe5eC9Sa2XIgH3gvPiWKiEg8jxjeALLMbIiZpQJXA4ubzPMMcKGZpZhZN2AU8E4caxQR6fTidsTg7nVm9l3gRSAZWOju682sNDr9QXd/x8xeANYCDcDD7r4uXjWKiEh8TyXh7kuAJU3GPdhk+F7g3njWJSIif6M7n0VEJEDBICIiAQoGEREJUDCIiEiAgkFERAIUDCIiEqBgEBGRAAUDMG3aNNLT08nNzT067oknniAnJ4ekpCTKyspCrC7xNNevmTNnMmzYMPLy8rjiiivYtWtXeAUmmOb6dcstt5CXl0dBQQHjx4+nqqrpt8N0Xs3164hf/OIXmBk7duwIobLE1Fy/fvzjH9O/f38KCgooKChgyZIlx1nCsRQMwHXXXccLL7wQGJebm8tTTz3FmDFjQqoqcTXXr5KSEtatW8fatWsZOnQod911V0jVJZ7m+jVz5kzWrl1LeXk5X/nKV7j99ttDqi7xNNcvgIqKCl566SUGDRoUQlWJq6V+ff/736e8vJzy8nIuu+yyVi1TwQCMGTOG3r17B8ZlZ2dz3nnnhVRRYmuuX+PHjyclJXIj/ejRo6msrAyjtITUXL969ux59O/a2lrMmvvy4c6puX5BZEd3zz33qFdNtNSvk6FgkDa3cOFCJk6cGHYZCe9HP/oRAwcO5LHHHtMRwwksXryY/v37k5+fH3YpHcYDDzxAXl4e06ZNY+fOna167nGDwcxGn1RlCaz3wYPMLS+HrVvDLqXjqK6GsWOP27M777yTlJQUpkyZEsfCElPfhgae27evxX7deeedVFRUMGXKFB544IE4V5eAWti+9u/fz5133qnwbOJ4+7AZM2bwwQcfUF5eTkZGBv/+7//eqmWf6IhhuZndYWZx/bK9eLj2o4/I270btLHF7o474NVXW+zZo48+yrPPPstjjz2mw31g5sGDjK6vP+E2NnnyZJ588sk4VZXAWti+PvjgAzZv3kx+fj6ZmZlUVlZSWFjI1k7+pu54+7C+ffuSnJxMUlIS//zP/8yqVatatWxzb/ojao0mml0CPAJ8CnzL3d9uXentr7i42Fv1qaGuXeHAgWNGbwDOT0qiR48egfH79u0jLS3t6Pnzzmjr7t2kNRn3IZGf2dtyxhkAHD58mAMHDnD66aeTlNS5z1DG0q/6+nqSk5MBOHjwIHV1dZx++unxLDNhxNKvxvbs2UP37t077XbWUr++Ysa6hgYAqqurycjIAGDOnDm8/vrrLFq0KPAcM1vt7sXNvcZxO+vuS4HhwFtAmZn94HOsR2LZtAkmT6Y2Orgf+D9m5JjR0NDAnj17OHToEIcPH2bPnj3U19ezf/9+amtrj7fUU1p+jx784bTT2B8d/iYwAngbjvbrwIEDuDu1tbXs3buXzz77LLyCQ3akX0e2mJb6tXfvXvbu3UtdXR1du3YNr+CQxbJ9yd803b4mJSfzxbQ0NiQnM2DAAB555BFmzZrF8OHDycvL45VXXmHOnDmtexF3j+kBXAnUAfuAPY0fsS6jPR5FRUXeaqWlXg9+ICnJPSnJfcaM1i+jsyktjfQqLU09i0VpqdeB7wf1KxbavlqnDfZhQJm3sF+N6fyImRUDPwXeB34RDYiOq6aGxf368d8ZGTwyalTkopccX00NlJbC9Onw0EPq2YnU1PDb1FR+l5rKq9deq36diLav1mnnfdiJrjGkALcBNwO/Ama7+7En6EPU6msMUePGjQNg2bJlbVuQSFSvXr0AdBe4tIuT3Ycd7xrDiY4Y3gB6AxM9cr1BREROcSe6rL8OyFMoiIh0Hsc9YnD3a+JViIiIJIbO+UFgERFpkYJBREQCFAwiIhKgYBARkQAFg4iIBCgYREQkQMEgIiIBCgYREQlQMIiISICCQUREAuIaDGY2wcw2mNlGM5t9nPkuMLN6M7sqnvWJiEgcg8HMkoH5wETgfGCSmZ3fwnw/B16MV20iIvI38TxiGAlsdPdN7n4IWARc3sx8/wo8CWyLY20iIhIVz2DoD1Q0Gq6MjjvKzPoDVwAPHm9BZjbdzMrMrGz79u1tXqiISGcWz2CwZsY1/fm4ucDN7l5/vAW5+0PuXuzuxX369Gmr+kREhBP/gltbqgQGNhoeAFQ1macYWGRmAGcDl5lZnbs/HZcKRUQkrsHwBpBlZkOALcDVwOTGM7j7kCN/m9nvgGcVCiIi8RW3YHD3OjP7LpFPGyUDC919vZmVRqcf97qCiIjERzyPGHD3JcCSJuOaDQR3vy4eNYmISJDufBYRkQAFg4iIBCgYREQkQMEgIiIBCgYREQlQMIiISICCQUREAhQMIiISoGAQEZEABYOIiAQoGEREJEDBICIiAQoGEREJUDCIiEiAgkFERAIUDCIiEqBgAHbt2sVVV13FsGHDyM7O5i9/+UvYJSWsDRs2UFBQcPTRs2dP5s6dG3ZZCW3OnDnk5OSQm5vLpEmTOHDgQNglJbR58+aRm5tLTk6Otq0YTJs2jfT0dHJzc4+O+/TTTykpKSErK4uSkhJ27tzZqmUqGIDvfe97TJgwgXfffZc1a9aQnZ0ddkkJ67zzzqO8vJzy8nJWr15Nt27duOKKK8IuK2E1NDRw//33U1ZWxrp166ivr2fRokVhl5Ww1q1bx29+8xtWrVrFmjVrePbZZ3n//ffDLiuhXXfddbzwwguBcXfffTeXXHIJ77//Ppdccgl33313q5bZ6YNhz549LF++nBtuuAGA1NRUevXqFW5RHcTSpUs599xzGTx4cNilJLS6ujo+++wz6urq2L9/P/369Qu7pIT1zjvvMHr0aLp160ZKSgpjx47lj3/8Y9hlJbQxY8bQu3fvwLhnnnmGqVOnAjB16lSefvrpVi2z0wfDpk2b6NOnD9dffz0jRozgxhtvpLa2NuyyOoRFixYxadKksMtIaElJSdx0000MGjSIjIwMzjjjDMaPHx92WQkrNzeX5cuX88knn7B//36WLFlCRUVF2GV1ODU1NWRkZACQkZHBtm3bWvX8ThsMvQ8eZG55OXVbt/Lmm28yY8YM3nrrLU4//fRWH3Z1GtXVMHYsbN3KoUOHWLx4Md/85jfDriph9W1o4PG9e3nmiSfYvHkzVVVV1NbW8vvf/z7s0hJTdTXZpaXcXFpKSUkJEyZMID8/n5SUlLArS0hH9mFs3drmy+60wXDtRx+Rt3s3AxYtYsCAAYwaNQqAq666ijfffDPk6hLUHXfAq6/C7bfz/PPPU1hYSN++fcOuKmHNPHiQPQ0NDPnkE/r06cNpp53GN77xDV577bWwS0tM0e3rho0befPNN1m+fDm9e/cmKysr7MoS0pF9GLfffsy0vn37Ul1dDUB1dTXp6emtWnbnC4auXcGMr1dXkwT83aOPMvCjj9jQpQsQOW9+/vnnh1tjoon2jAULoKEBFizg8a9/nUmvvBJ2ZYkp2q8bDx0iE1j57rvsN8PT0li6dKk+3NBUk+1r24IFYMbHaWk89dRTOl3ZVJN9GAsWwJAh8PbbR2f52te+xqOPPgrAo48+yuWXX96ql+h8wbBpE0yezIGk6Kp368YvJ05kyrBh5OXlUV5ezg9/+MNwa0w00Z7RrRsA+7t25aXUVL6xfn3IhSWoaL/2A6OAq1JSKOzZk+GZmTQ0NDB9+vSwK0wsTbavK5OSOL9nT7567rnMnz+fM888M+QCE0yTfdik5GS+mJbGhuRkBgwYwCOPPMLs2bN56aWXyMrK4qWXXmL27NmteonOd/IuIwN69iS1oYGDSUl0OXCAgsxMypYsCbuyxBXtGQcOQFoa3Q4e5JNvfxuGDg27ssQU7VcX4DPgJw0N/GTKFPjVr8KuLDE12b5WHDoE6lfLmuzDHneH668/pl9Lly793C/R+Y4YAGpqWNyvH98ZMQJKS9vl4s0pp6Ym0quVK9WzWNTU8NvUVEq6d1e/YqHtq3XaeR9m7t6mC4y34uJiLysra/Xzxo0bB8CyZcvatiCRqCP3w+zatSvUOuTUdLL7MDNb7e7FzU3rnEcMIiLSIgWDiIgEKBhERCQgrsFgZhPMbIOZbTSzYz4/ZWZTzGxt9PGameXHsz4REYljMJhZMjAfmAicD0wys6Z3km0Gxrp7HnAH8FC86hMRkYh4HjGMBDa6+yZ3PwQsAgK347n7a+5+5IvDVwID4lifiIgQ32DoDzT+msTK6LiW3AA839wEM5tuZmVmVrZ9+/Y2LFFEROIZDNbMuGZvojCzi4gEw83NTXf3h9y92N2L+/Tp04YliohIPL8SoxIY2Gh4AFDVdCYzywMeBia6+ydxqk1ERKLiecTwBpBlZkPMLBW4GljceAYzGwQ8BVzj7u/FsTYREYmK2xGDu9eZ2XeBF4FkYKG7rzez0uj0B4FbgbOAX5kZQF1Lt2yLiEj7iOu3q7r7EmBJk3EPNvr7RuDGeNYkIiJBuvNZREQCFAwiIhKgYBARkQAFg4iIBCgYREQkQMEgIiIBCgYREQlQMIiISICCQUREAhQMIiISoGAQEZEABYOIiAQoGEREJEDBICIiAQoGEREJUDCIiEiAgkFERAIUDCIiEqBgEBGRAAWDiIgEKBhERCRAwSAiIgEKBhERCVAwiIhIgIJBREQCFAwiIhKgYBARkQAFAzBt2jTS09PJzc0Nu5SEV1FRwUUXXUR2djY5OTnMmzcv7JISmrszcuRI8vPzycnJ4bbbbgu7pA6hvr6eESNG8JWvfCXsUjqEzMxMhg8fTkFBAcXFxSe9PAUDcN111/HCCy+EXUaHkJKSwn333cc777zDypUrmT9/Pm+//XbYZSW0l19+mTVr1lBeXs4LL7zAypUrwy4p4c2bN4/s7Oywy+hQXnnlFcrLyykrKzvpZSkYgDFjxtC7d++wy+gQMjIyKCwsBKBHjx5kZ2ezZcuWkKtKXGZG9+7dATh8+DCHDx/GzEKuKrFVVlby3HPPceONN4ZdSqelYJDP7cMPP+Stt95i1KhRYZeS0Orr6ykoKCA9PZ2SkhL16wT+7d/+jXvuuYekJO2eYmVmjB8/nqKiIh566KGTXl5cO29mE8xsg5ltNLPZzUw3M7s/On2tmRW2Vy29Dx5kbnk5bN3aXi9x6qmuhrFjYetW9u3bx5VXXsncuXPp2bNn2JUlpL4NDTy3bx/J27dTXl5OZWUlq1atYt26dWGXlpiqq3k2N5f07t0pKioKu5qE13gf9uc//5k333yT559/nvnz57N8+fKTWnbcgsHMkoH5wETgfGCSmZ3fZLaJQFb0MR1Y0F71XPvRR+Tt3g23395eL3HqueMOePVVDt92G1deeSVTpkzhG9/4RthVJayZBw8yur7+6DbWq1cvxo0bp+tZLbnjDv68fj2L/+u/yMzM5Oqrr+bll1/mW9/6VtiVJaTG+7B+/foBkJ6ezhVXXMGqVatOatnm7m1R44lfyOyLwI/d/dLo8P8FcPe7Gs3za2CZuz8eHd4AjHP36paWW1xc7K262NK1Kxw4cMzoDWYUdu3KBRdcEPuyOokXV6ygS0MDAA5MBXoDP09K4tILLwyztITUuF/bgdOAXsAuMwb36MGgQYM466yzQqwwsTTuV2P/Y8akM89k+PDhIVSVuJr2qxZoAHqkpVG7YwclJSXceuutTJgw4bjLMbPV7t7sR5jieSqpP1DRaLgyOq6182Bm082szMzKtm/f3roqNm2CyZOpjQ4eSErioi5dyE9J4bPPPuMvf/kL1dUt5lCnNGnkSF5KT+dAUhJ/Bv4fsDg5mb5paZSVlfHJJ5+EXWJCOdKvWqAaGAeck5zMoLQ0zjzzTIVCE423L4j8P/lSejp36FNJzWrar5q0NL7Uqxf555zDyJEj+fKXv3zCUDghd4/LA/gm8HCj4WuAXzaZ5zngS42GlwJFx1tuUVGRt1ppqXtSkntaWuS/M2a0fhmdjXrWOupX66hfrdMG/QLKvIX9ajyPGCqBgY2GBwBVn2Oek1dTA6WlsHJl5L+6AH1i6lnrqF+to361Tjv3K57XGFKA94BLgC3AG8Bkd1/faJ4vA98FLgNGAfe7+8jjLbfV1xhEROS41xhS4lWEu9eZ2XeBF4FkYKG7rzez0uj0B4ElREJhI7AfuD5e9YmISETcggHA3ZcQ2fk3Hvdgo78d+Jd41iQiIkG6tVBERAIUDCIiEqBgEBGRAAWDiIgExO3jqu3FzLYDH33Op58N7GjDcjoCrXPnoHXuHE5mnQe7e5/mJnT4YDgZZlbW0ud4T1Va585B69w5tNc661SSiIgEKBhERCSgswfDyf/UUcejde4ctM6dQ7usc6e+xiAiIsfq7EcMIiLShIJBREQCOkUwmNkEM9tgZhvNbHYz083M7o9OX2tmhWHU2ZZiWOcp0XVda2avmVl+GHW2pROtc6P5LjCzejO7Kp71tYdY1tnMxplZuZmtN7M/xbvGthbDtn2Gmf23ma2JrnOH/pZmM1toZtvMbF0L09t+/9XSL/icKg8iX/H9AXAOkAqsAc5vMs9lwPOAAaOB18OuOw7r/PfAmdG/J3aGdW4038tEvuX3qrDrjsO/cy/gbWBQdDg97LrjsM4/BH4e/bsP8CmQGnbtJ7HOY4BCYF0L09t8/9UZjhhGAhvdfZO7HwIWAZc3medy4D89YiXQy8wy4l1oGzrhOrv7a+6+Mzq4ksiv5XVksfw7A/wr8CSwLZ7FtZNY1nky8JS7fwzg7h19vWNZZwd6mJkB3YkEQ118y2w77r6cyDq0pM33X50hGPoDFY2GK6PjWjtPR9La9bmByDuOjuyE62xm/YErgAc5NcTy7zwUONPMlpnZajO7Nm7VtY9Y1vkBIJvIzwL/FfieuzfEp7xQtPn+K64/1BMSa2Zc08/oxjJPRxLz+pjZRUSC4UvtWlH7i2Wd5wI3u3t95M1khxfLOqcARUR+Urcr8BczW+nu77V3ce0klnW+FCgHLgbOBV4ysxXuvqedawtLm++/OkMwVAIDGw0PIPJOorXzdCQxrY+Z5QEPAxPd/ZM41dZeYlnnYmBRNBTOBi4zszp3fzouFba9WLftHe5eC9Sa2XIgn8jvr3dEsazz9cDdHjkBv9HMNgPDgFXxKTHu2nz/1RlOJb0BZJnZEDNLBa4GFjeZZzFwbfTq/mhgt7tXx7vQNnTCdTazQcBTwDUd+N1jYydcZ3cf4u6Z7p4J/BfwnQ4cChDbtv0McKGZpZhZN2AU8E6c62xLsazzx0SOkDCzvsB5wKa4Vhlfbb7/OuWPGNy9zsy+C7xI5BMNC919vZmVRqc/SOQTKpcBG4H9RN5xdFgxrvOtwFnAr6LvoOu8A38zZYzrfEqJZZ3d/R0zewFYCzQAD7t7sx977Ahi/He+A/idmf2VyGmWm929w34dt5k9DowDzjazSuA24DRov/2XvhJDREQCOsOpJBERaQUFg4iIBCgYREQkQMEgIiIBCgYREQlQMIiISICCQaSNmFmSmS03s6Y3E3aLfk30grBqE2kNBYNIG4l+Udt1wMVmNq3RpJ8TuZn0pjDqEmkt3eAm0said+HeAwwHvkDkLt1x7v5qqIWJxEjBINIOzOxFIt9mmgkscvdZ4VYkEjsFg0g7MLMhRH5p7AMg190PhlySSMx0jUGkfUwDPiPyFcjnhFyLSKvoiEGkjZnZBcBrwNeAGUBf4O/dvT7UwkRipCMGkTZkZmnAfwK/c/fngelELkDrGoN0GDpiEGlDZjYH+DqQ5+57o+OuBh4FijrybyFI56FgEGkjZjYGeBn4B3df1mTaH4hcaxjt7nUhlCcSMwWDiIgE6BqDiIgEKBhERCRAwSAiIgEKBhERCVAwiIhIgIJBREQCFAwiIhKgYBARkYD/D97KSXV1v1oBAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "mynewmesh.setmeshtype('quad9') # set the current mesh to 2nd order\n",
    "mynewmesh.update()\n",
    "mynewmesh.printnodes()\n",
    "mynewmesh.printelements()\n",
    "mynewmesh.plotmesh(withnode=True,withnodeid=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "id": "e0e09b2b-afb7-466e-9e5a-591f6c1289ee",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEKCAYAAAAFJbKyAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAdqklEQVR4nO2df7BcZZ3mP09IIqJRZJJAINHIjmsJDJSQDWH8lV0YwegMurtUCSwyOlWpsGZLy8IdXKuUkrVcZ2tnpsCdiQyyi1MzWFSNiKWhMDiTQpaS4sIECIPBwIDExHCVMQyrxIV8948+kVM3ffuec2+/7+nb7/OpOnX7/Oj3+T7ndvpNn9PP/SoiMMYYY5qyoOsCjDHGzC88cRhjjGmFJw5jjDGt8MRhjDGmFZ44jDHGtGJh1wXkYOnSpbF69equyzDGmHnD/fff/9OIWNZvXxETx+rVq5mYmOi6DGOMmTdIemq6fb5UZYwxphWeOIwxxrTCE4cxxphWeOIwxhjTCk8cxhhjWjEyE4ekVZL+TtKjkh6R9LE+x0jStZJ2S3pI0plJi9q3D971LvjJT5LKjIxul9r2PP66XWrb83CJiJFYgBXAmdXjJcBjwClTjtkA3A4IWAfc22Tss846K2bFFVdELFjQ+5mTrnS71Lbn8dftUtueWwNMxDTvqYoR/bPqkm4DvhQR22rbvgxsj4ibq/VdwPqI2DdorDVr1kSrHMcrXwkvvHDE5oMLFnD+O97RfJyW3PG97/GKQ4ey63apbc/jr9ultj3XOPpo+OUvG48j6f6IWNNv38hcqqojaTXwVuDeKbtOAp6ure+ptvUbY6OkCUkTk5OT7Qp44gm45BJeWNA7PS8sWMC25cu5+Oyz243TkovXrmXb8uXZdeva/7dat+f0ul29vnL7rWvbcz7Pvzi84Zhj4NJL4R//cXgi030U6WoBXg3cD/zbPvu+Dby9tv5d4KyZxpzVpapNm+IliBcWLMj7MXPTpp7e0Ufn/3i7aVO8CPELsOcMup147spvpW3P+TzP9f2LAZeqRuoTh6RFwN8AfxURX+9zyB5gVW19JbA3STH79/PNE0/kP771rbBpU74bW/v39/S+//28upX2/1q8mN959avtOYNuJ5678ltp23M+zynfv0bmHockATcBz0bEx6c55r3AZno3yc8Gro2ItTON3foeR8X69esB2L59e+vnzleOPfZYAH7+8593WkdOSvNcml8o0/Nc378G3eMYpT9y+DbgMuBhSTuqbf8FeD1ARGwBttKbNHYDvwA+nL9MY4wpm5GZOCLibnpfsx10TAAfzVORMcaYfozUPQ5jjDGjjyeOARx38CB/umPHeCU+Z+D4Q4f49vPP23MOOvLcmV+w54wkff+a7utW47TMNjl+64oV8RLMu8TnXPiLxYvjRXvOQ0eeO/MbYc8Zmev7F/MxOT5MhpUcb5u8bE1Xul1q2/P463apbc+z1p13yfHOmZIcT5K8HKDLMcfk1a1pJ02bDtC153y62f3WtO05n+eU71+eOPqxYgW85jUsPnSIgwsW9Gbv17wGTjghiy4vvND730Eu3Zr2K4Bfgj1n0M3uuSu/NW17zuc55fuXJ47pcHLcnhPrOkVtzyl1i0iOp8TJ8eaUmLAtzXNpfqFMzymT4/7EYYwxphWeOIwxxrTCE4cxxphWeOIYgJPjGbHnbJSYoi7Rs5PjTo5no7QUdUR5nktMUZfo2cnxOTKs5PgLwAmvfe3wCpvCTw4c4Og+21Prdqltz+Ov26W2PddwcjwxVfLycJ/iXwC3LFrE6UuWJJU9Y8kSblm06Ncp11y6dW17Hl/PXfmta9tzfs9F9BxPsbjneHPtovpvV9pFeS60/3aJnkvqOX6jpGck7Zxm/3pJByTtqJbPJCvGyXF7TqzrFLU9p9QtJjku6Z3A88BXI+K0PvvXA1dGxPvajOvkeHNKTNiW5rk0v1Cm52KS4xFxF/Bs13UYY4yZnpGaOBpyjqQHJd0u6dTpDpK0UdKEpInJycmc9RljzFgz3yaOB4A3RMQZwHXAN6Y7MCKuj4g1EbFm2bJlueozxpixZ15NHBHxXEQ8Xz3eCiyStDSVnpPjGbHnbJSYoi7Rc1HJcWA1sHOafSfw8g39tcCPDq8PWpwcb05pKeqI8jyXmKIu0XMxyXFJNwPrgaXAfuCzwCKAiNgiaTNwBfAivYZen4iIe2Ya1z3HR1jbnsdft0tte5617nz6VtXFEbEiIhZFxMqI+EpEbImILdX+L0XEqRFxRkSsazJpzAr3HLfnDLruv23PKXXdczw37jluzxl03X/bnlPquud4Fzg5bs+JdZ2itueUusUkx1Ph5HhzSkzYlua5NL9QpudikuPGGGNGH08cxhhjWuGJwxhjTCs8cQzAyfGM2HM2SkxRl+i5qOR4isXJ8eaUlqKOKM9ziSnqEj0XkxxPhXuOj662PY+/bpfa9lxjXJPjI4N7jttzBt2ue1GX2H+7RM/uOZ7zUpV7jttzYl3337bnlLrF9BwfKZwct+fEuk5R23NKXSfH54iT480pMWFbmufS/EKZnp0cN8YYMzJ44jDGGNMKTxzGGGNaMVITh6QbJT0jaec0+yXpWkm7JT0k6cyU9Tg5nhF7zkaJKeoSPReTHAfeCZzJ9D3HNwC3AwLWAfc2GdfJ8eaUlqKOKM9ziSnqEj0XlRyXtBr4VkSc1mffl4HtEXFztb4LWB8R+waN6Z7jI6xtz+Ov26W2Pc9ad5y+VXUS8HRtfU+17QgkbZQ0IWlicnKynYp7jttzBl3337bnlLruOf4y6rOt70emiLg+ItZExJply5a1U3HPcXvOoOv+2/acUtc9x19mD7Cqtr4S2JtEyclxe06s6xS1PafULSo5PsM9jvcCm+ndJD8buDYi1s40ppPjzSkxYVua59L8QpmeUybHF862qBRIuhlYDyyVtAf4LLAIICK2AFvpTRq76f2xyQ93U6kxxpTLSE0cEXHxDPsD+GimcowxxvRhvt3jMMYY0zGeOAbg5HhG7DkbJaaoS/RcTHI81eLkeHNKS1FHlOe5xBR1iZ6LSo6nwD3HR1fbnsdft0tte65RcHI8D+45bs8ZdLvuRV1i/+0SPbvneM5LVe45bs+Jdd1/255T6rrneBc4OW7PiXWdorbnlLpFJcdT4OR4c0pM2JbmuTS/UKZn9xw3xhgzMnjiMMYY0wpPHANwADAj9pyNEsNwJXp2ALCLb1WFA4BZsedslBiGK9GzA4BzxK1jR1jbnsdft0tte561rm+Ot8WtY+05g67bqNpzSl23js2NW8facwZdt1G155S6xbSOlXSBpF2Sdku6qs/+9ZIOSNpRLZ9JVowDgPacWNdhOHtOqVtEAFDSUcBjwO/Q6y1+H3BxRPxD7Zj1wJUR8b42YzsA2JwSg1KleS7NL5TpuZQA4Fpgd0Q8ERG/Ar4GXNhxTcYYY6YwShPHScDTtfU91bapnCPpQUm3Szp1usEkbZQ0IWlicnJy2LUaY0yxjNLEoT7bpl5HewB4Q0ScAVwHfGO6wSLi+ohYExFrli1bNrwqjTGmcEZp4tgDrKqtrwT21g+IiOci4vnq8VZgkaSlqQpycjwj9pyNElPUJXouIjkOLASeAN4ILAYeBE6dcswJvHxDfy3wo8PrgxYnx5tTWoo6ojzPJaaoS/RcTHJc0gbgT4GjgBsj4vOSNgFExBZJm4ErgBfpfS37ExFxz0zjunXs6Grb8/jrdqltzzXGNTkeEVsj4l9GxL+IiM9X27ZExJbq8Zci4tSIOCMi1jWZNGaFW8facwbdrluKlthGtUTPbh2b848cunWsPSfWdRtVe06p69axXeDkuD0n1nWK2p5T6haRHE+Jk+PNKTFhW5rn0vxCmZ5LSY4bY4yZB3jiMMYY0wpPHMYYY1rhiWMATo5nxJ6zUWKKukTPRSTHUy5OjjentBR1RHmeS0xRl+i5mOR4KtxzfIS17Xn8dbvUtudZ6/pbVW1xz3F7zqDr/tv2nFLXPcdz457j9pxB1/237TmlbjE9x0cKJ8ftObGuU9T2nFLXyfE54uR4c0pM2JbmuTS/UKZnJ8eNMcaMDJ44jDHGtMIThzHGmFYMnDgkrctVSKV3gaRdknZLuqrPfkm6ttr/kKQzU9bj5HhG7DkbJaaoS/TcWXIc+BVwDbBw0HHDWOi1i30cOJmXe46fMuWYDcDtgIB1wL1NxnZyvDmlpagjyvNcYoq6RM+dJcclnQt8BXgW+A8R8Q/Dn7p+rXUOcHVEnF+tfwogIr5QO+bLwPaIuLla3wWsj4h9g8Z2z/HR1bbn8dftUtuea+RKjkfEd4HfAv4emJD0icaq7TkJeLq2vqfa1vYYACRtlDQhaWJycrJdJe45bs8ZdLvuRV1i/+0SPXfacxz4d8CLwPPAc/Wl6RgzjH8RcENt/TLguinHfBt4e239u8BZM43tnuPNtYvqv11pF+W50P7bJXruvOe4pDXAfwV+CHwM+E9TlmGwB1hVW18J7J3FMcPByXF7TqzrFLU9p9TtLDkuaSHwWeAPgT8DroqIIy/+D6OQntZjwLnAj4H7gEsi4pHaMe8FNtO7SX42cG1ErJ1pbCfHm1NiwrY0z6X5hTI9p0yOL5zhufcBxwHvid79jmRExIuSNgN30PuG1Y0R8YikTdX+LcBWepPGbnqXDD+csiZjjDFHMtPEsRPYHBEHchQTEVvpTQ71bVtqjwP4aI5ajDHG9GfgxBERl+UqxBhjzPzAf3JkAE6OZ8Ses1FiirpEz+45PsfFyfHmlJaijijPc4kp6hI9u+f4HHHP8RHWtufx1+1S255nret+HG1xz3F7zqDr/tv2nFLXPcdz457j9pxB1/237TmlrnuOd4GT4/acWNcpantOqeue43PEyfHmlJiwLc1zaX6hTM/uOW6MMWZk8MRhjDGmFZ44jDHGtMITxwCcHM+IPWejxBR1iZ6dHHdyPBulpagjyvNcYoq6RM9Ojs8R9xwfXW17Hn/dLrXtuYaT44lxz3F7zqDbdS/qEvtvl+i5057j83lxz/Hm2kX13660i/JcaP/tEj133nM8NZKOk7RN0g+rn6+b5rgnJT0saYek9om+Njg5bs+JdZ2itueUumOfHJf0R8CzEfHfJF0FvC4i/rDPcU8CayLip23Gd3K8OSUmbEvzXJpfKNNzCcnxC4Gbqsc3Ae/vrhRjjDGDGJWJ4/iI2AdQ/Vw+zXEBfEfS/ZI2DhpQ0kZJE5ImJicnh1yuMcaUy8Ce48NE0p1Av7/r++kWw7wtIvZKWg5sk/SDiLir34ERcT1wPfQuVbUu2BhjTF+yfeKIiPMi4rQ+y23AfkkrAKqfz0wzxt7q5zPArcDalDU7OZ4Re85GiSnqEj2nfP8alUtV3wQurx5fDtw29QBJr5K05PBj4N3AzpRFfeippzj9wAH43OdSyhzJNdfA3Xfn1wU+efAg6156yZ5z0JHnzvyCPWck5fvXqHyr6jeAW4DXAz8CLoqIZyWdCNwQERsknUzvUwb0LrH9dUR8vsn47jk+wtr2PP66XWrb86x1R/5bVRHxs4g4NyLeVP18ttq+NyI2VI+fiIgzquXUppPGrHDPcXvOoOv+2/acUtc9x3PjnuP2nEHX/bftOaWue453gZPj9pxY1ylqe06pO/bJ8dQ4Od6cEhO2pXkuzS+U6bmE5Lgxxph5gicOY4wxrfDEYYwxphWeOAbg5HhG7DkbJaaoS/TsnuNdNHIK9xzPij1no8T+2yV6ds/xOeKe46Orbc/jr9ultj3XGLfk+MjhnuP2nEG3617UJfbfLtGze47nvFTlnuP2nFjX/bftOaXu2PccH0mcHLfnxLpOUdtzSl0nx+eIk+PNKTFhW5rn0vxCmZ6dHDfGGDMyeOIwxhjTipGYOCRdJOkRSYck9f1oVB13gaRdknZLuip1XQ4AZsSes1FiGK5Ez2MfAATeArwZ2A6smeaYo4DHgZOBxcCDwClNxncAsDmlheEiyvNcYhiuRM/FBAAlbQeujIgj7mRLOge4OiLOr9Y/BRARX5hpXLeOHWFtex5/3S617XnWuuNyc/wk4Ona+p5qW18kbZQ0IWlicnKynZJbx9pzBl23UbXnlLpj0TpW0p2SdvZZLmw6RJ9t035ciojrI2JNRKxZtmxZu2LdOtaeM+i6jao9p9Qdi9axEXFeRJzWZ7mt4RB7gFW19ZXA3uFXWuEAoD0n1nUYzp5T6hYTAJzhHsdC4DHgXODHwH3AJRHxyEzjOgDYnBKDUqV5Ls0vlOl57AOAkj4gaQ9wDvBtSXdU20+UtBUgIl4ENgN3AI8CtzSZNIwxxgyXhV0XABARtwK39tm+F9hQW98KbM1YmjHGmCmMxCcOY4wx8wdPHANwcjwj9pyNElPUJXoe++R46sXJ8eaUlqKOKM9ziSnqEj0XkxxPhVvHjq62PY+/bpfa9lyj0OR4Ptw61p4z6HbdUrTENqolenbr2JyXqtw61p4T67qNqj2n1HXr2C5wctyeE+s6RW3PKXWLSY6nwsnx5pSYsC3Nc2l+oUzPY58cN8YYM3/wxGGMMaYVnjiMMca0whPHAJwcz4g9Z6PEFHWJnp0c7+LruOHkeFbsORslpqhL9Ozk+Bxxz/ER1rbn8dftUtueZ63rb1W1xT3H7TmDrvtv23NK3bHoOT6vcM9xe86g6/7b9pxSdyx6jg9C0kWSHpF0SFLfj0bVcU9KeljSDkntE31tcHLcnhPrOkVtzyl1xz45LuktwCHgy0zTc7w67klgTUT8tM34To43p8SEbWmeS/MLZXpOmRwfldaxjwJI6roUY4wxMzASl6paEMB3JN0vaeOgAyVtlDQhaWJycjJTecYYM/5k+8Qh6U6g392ZT0fEbQ2HeVtE7JW0HNgm6QcRcVe/AyPieuB66F2qmlXRxhhjjiDbJ46IOC8iTuuzNJ00iIi91c9ngFuBtanqBSfHs2LP2SgxRV2i55TvX/PmUpWkV0lacvgx8G5gZ0rNDz31FKcfOACf+1xKmSO55hq4++78usAnDx5k3Usv2XMOOvLcmV+w54ykfP8alW9VfQC4DlgG/BzYERHnSzoRuCEiNkg6md6nDOhdYvvriPh8k/Hdc3x0te15/HW71LbnGuOWHI+IWyNiZUS8IiKOj4jzq+17I2JD9fiJiDijWk5tOmnMCvcct+cMul33oi6x/3aJnt1zPOcfOXTPcXtOrOv+2/acUtc9x7vAyXF7TqzrFLU9p9Qd++R4apwcb06JCdvSPJfmF8r07J7jxhhjRgZPHMYYY1rhicMYY0wrPHEMwMnxjNhzNkpMUZfo2T3Hu/g6brjneFbsORsl9t8u0bN7js8R9xwfYW17Hn/dLrXteda6/lZVW9xz3J4z6Lr/tj2n1HXP8dy457g9Z9B1/217Tqk79j3HRxInx+05sa5T1PacUtfJ8Tni5HhzSkzYlua5NL9Qpmcnx40xxowMnjiMMca0whOHMcaYVozExCHpv0v6gaSHJN0q6dhpjrtA0i5JuyVdlbouJ8czYs/ZKDFFXaLnsU+O0+sfvrB6/EXgi32OOQp4HDgZWAw8CJzSZHwnx5tTWoo6ojzPJaaoS/RcVHK86j/+7yPi0inbzwGujqqtrKRPAUTEF2Yac1jJ8YMLFnD+O97RfJyW3PG97/GKQ4ey63apbc/jr9ultj3XGPPk+EeA2/tsPwl4ura+p9rWF0kbJU1ImpicnGxXwZSe4y8sWMC25cu5+Oyz243TkovXrmXb8uW/Tnzm0q1r2/P4eu7Kb13bnvO/vuZ1z3HgTmBnn+XC2jGfBm6lypdMef5FwA219cuA65poz7bneFe9grvsv23PY+7Z59qeG8Io9ByPiPMi4rQ+y20Aki4H3gdcWhU9lT3Aqtr6SmBvsoI7THx2maK25zH37HNtz0NgJO5xSLoA+GPgXRHR97qSpIXAY8C5wI+B+4BLIuKRmcafbXLcGGNKZT7c4/gSsATYJmmHpC0Akk6UtBUgIl4ENgN3AI8CtzSZNIwxxgyXhV0XABARvznN9r3Ahtr6VmBrrrqMMcYcyah84jDGGDNP8MRhjDGmFZ44jDHGtMIThzHGmFaMxNdxUyNpEnhqlk9fCvx0iOUMC9fVDtfVDtfVjnGs6w0RsazfjiImjrkgaWK67zJ3ietqh+tqh+tqR2l1+VKVMcaYVnjiMMYY0wpPHDNzfdcFTIPraofraofrakdRdfkehzHGmFb4E4cxxphWeOIwxhjTCk8cgKSLJD0i6ZCkab+6JukCSbsk7ZZ0VW37cZK2Sfph9fN1Q6prxnElvbn6i8KHl+ckfbzad7WkH9f2bThCJFFd1XFPSnq40p5o+/wUdUlaJenvJD1a/c4/Vts3tPM13Wultl+Srq32PyTpzKbPnQsN6rq0quchSfdIOqO2r+/vM2Nt6yUdqP1+PtP0uYnr+mStpp2SXpJ0XLUvyTmTdKOkZyTtnGZ/2tfXdB2eSlqAtwBvBrYDa6Y55ijgceBkYDHwIHBKte+PgKuqx1cBXxxSXa3GrWr8Cb3gDsDVwJUJzlejuoAngaVz9TXMuoAVwJnV4yX0erwc/j0O5XwNeq3UjtlAr0WygHXAvU2fm7iu3wZeVz1+z+G6Bv0+M9a2HvjWbJ6bsq4px/8u8LepzxnwTuBMYOc0+5O+vvyJA4iIRyNi1wyHrQV2R8QTEfEr4GvAhdW+C4Gbqsc3Ae8fUmltxz0XeDwiZpuSb8pc/XZ2viJiX0Q8UD3+Z3q9XabtXT9LBr1W6rV+NXp8HzhW0oqGz01WV0TcExH/VK1+n16nzRzMxXen52wKFwM3D0l7WiLiLuDZAYckfX154mjOScDTtfU9vPyGc3xE7IPeGxOwfEiabcf9IEe+aDdXH1VvHNYloRZ1BfAdSfdL2jiL56eqCwBJq4G3AvfWNg/jfA16rcx0TJPnzpa2Y/8Bvf+1Hma632fO2s6R9KCk2yWd2vK5KetC0jHABcDf1DanPGeDSPr6GolGTjmQdCdwQp9dn46q7/lMQ/TZNufvMg+qq+U4i4HfAz5V2/znwDX06rwG+B/ARzLW9baI2CtpOb3ujj+o/qc0a4Z4vl5N7x/4xyPiuWrzrM/X1OH7bJv6WpnumCSvsxk0jzxQ+tf0Jo631zYP/ffZsrYH6F2Gfb66//QN4E0Nn5uyrsP8LvB/IqL+SSDlORtE0tdXMRNHRJw3xyH2AKtq6yuBvdXj/ZJWRMS+6uPgM8OoS1Kbcd8DPBAR+2tj//qxpL8AvpWzruh1cCQinpF0K72PyXfR8fmStIjepPFXEfH12tizPl9TGPRamemYxQ2eO1ua1IWk04EbgPdExM8Obx/w+8xSW22CJyK2SvozSUubPDdlXTWO+MSf+JwNIunry5eqmnMf8CZJb6z+d/9B4JvVvm8Cl1ePLweafIJpQptxj7i2Wr15HuYDQN9vYKSoS9KrJC05/Bh4d02/s/MlScBXgEcj4o+n7BvW+Rr0WqnX+qHq2y/rgAPV5bUmz50tM44t6fXA14HLIuKx2vZBv89ctZ1Q/f6QtJbe+9fPmjw3ZV1VPa8F3kXtNZfhnA0i7etr2Hf75+NC701iD3AQ2A/cUW0/EdhaO24DvW/hPE7vEtfh7b8BfBf4YfXzuCHV1XfcPnUdQ+8f0GunPP8vgYeBh6oXx4pcddH71saD1fLIqJwvepdeojonO6plw7DPV7/XCrAJ2FQ9FvA/q/0PU/s233SvsyGdo5nqugH4p9q5mZjp95mxts2V9oP0btz/9iics2r994GvTXlesnNG7z+J+4D/R++96w9yvr78J0eMMca0wpeqjDHGtMIThzHGmFZ44jDGGNMKTxzGGGNa4YnDGGNMKzxxGGOMaYUnDmMyImmBpLskTQ23HVP9qes/76o2Y5riicOYjETEIXphsX8jqf53sL5I708AXdlFXca0wQFAYzpA0iZ6/UN+C/hN4A5gfUTc3WlhxjTAE4cxHSHpDuCVwGp6f67iP3dbkTHN8MRhTEdIeiO9vxf0OHBaRBzsuCRjGuF7HMZ0x0eAX9L709Ynd1yLMY3xJw5jOkDSvwLuodd86wrgeHp/7fWlTgszpgH+xGFMZiQdDXwV+N8RcTuwkd4Nct/jMPMCf+IwJjOS/gR4P3B6RPxzte2DwE3AWRGRq9mPMbPCE4cxGZH0TuBvgfMiYvuUfbfQu9exLiJe7KA8YxrhicMYY0wrfI/DGGNMKzxxGGOMaYUnDmOMMa3wxGGMMaYVnjiMMca0whOHMcaYVnjiMMYY0wpPHMYYY1rx/wFV9UosOXycIAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "mynewmesh.setnx(5) # now, it is 5x2 mesh\n",
    "mynewmesh.setny(10)\n",
    "mynewmesh.setdomainsize(xmin=-1,xmax=1,ymin=-2,ymax=2)\n",
    "mynewmesh.update()\n",
    "mynewmesh.plotmesh(withnode=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "c3ac4b2b-027a-4e25-a044-6db5335a6dcf",
   "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.9.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}