{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "from pystencils.session import *"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Demo: Working with derivatives\n",
    "\n",
    "\n",
    "## Overview\n",
    "This notebook demonstrates how to formulate continuous differential operators in *pystencils* and automatically derive finite difference stencils from them.\n",
    "\n",
    "Instead of using the built-in derivatives in *sympy*, *pystencils* comes with its own derivative objects. They represent spatial derivatives of pystencils fields."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAD4AAAAYCAYAAACiNE5vAAAACXBIWXMAAA7EAAAOxAGVKw4bAAADhklEQVRYCd2Y600bQRCAjZUCrKSCUEKADqADp4RAB0H8gn8IOoBUwKMD6ICEDqADEB2Q71vfnvbWd4f3MLKSkcazj3nuzO4crL2+vo5yODo6mrC2Da6Dt8zvc55/fT7OAyBIAz4ADfYc3GTtCvpfwVqa8SrTB9D9NErmP5lvQb+n60PG6DhD7qWSPWYex0PUDZbJM26md3HGUk/BzE9Z/5Yulo6RP0HmM3gB7oJW10ogD/wBLy7zLFRzM/NeRw3WoA3+N3gLrgQ+pVYJ0Mx2wTMbW12bb62j20OzknwsPcSVBa2vjcBdwCmzsgHqnHDH2jXUF/7RhYGwg9xjFfRAFcsTyx+3m8q5vWgCR82UQfsonTOv9yJPH4XfgzRo9Vg1dovnUj3ILBXqwHHEoEdQnWwAa39Y8GHbY9x3HRpy6QQ5Pxj2oafp+qrG4XHDmfjCNtpY4pR3U7ickbJf9FsxQue9hseuEflm3Ev8RbctuYZ4xw34hU3LsA106Jr9cO+htiXhCfzCvOvAAhM/oQ126Wfdg/f+128I41Iboz4Z9k7BMzBc1TEDs2lgrdlg37speMdHzHXyCRrL9oZxuCbud4DdoA4q5UFW+zvQ2j7jYhsLylzBFw50jFF7qnA3I3O/ZvM+cUzB1EnH2+x7eF1gxruqSX3hUBPhITbelKli0NeJGTcToYQTw2HI3pSBGQ+fqswNzgzl2VM+lDO0DTZZ7DpYHUkPsthGoV8hUePKS7PaaFMoM5Bf4AbjGKhBt4FtKlZOYz9xai7j1V5+6MU2MFgiYwJ2wuOGA/Zne6t/hRlEhK+s5Y7FvZQadJfxUAnoqbOaCJrd1F6yNTfsszHHXC20yRjPenzVRzjm15nYB12HYNCxKnJ5vwu69CqX68znUV+fjRKZcNCx1KPyXsrhxPfATOVQlzJ89mSrR/B+X4TR/I/6DKiGRW3UAgwKZfT9sSjwytgxNLY4jTq2x6cZ94WdsBaCgrZmnHUPq+0Qe22oF9RGCr0yCaP2HooDx6CfnH60nIB+Dfni/wBT0Cl7uzR0hHQzG/vh0gh+ARvy78Jn1wmwgExknV09/wOzSjw8PNwGT4b4gNy0RA7+CXijTHHG47Eti5IpX3tLt5H1t/TD7xVr6xR9olZgaNsrD1wvCUJn/Muv8dC51wPr8He95nNi8PoZ7Cfr7C0qKZWP5rUUP8pGrvsv7qHl6SwSHmoAAAAASUVORK5CYII=\n",
      "text/latex": [
       "$\\displaystyle {\\partial_{0} {f}_{(0,0)}}$"
      ],
      "text/plain": [
       "D(f[0,0])"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "f = ps.fields(\"f: [2D]\")\n",
    "first_derivative_x = ps.fd.diff(f, 0)\n",
    "first_derivative_x"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This object is the derivative of the field $f$ with respect to the first spatial coordinate $x$. To get a finite difference approximation a discretization strategy is required:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAIIAAAAsCAYAAACpFWBjAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAFqUlEQVR4Ae2c7VHcMBCG7zIpgJAOoIMkdAAdQFIB0EGY/IJ/DOmAdMBAB5AKQugAOoChA/I+iuWzJdln6c4cp/POCH2tZOnd1erzGL+8vIza6OTk5Fz5zwXPqeI23FZsZfJyweddm8TUyTPlr8tdyB3IbcsNVCCQEz6tiqD+InyUAGW4lbuRG2iCQDb4jJumBmk7o/9a7oPCw3QwEb4J5YZPm0XYUY8fBiVwNGASzQofzyJI8Jg7OolFeJK7w1f6ofyVp1zx8RTBSlodZjtxJP+nTRv8CQK54ROcGtTJjaLLjYtD8exW+CYIzTGk+r/Psbq5VVXp90Lxie1QG55Bi6ACu/rIpfxx6GNKZ/pg/VACofCa0r7KHSr8OVTOTRMf21PoUe6j4kcmVvmjtHO5ztOSeGnHbzn8rrSnckyBnUi80fh0qthhKvoyV0xVZxDP9863bXRLgQcbqfpF43bk79l0hT8pbM8YOglAZVCmR/lm6pG/LXctx/qkSijkmZynJFUmGxbfs8KdFNGWSfCj8Emof6R+9IVpEM/g1KCG04imEcIo5rSxJDX6Tg6BBpWnZKwHqKdqUQijDHZaMtyK2/ROClb/RG+xKHxSWtEXpk14NinCFzX+T0MHEFYpwAae1mSVR9gI1lUcRjMgu2SUwU1cYLxXfFL6FYmph6enCJUKPYtQ5CGsWalpdD+p4vVA5SilO2UE2PpPeiV8UjoSg6mHp6cIaoEZkepwaNQzkhFWX4QShDqE8vHtt0CLxCel/yFMPTxDi0VG3lXDFxESlcxKTXVQvztd8K0+lS+2L9H4aFChPHaH1PY9dlyh/reVsXkxmHp4GkXQx9kOfZPPToD579TW7vg0EmHNRHRWjoYzyu+cytw42fClAuRUHx+dFR+Vp0+9Tm2RmHp4vitgQVvXVJkRsvygRSg6RCVRRL1y7ohA2eyWc6R8wlfyQwLnm/dRH50vc6/4pDRVOM2CqYdnVRG4aaTD5flAQwMZzTVlIC5H2R9yJqx49VQQ/gOlYXkMKcx2k0Mkzgjg5bv7JtP/02aOfe75p9C3ZHxSmiNM+sTUx5Nr6Bh3fHy8LXcWU8byqtyuDXf1VWZN7ror/6L51NZkfFLaHotpE57WInRWWmkquwnMUs0qTKtA/Jj+0E5kWlFGY+cj5mmV9Z2fik9KuxIxDeIZrQg0WA1AMKxwYxaOmDoWiJ1J/BxDcyQaWjd0rue1GdXeFHxSmhmFaSueKebIlsHM2HAfft/199Hmap1vrf1t7QnePqao5lBmuRFImhqWu8tD60MIjGUu2n/YECo1pGWHwDA1ZCfStA4NU0MabtmVGhQhO5GmdWhQhDTcsis1KEJ2Ik3rkLmGTiu6HKV0mrahltqHr1yxcxfP7zVq192Kc0p6KWd4FN9UeGUo611DoQQ83y7fAijMWTu3nbzE9u4+lPZXebfyl+Z+Yx7amvvU4F2wSMBYB+48GP0h4jURV84rRbkrAjee9xK+ezkWvEEVH/yQZyn+J+f7N3dFQKD2WVxIiq6CMIW08YfqyCIt6zVCk4Q08nn2xhVu7Sd9ipv1gfKYGnj9y4KRqWKWR6Uq/vYp+12DKwIJG8FWdxKGRelYB/Igns+ZNxDyWUuw1pj2hI9yS0u5Tw0hwSBYHsm6P/e364N9qwSVwu4UUsnKI7hSiiAB85tN1gCh0c36gN9w1s4XlIaVcNOUlBetjCJIwDx7W5dfnik4osQi1HYL4mUKwV04vNlFV0IRJFCe0W/KLy0BQi4EPZKP6Ufg7vkBh0rPyh8sgoBYapIQMe1b8u0xs+0PyvFURPhnFCPx1CyCkrAiv4o8fpeBsmRJWe8aCsGxOLxRuPY/HZTGz/vtgpHpwlUCBI6lsNMCFsUoBRm5UdaKIGFh6hnFjGyXquYegYeOnLEi/CYUpXEtilvfUsf/AagBNxLfuoCWAAAAAElFTkSuQmCC\n",
      "text/latex": [
       "$\\displaystyle \\frac{{f}_{(1,0)} - {f}_{(-1,0)}}{2 h}$"
      ],
      "text/plain": [
       "f_E - f_W\n",
       "─────────\n",
       "   2⋅h   "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "discretize_2nd_order = ps.fd.Discretization2ndOrder(dx=sp.symbols(\"h\"))\n",
    "discretize_2nd_order(first_derivative_x)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Strictly speaking, derivative objects act on *field accesses*, not *fields*, that why there is a $(0,0)$ index at the field:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAD4AAAAYCAYAAACiNE5vAAAACXBIWXMAAA7EAAAOxAGVKw4bAAADhklEQVRYCd2Y600bQRCAjZUCrKSCUEKADqADp4RAB0H8gn8IOoBUwKMD6ICEDqADEB2Q71vfnvbWd4f3MLKSkcazj3nuzO4crL2+vo5yODo6mrC2Da6Dt8zvc55/fT7OAyBIAz4ADfYc3GTtCvpfwVqa8SrTB9D9NErmP5lvQb+n60PG6DhD7qWSPWYex0PUDZbJM26md3HGUk/BzE9Z/5Yulo6RP0HmM3gB7oJW10ogD/wBLy7zLFRzM/NeRw3WoA3+N3gLrgQ+pVYJ0Mx2wTMbW12bb62j20OzknwsPcSVBa2vjcBdwCmzsgHqnHDH2jXUF/7RhYGwg9xjFfRAFcsTyx+3m8q5vWgCR82UQfsonTOv9yJPH4XfgzRo9Vg1dovnUj3ILBXqwHHEoEdQnWwAa39Y8GHbY9x3HRpy6QQ5Pxj2oafp+qrG4XHDmfjCNtpY4pR3U7ickbJf9FsxQue9hseuEflm3Ev8RbctuYZ4xw34hU3LsA106Jr9cO+htiXhCfzCvOvAAhM/oQ126Wfdg/f+128I41Iboz4Z9k7BMzBc1TEDs2lgrdlg37speMdHzHXyCRrL9oZxuCbud4DdoA4q5UFW+zvQ2j7jYhsLylzBFw50jFF7qnA3I3O/ZvM+cUzB1EnH2+x7eF1gxruqSX3hUBPhITbelKli0NeJGTcToYQTw2HI3pSBGQ+fqswNzgzl2VM+lDO0DTZZ7DpYHUkPsthGoV8hUePKS7PaaFMoM5Bf4AbjGKhBt4FtKlZOYz9xai7j1V5+6MU2MFgiYwJ2wuOGA/Zne6t/hRlEhK+s5Y7FvZQadJfxUAnoqbOaCJrd1F6yNTfsszHHXC20yRjPenzVRzjm15nYB12HYNCxKnJ5vwu69CqX68znUV+fjRKZcNCx1KPyXsrhxPfATOVQlzJ89mSrR/B+X4TR/I/6DKiGRW3UAgwKZfT9sSjwytgxNLY4jTq2x6cZ94WdsBaCgrZmnHUPq+0Qe22oF9RGCr0yCaP2HooDx6CfnH60nIB+Dfni/wBT0Cl7uzR0hHQzG/vh0gh+ARvy78Jn1wmwgExknV09/wOzSjw8PNwGT4b4gNy0RA7+CXijTHHG47Eti5IpX3tLt5H1t/TD7xVr6xR9olZgaNsrD1wvCUJn/Muv8dC51wPr8He95nNi8PoZ7Cfr7C0qKZWP5rUUP8pGrvsv7qHl6SwSHmoAAAAASUVORK5CYII=\n",
      "text/latex": [
       "$\\displaystyle {\\partial_{0} {f}_{(0,0)}}$"
      ],
      "text/plain": [
       "D(f[0,0])"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "first_derivative_x"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Sometimes it might be useful to specify derivatives at an offset e.g."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAO8AAAAyCAYAAAC00IxIAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAMbUlEQVR4Ae2d7ZXUNhSGZ/dQwIZUEOiAJBWwdECSCiAdkMMv+MeBDiAVBNJBkgr46CB0wLIdbN5Hq6vIsuSRPeP5WEvnaG1LV1fSK90PyfLsydXV1WpdeP78+V+i+U3XT+toW35DoCEwHQHJ2JlKv1P8SfeXQ5xO1gmvGMDoL13fDDFqefMgINxfi7MN4ot1AzpPKw6T603FRv26J8SRu++Hxvt0aFhUEAafdW2COwTUTHnC/aVY31b8Q/Gx4rliC0LgJmOjvuHh/qb4cWiwi8IrBkyWO7rCpIX9IMAYILgI8AfFvxVbuEbgRmMjuftT3fxTVwxoNmTdZhW4I2qkHrP9OVuyJc6KgHDHyrLX8I3uzW2etc5jYb4kbNTXfzUuL3Xteb8ly8ukYX3VBHd/M/qBqmbJ0gS3PwZLwgbPF+E9S2HoWV4R4Y5A/E1K3J7nR8Djz+TE8l4osv65UPqvui46LBUb9Rsv+ENvDrDbHMdnz559VXwZp7X7Lka7wENjcKX4pFSX8h4q3inlb5o+VPemvDctr7ZdDbVPebNiM6X9Q+1dx09lzxWvFDvj3XGbJdlPJOGY5xeLVvF77rzGgT0HQnaDSvl4R5e6hmWN7s9IVxzcoXRc/Z+hMsp7pchrqoMKatNobKZ0YAibEr+hMsqbjKfKMg/wwHj7EMKtcHd9g2vGDldbZyXA7PiR93wrjQMD1glKQ7k+0PUny9A99PYaqbc2Mrr4WlnmnehYQlW9cRAddf+jWNUG3x4OI/T6Gbc1uR+FTVK26lHtOQg8k8ZiUBkPlPQleUF4fYPRamFSQNDCXhD4UbUGq5q0AO3bsYgaOyb/J10fJrTFx5oyovlb0W2W6OomTJGhMjzN90M0W8gbhc2U+tSPg8Azabt5YXhdr8iL3eanesYVG6MF4dHC9hFA85fG4VxjZAO5/Zr7HKnLrHo/d/cph4TNlN5PwtMrRsqGjctYeBmgXU6KKR1fSpkf1NH3aWc1gHhGay1gWm7DZ9rB7vehhEPCZgomm+DJgQ0OTjEPri2vHtBmrFN4v9vCHhHwA8NY5Cwvg3ax4+ahLNxk2XG9veoOEJteGysSNsHzrefvPCFb85pbNNryClAmGuUZYNZIuUnn62yXCgRQpCvhmBsLsGbwdxl2rSyG+jYaG+FImc4ubaGCX0Vb2mcoFJmUPBlPtY9lLeOPJ/TGhJcHMkY1XvQILWXZQGGy/ay0p4pt00tgTAzgybnWXGB8EOBdBpTyqHkxY+NGY6O5iDGh3KGETfHkjLtTYqe+R6wjSKwOAoVJxCsLvvN1x/h05fzle13xzbcSxOu1IjuebtdzK0wPjIn69lDRMGMs+BihF0TDRGTwRwWV4/VCjfXJ8aU+ztfuJajds2IzpVN7xjPMgVMaog4Qx2pXdqcf+/IxBggwgDvtEGeMvRcPJtxtRSYzW+Tm3uv2RgX6iYA5q6pryfLSaRRlR4B5VoQHY+Lu9cyBGwvQM1YPLaGijJEOWTujmfO6ETZTGlaBzT7xdIpUbTzHbbaJMFa7Qv9WTC5jgHj2aQgaWmKTgMA+UkSA8Qxy60AlH31wE1S94LpuyQENrwvCwQnhjeINz7rvBOUzDnydFAvvYBkYiN6UyVjF3ql/w4eNsJlSt/o9iI3y94mnecj3YuEdNUDqQO8TpQioC93zMn1yEH+En8nDJhgK4qYK7moNlh0MRQsenErCwlaPmWjBcyyGpig6bdjlg9o9NM86TRHtJGw6TCof9oinjfndU7UVq0awxOunir/qAK5YvCY1zY41d1q7gk2JBHetfRKXQUeYY3nZHR2DMcKOEqwKosXr4Tje6HlRVcFMRGrvFGymtGYveEZjeBvLaxOgemDpqZjwThjhCic+dM/pHwadMGnQfXkEF0vBp3Bs5LRP4gRCHIQLG4U2dnFW9l601RbMM+gtibKMDzBxLDZTunAAeJ6d6DMjXCM2N6p/sUENd4c5dEXIOkFpfNXCZhWWYeyECbxUlp+1ZIK6c5who900BBaOgGTiqyD4HFveKkhU0HZ9SwfQzRq40yCiRzkQvih+q+fixoqj0h/R2CZacY0mGlx0DuNPsvBWV+kqvk8Um+IoAdTS94kAe0pnCK9b82qi1rrNCB87yuy45QKC5z4rFA2C/sWEQFfcan5GtmexE0ZY7pXosnUoHb647EFwdT9FSaBoflbES+goIz277y91DcsC0W09iP/6H87eeq2N4aEhoHlwMrJNbs1bXUYVMNmdcOYKKZ91KoETVwQE6r670x/lsxuI8K7bKS1+9qWytCH9nnW0khAfFIS1F565MOp71hyDdWlqx9hBW8ey5S8EgVP1ExO88kKxrtvOSovofYEQq4wri5Ai5AhFsI6+zKWuzrL659yF/KzVVToKwZSDlSUtuNjUr2esPG3IBuXRTtzitH2BPuJTEu5A224aAjtG4ALhRZgIJpjXT5m/msxMdKPvUCiPNSiWzA4ZlCY8ymJdXT+IpqQgEMpYUDdREqpmbXCKIEeldpT6mCNvaQ2BbSGA/FwivGwkEWonIta1sw7UJMZS/q5Y8zvPVFysS7xMGHuW1+ddqnwcSrxqlETMp3SPEumt0b3gftW1+jejShW09IbABATcbnMqDIN8NFnfKPLe1b1/jYi/U1rMK76PyLKudJyPIliJV7CuUSaCjVDWhEElUcPA09AP6u0E+qqIJ+LOJPPcIWgPDYH5EDgT6wt2m00YmOxVQROVg/NDh+dXTGw/oZn4qRVNn+N6sXIl3jQ6FZL02XhBW1zPGlHF1fDpkap/dxXtVFkvvyU4Jcz42+tBlkPg2fuPk8KR8cIgOBqw1X0LCQIeJ1Kd22wTvGddknJTHl+okO3orlQx97xGsjodTz1XffYlYmfpXCH/x/NCgHPtH1ISMZuhe/h22psQ/6g2lBRIQrqsR+ECdhyf5VWcvY5jTD7qOcwLUNEzngyKG6xzXhdkLfw/z7+cCjCb4J33nNtASbzZzeVgBt/icoqLzaxHGd7sFjv3kzzRZi2v0mlrTkgHlYTKwZs6pgTqy35xBV/l2Z7BFN43vQyYp/sj7pyA0rGyucCyqf0cUw6Z6zSb/27NSxLaDndl60ET3FymId5OeEXA1XarS/S44533xHrmQIX7YF+FECZcrlhJ0OHHouGHApxi0JU0JhYWAH7UHQ6U6N4C1qAzAS1DV3i2U1gRIMkt2P4rjNKjt1hWvK10HM0aN8ubABk9mvB+uuUTsWgGXES3m1sN4pgz0KbNO0pBPDrPccuVR/+2/j2r+DbBjYHu3yOE94RTaVmB5xIH5zYP0Me0S73nANNKGAXLi5uCJsS9LAF9EGCpfZO+2VQ5lNNYjW6K4iD6fmyNEOYlL8reKNiSzbrmxkjl2ARkAxUPClrWy0P7DiJZTAjz+NR3+YO/7s36joFeA4kby4CmmnuIDS5atWISLUcuj+571iEADiFPuCKMuH4dT8mPJXks3zj9xitJaBgzlOjig/AAN+a82xNwwqtENCAg/aJ4FMEPbHVbRT/GNYfvW5UZa6mr27NgQjaqeOOQLjnMcDxSXmplmbAt/L+0dfs2ZnkB5q2iAXgUQGmQUTizhDl5z9LgI2AqTDmTzoZjzp1mvYvFTV1prHGadgS9naWJtifglFssvGhE1rxHJcCzQNSYbh0BzSuWIbd1ZQLmAvOu4+mIFjeR+EeuwALT2AsIH+UE4RVQAIdE57TiAnFqXd4WAppbTDpOo4W5hWASqUNX3GLu0/e77G1cKn/xllcYoPwIYfkXhPc63W0MPPZg+qR2aQhMR0BzCbeXU2idDSqlIdAXnjM/iLASTcfyKokJ6yar8piXTtihXWAAv87/zr4VgyBw2OFjZ++pYgp2TNruGwJrEfDCxnKM13vB3fMFz5Vmm1a40qngQoZFNpcZyx2sDplLCeo3SwoUV2fJkVpe8EBo+f2mtsMHGi1sggBuMJMOC5rGeLORuZY7Lslc/EVzkfm4ZGNC33vfBJxcXfV/QklA8Y0q2nLJgAmCFhoC+0VAMojVRQmmR0yv/z9vpnmPlLb0NUYGlpbUENg5Aiw3+A8ZsafiGpFzm1ciZHfvhWK6TnGF2p+GQENgfgQkh+w/8e47+5VdVnhplgqwmcA2PQxaaAg0BHaIgOQOd5lNvfB6La2+KLwQ+oIwYFu/hYZAQ2AHCEjeeL2G0Qw/m5yrdlB4KSBGfKTPRwAwbKEh0BCYEQHJGTvvvyve131vnRtX/R+hQaovjWjuKAAAAABJRU5ErkJggg==\n",
      "text/latex": [
       "$\\displaystyle \\left( {\\partial_{0} {f}_{(0,1)}}, \\  \\frac{{f}_{(1,1)} - {f}_{(-1,1)}}{2 h}\\right)$"
      ],
      "text/plain": [
       "⎛           f_NE - f_NW⎞\n",
       "⎜D(f[0,1]), ───────────⎟\n",
       "⎝               2⋅h    ⎠"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "derivative_offset = ps.fd.diff(f[0, 1], 0)\n",
    "derivative_offset, discretize_2nd_order(derivative_offset)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Another example with second order derivatives:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMIAAAAYCAYAAACoRCJ4AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAFqklEQVR4Ae2b7XEUORCGB5cDcEEE+DI4IAOTgSGDwxkcxS/7H2UyACLgIwNDBNyRwZEBlDPwvY8YzUra0aw09LLaXXdVW99vt7rVkmZmfefm5qZL6eLi4kh1J+Jj8SeVv6Z9asrWeDWy97Gvtb2t8Vr0yUGqlCZNALwQs/jfiB+q7oPSWWSNN0uJPRpkbW9rvFZdcSc8ETRpToIXSp+HCqv8t8qPlD4J61flrfFy8iTntdqu+/aXKvt8bshO1lvb2xpvzOit+C49ETgJnvUGCPXmZDhV/Z9hZUHeGm9JpHS6VOVd8TvxMzEn2r6Stb2t8SK/tOS7NBD+k6bvpWC0o/Zl6moXmTVeZMi+wOInCAiGf8SfxFtHsvGJmLn8Clnb2xovnVszvjsMNZMj2Plz9EMNj3KNY/XWeKkMFo/qjsQ80BOoWxkE/byYBzybrO1tjRdOrDXfRYGAolKQKH0gZmFBX1T3Uemx+BsVNWSNl8h+rPI3yfC6Js37V7S2tzVe4JGmfBcFgiZ91S+sM6+w6sIjuyoQrPECnQhWDMmJ8ENyeKtFOuit8t6Rtb2t8XCIMJv03fCM0E8aRaPFpDLXDV/3b+nqssYL5Qr7jZg3WFwlXpMXex3DrnuT1/zZxLrUDio35T/p06TvXCBIOaKU3TV6bYphe/J31/e+Yiq1xhuTJRlc1aDsc4H68KbL9/vZ2/CvsHmtvHGSHlvlv8AnzfjOX40IgGspmPuCzGL6qHZ3F1fKK0vou/ieymkAWeM5Yckf9yo3p7PqWRw8PwzXOeVX6Z2IcDtsdozwXok5kYpPI/pLCJtOSrz16jJYX1U/9Q2n1t5sbE/FZ8LleTAla7wUvznfHcoQGMUt9FRbymr3TsOBlFlg35W+6ssnyl+JubN3Sk3xwMwQb7CGRR726XV4rHRYPMpP6h2O9/nCMR/U71KcbgYeJkrVbzRoVH+qjsdKnV2jQRMF9a+1N4vQ+5SxEVnjReCLQnO+O5BubidS+mWhZ5TDwexI/hhjh/T5rq8nGPwVxBovUiYo4NDcCYaOLnCD/qv0DroO2ZVjgvkvLaoBZb2ZKntLX3xJsI1uIqq3xhubfXO+O5BRMMj1mLZqY5di93A7q8osdhyeGpHxTK6zxgMzQw9VnwteAjMM1pV6pzJK5hqMQZbfZYPq9Wdr7F2ijTVeRmZzvuNEgNj1oyNbBmFhvxU/6I2jrAsC0pT42OZ3Etqs8SJ50scv7KUToW9LA5vgHaNU77BPzRgC0l0NQ4DfmC+1d6lK1niD3FZ9d4iGUo5XWryHd+/jB6277r7q0kUVNA9ZgmBYONZ4g5RFxp8+w66/aHLPOyzwEor0LhmgPmNjsBHBuREysHektzVeBL64OTTlOxcIKKrJ8/UYnqJcUBAE0XXJGi9Rit03pyu6pHqmZQ+3pLdvGMHwTWNjSgPPY5inhfYulmuNFwhu0nf+ahTomc/KOP55Ymz3W7qm5JF+ttTgqS/fBDixIO6Y71xu+Q86slgHqpHjB1WOwR7RRuBxKlKCFd452gbfVQVC76GXSocHQ02SPN8Y5i6EUjze4BxJjlvkSkdPBNUTkGOBOikHXDEyQpocE3REHr/UnE2SzQ8Hp370OBu7gYHN+y76x5xSgwULhg9qf4ifq272blaCpz58ByAIkMd7+2zgqY2fG/CxKOqjsl/oS3qrjeeOz+K/lB+CbGqM+jpSn1F5vr2lVLoStGdiNjDmzKvU4buQ8lVUgqc+zfuu4z/Udo3Pz89PxJdz5qVxpzXj1P9IfFUz5rZvfs3Jlhvx3ZyrUdWOsYnO2oF4I8FVZ+yKlFVJ/dklx95mZMeogVOGHfaWDCywKd/tZCDgDxmUxcn1yD1TFPqInzgUX/HUlyOfn1hEV7BCWbfdMhaQPX+/73b9mObqsq45rhN7XTpvE+467Zti/w9+uITdWpTEXgAAAABJRU5ErkJggg==\n",
      "text/latex": [
       "$\\displaystyle {\\partial_{0} {\\partial_{0} {f}_{(0,0)}}} + {\\partial_{1} {\\partial_{1} {f}_{(0,0)}}}$"
      ],
      "text/plain": [
       "D(D(f[0,0])) + D(D(f[0,0]))"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "laplacian = ps.fd.diff(f, 0, 0) + ps.fd.diff(f, 1, 1)\n",
    "laplacian"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeYAAAAsCAYAAACwlxupAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAK40lEQVR4Ae2d7XHVOBSGCbMFBLaD0AEfHYQOAlRA6ABmf8E/hu0AOmCgg7AVLKQD6AA2HWTf11g3tq5kS7Z0LdmvZnRlS0ey9Eg+x5I/7tH19fWtFO7NmzcnKOdVW9ZDhL+4j/jLNm4wgNx7CFy1Qm+xb7YH820lUXz8PS02fjaHTkFfSA9khK6x7oe7Jja3/c0MT2lPxvcIX7T+AXLTIH/D/ulYSZB5B5m78B/hz+FH80BmM058/F0tNn42h05BX9AoSw9kAq+x7ge7NjZJDDNw0bC+6GIDKM6eOev91I33bNMY0yjTOH+F/wJfpUO7T+HZnpRuFXzEJuWQKLIs6YG2WzTW/eNTbPxsTMofZmNmyBnudwC/A99dgqaBPUPcCfwP1zEQz7zH8F/avNUa5bZ9bAt9ErcyPmKTZFQUW4j0wE3XaKzfsLC3xMYmYu2nmjHTmP5oDat1iGZ3yFA9HsnrKm9LceLj722x8bNZIkV6IB91jXU/29WxSTJjhkF+4mF2n/FI33sADHFcniVQXmX/wj6XvBn2lsQRt0knPv5uFxs/myVT0C/SA4k7QGPdD3TNbJIYZhc6QKNR5sMg5kntnhjSPyDiA0I+Fs6nsP/uCWx8R3z8A0Bs/GxKS0FfSQ/M6BSNdT+8NbNJtZTtoscZ8GfA8xpcpNFw03nvK0OmuUf9WyzPL47xMk/J80otgY/YzOtDO3epPO16JtyXHpgJU3rAD7AENv7a+VPG9MAR32OGEO8B/wM/dC/YPsoT5NtboqYQ4vlOMh/44lK11yH9DImfEB65hBDP5W7eu+4Zbuzz6U+6n/B/Yt85K28k2p+xPEhvXvPo5hnapjzSuQxvOz5ZTvfrd9D7vUQ+33JfT5A7kI3mgzzsw6fwfHWNr62NOsgN8kT6KtiMgnAIlMDTUa0sUW1bq9ADY2PWBSgkD2RWMdZD2mozQp5B3bFlNjarkP05PJulbBRwhQMFKfGxCqEsGtO7CAeNclvOI4S+p7U5SB6jnJ4ha8v/ibCZiSPk60kX8N7jIY11GsvDC4R38KNGnnWHnPNeOOJpTHlR4l0pYP5AF8UHx+SyoblYIL9RhzybYDMKwiFQCk9H1bJEob1V6IHAMdtjFJFnDXog5Jy2+YTojk2y6YEK3AnUHV6eSZeyURkapXsId8YU2zRSZsnabhYHg3PWjXjO4jgrtR3jdzNolM1tGmffMZh/NE+nnCCDxkIP4KL4oA2ckfOCwHmx46nvJth42j4YvTGegyxiEsEttx4YHbOO+gblQd2NPqlWD6DtQW3tMgoZ61tl0+UUuj2XZzLDjIrQiDxCaM84eZK6lnXZxofw/3LD4WhsdwaY6din8eUJYxseXunz+HsuMk9zUu4VslxEFJ/YaopNLLFh+cp5DjcuMBUMsuqBSMZNrSfkqVYPTGhrYM/uxMRmhyLJhpNnkqey28HAhzz4kRB7lksDu7es2xlAezPmNo3G1na+q1ga/ru2cLsfk4cXCVwS/+wp62DRE/nE1k9sYokNy1fJc7hJ4antmM2tB2IYm8rH5qlZD8S21TAKDcUmlFSYnJNnEsOM41/AczbLexu22zO8rUAzw8XJzCsG27Es3yzbluU+jbJvQLrkfXmukMBjl+BS8oltj4un2MRSvJEvnedNTedtLakHXIzHWuPLs8ax7mvrGCM7XWxsIvP2nTyTGGYY13sT6jY0Mz1Geayw7VxxlKG8vbxt8sbkibkYMOXnCqfwia3LZthgjN4HHN57G3N8mt03lsby1spzrF1B6QfSAzGMTb1j89SsB2LbahiFhptjk1l3OHkmMcyhPYoG8n7zM4R8OIz3T9968lIxHttpyGc++8lZ7aWVbu83yZF5WO5UpWyqwxPDd3IYGWc4l4+z0IHILbFBWzk+vE/uD2AKTlqAZ3DdShKcM84jGTfNnpCnWj0woa2xQ2NzbMA0p+5w8rwd2ysz5TljOUZDG6OL0HkvtwXBCrscjfmpSYAst/khk8agImT59sxoMI8pCyGP+b2zH72JY/M+O79qNsWl4OM9rth40UxKWJrnpEqXkWnuOB88n2f2CwnVrgem8AkdGWITSipMzsnzdljeZFI8IXkfiuHulSpP6Zwds9I9hzg+SMaPivCdY36xi+U87wgxzznSODtvXEAeIzq0fGxkcoaz+KCdfDWNZfwF32xjv/tVs82ymdJpK+c5BUmqPHPHeU4dwDbWrgem8BnTHabvN8fGNDwmDNAdpjg3T375q0T/+vXrU/h3U+uGvGcxeSF/DH8Rk2dJWdR1Mh+xSTvmS+DZjofzJcdkjmMfcpyz/jie9IDHJohNcr3hHWuHnjGbq4TREFccfFqby9J7s+axzMjD5W3X095DWXkV/2JIoKS0qXzEJm0vFsSTt4eaW0RpW7hsaQce52ys9IC/y8XGz2ZKipdnsYaZrcRJSUPJp2RjFQ6XZa5CSUGWr3nx82hzH/wKPWQSOdR3Ch+xSUJ/V8gmeO5au8DGIcY5m4XjSA94+ldsPGAmRo/yzLH8lLpMLqGkLrNbXu7yu8fKsZ2z/jnLzsHCLrO0+ueqD8o9g39pt39N+7nYGUa5yzfHyRXmrH/OsnPx6JZbWv3H6tP8u9REo69sIiAChRDAFTgfduTsfe8re4VUUdUQAREIJFD0UnZgGyQmAiIgAiIgAqshcIQp9fVqWqOGiEDlBDDjdf43uWkW0rP8DzjKlR4wkBWKwMIEtJS9cAfo8CKQggAMq5ayU4BUGSJQAAEtZRfQCaqCCIiACIiACBgCMsyGhEIREAEREAERKIDAQf/EooD2rr4KWNLkO9/8JCfdye/g1nPEB7/X3eZRIAIiUCkB6YFKO66ttgxz3f3nqj2/Ib77ghm2+bDQN/gpf83pKl9xIiAC5ROQHii/j7w11FK2F005Cbz6hb+A/w9+7N+v+Ace/CSpcfzsG99v5f8Ry4mACFRKAOew9EClfRdbbRnmWGILyOOEvILnv5Dwk6Fj3wDnbPnrAtXUIZclwFsVul2xbB9kPbr0QFa8RRWupeyiumO0Mpz1vh2Swslr/xc0DTX/QpN/9i23UgLo37ELtpW2fJPNkh5YebdrxlxJB0PxmuXpYAWMPDyB+X7rg0qaqWqKgAgMEJAeGICzoiQZ5no6s1nKxokZtFwJuRM0jfeXH4TmqQeFaioCmyUgPbCBrteXvyrpZBhXPlnNe8cX8Hfh+ZQ1Z8T8W8ze31W2RvkVwubp7Hb/FsKeHPLKiYAIVEQA57D0QEX9NbWqusc8ldwB8+FkPMbhaITp+BpEY2ARfuI+/BMm0CGOM2W+IkXDbPLQQL+ClxMBEaiUAM5n6YFK+y622jLMscSWkTf3l/mhEHvWy5O163hFzTiGO4d8zex5F6ENERCB2ghID9TWYxPrK8M8EdyBs/G+0iWMq/1kNWfEn7t1gcyd7r62RUAEVkNAemA1XTncED38NcynlFReKfeexoYB5pI1/cdSKql6iIAIZCUgPZAVbzmFyzCX0xfOmsAAc1maBpgPfXUdl6b54RF7Ft2V0bYIiMAKCEgPrKATI5ogwxwBayHRpzwuTszejBlR5/DNx0SQxs9w0njLiYAIrJOA9MA6+9XZKhlmJ5aiInlfyTbKrCBn0mYZ+x4Ms/1QGGXkREAE1kFAemAd/RjUCj38FYRpUSEaYL4WZTu+/vQMBvkUXq9C2XS0LwLrIiA9sK7+HGzN/xb7P/aQN5W+AAAAAElFTkSuQmCC\n",
      "text/latex": [
       "$\\displaystyle \\frac{- 2 {f}_{(0,0)} + {f}_{(1,0)} + {f}_{(-1,0)}}{h^{2}} + \\frac{- 2 {f}_{(0,0)} + {f}_{(0,1)} + {f}_{(0,-1)}}{h^{2}}$"
      ],
      "text/plain": [
       "-2⋅f_C + f_E + f_W   -2⋅f_C + f_N + f_S\n",
       "────────────────── + ──────────────────\n",
       "         2                    2        \n",
       "        h                    h         "
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "discretize_2nd_order(laplacian)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Working with derivative terms\n",
    "\n",
    "No automatic simplifications are done on derivative terms i.e. linearity relations or product rule are not applied automatically."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQUAAAAXCAYAAAAC/1q1AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAIOElEQVR4Ae2c/1XcOBDHHR4FkFDBcR2QpIJbOsjRQaADePkL/uORDrirgCQdJKkgCR0cHYRsB9z3o/Uotiyv7cXyYrzznlbSaPRjvh6NxjLJs/v7+yyk8/PzHfFmSntKX1S/CWW61jXGntJt135Dy6fQfWgdxjLfBuv1P6nYvtwKlyUhnME7JRzBP0qvxPuofGVS/xN13l95gIE6ptB9oKWPbpoN1ukfmTCu7DnxdpQ47I04rNmfnp4VIwU6qOWd8lMvoULe6bXyv4v8NmX1eSM5+pbGbNN3SBmtr3fdw/Vrjivx5jn/QnUrh6JPuj4E1k8awJbKCedfEsWuLdKnDL0s2p7KOIW5coKALIwUiBCO1GidkYEQfiN+xfO41pqffJyKk6kRXze7V91DZYTFpXgvlK6VjpSIyKZKSbGeKqgRve/E4+Bh37KnPymVHILqmWzzvbLjfL9m2zAL9J/KH9TIQJ6o5zwM2byOb19SYCNwOiYlrY11EQY5T7fiZH3rHi4DR/BWCcfwXemL0uhoJFiPDtfYgnvA+kZjtI3u2afs1+OSU2jYVHid17HFL+EdaszjJe19NeEFSStTAt39WjQ2Tov1cWmLwx2lQ8gVetRY52t8KtmDsW4LBPav9EvptOQUGEBMTrSXShgv9E08wg4uJ1p/PVAf7hJay0t27dSX7hFFDsS71fiGaURkWqyEWE8LyH61Zb8elpyCHtRnMTFef7qrPMsfINN32eRshMYTUWMTsvxkcBGb5rt4XV5RXMeH/mjOPnV3y8lxAwcihTvV+YpD7vF1ghP7SYS1u9AWlGZLHGTY7kNeKUf/ZKQ/hzwRx64SBzsX3HX7i/16sKUfR/mDypSXDFZ1BI33Ixdvk72SEO/pUdK43AHQTiTynqQyEQpOYlDS3DiETLnp6eZXfVXdrT8hGe90PJQryuEcCIrHJS4PLAlp7NInpySTtBxUa+kda42JzbgLbZWxIxwCtlp5n1b7ZLCW/tgdd4TsL77+kX6ozCEVI/bjnosUJGS34WzKGDE49IEfydvGxSvvqh773Egf7iHqCOP4pL48QKNDFS6sMkSu+VPo7peu8W2zR6OmfH6iMx+FqdwGXz8HhWV91IZR4JRKTq80wAAVzd871hqTCAGn99xUEA88sT/ngAp85p8E1ugsDIhSPeW4YIdcKv7pG34X2K8Lp6ACm5ovDHVhBYbNBkYGYH8qxyNnymdKn5VKC1ATt+xzZEKSLItizJIDEN8/2LAP9bxfzMsxF+0xo2+6gU2hO8sx4nMQa6tgKx6Ge6Dcn2gqt8XXxs9a9vkouUulmAP3Y1lBcjyjMWD9r9bpbLOwdoc5/AJvalib6mHO4cOejf2FMW0722oELLfpw97U1W6GgZFAnGJ/uZJ+1M6NOk4hNomJhTkRgd3Eh221dc0R2/SZ+JwWzO8cVe0AQYPkh9CdLzY+CgiWAJaGqzWtgm9jH+kK3jiFHaW5TVaXS+bRY601svl5hteBHtgsB1gR98lgDRbSnSjphfKm6B9xIw7X+ZZ+KEDfFlnll5OF0xajwnnwEIpg0wEjM+9MHSIUQbZEGgMeqRTalYSGq6TSvagBuFSihFxgBq4m3BFf161jH+YyJ2/TDpWnwJp7KyjEN3bJPSWswQRsKvtPPPccZDchZvRB/m5LjWzwOZyQ1MYJjBFZeBubhG44AHvo1CHGxYmUSGPaXJaH7YMZbULdizrxcCoOV3ODTYhBF3xtji59WEf4mmfjJM1TYp2P7davMnhgQ/7QmRrW+YPkkjt2bwA2/iDKZS1jD98SKUBEA6VwUQPui8f7Gn8WGUYGYpeIwULjxBMROseIz0TmaHy75iHEw8EMSSl0d+uXPhZZxbwybW11jeHbhFGsz1ydmHdd1DfWdvFd1AmbhVzbouh0nhrWXCyXXk1Vt69Qlb2X48Srxs02FQnjVfx39FyA7A/xMSSjYtl45DiE0HHwnsd3+QppTP7OmvdbnIB9tmSMwf+RkNaQQnfTGceaaY6YZ0bfEM+wbuMgG+JrbV36tN0YNnaved9YazzuDYh8sCWioF0ldOQLQxEX8CvWVa3U4UHIPgWswQBczDFwSIBNuKfF8kQU8fvPnNWZm1p/W+vFCgXJGNh45ptCE8VSXbLcQ2RK+5QDWficGo+CtJZedS8ohcHWYYrhYYCetI7W+Fqnjn14bnUGb0MmzfvGWuPhcL3TVZ2/T/D1XJmpYo3epTeAHI9KJtxctAWeW5XWZsaFRPAojjQIZT4JxYyNSKDVohajrfzLKRCeBCsPtqRjo+7CgT+OsQiJ+wQipgpJBkfpHkTQuHQO9ePrAbgWaWmfgiDzWWRWYHcqPgqshQGHDZ9ZvWOFJ01Idjo6xcTfYN38iDmknV2V/j+F5n4LCYFsRskfL3GZcSpedFOKz4UPrwsxp7EYcES/TbqrnU2Hrry3fVW97pNQprYoNuLX4qs2jP6r0luVfRSyrI9kHUkmOp+1P7Z8mU5qw/FyILlwWHWcA7hcq1z5NC1eVHfxJ4+1MOCw4A7CXUJvq9KZ1LlL6M/m4AGu5da7s3INHVrojpFhoOR1Fzo2CzJEUiU8l82hNk6958r5MuRpWR+E1O5OVOWjcc4NOlk0cCQ57hLQj8MpfHVAfWiD9QKH2C9YYocL4n9eSp3Ozs72lE5SzzPG8YXLFfh0WbvkZ0o7Hft0nqfL+GOQ3WBd3evC5ESpZH/JHcIYjGXda9RDuVRqvckle9Rlzcgrzbr0eaqywmGDdUMgsHEKDQANtTlkrK2dQtc1pRy761oeg3xKPFKOPRR2/wMdsuuwbt1IoQAAAABJRU5ErkJggg==\n",
      "text/latex": [
       "$\\displaystyle {\\partial_{0} (c + {\\partial_{0} {f}_{(0,0)}} + {\\partial_{0} {g}_{(0,0)}} + 5) }$"
      ],
      "text/plain": [
       "D(c + Diff(f_C, 0, -1) + Diff(g_C, 0, -1) + 5)"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "f, g = ps.fields(\"f, g :[2D]\")\n",
    "c = sp.symbols(\"c\")\n",
    "δ = ps.fd.diff\n",
    "\n",
    "expr = δ( δ(f, 0) +  δ(g, 0) + c + 5 , 0)\n",
    "expr"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This nested term can not be discretized automatically. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Only derivatives with field or field accesses as arguments can be discretized\n"
     ]
    }
   ],
   "source": [
    "try:\n",
    "    discretize_2nd_order(expr)\n",
    "except ValueError as e:\n",
    "    print(e)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Linearity\n",
    "The following function expands all derivatives exploiting linearity:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPkAAAAYCAYAAAA4V1u3AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAHIklEQVR4Ae2c61UeNxCG1xwKIHEFIR34UkFwByQdxHQQH37BPw7pAKcCjDvAqcCOO4g7sI87IO+z1shasVcY7X4bvjlHSBpJrzQzGl1Wn/3o5uamyun09HRPvAOFfYV3yn/M66w17y2bN95a9TrHuEvougTmHLqY0sdOXllC49zHCjj2a4Vn4l0pXj15y+aNt3oFFxSghK5LYBZUwZ2hH6U7uYRmBz9W/CpFVP4P5Z8r/jXlryntLZs3Xpsu1ceF+F9D2Znylm6r/r/lldB1CcxNNUC+k7ODvwwKSMfMjn4o/pOUubK0t2zeeA11StfnYvyocKnwUoET1kOlErougbmR9smd/F+N8o0mWGPHCHl4i0w09X+gwES/D3nL5o2Xy4a8ODiO/kHhncIqycF+JXRdAnMj7bObjkrGYMfuoi8qeN5VWJi/J3zCnclbNm+8VDCcQnnk5aMni+tqHTzIdS/7ldB1Ccwg68ZFDSdndBKeHeSpApMLei/eW8X7Cp9grJW8ZfPGS/T6QulPwjcbJEUPM1lC1yUwN9E6DSeX0Ndhch3ZYMVLj8qTnVztuVt+DnhM2g/izf4k5y2bNx76ESYLLA7OTv5FeV41iKM9lH9wJPlLzMsSmIcyDqddm+9sjvhP3wm5uD13rIegyEpxY0Ipz1HReP9Y/aFY7fYVuPdwEviToDQnBJx+VlLfGLRSbHLU/St/V9lc8UwZGs9rBV4wON5ekFZojNnqPpRY8rvruhAm87p+mRI+cx0Hx18Wf5GqnVyDYgdh92g8nSlvxKSD3nyLRv3FOG+FjbBGvylRG80YpWNv2bzxcvmFvx94nfdw1eGlw+rlEPfOC5sn08VJ43Cfl4Uw2cHR2S+mNPXDqRe/uTXfVTar/ey4jnN/Veddx2gmFA7LcbtSbLsxx5LHyjcWB+UvxKfNmUIk8X+ImZZEaMdikxNfmCuVt+1qH8XvWy1dZdMwvPFyWetnSsnUagvxmfjc1+PVSelee+QdkO9rozJOXpwk2vTdBlfzaKOEp/1K6LoE5l+SO/oHypAuajvCJ28k/uz221WnrDa1E9tA0ljlZjQMWIVBflbMkYT8gcK1AndJI3Zs+zJsvMFYGK2TSnxWSo7/dZ+DQKGC6rvK5o3XIQd3uujAaZ3Q/wvFcVFTmkkzZI8UphrZ5kr1zhVwilGkum72E5ar7RCgECbOzFgv6SMh/IaNM9oy9D+7/XY0kHqXVPw+GWCaxMjslnZ8ZNewdBX4ODoLBXkEJtw6plA+M7nKprF747Wpg0nTuouLj+7rxTZp2GuPpF6aHGyT2BVbLkEldF0C81lQTm4zNr3oJ6HOIvbbkTFZab6GQTQilbGDsiLVO4fyODJGj6uT0hDtmZyV6hiWxbAjqdxOBpFXKqG+vGVzxeuQm0nTteCymMaJo/SgPfI+Jrahr9nslY7V23Zgl8C0MQfsOqs0PoLe8o1uEfuxk0Ps1o2jlgaK03LXeJoIwODb6IuYtkpSzpNBPFLCgITDSkbdOclbNm+8qAvpx5w23xXQHWX5wjnWHrEPJaa0YbFJr2EpzhzpErr2xqw/Rgf7mE7wGyh+qF7SfruMRAPg6Sa+y8IL9JP4+cSysjTGwePkUZsjBe5zODXPaBDls/8jC43BWzZXvFoz3//YaSju1t+LKpx87ALZsEeC0Zdsa4Pt6XcR8rYdQnhjCo97Nwsh851F8bECdsp/zLSY/Won14AqDZCvgI0vgfAz6nJ4HLhxhBceK+ZGUAHZXHWVKInJ0mUDdJzrP88b1C17WEELhhW1tRm7qBiGe+xtOwbojSk8FuW4MCvP+3jM06cI/eb2yvN1xVC34U9W0IJhRZ3227EaY2IN3u6kbav7rSPmGMyRdVBGl0JGQvRX85ZtLJ7q8WZ6FUbHffyyY6ToHkNGGttHbKDExDbYuWuypbBD6aL2myjT0Fjr8jGYqvNEgVeIaBd4AiDkH0gXs98kJw/aOVMcP8ZIKNK8EXpMhtBFMxI2z3Fz/DTQW7YxeFxp9iRfPVEUt+7k4rOIti2uvX2Aq0AfKfW2SSrSn123Eva0pPqfw35jZZoy+CHMY4GlvoANuY+/CvaKfYX8IvZr/KcRcUQDiWTS8GOYnxUQitV69eQt2xCeynnnZnKgR+51nYulyvhay/eORh3lzYlv2UNl7Cp/K/yudFxA+tqobk2q09qflW9aPEamqWPuw1QZDn6kYHdx7MjOnh/V627Fb9XnQB/3tl/F/wyzDevQwcnJyYHC+V3spXaHU9qp/p7C9ZQ227r980j6XMR+dzmu16vS9s/8Ggg7BMfvtmNf54BUnx2ndXfpbPTthzfsUlty0sBS9ts6uZMB54LRRMHxOLJzNBxL/CR49HVKdblCcOxsXAvGdrat160B6XR++22PWP1HrE3VD8fpUmMriV1qzGvDLanjHPs/3HftqL6VlNMAAAAASUVORK5CYII=\n",
      "text/latex": [
       "$\\displaystyle {\\partial_{0} c} + {\\partial_{0} {\\partial_{0} {f}_{(0,0)}}} + {\\partial_{0} {\\partial_{0} {g}_{(0,0)}}}$"
      ],
      "text/plain": [
       "D(c) + D(D(f[0,0])) + D(D(g[0,0]))"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ps.fd.expand_diff_linear(expr)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The symbol $c$ that was included is interpreted as a function by default. \n",
    "We can control the simplification behaviour by specifying all functions or all constants:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMIAAAAYCAYAAACoRCJ4AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAGLUlEQVR4Ae2b7VEcORCG1xQBUOcIjsvAHxEczoALwSaDc/EL/rlwBrYjwDgD7Ah85wzOGdhFBtz76NQ6jRjtaIbenYVFVUJfrbdb3eqWRsCj6+vrRZlOT0/31HegvK/8We1vJc2YtjfeGN7bSOutb2+8TbTJTimUFo0DHCuz+d8rP1PfhcpJyRtvkhBbNMlb3954m2qKR/mJoEVzEhyrfJ0LrPafaj9X+UfeP1T3xqvxE593GruK42/Utnptyr3s99a3N94mK708ETgJXkUF5HJzMhyq/0ne2VD3xrvBUjKdqfMX5XPlV8qcaNuavPXtjbexdikd4R9J+lGbqxNRY5u+sZvMG69PkWx+nABn+Ev5s/KdS9LxgTJruU3y1rc33m3WttK5uzm6DEHkr6WfGnheG+zr98YrebB51LenzAc9jnonnSCui3WQJydvfXvjTV7YGiZ2HAF+WjxR6akyG4v0VX2fVO4rf6djTPLGK3i/UPu7eJisxfD2Nb317Y23qRbpOIIWfRk31pEJrL78yB7lCN54mUw4K07AifBTfHjVokxyq711yVvf3ngYRJiHKrhZ/KCtRJBljy27jQTCVf7YMfC4aATtbCa1uW5Y399GP1R64+X8hP1emRcsrhLvqCubjDnp1tS1foLYotSD2htjP8nCw0Z4lVT9reo4AXtq1Guk6N1TcAQJRYQlunaeTTNudnf9mPVVq954fYzEYz/2V78LRMNLl9H1wdyqT9g8K8+eJMfG2w9bSFHo63dTmPq4YbC3ghNbPyX0ymuznV2NcIArMa79BhmBPmk83MVV4tkkjrfHapcO5I0XmBU/wlNuTWb1szn4fkjXOdWH5C5YBINU5wjvrTInUvNpBL2YEHTKxKvXooL1Tf3Loqa3vr3xWNoH5bSH6NCagg3pp21J/Wu33a6Y4pFho5sgealxMxoGXEQhf6jkaKN9oHypzJ2dtisemJXEPTNt8pwmyvBCZdo8qqPcqtz5fKs3zrkQ3ZlyGQwMplOKrtdp1E/E3FcZ9NqZtKQheld9e+MhujDZ8Mh5TjtL7C0CcLJj5L922+1IkBCJVH7NBMyrGJiIZFcQIqTVF7EfZ8CZSN54/6He/IlyaycYMgbHzaYNyZ2RpurgnGz9GHqO5K1vbzx08iwqprQXwTPtpUgzi+12ZEi88SoK0Sk0RpTCa0NkjZsdgycPjhOYz8ZceONF/L4C5dacF8dMCm6Ru2Qwcg687OQsoVba9ta3N16++IgdulRnH6Gz8vtgFttxIpCI+p0jW4KysbnXPc0WgPB9iV+2WSRh3Buvw1PycPogSxlhFnEMx8xTq9xT5+CQ4WqYA6yx7q1vb7zwyBJtY2phb5HSA8yctttFEgnAc2R6j6cvpl/VX24qG8tLnCBtNm+8nFGs2+mTon5Gg5PgmC2pI3fLBNH0zUFHdjVshPEj89b3CvD4DiBQ8C1F0HisjI3KX4bOZrvgCBJoIQH5cu98vdNfpJpT4ASd65I3XiEHSq3JiiylnGXb4G7IbQM9GDbUN6fV8QzDvfTW9wrwCFopcAmf3x+kdlQIui1tVbYjaQi8nT1nAz0YNlS13Y5RtJQS3r4n+qLfjWvKEOYYPNHyrnwRMfk+OK/gIyMLTmkMH5s0cg76qBnFIIdKDF4z+tDcpvGRaxrEbMETzRNlXtaSTegTOLl80JjNdqMcIWrmjcr0YahFUed9eOpGaMXjNWHPFKqy90RQPw7Z56hL+YCrDI88LZ2TEcKPv9ScnMSbPxxcx58ZtK6pdS1DeMcCyvcLDsH3wWuttxM8Y3sW23X+Mad15RLYNgy/UPstLmpyNGvBEw2/B0CJ8OOuWXU8jfEScVTSqF2VW2NEqC/KL1VPTrZsjmhDEk0vPxvftLJlTWNkXoanMZzgSNm+DbAhJ0R5LQos1d+rywEet7bdgv9Qu2/55OTkQPlsyro073DMPNHvKV+OmfNAW99z0uUstptyNQqeu8k/YrThqtN3zFZFFz3RqzdSVSfpdNIYEe8hOWhgLtvdS0fAHlIom5PrEUdxa+JPHJqveKLlusYxX72mtTJ+oPtfA9Ln+m13349pri6rWuMqsVcl813CXaV+S+x/AQ3mxxrhRin1AAAAAElFTkSuQmCC\n",
      "text/latex": [
       "$\\displaystyle {\\partial_{0} {\\partial_{0} {f}_{(0,0)}}} + {\\partial_{0} {\\partial_{0} {g}_{(0,0)}}}$"
      ],
      "text/plain": [
       "D(D(f[0,0])) + D(D(g[0,0]))"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ps.fd.expand_diff_linear(expr, functions=(f[0,0], g[0, 0]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMIAAAAYCAYAAACoRCJ4AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAGLUlEQVR4Ae2b7VEcORCG1xQBUOcIjsvAHxEczoALwSaDc/EL/rlwBrYjwDgD7Ah85wzOGdhFBtz76NQ6jRjtaIbenYVFVUJfrbdb3eqWRsCj6+vrRZlOT0/31HegvK/8We1vJc2YtjfeGN7bSOutb2+8TbTJTimUFo0DHCuz+d8rP1PfhcpJyRtvkhBbNMlb3954m2qKR/mJoEVzEhyrfJ0LrPafaj9X+UfeP1T3xqvxE593GruK42/Utnptyr3s99a3N94mK708ETgJXkUF5HJzMhyq/0ne2VD3xrvBUjKdqfMX5XPlV8qcaNuavPXtjbexdikd4R9J+lGbqxNRY5u+sZvMG69PkWx+nABn+Ev5s/KdS9LxgTJruU3y1rc33m3WttK5uzm6DEHkr6WfGnheG+zr98YrebB51LenzAc9jnonnSCui3WQJydvfXvjTV7YGiZ2HAF+WjxR6akyG4v0VX2fVO4rf6djTPLGK3i/UPu7eJisxfD2Nb317Y23qRbpOIIWfRk31pEJrL78yB7lCN54mUw4K07AifBTfHjVokxyq711yVvf3ngYRJiHKrhZ/KCtRJBljy27jQTCVf7YMfC4aATtbCa1uW5Y399GP1R64+X8hP1emRcsrhLvqCubjDnp1tS1foLYotSD2htjP8nCw0Z4lVT9reo4AXtq1Guk6N1TcAQJRYQlunaeTTNudnf9mPVVq954fYzEYz/2V78LRMNLl9H1wdyqT9g8K8+eJMfG2w9bSFHo63dTmPq4YbC3ghNbPyX0ymuznV2NcIArMa79BhmBPmk83MVV4tkkjrfHapcO5I0XmBU/wlNuTWb1szn4fkjXOdWH5C5YBINU5wjvrTInUvNpBL2YEHTKxKvXooL1Tf3Loqa3vr3xWNoH5bSH6NCagg3pp21J/Wu33a6Y4pFho5sgealxMxoGXEQhf6jkaKN9oHypzJ2dtisemJXEPTNt8pwmyvBCZdo8qqPcqtz5fKs3zrkQ3ZlyGQwMplOKrtdp1E/E3FcZ9NqZtKQheld9e+MhujDZ8Mh5TjtL7C0CcLJj5L922+1IkBCJVH7NBMyrGJiIZFcQIqTVF7EfZ8CZSN54/6He/IlyaycYMgbHzaYNyZ2RpurgnGz9GHqO5K1vbzx08iwqprQXwTPtpUgzi+12ZEi88SoK0Sk0RpTCa0NkjZsdgycPjhOYz8ZceONF/L4C5dacF8dMCm6Ru2Qwcg687OQsoVba9ta3N16++IgdulRnH6Gz8vtgFttxIpCI+p0jW4KysbnXPc0WgPB9iV+2WSRh3Buvw1PycPogSxlhFnEMx8xTq9xT5+CQ4WqYA6yx7q1vb7zwyBJtY2phb5HSA8yctttFEgnAc2R6j6cvpl/VX24qG8tLnCBtNm+8nFGs2+mTon5Gg5PgmC2pI3fLBNH0zUFHdjVshPEj89b3CvD4DiBQ8C1F0HisjI3KX4bOZrvgCBJoIQH5cu98vdNfpJpT4ASd65I3XiEHSq3JiiylnGXb4G7IbQM9GDbUN6fV8QzDvfTW9wrwCFopcAmf3x+kdlQIui1tVbYjaQi8nT1nAz0YNlS13Y5RtJQS3r4n+qLfjWvKEOYYPNHyrnwRMfk+OK/gIyMLTmkMH5s0cg76qBnFIIdKDF4z+tDcpvGRaxrEbMETzRNlXtaSTegTOLl80JjNdqMcIWrmjcr0YahFUed9eOpGaMXjNWHPFKqy90RQPw7Z56hL+YCrDI88LZ2TEcKPv9ScnMSbPxxcx58ZtK6pdS1DeMcCyvcLDsH3wWuttxM8Y3sW23X+Mad15RLYNgy/UPstLmpyNGvBEw2/B0CJ8OOuWXU8jfEScVTSqF2VW2NEqC/KL1VPTrZsjmhDEk0vPxvftLJlTWNkXoanMZzgSNm+DbAhJ0R5LQos1d+rywEet7bdgv9Qu2/55OTkQPlsyro073DMPNHvKV+OmfNAW99z0uUstptyNQqeu8k/YrThqtN3zFZFFz3RqzdSVSfpdNIYEe8hOWhgLtvdS0fAHlIom5PrEUdxa+JPHJqveKLlusYxX72mtTJ+oPtfA9Ln+m13349pri6rWuMqsVcl813CXaV+S+x/AQ3mxxrhRin1AAAAAElFTkSuQmCC\n",
      "text/latex": [
       "$\\displaystyle {\\partial_{0} {\\partial_{0} {f}_{(0,0)}}} + {\\partial_{0} {\\partial_{0} {g}_{(0,0)}}}$"
      ],
      "text/plain": [
       "D(D(f[0,0])) + D(D(g[0,0]))"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ps.fd.expand_diff_linear(expr, constants=[c])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The expanded term can then be discretized:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAeUAAAAsCAYAAABboKCqAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAL5ElEQVR4Ae2d7ZWWOBiG1WMB49jB2IEfFezYgR8VOHagx1/6z+N2oFuBqx2MW8GqHWgHjtOBe99sQAgJkBDeBLhzTl4g37mSPA8Jgffqr1+/rqQwr169OkE6z01ad3G84DXcvxq3wQPCvUWASxPoNa7r88F4e/EUH39Li42fzaF90Baz5MChy6v8RKA0AtdSFMgMxLc4PjX2DtKlMv6C69OxPBDmDcIcw76HPYMdjYMwuzHi429qsfGzObQP2oIKOVoOHLq8yk8ESiSQRCmjYlSqT9sVxADlrJmz3Q9td885FTEVMhXzZ9hPsKs0qPcpLOuT0myCj9ik7BJFpjVXDhRZqZhCLdTXY4pSVBxxGW+O6+NBJoXgzPYbgN+AbS87U7k+gNsJ7HdXSnBn3CPYTybuahWyqR/rQpvEbIyP2CTpFcUmEi0Hiq1RfMGS9vX4YhQXU1xGmiTVTJmK9LtRqq4sh5TU/ZG4rvT25CY+/tYWGz+bHD5z5ECO8ipPESiOQJKZMpTxQ0/NbtMd/r3NXnDjkiyFKu+uL3DNZW4eO8vgcNulER9/s4uNn01OH7RLjBx4gDLfg/1hyv4RRz4Cemeud3sAA7HxtP6W2SRRyi5ugEaFfAJb78juBDOD7h2O3P7N3dZ/dgLs/EJ8/B1AbPxsSvNBW3nlAPz4DJoKmBtDr+BIefEFlvtKdq2UxQY9wGO2zuaap94pnDnz/QiAXmVrBiHz8j5HRpjqmXSKAvnSQB7PfH453UvgIzZpe0CpPNPWspOaUw6AA2eBHHd/1KHhxn0nfNR1XrvVR4aHpdJexCDtYmQA64pKio2jpUtj4yjiqNNYX7vK95QRiAPhH9ihZ792Zg8Rr7cszUBw5zvH3NzF5WmvgT873wccr7oCwZ1L3HxW3VHauOYdNg2XvG7i2jkbr0KYn7E48K9e5WjHGTpnePhz6d023EFOc/H/ofP7FfF8S3ydgLxA2GA+iMM2fATL19OqGQjTGjIIN8gT/ptgM8TA51cCT1/ZUrubuh5EDiCvnyg/N3c24wHnnFVzpnwL583GUJz35ADcBvusi81YHPgH9XPmwTg4JJUDSHMum83KgLlsXP1izA15HpRntXyNTC9RsEkCfEIFOICOkeagQjbp8FlSM/jaaRsQ93FsBi39cc30f+BYzcBxPIU9h/XmB78pcXhz8AZ2VMGbcjiffSM+FSlvSLwrBIw/0QTxQZ4UarWAYEcaNYizCzajIBwBSuHpKNoiTqjvQeSA4cr++d6qCPvuJfwbmYBzhuvIAbhN6bOdpCfGCZIBzADpJpUDSI9jeA6bzcqAuWw6HWLihcnzoDI16fI1KkCFxLvcRpHinArKt+zEDuScbcOdd8K8C7UN3ZuZM9LmORWzLw/GH43TSocDohQTxAd14EycNwONUJtQkV2wmcChF2RnPHv1j3UAtzE5cNekbY993lg3Y9uEccmB0T7rKPtoHJS7liU5ZcAsNhvvs7PYOPrEqFMOnsmUMgpPBXIPR3umyQHqWsolEEL+lycOQ0XbGaC4puLlgLGVziXcmH/PBMapBmUvkXwOQXxCiyk2ocSGw6+c53DlJvqCwWQ5gLDNOMY5xzVnJPbz5I4cCGRclTowThEyIIbNxCbqBBObDo7ZFyl4Xp9dCiRgCsINHXxGZM9uOah6S7kmDgeifbd8xfhR0dqG4V2GSv/Y5QG3kDi8QeDdOl/LyGoi+YSWWWxCiQ2HXyXP4SpN9zV9dooc+Bup8vlt+6NCf5mc6FcZjxwIYVwnFRIntwyYw6aub8hRbEJojYedzTOJUkY5eXfLWSyf9dimp3RNgGpmi4HXmQ0bP6blm13b6fOaCtkHwxXeF+cSHsy7BJOST2h9XDzFJpTi7/Cl8/xd0nlnk+QAxjyfG/Pml3s4qARvwnK82x8gCpEDLsZIctC44mTt5wuxGYTg8RQbD5hI58k8kyhldKRbEQUdmpFSwXJw2MblxjAM3yyFWZFC4oTcCFjZJL+M4RNaiN2wQR+9DTh8rjhmuGvd15fG4q6V51i9JvmHyAGE5c14c0OOa+66bq5NhhzXNlP72gTdjgyYwaZmEXIM4ZldPs5hg7irkAFJlPLUHgAofL78GEduBOPz0teeuBSKHJAdg3j1nTTvoL92PPvXlXdgHKYbK5Dr4rCT+zp6HcZ5nMvHmeiA457YoK7sL7zRWcxk4LlYXZZKGIwoGF/APsF5NU6MG92fWPn25EAg4yq5wDgpZADzDZYDhkM0G4vdpMsMbIK5sCIp2CCNVciAa5NaLl0gzlSOAOeISeLofHZr4HFwuAwV+WntgbA850dKKmWKI9O3Z0SDceq0cGSe31rXwafIm8/VY79GlIKPt8xi40UT5ZGbZ1Sh80ei0mmPX8oCPk/u/fc6+FKIuuTA4HguoV1Qhhg5kIINkPlNbjaRXFihxdn4qfl9luCZQynzuROVT/PalKfKnBX3BiTcuGmMHwzh8yh+9YbptO+wGecMfpyVV2ZCnDro0JJxHWbJI7lE80E9uXGGabADV+e4bn+paLdsYhpt4zxjkKSI8xaJcJmaY5R9lZYKmePaZXpywITdogxIwWarMmA2G1fnGnPLIgP4Ra8S7cuXL09h38SWDXEfhMRF+CPY85A4OcOirNF8xCZtny+Bp+kPZzn75BJ5q5/7+6rYLMMmph+nlAGHnimP3Zg0/rhD4d00l6J7s+UmkOcEcbg8Zm8a8YRunHnH/rS5Kvwklo/YpG3YgnhyGbh6LJS2hnlTUz/38xeb9Gz8Kfp9UsuAYpUyEaCyVJLcDRsqbLiEc+nH2PVBWL7KxU/szd3k1U144SuUN4aP2KRtl13wTIssLDX1cz8vsUnOxp+g3yetDIiZqh86DpYGjpbMc+n0lyw7016y/EumvTSXpdnElH8pnkj3AeyzmDKtJc5S7ErsJ6FtIjaDS9mL6o/UbVX9S5T/BkA+IiACayCAGRM3NvKO3bdhag3VUBlFYPcEil6+3n3rCIAIiIAIiMCuCFzFssevXdVYlRWBgglgpuv8b/G6yPDnqyHNe761O47H5vyi5Vafjv6PN9KVHKhp6SgCGQlo+TojfGUtAqkIQKlq+ToVTKUjAhkJaPk6I3xlLQIiIAIiIAJtAlLKbRo6FwEREAEREIGMBK5nzFtZL0AAy5h8p/uFSfrEHJuP/y+QpZIUAREojIDkQGENElAcKeUAWCsJym+CN18mwzk3BvFv8WL+XnMlVVYxRUAELAKSAxaQtVxq+XoFLcW7Xthz2J+wY/9ixQ/9t3fn8vOhfH+Vf40nIwIisFICGMOSAyttu5BiSymH0MoUFoPxEpb/YMXPgI5905uz5M+Ziqps8xHgZ2Unf1o2XzGVcywByYFYcuuKp+XrdbUXZ7uvh4qMgWv/lzOVNP/+jv9NK7NRAmjfsZu1jdZ8l9WSHNhws2umvJLGhdCtl6QnC1/E4eDl+6t3VlJNFVMERGCAgOTAAJyNeEkpr6chq+VrDMpJS5QId4Kq8Xnynalx1oNCJRWB3RKQHNh40+uLXitpYChW7qDms+Jz2GNY7qbmTJh/bdn5y0mjkJ/jWO3CNtdXcOyEQ1wZERCBFRHAGJYcWFF7xRRVz5RjqB04DgYi3z2mAqbhqw6VcsXxA69hH9KDBm6cIfM1KCrlOg6V83NYGREQgZUSwHiWHFhp24UUW0o5hFa+sPXzZH4ExJ7tcqC2De+k6cZjYxCvmjU3DjoRARFYGwHJgbW1WER5pZQjoGWIwudI/Kcfewc1Z8If2+VBmBvta52LgAhshoDkwGaa0l8RbfTysynJh3fInV3XUL5cpqZ9X1JBVRYREIHFCEgOLIa2nISllMtpC2dJoHy5FE3lyw1ebcPlaH5UxJ49t8PoXAREYAMEJAc20IgTqyClPBFUxmCPmDcGZWemDKcz2OpDIfDjpzWpuGVEQAS2SUByYJvt2quVlHIPSXEOfI5kK2QWkjPoeun6FpSyvQGMYWREQAS2QUByYBvtOFoLbfQaRZQ9AJUvX32yDV9xegxlfAqr151sOroWgW0RkBzYVnt6a/MfM2R8O7RQA3cAAAAASUVORK5CYII=\n",
      "text/latex": [
       "$\\displaystyle \\frac{- 2 {f}_{(0,0)} + {f}_{(1,0)} + {f}_{(-1,0)}}{h^{2}} + \\frac{- 2 {g}_{(0,0)} + {g}_{(1,0)} + {g}_{(-1,0)}}{h^{2}}$"
      ],
      "text/plain": [
       "-2⋅f_C + f_E + f_W   -2⋅g_C + g_E + g_W\n",
       "────────────────── + ──────────────────\n",
       "         2                    2        \n",
       "        h                    h         "
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "discretize_2nd_order(ps.fd.expand_diff_linear(expr, constants=[c]))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Product rule\n",
    "\n",
    "The next cells show how to apply product rule and its reverse:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAH0AAAAXCAYAAAAm70AZAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAFwElEQVRoBe2a63EUORRGBxcBDEsETAg8MrAzYMlgxxmY8i/7H2UyMEQAJgPIANYZLBngdQbeczQtoVa/24zNVs2t0uh179Wn+1Cr235wc3OzKOn09HTJ2D5lRflC/7LkmdpHx4ryfarcjv92Fmiz+16pEiadfUzR0e8ozxm7oJ5NyB8h/HS2gp3gbSxgsmn/RA/yTGfSDD+mfp04aFRCL6j/zMfHtJF5CZ+yNZ1dsvCdM3ddzb+hH9tdIrvxAQtgQ51+TW0SL8pMN8PXTOr8nGR+yfikbK30NIIoV5y34T+j/wflA2VN8dTZ0S0tgF3fouKw8kfD6f8w+ZHJWnZVfcemOkEnmrljSUfrcB3/jfKFsqNfYwH9oD8WteO9TzeONyAuqUcf8fD+S3nUpzfOwWdAfaY8ol0Lusizq29nAf2BhicPSzVMmG3PKNHwXxn7RN+b/OjbNzI+y0fzw3sgP3JxXbo7+sUW0B+vak7H4Gaahj+Mi9HepxgI0lQnDh7PlW4dbqZf0fdNwTphoP/bEfjCBRVgPypwJoa2Cpel3w7wBpD+OEgXOcDq8AV1zdj0ZYxjf8szkp7D5yOhl9D/juIjY0k5t02J6yVZxrxIetpsjdBfe7XpWgg+n43hgkrbS5IO1zaNRx/zW8U9FjPYJP2xCpmOYLwpe6y3kQ6RPvoDf7gQ0DTKH9Nvex1T5ooySMhHZ7aeDMyLzxMonTS0x2CorT0kw/xbioHXCLqoiDkz3OBIdxXGxOZ+Q+JkvFvHzbqDmCMeav2xipmu03yP6/ryplM+MS+PG/lB/ZpilH+mrm2WMckb+HVoDf+EV0H0NNZnTGMeUKeAoD0WQ1p5gswFvDGgknzWeE872CKOwR/wO56N3SXuIcwRlkmz3AOw4HRqMmrksGbeZ60UX700SOJl3rbPMnXMpRcIpiwulLheXDtOzcEwSibbj3apEXM61/EPtYnNfcSEyPdwZ7j7MBc4QyLuMWhD+rqpGr+eAr6q+Q1ex7rpfHMKXFNitNuXPEoahgszzR9lG1lesRlQeZBNwRBUTMStTAjkav288p4ilVi9iCaMgWOTCGlsBobFRJkuzBWcUOmPKzNdB+q0BjHn88tMjxeULifq4Bg8UY96x2a/xmwEXbXpEtsUDBHLVBmx6MhWqmwW5mirWxulR9w94e7FXG1EH3030yWzuXZ5AbjZ5/PrWb5J+m2kstKwZoPHdi9VBlK2zB7lDBoDagy1YRiS65K5RrAtYONFNp/TRlKY2zTvBXcX5gpSqLyoX8bbu69N6R0543rCuMoi5e04Zr2klEe+z70LJwfI4FqwTjoKM371lmuW/cjehiHOTZVpDTRtQfEEOKM2sx5T5C0/KomlXLPswxJI3tJ21VRDRxxvk2nFHAWq2hPpMDjdATbhzTPdPh0rCZ64OSP9spiv9eH1HrCgPLVd8OZdjdi1rsZwg4nQNRpDFJoh4/5aHYEugzMFKH3fz1O/WvM+cHdiFhM4w+kk/r0K5JTqDcxGTCCU2PYVps1I3mBrjw2F4PWDRTwFfJ6Xt2HZ5DNYAtgw8POnFwNyS4pr59QrkzPSds3ahyX0Gby+GqUgdAw+S+3tgvH7wN3ADK6cfIQHm0x2Ohvy3dwPMh5xfqTwkvcXpUHM+0nSP+KXjnNxHRMMSN2V6eo0s2vy9IcwyL+Gz4tooBEykdW67fQ5ZjwPdrH7PPd7RdtJdte42zADL2W5fgifiB+G0Yk/CBs1Y8mgMKsFFSk4nY51fDOIc2Utj6dFbc0+DJUT/GtdcrpK+2Scl+CJgVieXDGbDSaf5fLp8PJoV410Z7h7MG+QbE6inyeu/zmz7XJycrKiHM1dB9lzdUyRh3+fspwiI++ctbrWmKNrDu6+dZg7otRst3WHdxlk6jjAzyijnQjvesYaa+T2p8r18W8bt/ucivl/43QNy+ZGO73PEV1z29K/Lb1zbfIfzLr+GVLUm+IAAAAASUVORK5CYII=\n",
      "text/latex": [
       "$\\displaystyle {\\partial_{0} ({f}_{(0,0)} {g}_{(0,0)}) }$"
      ],
      "text/plain": [
       "D(f_C*g_C)"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "expr = δ(f[0, 0] * g[0, 0], 0 )\n",
    "expr"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAPQAAAAYCAYAAADNqdAHAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAGe0lEQVR4Ae2c7XEUORCGhy0CcEEE+DLgIwJMBiYE7AyO8i/7H2UyACIwJgOTAXfO4MgAlzPwvY9WEhrtaD6W1a5mZ7pKlkZqtVqtftUz0sKjh4eHah/o4uLiQPM4UjpU+q7n232Y1zyHzVhgKv6x2Iy5ditFiwWQz5QA8Well6q7Vj7TbIFqSv7xqC1CyxCf5A/31ic+6NmVi3ET6URkPlP+PlRKz3/r+ZXyt2H9XJ6WBbbhHxqjGJwkI7SUvNTSP1G6UjpRIgqWSETmE7twoX5E6mPVPw8r5/LkLJDVP0rDSRLQWnZADJgB9T9K35VKpP+k1FcZtvb2YJ+pK3Uj2qotZY8jJdZ0apTbP4rCyeOm1WXxVc+rLIdLgKJUMFfSj0icojs1vEo1Tqye9SRNinL6h2QXh5NGQGvF3yj9lMKAuXiSnuySL5Scvj9U903Ph0o/lWaasAUy+kdxOKkB2k4cJdl57vTMSTH5aan+IN1upBubj9dR5fD1cjCg1f9YMonsv+y82RyQ2fY2YFn3OxubbaRvDv8ggBSJkxjQOOxnGYHLaU61P5bsnnaxKuUezOirZz4V/rW6u7zXVNSPw0DAS8SvlBPlkcE5wqQBPTbbSF/AzBpuzD+svGJxskDBkDR5HBhKfjeLh9Njx7fk3uBfyebKqZXEwy7Jm0Ttuiro5L4XvwZ1rUXmJQbGfu0YVUeER5ZxDldPDr/STu0Q6pOzzFwlfzS2kb4b94/QvsG67wwn6CM9alipRWirsLnmEWPjL61Uj6F4xfWvsioT1SBeUZ/qOQUyw8Sftj5q+6j0Sam2s/rOywJj3IunUU+1AbRvajff1cr76Pgl7KNypX7u2ovXbk+qL8UOXicK0os7UTa6mLitoL3Jpreq77qvH5ttcvhHaNOd48SuZw0ri1BDW+bb0YM1bNeiE6neKPe7kso49i/l75V4Rb9RvhLNIjl9+lxLjgNh2L1SPXoAWK9HyKB259A4N/yd44mHBULuFX0CQhYbh7eJykXYIdDRF6XbqdJfcRLDOyU+o1baVNcKZrWPyjbSd+P+4Q38u1AKTtDIY6UJ0CxeKuoBMAMSpFiizgNLxqR8pBzApaizTyCHxYnpia34ETfYZ3Znoo7Tq3M88b+0feO5c/jh5FiWqhQ7OH1y52OzTQ7/iG1cBE5QKsRKE6BZvBRQAKp3bpUBLYDz0Utl6F6JCa/QwD5mc4iFSAbjMcYKqY1vPaKqiToDx6usbCNXZeaGrPiNowg7GCW3+GcstrF6ZvGPwNwl4QS1DFZqgJYhHEDjKFXZtthIOHwT3anS7ZJx+5A+bCxEyCYiCte+B6Ujmwjfei9UdptM3/HM4ZmdpxsPWZA/WCvQDksN8/4do2027R/ewtYH8KtScIJuBiuPvZbLgomqUthH4aAdsAPUPgSYU0BK9W/qcy9mxl0h6cj1mr8rDxieqZ5+XVQbjz5KbB6XyjHOUyXmG//Apig7SL/sNEbbSOeN+kdk5NJwgnoGKzGgcejaaW4wEQAaAyV+duzwugjp6lw+pE/rBqJFQ9eUvoPHkzw2Mr+Z6Zn7Z/9sBTK3eA7xs2Vt/BxxbUP6tNrBCcyZj9E2m/aPwL6l4QTVjI8sNGnuUq+tsnwXXNlynAFQnNmT+rlv2aYoeusZg8LAPshNbQyB1HSxz3jiea7ESaGfH3WSSooPAUdpB82DDSS1iaipmfbdNppfLx8WX8k4YfEMVhYqcGJ7IIWNMytvjHiqB6BNwP2geg6ODImPMve/BojKkc0YIbX2CRgZj38t86fUNd6ZBgjngC34fuYqrrYx2efR2UF68+u5dX7ptve20Tp3+YdYisYJ+hmsLKyinOICutb7SLXzPVlzZj1z98yPSfj25FcryHin5Aj+E7Vx+myoRx/H2vZq43g68x7jEYV5tUZP7EBy9+oqrtAo7bAyi34Ve2+bHv6BpfCJUnGCfkus8D+W9E3n5+dHSpd9+UM+9TsOn7vK4j9Quuni20W79JrtkPCb2TYP1TZtgP9rPI8VInRv0k5GFOMVuhaluwSIn9fZ+HCpqxs74mkX0y7aZzukrT7bxv/QY1s4YTE8VgYBmp5aMEDGzwvNNzd1PehQ/L0PZMTLTzU5pPqjA7Eeeq3NIt1mOySsN9tmOzjB/LJ1HSvrvrIS5tft29Uvp+yusYe259Q1p+yh81yHP6f+OWWvM9dUn9x6xvL/B4h2hVR2GmXBAAAAAElFTkSuQmCC\n",
      "text/latex": [
       "$\\displaystyle {f}_{(0,0)} {\\partial_{0} {g}_{(0,0)}} + {g}_{(0,0)} {\\partial_{0} {f}_{(0,0)}}$"
      ],
      "text/plain": [
       "f_C⋅D(g[0,0]) + g_C⋅D(f[0,0])"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "expanded_expr = ps.fd.expand_diff_products(expr)\n",
    "expanded_expr"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAH0AAAAXCAYAAAAm70AZAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAFwElEQVRoBe2a63EUORRGBxcBDEsETAg8MrAzYMlgxxmY8i/7H2UyMEQAJgPIANYZLBngdQbeczQtoVa/24zNVs2t0uh179Wn+1Cr235wc3OzKOn09HTJ2D5lRflC/7LkmdpHx4ryfarcjv92Fmiz+16pEiadfUzR0e8ozxm7oJ5NyB8h/HS2gp3gbSxgsmn/RA/yTGfSDD+mfp04aFRCL6j/zMfHtJF5CZ+yNZ1dsvCdM3ddzb+hH9tdIrvxAQtgQ51+TW0SL8pMN8PXTOr8nGR+yfikbK30NIIoV5y34T+j/wflA2VN8dTZ0S0tgF3fouKw8kfD6f8w+ZHJWnZVfcemOkEnmrljSUfrcB3/jfKFsqNfYwH9oD8WteO9TzeONyAuqUcf8fD+S3nUpzfOwWdAfaY8ol0Lusizq29nAf2BhicPSzVMmG3PKNHwXxn7RN+b/OjbNzI+y0fzw3sgP3JxXbo7+sUW0B+vak7H4Gaahj+Mi9HepxgI0lQnDh7PlW4dbqZf0fdNwTphoP/bEfjCBRVgPypwJoa2Cpel3w7wBpD+OEgXOcDq8AV1zdj0ZYxjf8szkp7D5yOhl9D/juIjY0k5t02J6yVZxrxIetpsjdBfe7XpWgg+n43hgkrbS5IO1zaNRx/zW8U9FjPYJP2xCpmOYLwpe6y3kQ6RPvoDf7gQ0DTKH9Nvex1T5ooySMhHZ7aeDMyLzxMonTS0x2CorT0kw/xbioHXCLqoiDkz3OBIdxXGxOZ+Q+JkvFvHzbqDmCMeav2xipmu03yP6/ryplM+MS+PG/lB/ZpilH+mrm2WMckb+HVoDf+EV0H0NNZnTGMeUKeAoD0WQ1p5gswFvDGgknzWeE872CKOwR/wO56N3SXuIcwRlkmz3AOw4HRqMmrksGbeZ60UX700SOJl3rbPMnXMpRcIpiwulLheXDtOzcEwSibbj3apEXM61/EPtYnNfcSEyPdwZ7j7MBc4QyLuMWhD+rqpGr+eAr6q+Q1ex7rpfHMKXFNitNuXPEoahgszzR9lG1lesRlQeZBNwRBUTMStTAjkav288p4ilVi9iCaMgWOTCGlsBobFRJkuzBWcUOmPKzNdB+q0BjHn88tMjxeULifq4Bg8UY96x2a/xmwEXbXpEtsUDBHLVBmx6MhWqmwW5mirWxulR9w94e7FXG1EH3030yWzuXZ5AbjZ5/PrWb5J+m2kstKwZoPHdi9VBlK2zB7lDBoDagy1YRiS65K5RrAtYONFNp/TRlKY2zTvBXcX5gpSqLyoX8bbu69N6R0543rCuMoi5e04Zr2klEe+z70LJwfI4FqwTjoKM371lmuW/cjehiHOTZVpDTRtQfEEOKM2sx5T5C0/KomlXLPswxJI3tJ21VRDRxxvk2nFHAWq2hPpMDjdATbhzTPdPh0rCZ64OSP9spiv9eH1HrCgPLVd8OZdjdi1rsZwg4nQNRpDFJoh4/5aHYEugzMFKH3fz1O/WvM+cHdiFhM4w+kk/r0K5JTqDcxGTCCU2PYVps1I3mBrjw2F4PWDRTwFfJ6Xt2HZ5DNYAtgw8POnFwNyS4pr59QrkzPSds3ahyX0Gby+GqUgdAw+S+3tgvH7wN3ADK6cfIQHm0x2Ohvy3dwPMh5xfqTwkvcXpUHM+0nSP+KXjnNxHRMMSN2V6eo0s2vy9IcwyL+Gz4tooBEykdW67fQ5ZjwPdrH7PPd7RdtJdte42zADL2W5fgifiB+G0Yk/CBs1Y8mgMKsFFSk4nY51fDOIc2Utj6dFbc0+DJUT/GtdcrpK+2Scl+CJgVieXDGbDSaf5fLp8PJoV410Z7h7MG+QbE6inyeu/zmz7XJycrKiHM1dB9lzdUyRh3+fspwiI++ctbrWmKNrDu6+dZg7otRst3WHdxlk6jjAzyijnQjvesYaa+T2p8r18W8bt/ucivl/43QNy+ZGO73PEV1z29K/Lb1zbfIfzLr+GVLUm+IAAAAASUVORK5CYII=\n",
      "text/latex": [
       "$\\displaystyle {\\partial_{0} ({f}_{(0,0)} {g}_{(0,0)}) }$"
      ],
      "text/plain": [
       "D(f_C*g_C)"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "recombined_expr = ps.fd.combine_diff_products(expanded_expr)\n",
    "recombined_expr"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "assert recombined_expr == expr"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Evaluate derivatives\n",
    "\n",
    "Arguments of derivative need not be to be fields, only when trying to discretize them. \n",
    "The next cells show how to transform them to *sympy* derivatives and evaluate them."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAGYAAAAYCAYAAAAI94jTAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAFG0lEQVRoBeWZ23EURxRAF5UCkFEEXjIQKAJEBrJCQBng4kv6U+EMwBEYkwE4AkAZSBkgyEA+pzW3q6d3ZndmVC4N61vV1a/76vvqnt1Ht7e3i7nD+fn5Eh2PaDc0x4e0C9Yv6bcSdn+SU71Bz2sc8U596V/R/UP7xfk2ws/imIvK+PvMzZ6thU7HEJF7nNjSYdn4xPzeJQMeS9r1FEt2yD+Gj1m0FdBlm536ZCDpkNc0nWHpeMba3/STAXpLz8FkBg2hfGhfmb6lT2XtvjxnQm/QaqMMj8rLn00z5TX97xmDQUN0SP9buT5kDI3RLW3NU1k6/Bnthv0n9BsBvKD7yPiPjQT/AQJyrSRxnqS/c9ZblaXQdeMZwdUxP+hTwNUZY6a8bBgyzCDyMeujor7hs+JoubKnEi8YWt4+uTYEpANPo7xhPEqfIfw34SBTp5ixp017ylyHfGVutcnAfPAZwTXI5GngLWrHXLH2XoZuBjRz11qCY39N7z3wds2+Wxr3Yx+OitK+00onhH5G4miA1xHt5WjCOwLPdFrSwstAUae+kr/2jAUvbZXuzpZjEPCO1hJaEN0w9vthDJzIs4+AvXB0b8aA44GVXT4cpHP9PW0K7EFkmwLKvkKvmt4zGERmVIYhZwxkcLWVNtvbjcXoWTSSTE8PLnxm7QO9AkvjuNcL0ByzuQnfUma6h6w+ft5tlsRvDYIB8nQAXR+/+6zrgIM1smuHDT1j6KTNTlqOQZglxQ+5nDWMy7TfZOhgbq9CvZnQIBp9GQdZOlOjezgv0+Qwemu47cEBXQySLkilttG13B90xoJAe7zIpQyGqc7TZ6eIzFzEWPOpOhSs/95ZnQBfjZ9rL3OdIlhnT5rmfPaA7p7DihIvtaTzxDNqs2XKGBhYvvSsJawLIj1TTQc/Pu4sLfvMWwo1DKTxbugD5Ql+wKYx/QeahzQzcyaJNHPw0lf3+vk+5Yza7M4xDNKrAsZ95cJoULD3gU78FkrQW+r8prB0lfCYSSpF5WIxTrWXucp/gT45gl4d+gKkIN88hJfZF8YpCdRtwX5UgnLvkvW+clXipXEjw/LfRTPljAbl3i4Mjexk+BWpLLAfB4tnr9nyPHDZN+J1zNifXOSrEjrBi9xvg77ACHGjevh1GX7BumVTfesIH8vfIH0Mnzoog8+UM6aA3pFxw+VzcKt6s8ko0gE6UEfWjwAz44BWgikp7grAJ4LBn+41nqXAD7SaxwrtXBbQVec+oc+Zwlhna6MF/dQzSnezAwON/ENmNbCncL0ewiXqAp0QDo59+SYlY6HovdyFuEe+3E3v8JHb+t2o2ZtNh34G0MrPTKxpL20hTD2jdrzeTSzu7hgzI6d2I/xP1iwzdYY0ZLmTWe00y9JhxmgPTH2zsA6IkLPfRp/PDJ0NNjPcChLlPRQ8Yi1sOPWMlvbL5BiY+cXvD4kKDI8r7FfWSuOVY/cDdEoYNdb+YiC/LhD/IjaUQTMw/K3I52LeC5wZ9X5W6BzvlxoMxoCpZ7RCnS78dXlMOzs7+047KGmY39KW5Zpj1q5oLdwa5yHm6HRMe/UQstfJRKcl7UqcnXDviN5o1qsJiHDHPqXrjHHfF5yX+9zAzLfNDawa6Rux9X/MUC1xQiIG3w9M/0fJP5/UPMA19S1RXY6r0f+3c+xjefSTIT290x0z1hoQ69mh4IvOu6bvrT+Uz7bj+ZDI1WVSxoy1UBMN/tEWL5axLLYaH7v4edC6Dv4FMIw+Qo3ltiAAAAAASUVORK5CYII=\n",
      "text/latex": [
       "$\\displaystyle {\\partial_{0} (k^{3} + 2 k) }$"
      ],
      "text/plain": [
       "D(k**3 + 2*k)"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "k = sp.symbols(\"k\")\n",
    "expr = δ(k**3 + 2 * k, 0 )\n",
    "expr"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAEMAAAAVCAYAAAAdHVOZAAAACXBIWXMAAA7EAAAOxAGVKw4bAAADHklEQVRYCdWX21FbQQxADUMBhnQAHfDoADqAdBDoIJn88cdAB9BBJnQQlwDugJRA6MA5R1ntrG98CWPwcKMZWVppH5JW0l2vzWaz0arg/Px8zN5fy/7bhX5C/rSqM1+z78ZrFr9g7SWOn+U8+Gv4e3AnZUOi6ys25pQAHDZnXMJvI9ttZINhVx0Ms+JuMN7+w5C1VfaM7tlkhJlxDB1kmdSegYE2uKxvG59ja34CrcBY3XdwH3x8qWPMszSOwT3wXQAb9OlLOTzsd4x8qiwyg4EOdpudhuv0Cfpb6BwgsxHeQTOAc/p2wByNsHm617t8SdIG6FHaBm+mfgaP4CfrRXEKtdkZgITMiPw0pjypN/0jB32UPeM2oB745LjI+pb0yll3CGrrMqDjcxfHXmaJl+OljzIYponCemtMrLwTW9CoMs6AterKM89AmBHX8LsivAY8gsvAmEXiMqDND9jQXa8PY22NngGjYLM9AVlmic50wVTzlnsDVhZYSh4urcC6uRuqitUy+uiF9Nk8rg20tYMFRjHSCv6m1RVefc0K5hi4A1DHbUhxIHQuwOjeDbDlpOfwePOgn84FA4EKHdUxS+evNwJzdNh5kTGM2wwyA8RFAUQ8LCj+Rk/TsuwZYSXKKXgFGsVv4D18Ohtz+DFYwgRd8FC/Ngbpp3LwfwEb5y32X2nws48uJv1ijk5uwmfqmxEfQRuhn9b4RsO/GbCnZ2TQ2323ymBRA/Yi+0qh3SP4coZftvqp3VCDoNZNZ5VlolFivjXkzQAfT3tuCr5pQNhvYYNFbpbqQNwk/FLAej/PW9AaCDfKMolaR2kW9ELRW2MX8BpsmllKEczehQNSYKsB3YHWLIKPt08GwxKwdqIUGtv3C599wPIQcpwN1gCNWO9rbrCAfV7aATSf5GmrAXqMMoHpKnXMCWbKGXwGybSyNnPMMMCyET78IcP7xWYvzEy28XffTr5sryIYMDeggnaSi+PN3rhmcC5yzHwfXgbSgD20upwzIOpfB31a9JyPnvfs12RAjoQpBNxsfXUD7fNrvU8xULnl2S3RNzP1N/jzMT/6sRbgAAAAAElFTkSuQmCC\n",
      "text/latex": [
       "$\\displaystyle 3 k^{2} + 2$"
      ],
      "text/plain": [
       "   2    \n",
       "3⋅k  + 2"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ps.fd.evaluate_diffs(expr, var=k)"
   ]
  }
 ],
 "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": 2
}