12_two_dimensional_filtration.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import pandas as pd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "mesh_x_count = 4  # количество ячеек по пространству в направлении оси Ox\n",
    "mesh_y_count = mesh_x_count  # количество ячеек по пространству в направлении оси Oy\n",
    "mesh_t_count = 5  # количество шагов по времени"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "dt = 10  # сут, временной шаг\n",
    "dt *= 24 * 60 * 60  # с, временной шаг\n",
    "\n",
    "m = 0.18  # пористость\n",
    "\n",
    "mu = 4  # мПа*с, вязкость жидкости\n",
    "mu *= 1e-3  # Па*с, вязкость жидкости\n",
    "\n",
    "k = np.array([\n",
    "    [30.0, 30.0, 30.0, 30.0],\n",
    "    [30.0, 40.0, 40.0, 40.0],\n",
    "    [30.0, 40.0, 40.0, 50.0],\n",
    "    [30.0, 40.0, 50.0, 40.0]\n",
    "])  # мД, проницаемость\n",
    "k *= 1e-3 * (1e-6) ** 2  # м^2, проницаемость\n",
    "\n",
    "dx = 100  # м, ширина ячейки по оси Ox\n",
    "dy = dx  # м, ширина ячейки по оси Oy\n",
    "\n",
    "h = 50  # м, толщина пласта\n",
    "\n",
    "b_liq = 1  # м^3/м^3, объёмный коэффициент жидкости\n",
    "\n",
    "compr = 2.2e-9  # Па^(-1), общая сжимаемость\n",
    "\n",
    "p_initial = 20  # МПа, начальное пластовое давление\n",
    "p_initial *= 1e6  # Па, начальное пластовое давление\n",
    "\n",
    "inj_cells = np.array([13])  # массив номеров ячеек с нагнетательными скважинами\n",
    "q_inj = np.array([40])  # м^3/сут, массив расходов нагнетательных скважин в поверхностных условиях соответственно\n",
    "q_inj = q_inj / (24 * 60 * 60)  # м^3/с, расходы нагнетательных скважин в поверхностных условиях\n",
    "\n",
    "prod_cells = np.array([4])  # массив номеров ячеек с добывающими скважинами\n",
    "q_prod = np.array([40])  # м^3/сут, массив дебитов добывающей скважины в поверхностных условиях соответственно\n",
    "q_prod = q_prod / (24 * 60 * 60)  # м^3/с, дебиты добывающих скважин в поверхностных условиях"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "# отображать ли матрицы в Jupyter-тетради\n",
    "display_matrices = True\n",
    "# настройки отображения чисел в датафрейме\n",
    "pd.options.display.float_format = '{:.2e}'.format"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Матрица проницаемости:\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "      <th>5</th>\n",
       "      <th>6</th>\n",
       "      <th>7</th>\n",
       "      <th>8</th>\n",
       "      <th>9</th>\n",
       "      <th>10</th>\n",
       "      <th>11</th>\n",
       "      <th>12</th>\n",
       "      <th>13</th>\n",
       "      <th>14</th>\n",
       "      <th>15</th>\n",
       "      <th>16</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.43e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>0</td>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.43e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.43e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>3.43e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>0</td>\n",
       "      <td>3.43e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.43e-14</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.43e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.43e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.44e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>3.43e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.43e-14</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>4.44e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.44e-14</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.44e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.44e-14</td>\n",
       "      <td>5.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.44e-14</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.00e-14</td>\n",
       "      <td>3.43e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>3.43e-14</td>\n",
       "      <td>4.00e-14</td>\n",
       "      <td>4.44e-14</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.44e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.44e-14</td>\n",
       "      <td>5.00e-14</td>\n",
       "      <td>4.44e-14</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.44e-14</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>4.44e-14</td>\n",
       "      <td>4.00e-14</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "         1        2        3        4        5        6        7        8   \\\n",
       "1  3.00e-14 3.00e-14        0        0 3.00e-14        0        0        0   \n",
       "2  3.00e-14 3.00e-14 3.00e-14        0        0 3.43e-14        0        0   \n",
       "3         0 3.00e-14 3.00e-14 3.00e-14        0        0 3.43e-14        0   \n",
       "4         0        0 3.00e-14 3.00e-14        0        0        0 3.43e-14   \n",
       "5  3.00e-14        0        0        0 3.00e-14 3.43e-14        0        0   \n",
       "6         0 3.43e-14        0        0 3.43e-14 4.00e-14 4.00e-14        0   \n",
       "7         0        0 3.43e-14        0        0 4.00e-14 4.00e-14 4.00e-14   \n",
       "8         0        0        0 3.43e-14        0        0 4.00e-14 4.00e-14   \n",
       "9         0        0        0        0 3.00e-14        0        0        0   \n",
       "10        0        0        0        0        0 4.00e-14        0        0   \n",
       "11        0        0        0        0        0        0 4.00e-14        0   \n",
       "12        0        0        0        0        0        0        0 4.44e-14   \n",
       "13        0        0        0        0        0        0        0        0   \n",
       "14        0        0        0        0        0        0        0        0   \n",
       "15        0        0        0        0        0        0        0        0   \n",
       "16        0        0        0        0        0        0        0        0   \n",
       "\n",
       "         9        10       11       12       13       14       15       16  \n",
       "1         0        0        0        0        0        0        0        0  \n",
       "2         0        0        0        0        0        0        0        0  \n",
       "3         0        0        0        0        0        0        0        0  \n",
       "4         0        0        0        0        0        0        0        0  \n",
       "5  3.00e-14        0        0        0        0        0        0        0  \n",
       "6         0 4.00e-14        0        0        0        0        0        0  \n",
       "7         0        0 4.00e-14        0        0        0        0        0  \n",
       "8         0        0        0 4.44e-14        0        0        0        0  \n",
       "9  3.00e-14 3.43e-14        0        0 3.00e-14        0        0        0  \n",
       "10 3.43e-14 4.00e-14 4.00e-14        0        0 4.00e-14        0        0  \n",
       "11        0 4.00e-14 4.00e-14 4.44e-14        0        0 4.44e-14        0  \n",
       "12        0        0 4.44e-14 5.00e-14        0        0        0 4.44e-14  \n",
       "13 3.00e-14        0        0        0 3.00e-14 3.43e-14        0        0  \n",
       "14        0 4.00e-14        0        0 3.43e-14 4.00e-14 4.44e-14        0  \n",
       "15        0        0 4.44e-14        0        0 4.44e-14 5.00e-14 4.44e-14  \n",
       "16        0        0        0 4.44e-14        0        0 4.44e-14 4.00e-14  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# матрица проницаемости; заполняем её значениями\n",
    "# заполнение значениями выполнено в общем виде, чтобы была возможность экспериментировать с пластами произвольных размеров\n",
    "# текущее ограничение: длина пласта должна быть равна его ширине\n",
    "k_matrix = np.full((mesh_x_count ** 2, mesh_y_count ** 2), 0.0)\n",
    "k_matrix += np.diag(k.flatten())\n",
    "\n",
    "for i in range(mesh_x_count * mesh_y_count):\n",
    "    if not (i + 1) % mesh_x_count == 0:\n",
    "        k_matrix[i, i + 1] = 2 * dx / (dx / k[i // mesh_x_count, i % mesh_x_count + 1] + dx / k[i // mesh_x_count, i % mesh_x_count])\n",
    "\n",
    "for i in range(1, mesh_x_count * mesh_y_count):\n",
    "    if not i % mesh_x_count == 0:\n",
    "        k_matrix[i, i - 1] = 2 * dx / (dx / k[i // mesh_x_count, i % mesh_x_count - 1] + dx / k[i // mesh_x_count, i % mesh_x_count])\n",
    "\n",
    "for i in range(mesh_x_count * mesh_y_count - mesh_x_count):\n",
    "        k_matrix[i, i + mesh_x_count] = 2 * dx / (dx / k[i // mesh_x_count + 1, i % mesh_x_count] + dx / k[i // mesh_x_count, i % mesh_x_count])\n",
    "\n",
    "for i in range(mesh_x_count, mesh_x_count * mesh_y_count):\n",
    "        k_matrix[i, i - mesh_x_count] = 2 * dx / (dx / k[i // mesh_x_count - 1, i % mesh_x_count] + dx / k[i // mesh_x_count, i % mesh_x_count])\n",
    "\n",
    "if display_matrices:\n",
    "    print('Матрица проницаемости:')\n",
    "    df = pd.DataFrame(k_matrix)\n",
    "    df.index += 1\n",
    "    df.columns += 1\n",
    "    for i in range(1, df.shape[0] + 1):\n",
    "        for j in range(1, df.shape[1] + 1):\n",
    "            if df.at[i,j] == 0:\n",
    "                df.loc[i,j] = '0'\n",
    "    display(df)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Матрица проводимости:\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "      <th>5</th>\n",
       "      <th>6</th>\n",
       "      <th>7</th>\n",
       "      <th>8</th>\n",
       "      <th>9</th>\n",
       "      <th>10</th>\n",
       "      <th>11</th>\n",
       "      <th>12</th>\n",
       "      <th>13</th>\n",
       "      <th>14</th>\n",
       "      <th>15</th>\n",
       "      <th>16</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "         1        2        3        4        5        6        7        8   \\\n",
       "1         0 7.50e-16        0        0 7.50e-16        0        0        0   \n",
       "2  7.50e-16        0 7.50e-16        0        0 8.57e-16        0        0   \n",
       "3         0 7.50e-16        0 7.50e-16        0        0 8.57e-16        0   \n",
       "4         0        0 7.50e-16        0        0        0        0 8.57e-16   \n",
       "5  7.50e-16        0        0        0        0 8.57e-16        0        0   \n",
       "6         0 8.57e-16        0        0 8.57e-16        0 1.00e-15        0   \n",
       "7         0        0 8.57e-16        0        0 1.00e-15        0 1.00e-15   \n",
       "8         0        0        0 8.57e-16        0        0 1.00e-15        0   \n",
       "9         0        0        0        0 7.50e-16        0        0        0   \n",
       "10        0        0        0        0        0 1.00e-15        0        0   \n",
       "11        0        0        0        0        0        0 1.00e-15        0   \n",
       "12        0        0        0        0        0        0        0 1.11e-15   \n",
       "13        0        0        0        0        0        0        0        0   \n",
       "14        0        0        0        0        0        0        0        0   \n",
       "15        0        0        0        0        0        0        0        0   \n",
       "16        0        0        0        0        0        0        0        0   \n",
       "\n",
       "         9        10       11       12       13       14       15       16  \n",
       "1         0        0        0        0        0        0        0        0  \n",
       "2         0        0        0        0        0        0        0        0  \n",
       "3         0        0        0        0        0        0        0        0  \n",
       "4         0        0        0        0        0        0        0        0  \n",
       "5  7.50e-16        0        0        0        0        0        0        0  \n",
       "6         0 1.00e-15        0        0        0        0        0        0  \n",
       "7         0        0 1.00e-15        0        0        0        0        0  \n",
       "8         0        0        0 1.11e-15        0        0        0        0  \n",
       "9         0 8.57e-16        0        0 7.50e-16        0        0        0  \n",
       "10 8.57e-16        0 1.00e-15        0        0 1.00e-15        0        0  \n",
       "11        0 1.00e-15        0 1.11e-15        0        0 1.11e-15        0  \n",
       "12        0        0 1.11e-15        0        0        0        0 1.11e-15  \n",
       "13 7.50e-16        0        0        0        0 8.57e-16        0        0  \n",
       "14        0 1.00e-15        0        0 8.57e-16        0 1.11e-15        0  \n",
       "15        0        0 1.11e-15        0        0 1.11e-15        0 1.11e-15  \n",
       "16        0        0        0 1.11e-15        0        0 1.11e-15        0  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# матрица проводимости\n",
    "trans_matrix = k_matrix / (mu * b_liq * dx**2)\n",
    "trans_matrix -= np.diag(np.diag(trans_matrix))\n",
    "\n",
    "if display_matrices:\n",
    "    print('Матрица проводимости:')\n",
    "    df = pd.DataFrame(trans_matrix)\n",
    "    df.index += 1\n",
    "    df.columns += 1\n",
    "    for i in range(1, df.shape[0] + 1):\n",
    "        for j in range(1, df.shape[1] + 1):\n",
    "            if df.at[i,j] == 0:\n",
    "                df.loc[i,j] = '0'\n",
    "    display(df)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Матрица СЛАУ:\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "      <th>5</th>\n",
       "      <th>6</th>\n",
       "      <th>7</th>\n",
       "      <th>8</th>\n",
       "      <th>9</th>\n",
       "      <th>10</th>\n",
       "      <th>11</th>\n",
       "      <th>12</th>\n",
       "      <th>13</th>\n",
       "      <th>14</th>\n",
       "      <th>15</th>\n",
       "      <th>16</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>-1.96e-15</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>-2.82e-15</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>-2.82e-15</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>-2.07e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>-2.82e-15</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>-4.17e-15</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>-4.32e-15</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>-3.43e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>-2.82e-15</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>-4.32e-15</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>-4.68e-15</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>-3.79e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>7.50e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>-2.07e-15</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.00e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>8.57e-16</td>\n",
       "      <td>-3.43e-15</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>-3.79e-15</td>\n",
       "      <td>1.11e-15</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>1.11e-15</td>\n",
       "      <td>-2.68e-15</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "          1         2         3         4         5         6         7   \\\n",
       "1  -1.96e-15  7.50e-16         0         0  7.50e-16         0         0   \n",
       "2   7.50e-16 -2.82e-15  7.50e-16         0         0  8.57e-16         0   \n",
       "3          0  7.50e-16 -2.82e-15  7.50e-16         0         0  8.57e-16   \n",
       "4          0         0  7.50e-16 -2.07e-15         0         0         0   \n",
       "5   7.50e-16         0         0         0 -2.82e-15  8.57e-16         0   \n",
       "6          0  8.57e-16         0         0  8.57e-16 -4.17e-15  1.00e-15   \n",
       "7          0         0  8.57e-16         0         0  1.00e-15 -4.32e-15   \n",
       "8          0         0         0  8.57e-16         0         0  1.00e-15   \n",
       "9          0         0         0         0  7.50e-16         0         0   \n",
       "10         0         0         0         0         0  1.00e-15         0   \n",
       "11         0         0         0         0         0         0  1.00e-15   \n",
       "12         0         0         0         0         0         0         0   \n",
       "13         0         0         0         0         0         0         0   \n",
       "14         0         0         0         0         0         0         0   \n",
       "15         0         0         0         0         0         0         0   \n",
       "16         0         0         0         0         0         0         0   \n",
       "\n",
       "          8         9         10        11        12        13        14  \\\n",
       "1          0         0         0         0         0         0         0   \n",
       "2          0         0         0         0         0         0         0   \n",
       "3          0         0         0         0         0         0         0   \n",
       "4   8.57e-16         0         0         0         0         0         0   \n",
       "5          0  7.50e-16         0         0         0         0         0   \n",
       "6          0         0  1.00e-15         0         0         0         0   \n",
       "7   1.00e-15         0         0  1.00e-15         0         0         0   \n",
       "8  -3.43e-15         0         0         0  1.11e-15         0         0   \n",
       "9          0 -2.82e-15  8.57e-16         0         0  7.50e-16         0   \n",
       "10         0  8.57e-16 -4.32e-15  1.00e-15         0         0  1.00e-15   \n",
       "11         0         0  1.00e-15 -4.68e-15  1.11e-15         0         0   \n",
       "12  1.11e-15         0         0  1.11e-15 -3.79e-15         0         0   \n",
       "13         0  7.50e-16         0         0         0 -2.07e-15  8.57e-16   \n",
       "14         0         0  1.00e-15         0         0  8.57e-16 -3.43e-15   \n",
       "15         0         0         0  1.11e-15         0         0  1.11e-15   \n",
       "16         0         0         0         0  1.11e-15         0         0   \n",
       "\n",
       "          15        16  \n",
       "1          0         0  \n",
       "2          0         0  \n",
       "3          0         0  \n",
       "4          0         0  \n",
       "5          0         0  \n",
       "6          0         0  \n",
       "7          0         0  \n",
       "8          0         0  \n",
       "9          0         0  \n",
       "10         0         0  \n",
       "11  1.11e-15         0  \n",
       "12         0  1.11e-15  \n",
       "13         0         0  \n",
       "14  1.11e-15         0  \n",
       "15 -3.79e-15  1.11e-15  \n",
       "16  1.11e-15 -2.68e-15  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# оставшийся множитель перед производной по времени\n",
    "mult_before_time_derivative = m * compr / dt\n",
    "\n",
    "final_reservoir_matrix = trans_matrix.copy()\n",
    "\n",
    "for i in range(mesh_x_count * mesh_y_count):\n",
    "    final_reservoir_matrix[i, i] = -sum(final_reservoir_matrix[:, i]) - mult_before_time_derivative\n",
    "\n",
    "if display_matrices:\n",
    "    print('Матрица СЛАУ:')\n",
    "    df = pd.DataFrame(final_reservoir_matrix)\n",
    "    df.index += 1\n",
    "    df.columns += 1\n",
    "    for i in range(1, df.shape[0] + 1):\n",
    "        for j in range(1, df.shape[1] + 1):\n",
    "            if df.at[i,j] == 0:\n",
    "                df.loc[i,j] = '0'\n",
    "    display(df)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "# список, в который будут добавляться поля давлений на каждом временном шаге\n",
    "p = []\n",
    "p.append(np.full((mesh_x_count, mesh_y_count), p_initial))\n",
    "\n",
    "for time_step in range(mesh_t_count):\n",
    "    # вектор правых частей\n",
    "    d = -p[-1].flatten() * mult_before_time_derivative\n",
    "    # добавление скважин\n",
    "    for i, inj_cell in enumerate(inj_cells):\n",
    "        d[inj_cell - 1] -= q_inj[i] / (dx * dy * h)\n",
    "    for i, prod_cell in enumerate(prod_cells):\n",
    "        d[prod_cell - 1] += q_prod[i] / (dx * dy * h)\n",
    "    # решение СЛАУ\n",
    "    solution = np.matmul(np.linalg.inv(final_reservoir_matrix), d)\n",
    "    # добавление найденного поля давлений в список\n",
    "    p.append(solution.reshape(mesh_x_count, mesh_y_count))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Шаг по времени: 0 \n",
      "Распределение давления, МПа\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>20.000</td>\n",
       "      <td>20.000</td>\n",
       "      <td>20.000</td>\n",
       "      <td>20.000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>20.000</td>\n",
       "      <td>20.000</td>\n",
       "      <td>20.000</td>\n",
       "      <td>20.000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>20.000</td>\n",
       "      <td>20.000</td>\n",
       "      <td>20.000</td>\n",
       "      <td>20.000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>20.000</td>\n",
       "      <td>20.000</td>\n",
       "      <td>20.000</td>\n",
       "      <td>20.000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       1      2      3      4\n",
       "1 20.000 20.000 20.000 20.000\n",
       "2 20.000 20.000 20.000 20.000\n",
       "3 20.000 20.000 20.000 20.000\n",
       "4 20.000 20.000 20.000 20.000"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUAAAAEICAYAAAAumy2rAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAAAeWUlEQVR4nO3df7BfdX3n8eeLEKJLUHACGJKLONPUEW0J1kmxmW0j/grZVhwrHdgWqdpmdGAKLjpl0MXquI51XaoUxpgFCq6IqwvRVKMYFReZFSRkAiFeKFmGlruXmgJKCKKYm9f+cT5Xv/nm+733e/M9uff75bwenTP3nM/5nPP53Nj75vM5n/M5H9kmIqKJDpvrCkREzJUEwIhorATAiGisBMCIaKwEwIhorATAiGisBMCIaKwEwFkg6WFJr285XiVprC3PdZI+Ovu1i+lIWi7pbkk/Kz+Xz3Wdoh4JgBFTkHQE8FXg88AxwPXAV0t6DLkEwAEl6cuS/lXSk5Juk/SKlnPXSXpW0p6W7YUd7vHnkiba8u2TtKqcXyDpU5LGy/YpSQu61Kf1XrslfVfSknJulaQxSZdKeqy0eP+05doFkj4p6V8k/VjSOknPb7nWkj7dkv/kkvb5lrTTJP0fST+VdM/k71DOfU/SX7Qcv17Swy3Hv2qBS1pY6nB7j/9TrAIOBz5l+xe2rwAEnN7j9THAEgAH1zeAZcBxwFbghrbzn7C9sGV7sst9ftCaDxhvOfcB4DRgOXAKsAL44BR1+kG5x3HAL4D3tpx7MbAIWAKcB6yX9LJy7m+B3yzl/EbJc1nLtf8GnNESfP8CGJ08WQLt14GPAi8C3gfcJOnYKerazfuBX7YmSPqapEu65H8FcK/3nzN6b0mPIZcAOKBsX2v7Kdu/AP4GOKVTK69Pfwp8xPYu2/8GfBg4t4frDivb423p/7m0kv43VcD6E0kC/hJ4r+0nbD8FfAw4u+W6Z4FNwFtL1/IM4Cst5/8M2GR7k+19tjcDW4A1M/llJR0PvAu4vDXd9h/a/niXyxYC7f9xeRI4aiZlx2A6fK4rEAeSNA/4L8BZwLHAvnJqEQf+MfbjBOCfW47/uaR1c5qknwL/DngY+POWcz+x/XSHex1b8t9dxUKg6kLOa7v31cCngAngm1RBcdJLgLMk/VFL2nzg1pbjKyR9suwfDjzWof5/A/w98ETnX6+jPcAL2tJeADw1g3vEgEoLcDD9R+BM4PXAC4GTSrq6XXCQxqmCy6QT2b+L3O4O20cDz6MaFLiu5dwxko7scK/HgGeAV9g+umwvLF3pX7F9H1Wg/CBVMGz1CPA/Wq4/2vaRba22v5o8B7ylQ91/E3gTcMUUv18nO4DfVkv0Bn67pMeQSwCcPfMlPU/S84AjACaPS1pri+goqmdsj1MFhY8dojrdCHxQ0rGSFlE9l/v8NNcAmKql1v4M7sOSjpD074E/BL5sex/w34G/k3QcVM/0JL2pw30/Bnzbdntw+TzwR5LeJGle+TdbJWlpz79pFVg/YvuZGVwD8D2q3/WvymDOBSX9uzO8TwygBMDZs4mqJfQMcAvVQMAzLdufteT9HFUX8v8BPwLuOER1+ijVs7R7ge1Ugy1TvYv4Gkl7qLrhbwUuaDn3r8BPqFp9NwDvtn1/OffXwE7gDkm7gW8DL6ON7a/Z/k8d0h+hahFfSjVg8gjVYMZM/v/3cap/1wNI+oakSzuds/0sVYvy7cBPgXcCbynpMeSUD6JGv8orKZ+3PZMWWcScSwswIhqrrwAo6UWSNkt6sPw8pku+hyVtl7RN0pZ+yoyI5yZJI5JulTQqaYekC0v6f5V0v6R7JW2QdHSX61dLekDSzine69z/mn66wJI+ATxh++OlwGNs/3WHfA8Dr7bd6dWEiAgkLQYW294q6Sjgbqrnr0uB79reK+lvAdrjTHl17J+ANwBjwF3AObZ/NFWZ/XaBz6SaG0n5+ZY+7xcRDWX7Udtby/5TVLOBltj+lu29JdsdVAGx3Qpgp+2HygDVF6ni05T6fRH6eNuPTlZ+8jWHDgx8S5KBz9pe3+2GktYCawF0xBG/M//4breMiH7tfeIJJvY83df7pW967ZF+/ImJnvLefe8vdgA/b0la3ykeSDoJOBW4s+3UO4H/2eHWS6jeDpg0BvzudPWZNgBK+jbVPM92H5ju2hYrbY+XALlZ0v22b+uUsfxjrAdYcOKIT7j4ohkUExEzMf7fPtX3PR5/YoIf3nJiT3nnLX7w57ZfPVUeSQuBm4CLbO9uSf8AsJcD58VD50kC0z7fmzYA2n59t3PlqxqLS+tvMbCryz3Gy89dkjZQNVc7BsCIGC4G9v1qtmZ/JM2nCn432L65Jf08qpfrX+fOAxdjwEjL8VKmntUE9P8McCPVlz8oP7/ankHSkeWBJmWq1BuB+/osNyIGhDG/9ERP21TKdMNrgFHbl7ekr6Z6mf7Ntn/W5fK7gGWSXlo+qHE2VXyaUr8B8OPAGyQ9SDX68vFS4RMkbSp5jgdul3QP8EPg67a/2We5ETFA9vX4f9NYSfU1otPLK3PbJK0BrqSaHrq5pK2D/eNMGSS5gGqW1SjwpQ5TKg/Q1yCI7ceB13VIH6d8qsj2Q1TfmouI5yBjJmqYUWb7djo/y9vUIW2/OFOON3XL200+hxURfds3/XjDQEoAjIi+VJ8GSgCMiIZKCzAiGsnAL4f0q1IJgBHRF+N0gSOioQwTwxn/EgAjoj/VTJDhlAAYEX0SE7Wv1zU7EgAjoi/VIEgCYEQ0UPUeYAJgRDTUvrQAI6KJ0gKMiMYyYmJIF5hMAIyIvqULHBGNZMSznjfX1TgoCYAR0ZfqReh0gSOioTIIEhGNZIsJD2cLsJZaS1ot6QFJOyVd0uG8JF1Rzt8r6VV1lBsRg2Ef6mmbiqQRSbdKGpW0Q9KFJf2scrxPUtclNSU9LGl7WTdkSy/17rsFKGkecBXVokhjwF2SNtr+UUu2M4BlZftd4DP0sGhxRAy+ahCkls7kXuBi21vLSpJ3S9pMtYrkW4HP9nCP19p+rNcC66j1CmBnWfwISV8EzgRaA+CZwOfKep53SDp6cj3hGsqPiDlU1yBIiQePlv2nJI0CS2xvBqhWzaxXHV3gJcAjLcdjJW2meSJiSE1YPW29knQScCpw5wyqYeBbku6WtLaXC+poAXb6rdo/j9hLnipjVfG1APOOOaa/mkXEITfDmSCL2p7Prbe9vjWDpIXATcBFtnfPoCorbY9LOo5qDeH7bd821QV1BMAxYKTleCkwfhB5ACj/GOsBFpw4MqTfmY1oln29jwI/ZnuqgYz5VMHvBts3z6QOZZ1gbO+StIHq8dyUAbCOLvBdwDJJL5V0BHA2sLEtz0bg7WU0+DTgyTz/i3huqD6GcFhP21RUPeS7Bhi1fflM6iDpyDJwgqQjgTdSDZ5Mqe8WoO29ki4AbgHmAdfa3iHp3eX8OqrV2tcAO4GfAe/ot9yIGAxG/LKeqXArgXOB7ZK2lbRLgQXA3wPHAl+XtM32mySdAFxtew1wPLChDJQcDnzB9jenK7CWsWvbm6iCXGvaupZ9A+fXUVZEDBabWl6Etn07nccLADZ0yD9O1bCivIVyykzLzEyQiOjT9C85D6oEwIjoi6mnBTgXEgAjom/5IGpENJJRPogaEc1ULYs5nKFkOGsdEQMkC6NHREOZGc0EGSgJgBHRt7QAI6KRbKUFGBHNVA2CZFW4iGik4V0TJAEwIvpSDYLkGWBENFRmgkREI2UmSEQ0Wh2LIs2FBMCI6IsNv9yXABgRDVR1gRMAI6KhMhMkIhppmF+DqaXdKmm1pAck7ZR0SYfzqyQ9KWlb2S6ro9yIGARVF7iXbcq7SCOSbpU0KmmHpAtL+lnleJ+kqZbUnDIOddJ3C1DSPOAq4A1U6//eJWmj7R+1Zf2+7T/st7yIGDw1rQmyF7jY9tayxOXdkjZTLW/5VuCz3S6cQRzaTx1d4BXAzrIqE5K+CJwJTFlwRDw3VKPA/c8FLmuFP1r2n5I0CiyxvRmgLHnZzUHFoTq6wEuAR1qOx0pau9dIukfSNyS9otvNJK2VtEXSlok9T9dQvYg4lCZfhO5lAxZN/n2XbW2ne0o6CTgVuLPHavQah/ZTRwuwU1h22/FW4CW290haA3wFWNbpZrbXA+sBFpw40n6fiBhAM+gCP2a763M8AEkLgZuAi2zv7vG+vcShA9TRAhwDRlqOlwLj+9XC3m17T9nfBMyXtKiGsiNijk2OAvfYApySpPlUwe8G2zfPoBrTxqFO6giAdwHLJL1U0hHA2cDG1gySXqzSgZe0opT7eA1lR8QAqGkUWMA1wKjty2dYhWnjUCd9d4Ft75V0AXALMA+41vYOSe8u59cBbwPeI2kv8Axwtu10byOeA2yxt56ZICuBc4HtkraVtEuBBcDfA8cCX5e0zfabJJ0AXG17Tbc4NF2BtbwIXbq1m9rS1rXsXwlcWUdZETF46ngR2vbtdH6WB7ChQ/5xYE3L8QFxaDqZCRIRfRnmmSAJgBHRtwTAiGikfBA1Ihqtpqlwsy4BMCL6YsPefBA1IpoqXeCIaKQ8A4yIRnMCYEQ0VQZBIqKR7DwDjIjGEhMZBY6IpsozwIhopMwFjojmcvUccBglAEZE3zIKHBGN5AyCRESTpQscEY01rKPAtbRbJV0raZek+7qcl6QrJO2UdK+kV9VRbkTMPbsKgL1sU5E0IulWSaOSdki6sKS/SNJmSQ+Wn8d0uf5hSdslbZO0pZe619Vxvw5YPcX5M6jWAV4GrAU+U1O5ETEAaloWcy9wse2XA6cB50s6GbgE+I7tZcB3ynE3r7W9fLq1hyfVEgBt3wY8MUWWM4HPuXIHcLSkxXWUHRFzz+5tm/oeftT21rL/FDAKLKGKH9eXbNcDb6mr3rM1dLMEeKTleKykHUDSWklbJG2Z2PP0rFQuIg6eEfv2HdbTBiya/Psu29pO95R0EnAqcCdwvO1HoQqSwHFdqwLfknR3t/u2m61BkE5t347/PbC9HlgPsODEkSEdW4polhn8oT42XfdU0kLgJuAi27ur9dJ7stL2uKTjgM2S7i+9065mqwU4Boy0HC8Fxmep7Ig4lGoaBAGQNJ8q+N1g++aS/OPJR2bl566O1ajWCcb2Lqp1hFdMV95sBcCNwNvLaPBpwJOTTdqIeA5wj9sUVDX1rgFGbV/ecmojcF7ZPw/4aodrj5R01OQ+8Eag41sprWrpAku6EVhF1b8fAz4EzAewvY5qtfY1wE7gZ8A76ig3IgZDTe8BrgTOBbZL2lbSLgU+DnxJ0ruAfwHOApB0AnC17TXA8cCG0l0+HPiC7W9OV2AtAdD2OdOcN3B+HWVFxGAxsG9f/wHQ9u10Hi8AeF2H/ONUDStsPwScMtMyMxMkIvpjYEhngiQARkTfMhc4IporATAimqm3V1wGUQJgRPQvLcCIaCSDaxgFngsJgBFRgwTAiGiqdIEjorESACOikfIidEQ0WV6EjojmyihwRDSV0gKMiEbq4Vt/gyoBMCL6pAyCRESDpQUYEY21b64rcHASACOiP0P8HmAtiyJJulbSLkkdFyGRtErSk5K2le2yOsqNiMEg97ZNeQ9pRNKtkkYl7ZB0YUl/kaTNkh4sP4/pcv1qSQ9I2inpkl7qXdeqcNcBq6fJ833by8v2kZrKjYhBUMOqcMBe4GLbLwdOA86XdDJwCfAd28uA75Tj/UiaB1wFnAGcDJxTrp1SLQGwLD78RB33iohmsv2o7a1l/ylgFFgCnAlcX7JdD7ylw+UrgJ22H7L9LPDFct2UZmtdYIDXSLpH0jckvaJbJklrJW2RtGViz9OzWL2IOFgz6AIvmvz7LtvajveTTgJOBe4Ejp9cR7z8PK7DJUuAR1qOx0ralGZrEGQr8BLbeyStAb4CLOuU0fZ6YD3AghNHhnRwPaJBzEymwj1m+9VTZZC0ELgJuMj27rLW73Q6ZZo2fsxKC9D2btt7yv4mYL6kRbNRdkTMgnqeASJpPlXwu8H2zSX5x5IWl/OLgV0dLh0DRlqOlwLj05U3KwFQ0otVwrikFaXcx2ej7Ig49GoaBRZwDTBq+/KWUxuB88r+ecBXO1x+F7BM0kslHQGcXa6bUi1dYEk3Aquo+vdjwIeA+QC21wFvA94jaS/wDHC2Pawf0ImIA9Tz17wSOBfYLmlbSbsU+DjwJUnvAv4FOAtA0gnA1bbX2N4r6QLgFmAecK3tHdMVWEsAtH3ONOevBK6so6yIGEA1BEDbt9N9cZHXdcg/DqxpOd4EbJpJmZkJEhF96aV7O6gSACOif/kgakQ0VVqAEdFcCYAR0Uh5BhgRjZYAGBFNpSH9IOpsfgwhImKgpAUYEf1LFzgiGimDIBHRaAmAEdFYCYAR0URieEeBEwAjoj95BhgRjZYAGBGNlQAYEU2VLnBENNeQBsC+p8JJGpF0q6RRSTskXdghjyRdIWmnpHslvarfciNiQLgaBe5lm46kayXtknRfS9opkn4gabukf5T0gi7XPlzybJO0pZeq1zEXeC9wse2XA6cB50s6uS3PGVTrAC8D1gKfqaHciBgUNS2LCVwHrG5Luxq4xPZvARuA909x/WttL59u7eFJfQdA24/a3lr2nwJGOXBF9jOBz7lyB3D05DqfETH86lgWE8D2bcATbckvA24r+5uBP66r3rV+DUbSScCpwJ1tp5YAj7Qcj3FgkJy8x1pJWyRtmdjzdJ3Vi4hDpfcW4KLJv++yre3h7vcBby77Z7H/AujttfiWpLt7vG99gyCSFlKt6H6R7d3tpztc0vG/B7bXA+sBFpw4MqSPViMapPfuLcBjvXZPW7wTuELSZVSLnT/bJd9K2+OSjgM2S7q/tCi7qmth9PlUwe8G2zd3yDLG/lF7KTBeR9kRMbfEoX0Nxvb9wBsBJP0m8B+65BsvP3dJ2gCs4Ndd547qGAUWcA0wavvyLtk2Am8vo8GnAU/afrTfsiNiMNT1DLDjvasWHZIOAz4IrOuQ50hJR03uUwXM+9rztaujBbgSOBfYLmlbSbsUOBHA9jqq1drXADuBnwHvqKHciBgUNbUAJd0IrKJ6VjgGfAhYKOn8kuVm4B9K3hOAq22vAY4HNlTtMQ4HvmD7m9OV13cAtH07nZ/xteYxcP5UeSJiiNUUAG2f0+XUpzvkHadqWGH7IeCUmZaXmSAR0Z98DSYiGi0BMCKaKh9EjYjGShc4IpppZi9CD5QEwIjoXwJgRDTRoZ4JciglAEZE37RvOCNgAmBE9CfPACOiydIFjojmSgCMiKZKCzAimisBMCIayZkKFxENlfcAI6LZPJwRMAEwIvqWFmBENNMQvwhdx6JII5JulTQqaYekCzvkWSXpSUnbynZZv+VGxODQvt62ae8jXStpl6T7WtJOkfQDSdsl/aOkF3S5drWkByTtlHRJL/WuY2H0vcDFtl8OnAacL+nkDvm+b3t52T5SQ7kRMSDqCoDAdcDqtrSrgUts/xawAXj/AeVL84CrgDOAk4FzusSh/fQdAG0/antr2X8KGAWW9HvfiBgSphoE6WWb7lbVQuZPtCW/jF+v77sZ+OMOl64Adtp+yPazwBeBM6crr44W4K9IOgk4Fbizw+nXSLpH0jckvWKKe6yVtEXSlok9T9dZvYg4RGawLvCiyb/vsq3t4fb3AW8u+2cBIx3yLAEeaTkeo4eGWG2DIJIWAjcBF9ne3XZ6K/AS23skrQG+AizrdB/b64H1AAtOHBnSR6sRDdP7X+pjtl89w7u/E7iijB1sBJ7tkKfT0rzT1qqWFqCk+VTB7wbbNx9QC3u37T1lfxMwX9KiOsqOiLk1+SJ0jy3AGbN9v+032v4d4Ebg/3bINsb+LcOlwPh0965jFFjANcCo7cu75HlxyYekFaXcx/stOyIGgI329bYdDEnHlZ+HAR8E1nXIdhewTNJLJR0BnE3VWpxSHV3glcC5wHZJ20rapcCJALbXAW8D3iNpL/AMcLY9pK+OR8SBavprlnQjsIrqWeEY8CFgoaTzS5abgX8oeU8Arra9xvZeSRcAtwDzgGtt75iuvL4DoO3b6dz/bs1zJXBlv2VFxGCqayaI7XO6nPp0h7zjwJqW403AppmUl5kgEdEfA1kTJCIaazjjXwJgRPQvH0OIiMbKspgR0UxD/DWYBMCI6Ev1IvRwRsAEwIjoX9YEiYimSgswIpopzwAjorkOfp7vXEsAjIj+pQscEY2UhdEjotHSAoyIxhrO+JcAGBH9077h7AMnAEZEf0xehI6IZhLOi9AR0WBDGgDrWBTpeZJ+WNb83SHpwx3ySNIVknZKulfSq/otNyIGSE0Lo0u6VtIuSfe1pC2XdIekbWUt4RVdrn1Y0vbJfL1Uu45lMX8BnG77FGA5sFrSaW15zqBaB3gZsBb4TA3lRsQgmHwG2Ms2veuA1W1pnwA+bHs5cFk57ua1tpf3uvZwHYsiGdhTDueXrT3Unwl8ruS9Q9LRkhbbfrTf8iNi7tU1Cmz7NkkntScDLyj7L6SH9X57VcszQEnzgLuB3wCusn1nW5YlwCMtx2MlLQEwYuj11r3tw0XALZI+SdVr/b3uFeFbkgx81vb66W5cRxcY2xOleboUWCHplW1ZOi2b2fFfTNLa0s/fMrHn6TqqFxGHkpnJM8BFk3/fZVvbQwnvAd5rewR4L3BNl3wrbb+K6pHb+ZJ+f7ob1zoKbPunkr5H1Ye/r+XUGDDScryULs3YErXXAyw4cWQ4h5Yimqb3HvBjvT6fa3EecGHZ/zJwdadMZZ1gbO+StAFYAdw21Y3rGAU+VtLRZf/5wOuB+9uybQTeXkaDTwOezPO/iOcO2T1tB2kc+IOyfzrw4AHlS0dKOmpyH3gj+zfCOqqjBbgYuL48BzwM+JLtr0l6N4DtdVSrta8BdgI/A95RQ7kRMShqegYo6UZgFVVXeQz4EPCXwKclHQ78nOpNEiSdAFxtew1wPLBBElRx7Qu2vzldeXWMAt8LnNohfV3LvoHz+y0rIgaQDRO1jQKf0+XU73TIO07VsML2Q8ApMy0vM0Eion9DOhMkATAi+pcAGBGNZCBrgkREMxk8nN/DSgCMiP6Y2gZBZlsCYET0L88AI6KxEgAjopkO+ccQDpkEwIjoj4EsihQRjZUWYEQ0U31T4WZbAmBE9MfgvAcYEY2VmSAR0Vh5BhgRjWRnFDgiGiwtwIhoJuOJibmuxEFJAIyI/uRzWBHRaEP6Gkwdq8I9T9IPJd0jaYekD3fIs0rSk5K2le2yfsuNiMFgwPvc0zYdSddK2iXpvpa05ZLuKLFji6QVXa5dLekBSTslXdJL3etYGP0XwOm2TwGWA6vL0pftvm97edk+UkO5ETEIXD6I2ss2veuo1hVv9Qngw7aXA5eV4/2UVSmvoloU/WTgHEknT1dYHavCGdhTDueXbTgfCETEQalrEMT2bZJOak8GXlD2X0i1TnC7FcDOsjockr4InAn8aKryankGWKLv3cBvAFfZvrNDttdIuoeq8u+zvaPLvdZS1v0E9jx80fseqKOOPVgEPDZLZc2m/F7DZzZ/t5f0e4On+Mkt3/b/WtRj9udJ2tJyvN72+mmuuQi4RdInqXqtv9chzxLgkZbjMeB3p6tMLQHQ9gSwXNLRVIsTv9J266rsW4GX2N4jaQ3wFWBZl3utB6b7B6mdpC22Xz3b5R5q+b2Gz7D9brbbu6x1ew/wXts3SfoT4Brg9W151Klq0924jmeAvy7N/inwPdr68LZ3295T9jcB8yX1+l+MiGi284Cby/6Xqbq77caAkZbjpXTuKu+njlHgY0vLD0nPp4rM97flebEklf0VpdzH+y07IhphHPiDsn868GCHPHcByyS9VNIRwNnAxuluXEcXeDFwfXkOeBjwJdtfk/RuANvrgLcB75G0F3gGOLsMngySWe92z5L8XsPnufy7TUnSjcAqYJGkMeBDwF8Cn5Z0OPBzyhiBpBOAq22vsb1X0gXALcA84Npu4wz7lTd4cSgiYnbU+gwwImKYJABGRGM1PgAezPSZYdBpStFzgaQRSbdKGi1TLy+c6zrVoZcppVG/Rj8DLAM3/wS8gWoY/S7gHNtTvj0+DCT9PtUMnc/ZfuVc16cukhYDi21vlXQU1Qv4bxn2/83KWxJHlndl5wO3AxfavmOOq/ac1vQW4K+mz9h+FpicPjP0bN8GPDHX9aib7Udtby37TwGjVLMAhpormVI6y5oeADtNnxn6P6amKHNGTwU6Tb0cOpLmSdoG7AI2d5lSGjVqegA8qOkzMfckLQRuAi6yvXuu61MH2xPliydLgRWSnjOPLgZV0wPgQU2fiblVnpHdBNxg++bp8g+bblNKo35ND4AHNX0m5k4ZLLgGGLV9+VzXpy69TCmN+jU6ANreC0xOnxmlmsY37fSZYVCmFP0AeJmkMUnvmus61WQlcC5wessXxtfMdaVqsBi4VdK9VP9h3mz7a3Ncp+e8Rr8GExHN1ugWYEQ0WwJgRDRWAmBENFYCYEQ0VgJgRDRWAmBENFYCYEQ01v8HJybH9YP6uCIAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Шаг по времени: 1 \n",
      "Распределение давления, МПа\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>20.000</td>\n",
       "      <td>19.947</td>\n",
       "      <td>19.800</td>\n",
       "      <td>19.399</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>20.053</td>\n",
       "      <td>20.000</td>\n",
       "      <td>19.915</td>\n",
       "      <td>19.807</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>20.200</td>\n",
       "      <td>20.085</td>\n",
       "      <td>20.000</td>\n",
       "      <td>19.943</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>20.601</td>\n",
       "      <td>20.193</td>\n",
       "      <td>20.057</td>\n",
       "      <td>20.000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       1      2      3      4\n",
       "1 20.000 19.947 19.800 19.399\n",
       "2 20.053 20.000 19.915 19.807\n",
       "3 20.200 20.085 20.000 19.943\n",
       "4 20.601 20.193 20.057 20.000"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUAAAAEICAYAAAAumy2rAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAAAdYUlEQVR4nO3dfbAddZ3n8feHEESeBCaAIQni1ETGp+JhqAhDrZUFVIiMWKNswY6ID7NZXXTARWcYxtHRmnWYKYtVBstsFhBYKCwHUCMTZSLCIrWAQAwPISBZx4E7iYYHIQQZYnI/+0f3dQ8n597Tl+577znpz4vquqe7f92/38kN3/y6f0+yTUREG+0y0wWIiJgpCYAR0VoJgBHRWgmAEdFaCYAR0VoJgBHRWgmAEdFaCYDTQNLPJJ3Ysb9Y0khXmisk/fX0ly76kbRc0iOSRiV9YKbLE81JAIzo7z7gvwCrZ7og0awEwAEl6R8k/VzSs5Juk/TGjnNXSNoqaUvH9qoe9/iApO1d6UYlLS7Pv0LSlyRtKLcvSXrFOOXpvNdmST+QNK88t1jSiKQLJD1Z1nj/qOPaV0j6oqTHJP1C0jJJr+y41pK+3JH+DeWxqzuOHSPp/0h6RtJ9Y9+hPHerpD/u2D9R0s869n9TA5e0V1mG26v+Lmx/xfbNwL9VvSaGQwLg4PousBA4kKLmcU3X+b+zvVfH9uw497mjMx2woePcXwDHAEcAhwOLgE9PUKY7ynscCLwIfKLj3KuBOcA84CxguaTDynN/C7yuzOd3yjSf6bj2CeDkjuD7x8C6sZNloP1H4K+B/YFPAtdLOmCCso7nU8CvOw9IulHS+S/jXjHkEgAHlO3LbT9n+0Xgr4DDe9Xyavoj4PO2N9l+AvgccGaF63Ypt6e6jv+l7Rdt/2+KgPUfJAn4T8AnbD9t+zngC8DpHddtBVYCfyhpN+Bk4Fsd598HrLS90vao7VXAPcCSyXxZSQcBHwYu6jxu+xTbF07mXrFz2HWmCxA7kjQL+G/AacABwGh5ag4wXk3v5TgY+JeO/X8pj43nGEnPAHsAPwM+0HHul7af73GvA8r09xaxEAABs7rufSnwJWA78D2KoDjmNcBpkv6g49hs4JaO/YslfbH8vCvwZI/y/xXw98DTvb9etE1qgIPpPwKnAicCrwIOLY9rvAtepg0UwWXMIbz0Ebnbnbb3BXYHrgau6Di3n6Q9e9zrSeAF4I229y23V5WP0r9h+0GKQPlpimDY6XHgf3Vcv6/tPbtqbX8ydg54d4+yvw54B3DxBN8vWiYBcPrMlrS7pN2B3QDG9stjnTWivSnesT1FERS+MEVluhb4tKQDJM2heC93dZ9rAExRU+t+B/c5SbtJ+nfAKcA/2B4F/ifw3yUdCMU7PUnv6HHfLwDft7226/jVwB9IeoekWeWf2WJJ8yt/0yKwft72C5O4hrK8u5W/I/H/f4/5f2cnkF/i9FlJURN6AbiJoiHghY7tfR1pr6J4hPxX4CHgzikq019TvEu7H3iAorFlor6Ix0raQvEY/ofAxzrO/Rz4JUWt7xrgI7YfLs/9GbAeuFPSZuD7wGF0sX2j7f/a4/jjFDXiCygaTB6naMyYzN/fpyj+XHcg6buSLpjg2n+i+B39PrC8/PzWSeQdA0qZEDXqKrukXG17MjWyiBmXGmBEtFatAChpf0mrJD1a/txvnHQ/k/SApDWS7qmTZ0TsnCQtkHSLpHWS1ko6pzxeNc7sK+k6SQ+X9zi2b551HoEl/R3wtO0Ly46k+9n+sx7pfgYcbbtX14SICCTNBebaXi1pb+Beihb9D1AtzlwJ/ND2pWV/0j1sPzNhnjUD4CPAYtsby8LfanuHl9sJgBExWZK+DVxSbhPGGUn7UIzZ/m1PIqjVDYDPlP2uxvZ/aXuH6qmkf6ZoITTwP2wvn+CeS4GlANptt9+bfdCBL7t8g0rbZ7oEU2eXrf3TDKNdn3y+f6Ih9G88z1a/WKt/6Tv+/Z5+6ulqf6nvvf/Fm2yf1C+dpEOB24A3AY/1izOSjqBooX+IYljnvcA5XZ3zd9B3JIik71OM8+z2F/2u7XCc7Q1lP7BVkh62fVuvhGVwXA7wikMW+ODzzp1ENsNh1+eb7s88OPZ6bKZLMDXmLL9jposwJe7yzbXv8dTT2/nRTYdUSjtr7qO/29UOsLy7QiRpL+B64FzbmztGEE1kV+Ao4OO27yon1zgf+Mt+F03I9onjnStn1ZjbUTXdNM49NpQ/N0n6JsWg+54BMCKGi4HR34zW7OtJ20ePd1LSbIrgd43tG8rDVeLMCDBi+65y/zqKADihut1gVlDM/EH589vdCSTtWb7QpBwq9XbgwZr5RsSAMObX3l5pm0g5ccZlwDrbnRNW9I0ztn8OPN4xA9EJFI/DE6obAC8E3ibpUeBt5T6SDpa0skxzEHC7pPuAHwH/aPt7NfONiAEyWvG/Po6jmI3o+LLL3BpJS6gWZwA+Dlwj6X6Kqdf6DiGtNRuM7acoIm338Q2UUxXZ/inFS8mI2AkZs72BEWW2b2f8CT8mjDPl/hpg3MfrXjIdVkTUNspwDqlNAIyIWoqpgRIAI6KlUgOMiFYy8OshnVUqATAiajHOI3BEtJRh+3DGvwTAiKinGAkynBIAI6Imsb3x9bqmRwJgRNRSNIIkAEZECxX9ABMAI6KlRlMDjIg2Sg0wIlrLiO1DusBkAmBE1JZH4IhoJSO2etZMF+NlSQCMiFqKjtB5BI6IlkojSES0ki22ezhrgI2UWtJJkh6RtL5cub37vCRdXJ6/X9JRTeQbEYNhFFXaBk3tGqCkWcBXKBYrGQHulrTCdueKTCcDC8vtLcBXy58RMeSKRpDhfJhsoga4CFhv+6e2twJfB07tSnMqcJULdwL7lut7RsSQG2sEqbINmiZKNA94vGN/pDw22TQRMaS2W5W2QdNEvbXXt+qeHrFKmiKhtBRYCjBrv/3qlSwiplzbR4KMAAs69ucDG15GGgBsLweWA7zikAVDOs9sRLuMtrgV+G5goaTXStoNOB1Y0ZVmBfD+sjX4GOBZ2xsbyDsiZlgxGcIulbZBU7sGaHubpI8BNwGzgMttr5X0kfL8MmAlxQru64FfAR+sm29EDAYjft3moXC2V1IEuc5jyzo+Gzi7ibwiYrDYDG1H6OHsvBMRA2QwOzlXMZxhOyIGhilqgFW2iUhaIOkWSeskrZV0Tnl8f0mrJD1a/hy3e4ikWZJ+LOnGKmVPAIyI2hpqBNkGnGf79cAxwNmS3gCcD9xseyFwc7k/nnOAdVXLnQAYEbUYMepq24T3sTfaXl1+fo4ikM2jGEl2ZZnsSuDdva6XNB94J3Bp1bLnHWBE1FIsi1k5lMyRdE/H/vKy7+9LSDoUOBK4CzhorNuc7Y2SDhzn3l8C/hTYu2phEgAjoqZJLYz+pO2jJ7ybtBdwPXCu7c1S/3tLOgXYZPteSYurFiYBMCJqMc2NBJE0myL4XWP7hvLwLyTNLWt/c4FNPS49DniXpCXA7sA+kq62/b6J8ss7wIiobXtZC+y3TURFVe8yYJ3tizpOrQDOKj+fBXy7+1rbf257vu1DKUaj/aBf8IMEwIioyRaj3qXS1sdxwJnA8ZLWlNsS4ELgbZIepZh39EIASQdLWjn+7frLI3BE1FI0gtQfCmf7dnrPHAVwQo/0GyiG2HYfvxW4tUqeCYARUdPwrgmSABgRtRSNIMM5FC4BMCJqG8SprqpIAIyIWsZGggyjBMCIqG0QFzyqIgEwImqx4dejCYAR0ULFI3ACYES01CTGAg+UBMCIqGWYu8E0Um+VdJKkRyStl7TDZIWSFkt6tmN4y2eayDciBkFjQ+GmXe0aoKRZwFcoxuiNAHdLWmH7oa6kP7R9St38ImLwDOuaIE08Ai8C1tv+KYCkr1PM4NodACdvFzO6x/batxk02xjOJQSreOGgwftXvglP/udjZ7oIU2Lb9XfWvkfRCjycf6eb+Ns6D3i8Y3+kPNbtWEn3SfqupDeOdzNJSyXdI+me7c8930DxImIqNTUl/kxoogbY61u5a3818BrbW8rpbb4FLOx1s3J67OUArzh0fvd9ImIADesjcBM1wBFgQcf+fGBDZwLbm21vKT+vBGZLmtNA3hExw8ZagYexBthEALwbWCjptZJ2o5iNdUVnAkmvLmd7RdKiMt+nGsg7IgZAa1uBbW+T9DHgJmAWcLnttZI+Up5fBrwX+KikbcALwOm283gbsROwxbYBDG5VNNIRunysXdl1bFnH50uAS5rIKyIGzyA+3laRkSARUcswjwRJAIyI2hIAI6KVMiFqRLTasPYDTACMiFps2JYJUSOirfIIHBGtlHeAEdFqTgCMiLYa1kaQ4XxzGREDw25mMgRJCyTdImmdpLWSzimP7y9plaRHy5/7Vb22nwTAiKhJbB/dpdLWxzbgPNuvB44Bzpb0BuB84GbbC4Gby/2q104oATAiarNVaZv4Ht5oe3X5+TlgHcXkyqcCV5bJrgTePYlrJ5R3gBFRyyTHAs+RdE/H/vJyEuSXkHQocCRwF3CQ7Y1QBDpJB06UQde1E0oAjIh6XLwHrOhJ20dPlEDSXsD1wLm2N5dTiVbSfW2/9HkEjojaRlGlrR9JsykC2DW2bygP/0LS3PL8XGDTJK6dUAJgRNTihhpBylnjLwPW2b6o49QK4Kzy81nAtydx7YQSACOiNrva1sdxwJnA8ZLWlNsS4ELgbZIepVh//EIASQdLWtnn2gnlHWBE1NbESBDbt9N7lUmAE3qk3wAsqXDtuBqpAUq6XNImSQ+Oc16SLpa0XtL9ko5qIt+ImHlF7a5+N5iZ0NQj8BXASROcP5liHeCFwFLgqw3lGxEDoM3LYmL7NuDpCZKcClzlwp3AvmOtOhEx/Bp6Bzjtpusd4Dzg8Y79kfLYxu6EkpZS1BKZtf++01G2iKjBiNEhnRB1ukrdq+7b898D28ttH2376Fl77znFxYqIJrjiNmimqwY4Aizo2J8PbJimvCNiKnl45wOcrhrgCuD9ZWvwMcCzY2P7ImInMKRVwEZqgJKuBRZTDHQeAT4LzAawvQxYSdFfZz3wK+CDTeQbEYNhWGuAjQRA22f0OW/g7CbyiojBYmB0tMUBMCJazECba4AR0W6D2MevigTAiKgvATAi2mkwx/lWkQAYEfWlBhgRrWRwWoEjor0SACOirfIIHBGtlQAYEa2UjtAR0WbpCB0R7ZVW4IhoK6UGGBGtNKBz/VWRABgRNSmNIBHRYqkBRkRrjc50AV6eBMCIqGeI+wE2siiSpMslbZL04DjnF0t6VtKacvtME/lGxGCQq20T3kNaIOkWSeskrZV0Tnl8f0mrJD1a/txvnOtPkvSIpPWSzq9S7qZWhbsCOKlPmh/aPqLcPt9QvhExCJpZFW4bcJ7t1wPHAGdLegNwPnCz7YXAzeX+S0iaBXwFOBl4A3BGee2EGgmAtm8Dnm7iXhHRTrY32l5dfn4OWAfMA04FriyTXQm8u8fli4D1tn9qeyvw9fK6CU3nO8BjJd1HsSD6J22v7ZVI0lJgKcBuB+7DQQt+OY1FnB6bNr1qposwZV4sVkPd6WzfbTjfcfUz2tCvaxIdoedIuqdjf7nt5TvcTzoUOBK4CzhobB1x2xslHdjjvvOAxzv2R4C39CvMdAXA1cBrbG+RtAT4FrCwV8LyD2M5wJ6vmzukjesRLWImMxTuSdtHT5RA0l7A9cC5tjdLle7dK1Hf+NHUO8AJ2d5se0v5eSUwW9Kc6cg7IqZBM+8AkTSbIvhdY/uG8vAvJM0tz88FNvW4dARY0LE/n+Jpc0LTEgAlvVplGJe0qMz3qenIOyKmXkOtwAIuA9bZvqjj1ArgrPLzWcC3e1x+N7BQ0msl7QacXl43oUYegSVdCyymeL4fAT4Lxcsg28uA9wIflbQNeAE43R7WCXQiYgfN/N98HHAm8ICkNeWxC4ALgW9I+jDwGHAagKSDgUttL7G9TdLHgJuAWcDl47UzdGokANo+o8/5S4BLmsgrIgZQAwHQ9u2Mv7jICT3SbwCWdOyvBFZOJs+MBImIWqo83g6qBMCIqC8TokZEW6UGGBHtlQAYEa2Ud4AR0WoJgBHRVhrSCVGnZSRIRMQgSg0wIurLI3BEtFIaQSKi1RIAI6K1EgAjoo3E8LYCJwBGRD15BxgRrZYAGBGtlQAYEW2VR+CIaK8hDYC1h8JJWiDpFknrJK2VdE6PNJJ0saT1ku6XdFTdfCNiQLhoBa6yDZomaoDbgPNsr5a0N3CvpFW2H+pIczLFOsALKRYr/ioVFi2OiCHR1hqg7Y22V5efnwPWUazS3ulU4CoX7gT2HVvnMyKGXxPLYs6ERmeDkXQocCRwV9epecDjHfsj7Bgkx+6xVNI9ku7Z9uyvmixeREyVhhZGn26NBUBJe1Gs6H6u7c3dp3tc0vOPw/Zy20fbPnrXV+3RVPEiYqpUDX4DGACbWhh9NkXwu8b2DT2SjAALOvbnAxuayDsiZpYYzMfbKppoBRZwGbDO9kXjJFsBvL9sDT4GeNb2xrp5R8RgGNZ3gE3UAI8DzgQekLSmPHYBcAiA7WUUq7UvAdYDvwI+2EC+ETEoBjC4VVE7ANq+nd7v+DrTGDi7bl4RMaCGNABmTZCIqKfi42+VR2BJl0vaJOnBjmOHS7pD0gOSviNpn3Gu/UQ5GONBSddK2r1ffgmAEVFfc63AVwAndR27FDjf9puBbwKf6r5I0jzgT4Cjbb8JmAWc3i+zBMCIqK2poXC2bwOe7jp8GHBb+XkV8J5xLt8VeKWkXYE9qNDTJAEwImqbxCPwnLGBDuW2tMLtHwTeVX4+jZd2qQPA9r8CXwQeAzZS9DT5p343TgCMiHom1xH6ybGBDuW2vEIOHwLOlnQvsDewtTuBpP0ohty+FjgY2FPS+/rdOAEwIuqbwpEgth+2/XbbvwdcC/zfHslOBP7Z9hO2fw3cAPx+v3snAEZELWMjQaaqI7SkA8ufuwCfBpb1SPYYcIykPcrBGSdQTMwyoQTAiKhNo6609b2PdC1wB3CYpBFJHwbOkPQT4GGKho2vlWkPlrQSwPZdwHXAauABitjW9/E6M0JHRD0NTnRg+4xxTn25R9oNFCPMxvY/C3x2MvklAEZEbYM4zreKBMCIqC8BMCLaKjXAiGivBMCIaCUP5opvVSQARkQtwzwjdAJgRNTn4YyACYARUVtqgBHRTgO64lsVTSyKtEDSLZLWlbOxntMjzWJJz0paU26fqZtvRAyOpuYDnG5N1AC3AefZXi1pb+BeSatsP9SV7oe2T2kgv4gYMIMY3KpoYlGkjRQTEGL7OUnrgHlAdwCMiJ2RSSMIgKRDgSOBu3qcPlbSfRSzOXzS9tpx7rEUWApwyLxdufOI65os4kBY9OPTZroIU+aJrfvOdBGmxDZmzXQRpoQb+lrD2gjS2HRYkvYCrgfOtb256/Rq4DW2Dwf+HvjWePexvXxsttgDfmvn/EsXsdOZwglRp1IjAVDSbIrgd43tG7rP295se0v5eSUwW9KcJvKOiJk11ROiTqXaj8Dl7KuXAetsXzROmlcDv7BtSYsoAu9TdfOOiAHgapOdDqIm3gEeB5wJPCBpTXnsAuAQANvLgPcCH5W0DXgBON0e0remEbGjIf2/uYlW4NspasETpbkEuKRuXhExmAbx8baKjASJiHoMtPgROCLabjjjXwJgRNSXR+CIaK02twJHRJsNaCfnKhIAI6KWoiP0cEbABMCIqK+ts8FERKQGGBHtNMTvABubDSYi2qoYC1xl60fS5ZI2SXqw49jhku6Q9ICk70jaZ5xr95V0naSHyxnqj+2XXwJgRNRnV9v6uwI4qevYpcD5tt8MfBP41DjXfhn4nu3fBQ4H1vXLLAEwIupxc2uC2L4NeLrr8GHAbeXnVcB7uq8ra4VvpZiZCttbbT/TL78EwIior3oNcI6kezq2pRXu/iDwrvLzacCCHml+G3gC+JqkH0u6VNKe/W6cABgR9VWfEfrJsRnfy215hbt/CDhb0r3A3sDWHml2BY4Cvmr7SOB54Px+N04rcETUptGp6who+2Hg7QCSXge8s0eyEWDE9th6RNdRIQCmBhgR9ZiiI3SV7WWQdGD5cxfg08CyHYpg/xx4XNJh5aETqLAyZQJgRNQijFxt63sv6VrgDuAwSSOSPgycIeknwMMUq0p+rUx7sKSVHZd/HLhG0v3AEcAX+uWXR+CIqK+hkSC2zxjn1Jd7pN0ALOnYXwMcPZn8atcAJe0u6UeS7pO0VtLneqSRpIslrZd0v6Sj6uYbEQOkuX6A06qJGuCLwPG2t5TLY94u6bu27+xIczKwsNzeAny1/BkRw27sHeAQamJRJANbyt3Z5dYd6k8FrirT3lkOWZlre2Pd/CNi5k1lK/BUamph9FnlkpibgFUdTdFj5gGPd+yPlMciYuhVfPwdwEfgRgKg7e22jwDmA4skvakrSa9lM3v+aUhaOtZL/ImntjdRvIiYSqbdAXBMOfbuVnYczDzCS4evzKdozu51j+VjvcQP+K1ZTRYvIqbKFPYDnEpNtAIfIGnf8vMrgRMp+ut0WgG8v2wNPgZ4Nu//InYeTfUDnG5NtALPBa6UNIsioH7D9o2SPgJgexmwkqK/znrgV8AHG8g3IgbFAAa3KppoBb4fOLLH8WUdnw2cXTeviBhANmwfwOfbCjISJCLqa2sNMCIiATAi2slAhfU+BlECYETUZHDeAUZEG5k0gkREi+UdYES0VgJgRLTTYI7zrSIBMCLqMTCk02ElAEZEfakBRkQ7ZShcRLSVwekHGBGtlZEgEdFaeQcYEa1kpxU4IlosNcCIaCfj7cO5gFkCYETUM8TTYTW6KlxEtJRHq219SLpc0iZJD3YcO1zSHZIekPQdSftMcP0sST+WdGOVYjexKtzukn4k6T5JayV9rkeaxZKelbSm3D5TN9+IGAwGPOpKWwVXsOOyupcC59t+M/BN4FMTXH8OsK5q2ZuoAb4IHG/7cOAI4KRy6ctuP7R9RLl9voF8I2IQ2I3VAG3fBjzddfgw4Lby8yrgPb2ulTQfeCdFwKykdgB0YUu5O7vchvOFQES8LN6+vdL2Mj0IvKv8fBqwYJx0XwL+lEkswd5II0i5JvC9wO8AX7F9V49kx0q6D9gAfNL22nHutRRYWu5umTV3/SNNlLGCOcCT05PV30xPNoVp/F7Tamf9XjC93+01dW/wHL+86fu+bk7F5LtLuqdjf7nt5X2u+RBwcfnqbAWwtTuBpFOATbbvlbS4YlmQG+y/I2lfimf0j9vufIm5DzBqe4ukJcCXbS9sLOMGSLrH9tEzXY6m5XsNn535u1Uh6VDgRttv6nHudcDVthd1Hf8b4ExgG7A7sA9wg+33TZRXo63Atp8BbqXrJabtzWOPybZXArMlVf0XIyJaTNKB5c9dgE8Dy7rT2P5z2/NtHwqcDvygX/CDZlqBDyhrfkh6JXAi8HBXmldLUvl5UZnvU3Xzjoidi6RrgTuAwySNSPowcIakn1DElQ3A18q0B0taWSe/Jt4BzgWuLN8D7gJ8w/aNkj4CYHsZ8F7go5K2AS8Ap7vJZ+9m9HsPMazyvYbPzvzdJmT7jHFOfblH2g3Akh7Hb6V4Eu2r0XeAERHDJCNBIqK1EgAjorVaHwAlnSTpEUnrJZ0/0+VpSq8xlTsDSQsk3SJpXTn08pyZLlMTqgwpjea1+h1g2XDzE+BtwAhwN3CG7YdmtGANkPRWYAtwVa/+VMNK0lxgru3Vkvam6ID/7mH/nZW9JPYs+8rOBm4HzrF95wwXbafW9hrgImC97Z/a3gp8HTh1hsvUiHHGVA492xttry4/P0cx8H3ezJaqvgwpnRltD4DzgMc79kfYCf5naotyxMCRQK+hl0OnnMppDbAJWDXOkNJoUNsDoHocy7+6Q0DSXsD1wLm2N890eZpge7vtI4D5wCJJO82ri0HV9gA4wktnlphP0dM8Blj5jux64BrbN8x0eZo23pDSaF7bA+DdwEJJr5W0G8UYwhUzXKaYQNlYcBmwzvZFM12eplQZUhrNa3UAtL0N+BhwE8XL9G+MN03XsBlnTOXO4DiKWT+O75hhfIfhUENoLnCLpPsp/mFeZbvStO7x8rW6G0xEtFura4AR0W4JgBHRWgmAEdFaCYAR0VoJgBHRWgmAEdFaCYAR0Vr/D94WvMv482kDAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Шаг по времени: 2 \n",
      "Распределение давления, МПа\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>20.000</td>\n",
       "      <td>19.906</td>\n",
       "      <td>19.679</td>\n",
       "      <td>19.175</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>20.094</td>\n",
       "      <td>20.000</td>\n",
       "      <td>19.856</td>\n",
       "      <td>19.695</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>20.321</td>\n",
       "      <td>20.144</td>\n",
       "      <td>20.000</td>\n",
       "      <td>19.904</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>20.825</td>\n",
       "      <td>20.305</td>\n",
       "      <td>20.096</td>\n",
       "      <td>20.000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       1      2      3      4\n",
       "1 20.000 19.906 19.679 19.175\n",
       "2 20.094 20.000 19.856 19.695\n",
       "3 20.321 20.144 20.000 19.904\n",
       "4 20.825 20.305 20.096 20.000"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUAAAAEICAYAAAAumy2rAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAAAflklEQVR4nO3df7AdZZ3n8feHEED5IWgI5JcGy5gRtRDNYhh2raygQmSEHcWCHRGR2ZQWKrjoiAzKaM1YzJTFCINjNhuQsFA4qCiRiWJEWGQXEMgEEgxqBhGuiYSAEEAEkvvZP7qve3Jy7j3n5vS995z052V13f7x9PM8N5hvnu6nu7+yTUREHe020R2IiJgoCYARUVsJgBFRWwmAEVFbCYARUVsJgBFRWwmAEVFbCYDjQNJDko5p2F4gaaCpzBWS/nb8excjkfRaSddLekzSE5JulDR3ovsV1UgAjBjZ/sByYC5wEPBT4PqJ7FBUJwGwR0n6pqTfSnpK0q2SXt9w7ApJL0h6pmF5WYs6PiRpW1O5QUkLyuN7SvqKpA3l8hVJew7Tn8a6tkj6saQZ5bEFkgYknSdpczni/YuGc/eU9GVJD0t6VNJiSS9pONeSLm4of2i576qGffMl/V9JT0q6d+h3KI/dIukvG7aPkfRQw/YfR+CS9in7cFsn/x1s/9T2ZbafsP0i8I/AXEmv6OT86G0JgL3r+8AcYCqwCri66fg/2N6nYXlqmHpubywHbGg49tfAfOBNwGHAEcD5I/Tp9rKOqcDzwCcbjh0MTAFmAKcBSxouFf8eeG3ZzmvKMp9vOPcx4LiG4PuXwLqhg2Wg/Vfgb4GXA58Cvi3pwBH6OpxPAy827pB0g6RzOzz/bcBvbT++E21Hj0kA7FG2L7f9tO3ngb8BDms1yuvSXwBftL3J9mPAF4BTOzhvt3JpDgKfs/287f9NEbDeL0nAfwM+WY6inga+BJzccN4LwArgzyXtARwHfLfh+AeAFbZX2B60vRK4G1g4ml9W0kHAGcBFjfttH2/7wg7Onwl8Ffjvo2k3etfuE92B2JGkScDfAScBBwKD5aEpwHAjvZ0xHfh1w/avy33DmS/pSeClwEPAhxqO/c72sy3qOrAsf08RCwEQMKmp7qXAV4BtwA8oguKQVwEnSfqzhn2TgZsbti+R9OVyfXdgc4v+/w3wT8ATrX+94ZWjzR8C/2z7mtGeH70pI8De9F+BE4BjgJcBs8v9Gu6EnbSBIrgMeSXbXyI3u8P2/sBewFXAFQ3HDpC0d4u6NgPPAa+3vX+5vKy8lP4j22spAuX5FMGw0SPA/2o4f3/bezeN2j4xdAw4sUXfXwu8C7hkhN+vJUkHUAS/5bb/brTnR+9KABw/kyXtJWkvYA+Aoe1yX+OIaF+Ke2yPUwSFL41Rn64Bzpd0oKQpFPflrmpzDoApRmrN9+C+IGkPSf8JOB74pu1B4H8C/yhpKhT39CS9q0W9XwJ+ZPv+pv1XAX8m6V2SJpV/ZgvKS9JOnU9xuf/cKM5B0n7AjcD/sd3pfcLoEwmA42cFxUjoOYq/UDMatp+juM815EqKS8jfAD8D7hijPv0txb20+4A1FJMtIz2LeKSkZyguw/8c+FjDsd8Cv6MY9V0NfMT2A+WxzwDrgTskbQF+RPFYyXZs32B7h/trth+hGBGfRzFh8gjFZMZo/v/7OMWf6w4kfV/SecOc91+A/wCc3jSb/spRtB09SvkganSrfCTlKtujGZFFTLiMACOitrqaBZb0cuBfKG7SPwS83/bvWpR7CHia4r7RVtvzumk3IqIKXV0CS/oH4AnbF5YPkh5g+zMtyj0EzLPd6tGEiIgJ0e0l8AnAsnJ9Ga0fP4iI6EndjgCfLJ+7Gtr+ne0DWpT7FcUMoYH/YXvJCHUuAhYBaI893jJ56tSd7l+v0mD7Mv1q0vMT3YOxMenxZ9sX6kN/4Fle8PNdPV/6rv+8tx9/YltHZe+57/kbbR/bTXtVansPUNKPKN7zbPbXo2jnKNsbyufAVkp6wPatrQqWwXEJwJ6zZnn6OWePopn+MPnpXXfuab8Hd82nCg5YdvtEd2FM3Ombuq7j8Se28dMbO3sqaNK0X04Z7pikWRSPKh1M8fbTEtsXj2Ku4ZMU75Gb4rGu023/YaT+tP2baPsY229osVwPPCppWtn4NGDTMHVsKH9uAr5D8dJ9ROwCDAx2+L82tgLn2H4dxUc6zpR0KHAucJPtOcBN5fZ2yg9mfIJiruENFC8WnNxcrlm3Q5HlFF/+oPy5w3fSJO0tad+hdeCdwNou242IHmHMi97W0TJiPfZG26vK9acpvgg0g87nGnYHXiJpd4o3qEZ6rRPoPgBeCLxD0i+Bd5TbSJouaUVZ5iDgNkn3UnxM8l9t/6DLdiOih4xiBDhF0t0Ny6JW9UmaDRwO3AkcZHsjFEGS4nNs27H9G+DLwMPARuAp2z9s1++ungMsv4l2dIv9Gyg/VWT7QYpvzUXELsiYbZ1Ppm5u9xywpH2AbwNn297S8BWhkc45gGKkeAjwJPBNSR+wPeK77bvu3fiIGDeDuKOlHUmTKYLf1bavK3d3MtdwDPAr24+VX+6+DvjTdu0lAEZEV4pPA7mjZSTlx3MvA9bZbvxobdu5BopL3/mSXlrWczQNXxUfTgJgRHStohHgURRfJH+7pNXlspAO5hps3wl8i+KLRmsoYtuwzxsPyRehI6IrBl6s4KtStm9j+I/+jjjXUG5fAFwwmjYTACOiK+7g8rZXJQBGRHcM2/oz/iUARkR3ijdB+lMCYER0SWyrPF/X+EgAjIiuFJMgCYARUUPFc4AJgBFRU4MZAUZEHWUEGBG1ZcS2Pn2pLAEwIrqWS+CIqCUjXvCkie7GTkkAjIiuFA9C5xI4ImoqkyARUUu22Ob+HAFW0mtJx0r6uaT1klplbJKkS8rj90l6cxXtRkRvGEQdLb2m6xGgpEnAVyk+VDgA3CVpue2fNRQ7DphTLm8Fvlb+jIg+V0yC9OfFZBUjwCOA9bYftP0C8A2K5CSNTgCudOEOYP+hb/xHRH8bmgTpZOk1VfRoBvBIw/ZAuW+0ZSKiT22zOlpGImmWpJslrZN0v6Szyv0vl7RS0i/LnwcMc/7+kr4l6YGyjiPb9buKANjqt2r+PGInZYqC0qKhnKHbnn22685FxNgaehOkk6WNrcA5tl8HzAfOlHQocC5wk+05wE3ldisXAz+w/ScUqXjHJSnSADCrYXsmO2Zk76QMALaX2J5ne96kvfeuoHsRMdYGvVtHy0hsb7S9qlx/miKAzaC4hbasLLYMOLH5XEn7AW+jyCqH7RdsP9mu31UEwLuAOZIOkbQHcDJFGrtGy4EPlrPB8ymytm+soO2ImGDFxxA6HgFOGbrCK5dFreqUNBs4HLgTOGgoXpQ/p7Y45dXAY8DXJf2bpKWS2o6gup66sb1V0seAG4FJwOW275f0kfL4YmAFRfam9cDvgdO7bTcieoMRL3b+Ktxm2/NGKiBpH4rk6Gfb3lKk+W1rd+DNwMdt3ynpYopL5c+1O6lrtldQBLnGfYsb1g2cWUVbEdFbbCp7EFrSZIrgd7Xt68rdj0qaZntj+fTIphanDgADZX5gKHIED3ev8I96b146IvpMZw9Bt3sQWsVQ7zJgne2LGg4tB04r108Drm8+1/ZvgUckzS13HQ38rLlcs/58ejEieoapbAR4FHAqsEbS6nLfecCFwLWSzgAeBk4CkDQdWGp7KDn6x4Gry7mIB+ngVlsCYER0rYoPotq+jdaPzEExomsuv4FibmFoezUw4v3FZgmAEdEVo3wQNSLqqUiL2Z+hpD97HRE9JInRI6KmDG3f8uhVCYAR0bWMACOilmxlBBgR9VRMgiQrXETUUv/mBEkAjIiuFJMguQcYETVVxZsgEyEBMCK6kjdBIqLWejHhUScSACOiKza8OJgAGBE1VFwCJwBGRE3lTZCIqKV+fgymknGrpGMl/VzSekk7fIdf0gJJT0laXS6fr6LdiOgFqiQt5kToegQoaRLwVeAdFIlJ7pK03Hbz9/h/Yvv4btuLiN7TLt9HJyTNAq4EDgYGgSW2L5b0cuBfgNnAQ8D7bf9umDomAXcDv+kk3lRxCXwEsN72g2UHvkGRyLhtQpK2JpnB/bZ2XU2veVG77p2HZ6f33r/yVfDpR050F8bEtuvv6LqOYha4kneBtwLn2F4laV/gHkkrgQ8BN9m+sLzCPBf4zDB1nEWRUH2/Thqs4v+tM4BHGrYHyn3NjpR0r6TvS3r9cJVJWjSUNHnbM89W0L2IGEtDD0J3soxYj73R9qpy/WmKQDaDYkC1rCy2DDix1fmSZgLvBpZ22vcqhiKtfis3ba8CXmX7GUkLge8Cc1pVZnsJsARgz9kzm+uJiB40ikvgKZLubtheUv6d346k2cDhwJ3AQbY3QhEkJU0dpu6vAH8F7NtpZ6oIgAPArIbtmcCGxgK2tzSsr5D0z5Km2N5cQfsRMYFGOQu82faImdsk7UORHP1s21uKdMEjk3Q8sMn2PZIWdNqZKi6B7wLmSDqkzMd5MkUi48bOHVwmPUbSEWW7j1fQdkT0gKpmgSVNpgh+V9u+rtz9qKRp5fFpwKYWpx4FvEfSQ8A3gLdLuqpde10HQNtbgY8BN1Jcs19r+35JH5H0kbLY+4C1ku4FLgFOtp3L24hdgC22ereOlpGUg6TLgHW2L2o4tBw4rVw/Dbh+xz74s7Zn2p5NMQj7se0PtOt7JdORtlcAK5r2LW5YvxS4tIq2IqL3VPQg9FHAqcAaSavLfecBFwLXSjoDeBg4CUDSdGCp7YUt6urIrvs8RkSMi6reBLF9G60nVQGOblF+A7BD8LN9C3BLJ20mAEZE1/r1VbgEwIjoSj6IGhG1VsWrcBMhATAiumLD1nwQNSLqKpfAEVFLuQcYEbXmBMCIqKtMgkRELdm5BxgRtSW2ZRY4Iuoq9wAjopb6OStcAmBEdMfFfcB+lAAYEV3LLHBE1JIzCRIRdZZL4IiorX6dBa5k3CrpckmbJK0d5rgkXSJpvaT7JL25inYjYuLZRQDsZOk1VV24XwEcO8Lx4yjyAM8BFgFfq6jdiOgBVSRGlzRL0s2S1km6X9JZ5f6XS1op6ZflzwM6PbedSgKg7VuBJ0YocgJwpQt3APsPpbmLiP5nd7a0sRU4x/brgPnAmZIOBc4FbrI9B7ip3O703BGN19TNDOCRhu2Bct8OJC2SdLeku7c98+y4dC4idp4Rg4O7dbSMWI+90faqcv1pijS7MygGUMvKYsuAE0dx7ojGKwC2Gvu2/PfA9hLb82zPm7TP3mPcrYiogjtcgClDA5xyWdSqPkmzgcOBO4GDbG+EItABU0fqS9O5IxqvWeABYFbD9kxgwzi1HRFjyaOaBd5se95IBSTtA3wbONv2liJfemeaz21XfrxGgMuBD5azwfOBp4YiekTsAkYxBByJpMkUAexq29eVux8dmjMof24axbkjqmQEKOkaYAHF8HYAuACYDGB7MbCCIoHxeuD3wOlVtBsRvaGKR1xUDPUuA9bZvqjh0HLgNODC8uf1ozh3RJUEQNuntDlu4Mwq2oqI3mJgcLCSZ/yOAk4F1khaXe47jyLwXSvpDOBh4CQASdOBpbYXDneu7RUjNZg3QSKiOwYqGAHavo3WE6YAR7cov4HiyrLducNKAIyIruVd4IiorwTAiKin3nzPtxMJgBHRvYwAI6KWDK5mFnjcJQBGRAUSACOirnIJHBG1lQAYEbVU0YPQEyEBMCK6lgehI6K+MgscEXWljAAjopY6/NZfL0oAjIguKZMgEVFjGQFGRG0NTnQHdk4CYER0p4+fA6wkKZKkyyVtkrR2mOMLJD0laXW5fL6KdiOiN8idLSPWIc2SdLOkdZLul3RWuf/lklZK+mX584Bhzj9W0s8lrZfUKnn6DqrKCncFcGybMj+x/aZy+WJF7UZEL6gmK9xW4BzbrwPmA2dKOhQ4F7jJ9hzgpnJ7O5ImAV8FjgMOBU4pzx1RJQHQ9q3AE1XUFRH1ZHuj7VXl+tPAOmAGcAKwrCy2DDixxelHAOttP2j7BeAb5XkjGs97gEdKupciIfqnbN/fqlCZKX4RwJ5T9+XVs1umAO1rD/32FRPdhTHzB/ac6C6MiW17jVcK7fG1raL/XKN4EHqKpLsbtpfYXrJDfdJs4HDgTuCgoTzitjdKmtqi3hnAIw3bA8Bb23VmvALgKuBVtp+RtBD4LjCnVcHyD2MJwL5zD+7TyfWIGjGjeRVus+15IxWQtA9FgvOzbW8pUv621apQ2/gxLv+s2d5i+5lyfQUwWdKU8Wg7IsZBNfcAkTSZIvhdbfu6cvejkqaVx6cBrS4LB4BZDdszKa42RzQuAVDSwWXmdiQdUbb7+Hi0HRFjr6JZYAGXAetsX9RwaDlwWrl+GnB9i9PvAuZIOkTSHsDJ5XkjquQSWNI1wAKK6/sB4AJgMoDtxcD7gI9K2go8B5xs9+sHdCJiB9X8bT4KOBVYI2l1ue884ELgWklnAA8DJwFImg4stb3Q9lZJHwNuBCYBlw83z9CokgBo+5Q2xy8FLq2irYjoQRUEQNu3MXxykaNblN8ALGzYXgGsGE2beRMkIrrSyeVtr0oAjIju5YOoEVFXGQFGRH0lAEZELeUeYETUWgJgRNSV+vSDqLvmG94RER3ICDAiupdL4IiopUyCREStJQBGRG0lAEZEHYn+nQVOAIyI7uQeYETUWgJgRNRWAmBE1FW/XgLnTZCI6F51SZEul7RJ0tqGfYdJul3SGknfk7TfMOd+UtL9ktZKukbSXu3a6zoASpol6WZJ68rGz2pRRpIukbRe0n2S3txtuxHRI1zMAneydOAK4NimfUuBc22/EfgO8OnmkyTNAD4BzLP9Boq8ICe3a6yKEeBW4BzbrwPmA2dKOrSpzHEUeYDnUCQ9/1oF7UZEr6hoBGj7VuCJpt1zgVvL9ZXAe4c5fXfgJZJ2B17KeKTFtL3R9qpy/WlgHUWW9kYnAFe6cAew/1Cez4jof6NIizlF0t0Ny6IOql8LvKdcP4nt8/8CYPs3wJcpssZtBJ6y/cN2FVd6D1DSbOBw4M6mQzOARxq2B9gxSA7VsWjoD+fFJ39fZfciYqx0PgLcbHtew7Kkg9o/THFleQ+wL/BCcwFJB1AMtA4BpgN7S/pAu4orC4CS9qHI6H627S3Nh1uc0nJAbHvJ0B/O5P1fWlX3ImKsdBr8dnKm2PYDtt9p+y3ANcC/tyh2DPAr24/ZfhG4DvjTdnVXEgAlTaYIflfbvq5FkQG2H7bOpIPr84jofWJUl8Cjr1+aWv7cDTgfWNyi2MPAfEkvlSSKPMLr2tVdxSywgMuAdbYvGqbYcuCD5WzwfIrr843dth0RvaGqACjpGuB2YK6kAUlnAKdI+gXwAMXA6etl2emSVgDYvhP4FrAKWEMR29peXlfxIPRRwKnAGkmry33nAa8sO7aYIlv7QmA98Hvg9ArajYheUdGD0LZPGebQxS3KbqCIK0PbFwAXjKa9rgOg7dtofY+vsYyBM7ttKyJ6VJ++CZJX4SKiO/kaTETUWgJgRNRVPogaEbWVS+CIqKcuHnKeaAmAEdG9BMCIqKOhN0H6UQJgRHRNg/0ZARMAI6I7uQcYEXWWS+CIqK8EwIioq4wAI6K+EgAjopacV+EioqbyHGBE1Jv7MwImAEZE1/p1BFhpWsyIqKEKs8JJulzSJklrG/YdJul2SWskfU/SfsOcu7+kb0l6QNI6SUe2a6+KpEizJN1cNni/pLNalFkg6SlJq8vl8922GxG9Q4OdLR24Aji2ad9S4FzbbwS+A3x6mHMvBn5g+0+Aw+ggK1wVl8BbgXNsr5K0L3CPpJW2f9ZU7ie2j6+gvYjoMVXNAtu+VdLspt1zgVvL9ZXAjcDntmu/GBW+DfhQWc8LtEig3qzrEaDtjbZXletPU0TdGd3WGxF9whSTIJ0sMEXS3Q3Log5aWAu8p1w/ie1zjA95NfAY8HVJ/yZpqaS921Vc6SRIGbkPB+5scfhISfdS5PX8lO37h6ljEbAI4JUzduemQ5dX2cWesGDwxInuwpj59R+mTHQXxsSLu+h8oSdVU88oJkE22543yuo/DFxS3jpbTuuR3e7Am4GP275T0sXAuTSNFJtVNgkiaR/g28DZtrc0HV4FvMr2YcA/Ad8drh7bS2zPsz3vwFdU9F8nIsZWRZMgLau2H7D9TttvAa4B/r1FsQFgoEyQDkWS9De3q7uSAChpMkXwu9r2dc3HbW+x/Uy5vgKYLGnXHCpE1MzQg9CdLDtVvzS1/LkbcD6wuLmM7d8Cj0iaW+46Gmieh9hBFbPAAi4D1tm+aJgyB5flkHRE2e7j3bYdET3ARoOdLe1Iuga4HZgraUDSGcApkn4BPEBxC+3rZdnpklY0nP5x4GpJ9wFvAr7Urr0qbmwcBZwKrJG0utx3HvBKANuLgfcBH5W0FXgOONnu00fHI2JHFf1ttn3KMIcublF2A7CwYXs1MKr7i10HQNu3UYyCRypzKXBpt21FRG/q1zdBds2prYgYPwaSEyQiaqs/418CYER0L5fAEVFbSYsZEfWUtJgRUVfFg9D9GQETACOie8kJEhF1lRFgRNRT7gFGRH119p5vL0oAjIju5RI4ImopidEjotYyAoyI2urP+JcAGBHd02B/XgMnAEZEd0zfPghdWVKkiKgnYeTOlrZ1SZdL2iRpbcO+wyTdLmmNpO+VOYCHO39SmRbzhk76ngAYEd3rPC9wO1cAxzbtWwqca/uNwHeAT49w/lkUuck7UkVSpL0k/VTSvZLul/SFFmUk6RJJ6yXdJ6lturqI6CMVBUDbtwJPNO2eC9xarq8E3tvqXEkzgXdTBMyOVDECfB54e5nz903AsZLmN5U5DphTLouAr1XQbkT0gqF7gJ0sMEXS3Q3Log5aWAu8p1w/CZg1TLmvAH/FKO5IVpEUycAz5ebkcmkO9ScAV5Zl75C0v6Rptjd2235ETLxRzAJvtj2qzG3Ah4FLJH0eWA68sEP70vHAJtv3SFrQacWVzAJLmgTcA7wG+GpDdvYhM4BHGrYHyn0JgBF9r+P7eztXu/0A8E4ASa+luMxtdhTwHkkLgb2A/SRdZfsDI9VdySSI7W223wTMBI6Q9IamIq3SZrb8E5O0aGh4/Njj26roXkSMJVPlJMgOJE0tf+4GnA8s3qEL9mdtz7Q9GzgZ+HG74AcVzwLbfhK4hR1ncQbY/rp9JkWG91Z1LLE9z/a8A18xqcruRcRY6fwe4IgkXQPcDsyVNCDpDOAUSb8AHqCIG18vy06XtKKbbnd9CSzpQOBF209KeglwDPD3TcWWAx+T9A3grcBTuf8Xseuo6oOotk8Z5tDFLcpuABa22H8LxUCsrSruAU4DlpX3AXcDrrV9g6SPlJ1ZDKwoO7oe+D1wegXtRkSvqOvHEGzfBxzeYv/ihnUDZ3bbVkT0IBu29ee7cHkXOCK6V9cRYEREAmBE1JOB5ASJiHoyOPcAI6KOTCZBIqLGcg8wImorATAi6mlsP4YwlhIAI6I7BpIUKSJqKyPAiKinvAoXEXVlcJ4DjIjaypsgEVFbuQcYEbVk9+0scBKjR0T3KsoJIulySZskrW3Yd5ik2yWtkfQ9Sfu1OG+WpJslrSvzk5/VSbcTACOiS8bbtnW0dOAKdswptBQ41/Ybge8An25x3lbgHNuvA+YDZ0o6tF1jCYAR0Z2hz2F1srSryr4VeKJp91zg1nJ9JfDeFudttL2qXH8aWEeRendECYAR0T0PdrbAlKG0t+WyqIPa1wLvKddPYvsMkzuQNJsiTUdzfvIdVJEVbi+K6LxnWd+3bF/QVGYBcD3wq3LXdba/2G3bETHxDLjzx2A22543yiY+DFwi6fMUGSZfGK6gpH2AbwNn297SruIqZoGfB95u+xlJk4HbJH3f9h1N5X5i+/gK2ouIXuKx/SCq7QeAdwJIei3w7lblyvjzbeBq29d1UncVWeEMPFNuTi6X/nwoKCJ2SocTHDtF0lTbmyTtBpwPLG5RRsBlwDrbF3Vctyt4gLHMCXwP8Brgq7Y/03R8AUVkHqDI7P4p2/cPU9ciYOi+wFzg5113sDNTgM3j1NZ4yu/Vf8bzd3uV7QO7qUDSDyj63InNtptneRvrugZYUNb3KHABsA//P63udcBnbVvSdGCp7YWS/iPwE2ANMDQcPc/2ihH7XkUAbOj8/hTT1B+33fgcz37AYHmZvBC42PacyhqugKS7d+LeRM/L79V/duXfrddUOgts+0ngFpqe47G9xfYz5foKYLKkTv/FiIgYE10HQEkHliM/JL0EOAZ4oKnMweU1OpKOKNt9vNu2IyK6UcUs8DRgWXkfcDfgWts3SPoIgO3FwPuAj0raCjwHnOwqr72rsWSiOzBG8nv1n135d+spld4DjIjoJ3kTJCJqKwEwImqr9gFQ0rGSfi5pvaRzJ7o/VWn1WaFdwc5+9qjXSdpL0k8l3Vv+Xl+Y6D7VQa3vAZYTN78A3kHxkPZdwCm2fzahHauApLdRvKFzpe03THR/qiJpGjDN9ipJ+1I8gH9iv/83K5+S2LvxlVLgrBavlEaF6j4CPAJYb/tB2y8A3wBOmOA+VWKYzwr1vZ397FGvcyGvlI6zugfAGcAjDdsD7AJ/mepiNJ896geSJklaDWwCVtreJX6vXlb3AKgW+/Kvbh8Y7WeP+oHtbbbfBMwEjpC0y9y66FV1D4ADbP9xxZkUH2uIHrYznz3qJ8O9UhrVq3sAvAuYI+kQSXsAJ1N8cDF61M5+9qjXdfJKaVSv1gHQ9lbgY8CNFDfTrx3uM139pvys0O3AXEkDks6Y6D5V5CjgVODtklaXy8KJ7lQFpgE3S7qP4h/mlbZvmOA+7fJq/RhMRNRbrUeAEVFvCYARUVsJgBFRWwmAEVFbCYARUVsJgBFRWwmAEVFb/w/FVEGhmf+yeQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Шаг по времени: 3 \n",
      "Распределение давления, МПа\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>20.000</td>\n",
       "      <td>19.883</td>\n",
       "      <td>19.619</td>\n",
       "      <td>19.081</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>20.117</td>\n",
       "      <td>20.000</td>\n",
       "      <td>19.826</td>\n",
       "      <td>19.640</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>20.381</td>\n",
       "      <td>20.174</td>\n",
       "      <td>20.000</td>\n",
       "      <td>19.883</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>20.919</td>\n",
       "      <td>20.360</td>\n",
       "      <td>20.117</td>\n",
       "      <td>20.000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       1      2      3      4\n",
       "1 20.000 19.883 19.619 19.081\n",
       "2 20.117 20.000 19.826 19.640\n",
       "3 20.381 20.174 20.000 19.883\n",
       "4 20.919 20.360 20.117 20.000"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUYAAAEICAYAAAAjhV3sAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAAAfGklEQVR4nO3df7BfdX3n8eeLmIAF5IeBGEIU18aMtl2CzUQcdnf4IQqpFW2lQ9pGqrQZOrALLjoy1LXqdDuso1ipDmmEFFhZXB1AKY3SlKWL7AIKGQjEgKQulpiUSFBCFInJfe0f53Ph3G++997v5Zzce7/fvB7Omfs953w+5/P5Xsz7fj7n8znnI9tERMRLDpjqCkRETDcJjBERHRIYIyI6JDBGRHRIYIyI6JDAGBHRIYExIqJDAuMkkPSEpLfX9k+WtLkjzbWS/mLyaxdjkTRb0v+RtF3STyXdI+mkqa5X7FuvmOoKRExzO4EPAo8DBs4C/k7S0bZ3T2nNYp9Ji3GakvQ1Sf8q6VlJd0n6tdq5ayXtkrSzth3W5Rp/JGlPR7ohSSeX8wdK+itJW8r2V5IOHKU+9WvtkPS/JM0r506WtFnSZZKeLi3kP6jlPVDSZyT9i6SnJK2U9MpaXkv6fC39m8uxL9eOnSjp/5ZW20PD36Gc+ydJf1zbf7ukJ2r7L7bYJR1S6nB3L/8dbP/C9mO2hwABe4AjgCN7yR/9KYFx+vomsAA4GlgH3NBx/tO2D6ltz45ynXvq6YAttXN/BpwILAKOB5YAHxujTveUaxwNvAB8qHbuNcBsYB5wLrBK0sJy7r8Bbyzl/GpJ8/Fa3h8DZ9aC8h8DG4dPlgD898BfUAWkDwM3STpqjLqO5iPAL+sHJN0m6dKxMklaD/wCuBW42va2l1F29IkExmnK9mrbz9l+AfgEcHy3VmFDfwB8yvY22z8GPgks7yHfAWXb3nH8v9h+wfb/pgpkvydJwJ8AH7L9jO3ngL8Ezqnl2wWsAX5H0izgTODrtfN/CKyxvcb2kO21wP3A0ol8WUlzgPOAK+rHbb/L9uVj5bX9b4FXAb8P9NTajP6Ve4zTkKQZwH8FzgaOAobKqdnAaC3Dl+MY4Ie1/R+WY6M5UdJPgV8BngD+qHbuJ7Z/1uVaR5X0D1QxEqi6pDM6rn018FdUXdVvUQXLYa8Dzpb027VjM4E7a/tXSvpM+fwK4Oku9f8E8NfAM92/3ths/wK4UdJGSQ/afujlXCemv7QYp6ffp7rJ/3bgMOC4clyjZXiZtlAFnWGvZWRXu9O9tg8HDgK+DFxbO3eEpIO7XOtp4Hng12wfXrbDSpf8RbYfoQqgH6MKknVPAv+9lv9w2wd3tPL+0/A54D1d6v5G4J3AlWN8v17NBP5NC9eJaSqBcfLMlHSQpIOAWQDD++VYvQV1KNU9vO1UweIv91GdbgQ+JukoSbOp7vt9eZw8UI3O7qFqDdZ9UtIsSf8eeBfwtTJo8SXgc5KOhuqeoaR3drnuXwL/aHtDx/EvA78t6Z2SZpTf2cmSju35m1YB91O2n59AnuFBn39XvtcrJX0UmAPcN5HrRH9JYJw8a6haTs8Dt1MNQDxf2/6wlvZ6qq7oj4DvAffuozr9BdW9uvXAw1SDPGPNpXybpJ1U3fnfAS6snftX4CdUrcQbgPNtP1rOfRTYBNwraQfwj8BCOti+zfZ/7nL8SaoW9GVUAzVPUg2iTOT/v9upfq97kfRNSZeNku9A4Isl/4+o7mv+lu2xWtbR55QX1UZTZerMl21PpAUXMW2lxRgR0aHRqLSkI4H/STU48ATwe7Z/0iXdE8BzVPeldtte3KTciIh9qVFXWtKngWdsX14myB5h+6Nd0j0BLLbdbQpFRMS00rQrfRZwXfl8Hd2nSURE9JWmLcaflnljw/s/sX1El3T/j2rE0sDf2F41xjVXACsANGvWb86cc/TLrt90pQF+9cCMF6a6BvvGjGd+Nn6iPvQLfsYuv9Bofuw7TznY25/Z01PaB9a/cLvtM5qUNxnGvcco6R+pnoPt9GcTKOck21vKPLa1kh61fVe3hCVorgI48LXzfcwlF0+gmP4w66eDO+Z12Kah8RP1ocNu2FczpqbWfb6j8TW2P7OH79z+2p7Szpj7+OzGBU6CcQOj7bePdq68pWSu7a2S5gJdH6wfnvNle5ukW6heVtA1MEZEfzEwxGD9QWzadLmV6k0qlJ/f6Ewg6WBJhw5/Bt4BPNKw3IiYJoz5pff0tPWLpoHxcuB0SY8Dp5d9JB0jaU1JMwe4W9JDwHeAv7f9rYblRsQ0MtTj//pFo3mMtrcDp3U5voXySijbP6B6119EDCBj9gzYE3R57VhENDZEAmNExIuqVy0lMEZEjJAWY0REjYFf5h5jRMRLjNOVjogYwbBnsOJiAmNENFM9+TJYEhgjoiGxp/V12qZWAmNENFINviQwRkS8qJrHmMAYETHC0IC1GAf3xYARMSmGW4y9bGORNF/SnZI2Stog6aJy/EhJayU9Xn52exn2QkkP1rYdki4u5z4h6Ue1c0vH+05pMUZEI0bsaaeNtRu4xPa68qrCByStBf4IuKO2ttSlVGuVv1QH+zFgEYCkGVRrgN9SS/I525/ptSJpMUZEY0NWT9tYbG+1va58fg7YCMxj4mtLnQb8s+0fvtzvk8AYEY0YscszetqA2ZLur20rul1T0nHACcB9wBzbW6EKnsB4C0GdA9zYcexCSeslre7WFe+UwBgRjVQTvA/oaQOetr24tu21MJ6kQ4CbgItt75hIXSTNAt4NfK12+CrgDVRd7a3AZ8e7TgJjRDTWxuALgKSZVEHxBts3l8NPlTWlGGttqeJMYJ3tp4YP2H7K9h7bQ8CXqNacGlMCY0Q0Yos9PqCnbSySBFwDbLR9Re3UuGtL1Syjoxs9HFSL99LDmlOtBEZJZ0h6TNKmMmrUeV6Srizn10t6SxvlRsT0MIR62sZxErAcOLVjak0va0sh6VfK+Zs7rvtpSQ9LWg+cAnxovIo0nq5Thsa/WCq0GfiupFttf6+W7ExgQdneStXnf2vTsiNi6lWDL81n/tm+G0aNnmOuLVX2fw68uku65ROtSxstxiXAJts/sL0L+ArV8HrdWcD1rtwLHN7RvI2IPjXBwZe+0EZN5wFP1vY3l2MTTRMRfWqP1dPWL9p48qXbt+18bWUvaaqE1bymFQAzjhh3ulFETLEWn3yZNtoIjJuB+bX9Y4EtLyMNAGVe0yqAA187f8DeCxwxmIbGGXHuN218m+8CCyS9vkyuPIdqeL3uVuD9ZXT6RODZ4ZnsEdHfqpdIHNDT1i8atxht75Z0IXA7MANYbXuDpPPL+ZXAGqrRo03Az4EPNC03IqYHI35ZPe43MFp5u47tNVTBr35sZe2zgQvaKCsiphebcSdv95u8diwiGupp8nZfSWCMiEZMWowREXvpp4GVXiQwRkQjZvyX0PabBMaIaKRaPnWwQslgfZuImAK9vWuxnyQwRkQjZvCefElgjIjG0mKMiKixlRZjRERdNfiSRwIjImqUCd4REXXV4EvuMUZEjJAnXyIiagbxyZfBCvMRMSXaWAxL0nxJd0raKGmDpIvK8SMlrZX0ePnZdc0TSU+UZVIflHR/7XhP+esSGCOiERt+OXRAT9s4dgOX2H4TcCJwgaQ3A5cCd9heANxR9kdziu1FthfXjk0kP5DAGBENVV3pA3raxryOvdX2uvL5OWAj1WqiZwHXlWTXAe+ZYBUnnD+BMSIa21Oelx5vA2ZLur+2reh2PUnHAScA9wFzhteIKj+PHqUaBv5B0gMd1+01/4sy+BIRjUxwus7THd3cvUg6BLgJuNj2Dqnna59ke4uko4G1kh61fVevmetaaTFKOkPSY5I2Sdqr/y7pZEnPlpuiD0r6eBvlRsR00E5XGkDSTKqgeIPtm8vhpyTNLefnAtu65bW9pfzcBtwCLJlI/rrGgVHSDOCLwJnAm4Fl5YZpp2+Xm6KLbH+qabkRMX0MlXVfxtvGoqppeA2w0fYVtVO3AueWz+cC3+iS92BJhw5/Bt4BPNJr/k5tdKWXAJts/6BU6itUNzu/1/jKM8zQYbsbX2a62TXAdzB2zh/M29Z+/9umugr7xJ7b7m18jWpUupVnpU8ClgMPS3qwHLsMuBz4qqTzgH8BzgaQdAxwte2lwBzgltLtfgXwP2x/q1yja/6xtPEvdB7wZG1/M/DWLuneJukhYAvwYdsbul2s3DRdATDj1Ye3UL2I2JfamuBt+24YtVl5Wpf0W6jWq6c0zI4f5brbu+UfSxuBsdsXccf+OuB1tndKWgp8HVjQ7WK2VwGrAA58/bGd14mIaWjQlk9to9+zGZhf2z+WqlX4Its7bO8sn9cAMyXNbqHsiJhiw6PSvWz9oo3A+F1ggaTXS5oFnEN1s/NFkl5TbqwiaUkpd3sLZUfENNDWqPR00bgrbXu3pAuB24EZwGrbGySdX86vBN4H/Kmk3cDzwDm2002OGAC22N1HQa8XrQyPlu7xmo5jK2ufvwB8oY2yImL66aduci8Gd95IREyKvKg2IqKLBMaIiJpBfFFtAmNENDZo8xgTGCOiERt2j/8S2r6SwBgRjaUrHRFRk3uMERFdOIExImKkDL5ERNTYuccYEdFB7MmodETESLnHGBFRk2elIyI6ubrPOEgSGCOisUEblR6sO6YRMelcBl962cYiab6kOyVtlLRB0kXl+JGS1kp6vPw8ote85dwnJP2otq790vG+UwJjRDRm97aNYzdwie03AScCF5Q16i8F7rC9ALij7Pead9jnauvar+mSf4QExohozFZP29jX8Fbb68rn54CNVMsznwVcV5JdB7xnAnlfllYCo6TVkrZJemSU85J0paRNktZLeksb5UbE1Ktag80DY52k44ATgPuAOba3VmV5K3D0BPIOu7DEntXduuKd2moxXgucMcb5M6nWkV4ArACuaqnciJgGJrB86mxJ99e2FZ3XknQIcBNwse0dE6nHKHmvAt4ALAK2Ap8d7zptLYZ1V4nSozkLuL6sDHivpMMlzR3+KxAR/W0C03Wetr14tJOSZlIFthts31wOPzUcLyTNBbZNIC+2n6ql+RJw23iVnKx7jPOAJ2v7mxml/y9pxfBfkz3P/WxSKhcRL58RQ0MH9LSNpaw9fw2w0fYVtVO3AueWz+cC35hAXkowHfZeoOstv7rJCozdbi50/Rtje5XtxbYXzzj04H1crYhog3vcxnESsBw4tWNqzeXA6ZIeB04v+0g6RtKacfICfFrSw5LWA6cAHxqvIpM1wXszML+2fyywZZLKjoh9ye08K237bro3ogBO65J+C7B0vLy2l0+0LpPVYrwVeH8ZnT4ReDb3FyMGSEtNxumilRajpBuBk6lGnDYDfw7MBLC9ElhDFdk3AT8HPtBGuRExPeTtOl3YXjbOeQMXtFFWREwvBoaGEhgjIl5iIC3GiIiR8tqxiIhOCYwREXUTew66HyQwRkRzaTFGRNQYnFHpiIhOCYwRESOlKx0R0SGBMSKiJhO8IyL2lgneERGdMiodETGS0mKMiKjps3ct9iKBMSIaUgZfIiL2khZjRESHoamuQLsSGCOimQGcx9jKYliSVkvaJqnreq2STpb0bG1Zw4+3UW5ETA9yb9uY15DmS7pT0kZJGyRdVI4fKWmtpMfLzyNGyX+GpMckbZJ0ae14T/nr2lol8FrgjHHSfNv2orJ9qqVyI2I6aGeVwN3AJbbfBJwIXCDpzcClwB22FwB3lP0RJM0AvgicCbwZWFby0kv+Tq0ERtt3Ac+0ca2I2D/Z3mp7Xfn8HLARmAecBVxXkl0HvKdL9iXAJts/sL0L+ErJR4/5R5jMe4xvk/QQsAX4sO0N3RJJWgGsADhozqG86fVbJrGKk+P7W4+e6irsM89z0FRXYZ/Yc+BkLcE+ufbMauc6E5jgPVvS/bX9VbZX7XU96TjgBOA+YM7wOvS2t0rq9g9oHvBkbX8z8NbyuZf8I0xWYFwHvM72TklLga8DC7olLL+kVQCHLZwzYJMAIgaQmcgjgU/bXjxWAkmHADcBF9veIfV07W6JXnb8mJQ/g7Z32N5ZPq8BZkqaPRllR8QkaOceI5JmUgXFG2zfXA4/JWluOT8X2NYl62Zgfm3/WKreaa/5R5iUwCjpNSphX9KSUu72ySg7Iva9lkalBVwDbLR9Re3UrcC55fO5wDe6ZP8usEDS6yXNAs4p+XrNP0IrXWlJNwInU90/2Az8OTATwPZK4H3An0raDTwPnGMP2ouKIvZj7fxrPglYDjws6cFy7DLgcuCrks4D/gU4G0DSMcDVtpfa3i3pQuB2YAawujaO0TX/WFoJjLaXjXP+C8AX2igrIqahFgKj7bsZffGY07qk3wIsre2vAdZ0Sbe9W/6x5MmXiGikl25yv0lgjIjm8qLaiIiR0mKMiOiUwBgRUZN7jBERXSQwRkSMpAF7Ue1gPhkfEdFAWowR0Vy60hERNRl8iYjoIoExIqJDAmNExEvE4I1KJzBGRDO5xxgR0UUCY0REhwTGiIiR0pWOiOg0YIGx8SOBkuZLulPSRkkbJF3UJY0kXSlpk6T1kt7StNyImCZcjUr3svWLNlqMu4FLbK+TdCjwgKS1tr9XS3Mm1TrSC6gWwb6KlxbDjoh+lxbjSLa32l5XPj8HbATmdSQ7C7jelXuBw4fXeY2I/tfG8qkAklZL2ibpkdqx4yXdI+lhSX8n6VVd8i2U9GBt2yHp4nLuE5J+VDu3tDN/p1bfriPpOOAE4L6OU/OAJ2v7m9k7eA5fY4Wk+yXdv+vZ59usXkTsK+5xG9+1wBkdx64GLrX9G8AtwEf2Kt5+zPYi24uA3wR+XtIO+9zw+bKa4JhaC4ySDgFuAi62vaPzdJcsXX9NtlfZXmx78azDXtlW9SJiX+k1KPYQGG3fBTzTcXghcFf5vBb43XEucxrwz7Z/2FP9u2glMEqaSRUUb7B9c5ckm4H5tf1jgS1tlB0RU0u015UexSPAu8vnsxkZS7o5B7ix49iFZeB3taQjxiuwjVFpAdcAG21fMUqyW4H3l9HpE4FnbW9tWnZETA8TCIyzh2+VlW1FD5f/IHCBpAeAQ4Fdo9ZDmkUVRL9WO3wV8AZgEbAV+Ox4BbYxKn0SsBx4WNKD5dhlwGsBbK8E1gBLgU1Uff8PtFBuREwXvbcGn7a9eEKXth8F3gEg6Y3Ab42R/Exgne2navlf/CzpS8Bt45XZODDavpvu9xDraQxc0LSsiJim9uF0HUlH294m6QDgY8DKMZIvo6MbLWlurYf6Xqqu+Ziy5ktENNNjN7rH6To3AvcACyVtlnQesEzS94FHqcYm/rakPUbSmlreXwFOBzrHOT5dpvqsB04BPjRePfJIYEQ011KL0fayUU59vkvaLVS36Ib3fw68uku65ROtRwJjRDTWT4/79SKBMSIay9t1IiLqen+qpW8kMEZEcwmMEREvGX7yZZAkMEZEYxoarMiYwBgRzeQeY0TE3tKVjojolMAYETFSWowREZ0SGCMiapxHAiMiRsg8xoiIbjxYkTGBMSIaS4sxIqJuACd4t7EY1nxJd0raKGmDpIu6pDlZ0rO1Ba8/3rTciJg+NNTb1i/aaDHuBi6xvU7SocADktba/l5Hum/bflcL5UXENNNPQa8XbSyGtZVqSUJsPydpIzAP6AyMETGITAZfxiLpOOAE4L4up98m6SGqxWw+bHvDKNdYAawAeO28V7Bm4Zpuyfra6UO/PdVV2Gc2/WLOVFdhn9g1oLfj3dLXGrTBl9ZWCZR0CHATcLHtHR2n1wGvs3088NfA10e7ju1VthfbXnzUq2e0Vb2I2Jfc49YnWgmMkmZSBcUbbHcuXYjtHbZ3ls9rgJmSZrdRdkRMreEJ3i0tn7pa0jZJj9SOHS/pnrIE6t9JetUoeZ8oaR6UdH/t+JGS1kp6vPw8Yrx6tDEqLeAaYKPtK0ZJ85qSDklLSrnbm5YdEdOAjYZ623pwLXBGx7GrgUtt/wZwC/CRMfKfYnuR7cW1Y5cCd9heANxR9sfURovxJGA5cGptOs5SSedLOr+keR/wSLnHeCVwjj1gd2sj9mctdaVt3wU803F4IXBX+bwW+N0J1u4s4Lry+TrgPeNlaGNU+m6q1vRYab4AfKFpWRExPU1g8GV2vZsLrLK9apw8jwDvBr4BnA3MHyWdgX+QZOBvatedU2bPYHurpKPHq+RgDrVFxOQx0PuaL093dHN78UHgyvJgyK3ArlHSnWR7Swl8ayU9WlqgE9baqHRE7Mf24ai07Udtv8P2bwI3Av88Srot5ec2qnuRS8qppyTNBSg/t41XZgJjRDTW1qh012uXrq+kA4CPASu7pDm4PHmHpIOBd1B1waFqZZ5bPp9L1SUfUwJjRDTW1qi0pBuBe4CFkjZLOg9YJun7wKNUD4j8bUl7jKThJ0DmAHeXAd7vAH9v+1vl3OXA6ZIeB04v+2PKPcaIaKbFydu2l41y6vNd0m4BlpbPPwCOH+Wa24HTJlKPBMaIaKSa4D1Ys+8SGCOiubxdJyJipLQYIyLq+uwFEb1IYIyIhnp+DrpvJDBGRHPpSkdE1DhLG0RE7C0txoiIDoMVFxMYI6I5DQ1WXzqBMSKaMZngHRFRJ5wJ3hERexmwwNjGYlgHSfqOpIckbZD0yS5pJOlKSZskrZf0lqblRsQ0Yve29Yk2WowvAKfa3lmWUb1b0jdt31tLcyawoGxvBa4qPyOi3+Ue497Kan87y+7MsnX+aTgLuL6kvVfS4ZLmDi9QExH9bdBGpVt5g7ekGZIepFpLYa3t+zqSzAOerO1vLsciou/12I3uo650K4HR9h7bi4BjgSWSfr0jSbflVbv+liStkHS/pPt/vH1PG9WLiH3JJDCOxfZPgX8Czug4tZmRa8EeS7V2Q7drrLK92Pbio149o83qRcS+MtTj1ifaGJU+StLh5fMrgbdTLVpTdyvw/jI6fSLwbO4vRgwO2T1t/aKNUem5wHWSZlAF2q/avk3S+QC2VwJrqBat2QT8HPhAC+VGxHTRR0GvF22MSq8HTuhyfGXts4ELmpYVEdOQDXva6SdLWg28C9hm+9fLseOp1pI+BHgC+APbOzryzQeuB15D1WlfZfvz5dwngD8BflySX2Z7DWPIutIR0Vx7gy/XsvcYxdXApbZ/A7gF+EiXfLuBS2y/CTgRuEDSm2vnP2d7UdnGDIqQwBgRbWgpMNq+C3im4/BC4K7yeS3wu13ybbW9rnx+DthIgymBCYwR0YyBIfe2wezh6XhlW9FDCY8A7y6fz2bkDJe9SDqO6vZefT71heVx5NWSjhivwATGiGjI4KHeNnh6eDpe2Vb1UMAHqbrGDwCHArtGSyjpEOAm4OLafcirgDcAi4CtwGfHKzBv14mIZkxrgy9dL28/CrwDQNIbgd/qlq68q+Em4AbbN9fyP1VL8yXgtvHKTIsxIprbh0++SDq6/DwA+BjVCHVnGgHXABttX9Fxbm5t971UXfMxJTBGRHMtBUZJNwL3AAslbZZ0HrBM0vepHhzZAvxtSXuMpOER5pOA5cCpkh4s29Jy7tOSHpa0HjgF+NB49UhXOiIaau85aNvLRjn1+S5pt1A9OILtu+n+TgZsL59oPRIYI6IZAwP22rEExohoLo8ERkTUtfdI4HSRwBgRzRjsBMaIiJGG0pWOiBgp9xgjImrsjEpHROwlLcaIiDrjPYO1cF0CY0Q0M/zasQGSwBgRzQ3YdJ02Vgk8SNJ3JD0kaYOkT3ZJc7KkZ2sPd3+8abkRMT0Y8JB72vpFGy3GF4BTbe8s70O7W9I3bd/bke7btt/VQnkRMZ3YA9dibGOVQAM7y+7MsvXPn4aIaGzQBl/kFobZy5rSDwC/CnzR9kc7zp9M9WbdzVTvU/uw7Q2jXGsFMLwOxELgscYV7M1s4OlJKmsy5Xv1n8n8bq+zfVSTC0j6FlWde/G07c5VAKedVgLjixeTDqda3vA/2n6kdvxVwFDpbi8FPm97QWsFt0DS/bYXT3U92pbv1X8G+bv1i1bf4G37p8A/0bEurO0dtneWz2uAmZJ6/QsTETGp2hiVPqq0FJH0SuDtVK8gr6d5TVmTAUlLSrnbm5YdEbEvtDEqPRe4rtxnPAD4qu3bJJ0PYHsl8D7gTyXtBp4HznGbffh29LKMYz/K9+o/g/zd+kKr9xgjIgZBVgmMiOiQwBgR0WG/D4ySzpD0mKRNki6d6vq0RdJqSdskjbu4eD+RNF/SnZI2lkdQL5rqOrWhl0drY/Ls1/cYy4DR94HTqSaffxdYZvt7U1qxFkj6D1RPJF1v+9enuj5tkTQXmGt7naRDqR4seE+//zcrszYOrj9aC1zU5dHamAT7e4txCbDJ9g9s7wK+Apw1xXVqhe27gGemuh5ts73V9rry+TlgIzBvamvVnCt5tHaa2N8D4zzgydr+ZgbgH9n+QtJxwAnAfVNclVZImiHpQWAbsNb2QHyvfrS/B0Z1OZa/0n1A0iFUz99fbHvHVNenDbb32F4EHAsskTQwt0D6zf4eGDcD82v7x1K95CKmsXIP7ibgBts3T3V92jbao7Uxefb3wPhdYIGk10uaBZwD3DrFdYoxlEGKa4CNtq+Y6vq0pZdHa2Py7NeB0fZu4ELgdqqb+F8d7XVo/UbSjcA9wEJJmyWdN9V1aslJwHLg1Nob4ZdOdaVaMBe4U9J6qj/Ya23fNsV12m/t19N1IiK62a9bjBER3SQwRkR0SGCMiOiQwBgR0SGBMSKiQwJjRESHBMaIiA7/H5oVhcm12pblAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Шаг по времени: 4 \n",
      "Распределение давления, МПа\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>20.000</td>\n",
       "      <td>19.872</td>\n",
       "      <td>19.590</td>\n",
       "      <td>19.039</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>20.128</td>\n",
       "      <td>20.000</td>\n",
       "      <td>19.811</td>\n",
       "      <td>19.615</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>20.410</td>\n",
       "      <td>20.189</td>\n",
       "      <td>20.000</td>\n",
       "      <td>19.873</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>20.961</td>\n",
       "      <td>20.385</td>\n",
       "      <td>20.127</td>\n",
       "      <td>20.000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       1      2      3      4\n",
       "1 20.000 19.872 19.590 19.039\n",
       "2 20.128 20.000 19.811 19.615\n",
       "3 20.410 20.189 20.000 19.873\n",
       "4 20.961 20.385 20.127 20.000"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUYAAAEICAYAAAAjhV3sAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAAAe3UlEQVR4nO3df7ReVX3n8feHmIgFFDAQQwjiaiMjtotgs1Ic1sxCEYXUirbSBdNG2tJJtdABB11lGMdq27GMq2KlWDIRUnCkWF2AUozSlMGJrAEEssKPGJCUQYmJpAlCglJi7v3MH2dfPPfJc+99LufcX08+L9ZZz/mx99n7ueF+7z57n3O2bBMRET9zwFRXICJiuklgjIjokMAYEdEhgTEiokMCY0REhwTGiIgOCYwRER0SGCeBpCckva22fYqkLR1prpX055Nfu+iVpHMlWdLvT3VdYmIlMEb0QNJhwH8BNk51XWLiJTBOU5K+LOmHkp6VtE7SG2vHrpW0R9JzteVVXc7xO5IGOtINSjqlHH+5pL+StLUsfyXp5SPUp36uXZL+t6QF5dgpkrZIulTSjtJC/q1a3pdL+ktJ35f0lKSVkl5Ry2tJn6mlP77s+0Jt30mS/q+kZyQ9MPQdyrFv1ltxkt4m6Yna9ostdkkHlzrcOY5/DoC/AK4AdowzX8xACYzT19eBRcCRwHrg+o7jn7R9cG15doTz3FVPB2ytHfuvwEnAYuAEYCnwkVHqdFc5x5HAC8AHa8deA8wFFgDnAqskHVeO/Q/g9aWcXyhpPlrL+y/AGbWg/PvApqGDJQB/Dfhz4HDgQ8CNko4Ypa4j+TDw0/oOSbdKumSkDJKWAkuAlS+hvJiBEhinKdurbe+2/QLwMeCEbq3Chn4L+FPb223/C/BxYHkP+Q4oy86O/f/N9gu2/w9VIPtNSQL+I/BB20/b3g18Aji7lm8PsAb4dUlzgDOAr9SO/zawxvYa24O21wL3AcvG82UlzQPOAy6v77f9TtuXjZBnFvA3wB/ZHhxPeTFzvWyqKxD7Kr+M/x04CzgCGPqFnAuM1DJ8KY4Cvlfb/l7ZN5KTJD0D/BzwBPA7tWM/sv3jLuc6oqS/v4qRAAiY1XHuq4G/AgaAb1AFyyGvBc6S9Gu1fbOBO2rbV0j6y7L+Mrpf8n4M+Gvg6e5fr6s/BB60fdc48sQMlxbj9PQfgDOBtwGvAo4t+zVShpdoK1XQGXIMwy+1O91t+1DgQOALwLW1Y4dJOqjLuXYAzwNvtH1oWV5VLslfZPthqgD6EaogWfck8L9q+Q+1fVBHK+8/DR0D3t2l7q8H3kHVTzgepwLvKf29PwT+LfApSVeO8zwxgyQwTp7Zkg6UdCAwB2Bou+yrt6AOoerD20kVLD4xQXW6AfiIpCMkzaXq9/vCGHkATNWy6+zj+7ikOZL+HfBO4Mvl8vNzwKclHQlVn6Gkd3Q57yeAf7LdOfL7BeDXJL1D0qzyMztF0tE9f9Mq4P6p7efHkQeqVvEbqPpHF1Ndwn+cqn82+lQC4+RZQ9Vyeh64jWoA4vna8tu1tJ+nuhT9AfAd4O4JqtOfU/2iPwg8RDXIM9q9lG+W9BzV5fyvAxfUjv0Q+BFVK/F64P22HynH/hjYDNwtaRfwT8BxdLB9q+3/3GX/k1Qt6EupBmqepBpEGc//vzupfq77kPR1SZd2O2b7Gds/HFqoLvF3jTLYFX1AeVFtNFVunfmC7fG04CKmrbQYIyI6NBqVlnQ48PdUgwNPAL9p+0dd0j0B7Kbql9pre0mTciMiJlKjS2lJnwSetn1ZuUH2MNt/3CXdE8AS23lqICKmvaaX0mcC15X16+h+m0RExIzStMX4TLlvbGj7R7YP65Lu/1GNWBr4n7ZXjXLOFcAKAM2Z88uz5x35kus3XWnvVNdg4sx6YaprMDFmPf3jsRPNQP/Kj9njFxrdH/uOtxzknU8P9JT2/gdfuM326U3Kmwxj9jFK+ieq52A7jec+rpNtby33sa2V9Ijtdd0SlqC5CuDlxyz0URdfNI5iZoY5P+rfMa9DN/fnU3Ov/LuJumNqat3j2xufY+fTA3z7tmN6Sjtr/mNzGxc4CcYMjLbfNtKx8paS+ba3SZoPbB/hHFvL53ZJN1O9rKBrYIyImcXAIP31B7Fp0+UWqjepUD6/2plA0kGSDhlaB94OPNyw3IiYJoz5qQd6WmaKpoHxMuA0SY8Bp5VtJB0laU1JMw+4U9IDwLeBr9n+RsNyI2IaGezxv5mi0X2MtndSPWTfuX8r5ZVQth+netdfRPQhYwb67Am6vHYsIhobJIExIuJF1auWEhgjIoZJizEiosbAT9PHGBHxM8a5lI6IGMYw0F9xMYExIpqpnnzpLwmMEdGQGGh9nraplcAYEY1Ugy8JjBERL6ruY0xgjIgYZjAtxoiIn0mLMSKigxEDfTbhaAJjRDSWS+mIiBoj9njWVFejVf3V/o2ISVfd4H1AT8toJC2UdIekTZI2Srqw7D9c0lpJj5XPbhPuHSdpQ23ZJemicuxjkn5QO7ZsrO+UFmNENNbS4Mte4GLb68t0KPdLWgv8DnB7bf76S4Bh89fbfhRYDCBpFvAD4OZakk/b/steK5LAGBGN2GLAzS8+bW8DtpX13ZI2AQuo5q8/pSS7DvgmHYGxw6nAP9v+3kutSyuX0pJOl/SopM0loncel6QryvEHJb2pjXIjYnoYRD0twFxJ99WWFd3OJ+lY4ETgHmBeCZpDwXOsyebPBm7o2HdBiT2ru12Kd2ocGEuz9bPAGcDxwDmSju9IdgawqCwrgKualhsR00M1+PKynhZgh+0ltWVV5/kkHQzcCFxke9d46iJpDvAu4Mu13VcBP091qb0N+NRY52mjxbgU2Gz7cdt7gC9SNX3rzgQ+78rdwKFlHuqImOHaGnwBkDSbKiheb/umsvupoXgx2vz1xRnAettPvVg/+ynbA7YHgc9RxaxRtREYFwBP1ra3lH3jTRMRM9SA1dMyGkkCrgE22b68dmjM+etrzqHjMrqjEfYeepjXvo3Bl27ftvO1lb2kqRJWfQ4rAGYdNmZXQERMsRaffDkZWA48JGlD2Xcp1Xz1X5J0HvB94Cyo5q8Hrra9rGz/HNX89n/Qcd5PSlpMFXOe6HJ8H20Exi3Awtr20cDWl5AGgNLnsArg5ccs7LP3Akf0p8F2RqXvpHsjCsaYv75s/wR4dZd0y8dblzbC/L3AIkmvKx2fZ1M1fetuAd5XRqdPAp4dGmWKiJmteonEAT0tM0XjFqPtvZIuAG4DZgGrbW+U9P5yfCWwhiqybwZ+Avxu03IjYnow4qd99khgKzd4215DFfzq+1bW1g2c30ZZETG92LRyg/d0kidfIqKhF2/e7hsJjBHRiEmLMSJiHzNpYKUXCYwR0YhRXlQbEVFXTZ/aX6Gkv75NREwBZTKsiIg6086TL9NJAmNENJYWY0REja20GCMi6qrBlzwSGBFR086cL9NJAmNENFINvqSPMSJimDz5EhFRkydfIiK66GWiq5kkgTEiGrHhp4P9FRj769tExKSrLqUP6GkZjaSFku6QtEnSRkkXlv2HS1or6bHy2XWWPElPSHpI0gZJ99X295S/LoExIhobKM9Lj7WMYS9wse03ACcB50s6HrgEuN32IuD2sj2St9hebHtJbd948gMJjBHR0NDtOr0so57H3mZ7fVnfDWyimn/+TOC6kuw64N3jrOK487cSGCWdLulRSZsl7RONJZ0i6dnSxN0g6aNtlBsR08G4LqXnSrqvtqzoekbpWOBE4B5g3tCsouXzyBEqYuAfJd3fcd5e87+o8eCLpFnAZ6kmut4C3CvpFtvf6Uj6LdvvbFpeREw/45jzZUfHZe4+JB0M3AhcZHuX1PO5T7a9VdKRwFpJj9he12vmujZGpZcCm20/DiDpi1RN187AOH6zzOChP218mulmD7OnugoTZvfC/uydGVz+5qmuwoQY+Nrdjc9RjUq386y0pNlUQfF62zeV3U9Jmm97m6T5wPbu9fDW8rld0s1UsWldr/nr2vi/eAHwZG17S9nX6c2SHpD0dUlvHOlkklYMNbMHdv+4hepFxEQausG7aR+jqqbhNcAm25fXDt0CnFvWzwW+2iXvQZIOGVoH3g483Gv+Tm20GLt9W3dsrwdea/s5ScuArwCLup3M9ipgFcDLX3d053kiYhpqafrUk4HlwEOSNpR9lwKXAV+SdB7wfeAsAElHAVfbXgbMA24ul90vA/7O9jfKObrmH00bgXELsLC2fTSwtZ7A9q7a+hpJfyNpru0dLZQfEVOorZdI2L6T7g0tgFO7pN8KLCvrjwMnjHDend3yj6aNS+l7gUWSXidpDnA2VdP1RZJeU5rJSFpayt3ZQtkRMQ20cYP3dNK4xWh7r6QLgNuAWcBq2xslvb8cXwm8F/iApL3A88DZtnOZHNEHbLF3BgW9XrTyrLTtNcCajn0ra+tXAle2UVZETD95u05ERE1eVBsR0UUCY0RETV5UGxHRRUv3MU4bCYwR0YgNe/vsRbUJjBHRWC6lIyJq0scYEdGFExgjIobL4EtERI2dPsaIiA5iIKPSERHDpY8xIqImz0pHRHRy1c/YTxIYI6KxjEpHRNQ4gy8REfvqt0vp/grzETElbPW0jEbSQkl3SNokaaOkC8v+wyWtlfRY+Tys17zl2Mck/UDShrIsG+v7tBIYJa2WtF3SwyMcl6QrJG2W9KCkN7VRbkRMPbudwAjsBS62/QbgJOB8SccDlwC3214E3F62e8075NO2F5dlTZf8w7TVYrwWOH2U42dQzSO9CFgBXNVSuRExDQxaPS2jsb3N9vqyvhvYBCwAzgSuK8muA949jrwvSSuB0fY64OlRkpwJfN6Vu4FDJc1vo+yImHp2bwswV9J9tWVFt/NJOhY4EbgHmGd7W1WOtwFHjlaXjrxDLihXq6u7XYp3mqw+xgXAk7XtLYwQzSWtGPqhDez+8aRULiJeOiMGBw/oaQF22F5SW1Z1nk/SwcCNwEW2d42nLiPkvQr4eWAxsA341FjnmazA2K0N3XUcy/aqoR/arEMOmuBqRUQb3OMyFkmzqQLb9bZvKrufGrrCLJ/bx5EX20/ZHrA9CHwOWDpWPSYrMG4BFta2jwa2TlLZETGRWhp8kSTgGmCT7ctrh24Bzi3r5wJfHUdeOrrt3gN0HSSum6zAeAvwvjI6fRLw7FCfQUT0gXaajCcDy4G3dtxacxlwmqTHgNPKNpKOkrRmjLwAn5T0kKQHgbcAHxyrIq3c4C3pBuAUqo7VLcCfALMBbK8E1gDLgM3AT4DfbaPciJge2ni7ju076d7tBnBql/RbqeLKqHltLx9vXVoJjLbPGeO4gfPbKCsiphcDg4N5Vjoi4mcM5LVjERHD9duz0gmMEdFcAmNERF1Pz0HPKAmMEdFcWowRETUGZ1Q6IqJTAmNExHC5lI6I6JDAGBFRkxu8IyL2lRu8IyI6ZVQ6ImI4pcUYEVHT6+u5Z5AExohoSBl8iYjYR1qMEREdBqe6Au1KYIyIZvrwPsZWJsMqk1hvl9R19i1Jp0h6tjZJzUfbKDcipge5t2XUc0gLJd0haZOkjZIuLPsPl7RW0mPl87AR8p8u6VFJmyVdUtvfU/66tmYJvBY4fYw037K9uCx/2lK5ETEdtDNL4F7gYttvAE4Czpd0PHAJcLvtRcDtZXsYSbOAzwJnAMcD55S89JK/UyuB0fY64Ok2zhUR+yfb22yvL+u7gU3AAuBM4LqS7Drg3V2yLwU2237c9h7giyUfPeYfZjL7GN8s6QFgK/Ah2xu7JZK0AlgBcOC8Q3jj67ZOYhUnx6Nb5011FSbM8xw41VWYEAMHTtYU7JNrYE475xnHDd5zJd1X215le9U+55OOBU4E7gHmDc1Db3ubpCO7nHcB8GRtewvwK2W9l/zDTFZgXA+81vZzZRLsrwCLuiUsP6RVAK/6N/P67CaAiD5kxvNI4A7bS0ZLIOlg4EbgItu7pJ7O3S3RS44fk/Jn0PYu28+V9TXAbElzJ6PsiJgE7fQxImk2VVC83vZNZfdTkuaX4/OB7V2ybgEW1raPpro67TX/MJMSGCW9RiXsS1payt05GWVHxMRraVRawDXAJtuX1w7dApxb1s8Fvtol+73AIkmvkzQHOLvk6zX/MK1cSku6ATiFqv9gC/AnwGwA2yuB9wIfkLQXeB442+63FxVF7Mfa+W0+GVgOPCRpQ9l3KXAZ8CVJ5wHfB84CkHQUcLXtZbb3SroAuA2YBayujWN0zT+aVgKj7XPGOH4lcGUbZUXENNRCYLR9JyNPHnNql/RbgWW17TXAmi7pdnbLP5o8+RIRjfRymTzTJDBGRHN5UW1ExHBpMUZEdEpgjIioSR9jREQXCYwREcOpz15U259PxkdENJAWY0Q0l0vpiIiaDL5ERHSRwBgR0SGBMSLiZ0T/jUonMEZEM+ljjIjoIoExIqJDAmNExHC5lI6I6NRngbHxI4GSFkq6Q9ImSRslXdgljSRdIWmzpAclvalpuRExTbgale5lmSnaeFZ6L3Cx7TcAJwHnSzq+I80ZVPNILwJWAFe1UG5ETBftTZ+6WtJ2SQ/X9p0g6S5JD0n6B0mv7JLvOEkbassuSReVYx+T9IPasWWd+Ts1Doy2t9leX9Z3A5uABR3JzgQ+78rdwKFD87xGxMzXxvSpxbXA6R37rgYusf1LwM3Ahzsz2X7U9mLbi4FfBn5S0g759NDxMmnWqFp9u46kY4ETgXs6Di0Anqxtb2Hf4Dl0jhWS7pN0355nnm+zehExUVpqMdpeBzzdsfs4YF1ZXwv8xhinORX4Z9vf663y+2otMEo6GLgRuMj2rs7DXbJ0/THZXmV7ie0lcw59RVvVi4iJ0mtQrH7j5w41fMqyoocSHgbeVdbPAhaOkf5s4IaOfReU8Y3Vkg4bq8BWAqOk2VRB8XrbN3VJsoXhX+ZoYGsbZUfE1BLjupTeMdTwKcuqHor4Paqxi/uBQ4A9I9ZFmkMVRL9c230V8PPAYmAb8KmxCmxjVFrANcAm25ePkOwW4H1ldPok4Fnb25qWHRHTQ4t9jPuw/Yjtt9v+ZaqW4D+PkvwMYL3tp2r5n7I9YHsQ+BywdKwy27iP8WRgOfCQpA1l36XAMaVSK4E1wDJgM1Wn6O+2UG5ETBcTeB+jpCNtb5d0APARYOUoyc+h4zJa0vxaQ+w9VJfmo2ocGG3fSfc+xHoaA+c3LSsipqmWAqOkG4BTqPoitwB/AhwsaSh+3AT8bUl7FHC17WVl++eA04A/6DjtJyUtLrV8osvxfeTJl4hopsW369g+Z4RDn+mSdivVlejQ9k+AV3dJt3y89UhgjIjm+uyRwATGiGhsJj3u14sExohoLG/XiYio6/GplpkkgTEimktgjIj4maEnX/pJAmNENKbB/oqMCYwR0Uz6GCMi9pVL6YiITgmMERHDpcUYEdEpgTEiosZ5JDAiYpjcxxgR0Y37KzImMEZEY2kxRkTU9eEN3m1MhrVQ0h2SNknaKOnCLmlOkfSspA1l+WjTciNi+tBgb8tM0UaLcS9wse31kg4B7pe01vZ3OtJ9y/Y7WygvIqaZmRT0etHGZFjbqOZqxfZuSZuABUBnYIyIfmQy+DIaSccCJwL3dDn8ZkkPAFuBD9neOMI5VgArAI5Z8DJuff3X26zitPCOgf5tOH/3X18z1VWYEHv6tDveLX2ttgZfJK0G3glst/2LZd8JVFOmHkw1y99v2d7VJe8TwG5gANhre0nZfzjw98CxJf9v2v7RaPVo3MdYq9TBwI3ARV0qvR54re0TgL8GvjLSeWyvsr3E9pIjXj2rrepFxERyj8vYrgVO79h3NXCJ7V8CbgY+PEr+t9hePBQUi0uA220vAm4v26NqJTBKmk0VFK+3fVPncdu7bD9X1tcAsyXNbaPsiJhaQzd497KMxfY64OmO3ccB68r6WuA3xlnFM4Hryvp1wLvHytDGqLSAa4BNti8fIc1rSjokLS3l7mxadkRMAzYa7G0B5kq6r7as6KGEh4F3lfWzgIUj1QT4R0n3d5x3XhkLGRoTOXKsAtvoYTgZWA48JGlD2XcpcEypyErgvcAHJO0FngfOtvustzZif9b7b/OOjsvcXvwecEW5ze8WYM8I6U62vVXSkcBaSY+UFui4tTEqfSdVa3q0NFcCVzYtKyKmp4l88sX2I8DbASS9HvjVEdJtLZ/bJd0MLKW6BH9K0nzb2yTNB7aPVWZrgy8RsZ8yMOjelpegtACRdADwEaoR6s40B5X7qJF0EFUgfbgcvgU4t6yfC3x1rDITGCOiuZZGpSXdANwFHCdpi6TzgHMkfRd4hOp2v78taY+StKZknQfcWW4J/DbwNdvfKMcuA06T9BhwWtkeVX/enBURk6qtS2nb54xw6DNd0m4FlpX1x4ETRjjnTuDU8dQjgTEiGsv0qRERdX34dp0ExohopLrBu78iYwJjRDSXt+tERAyXFmNERF36GCMiOjmj0hER+8ildEREjTO1QUTEvtJijIjo0F9xMYExIprTYH9dSycwRkQzJjd4R0TUCecG74iIffRZYGxjMqwDJX1b0gOSNkr6eJc0knSFpM2SHpT0pqblRsQ0Yve2zBBttBhfAN5q+7kyjeqdkr5u++5amjOARWX5FeCq8hkRM136GPdVZvt7rmzOLkvnn4Yzgc+XtHdLOnRocpqm5UfE1Ou3UelW5nyRNKtMnbodWGv7no4kC4Ana9tbyr6ImPF6vIyeQZfSrQRG2wO2FwNHA0sl/WJHkm7Tq3b9KUlaMTQZ97/sHGijehExkUxrgVHSaknbJT1c23eCpLskPSTpHyS9sku+hZLukLSpjHVcWDv2MUk/kLShLMvGqkerswTafgb4JnB6x6EtwMLa9tFUs311O8cq20tsLzni1bParF5ETJTBHpexXcu+8eNq4BLbvwTcDHy4S769wMW23wCcBJwv6fja8U/bXlyWNV3yD9PGqPQRkg4t668A3kY1zWHdLcD7yuj0ScCz6V+M6B+ye1rGYnsd8HTH7uOAdWV9LfAbXfJts72+rO8GNtGgu66NFuN84A5JDwL3UvUx3irp/ZLeX9KsAR4HNgOfA/6whXIjYrro/VJ67lBXWVlW9HD2h4F3lfWzGH71uQ9JxwInAvWxjgvKrYKrJR02VoFtjEo/WCrRuX9lbd3A+U3LiohpyIaBnkeld9heMs4Sfg+4QtJHqa4+94yUUNLBwI3ARbZ3ld1XAX9G1Rv6Z8CnyjlHlCdfIqK5CRxxtv0I8HYASa8HfrVbunIf9Y3A9bZvquV/qpbmc8CtY5XZ6uBLROynJvB2HUlHls8DgI8AK7ukEXANsMn25R3H5tc230N1aT6qBMaIaMbAoHtbxiDpBuAu4DhJWySdB5wj6btUg7pbgb8taY+SNDTCfDKwHHhrl9tyPllu9XkQeAvwwbHqkUvpiGjI4HaefLF9zgiHPtMl7VZgWVm/k+73S2N7+XjrkcAYEc2Y8Qy+zAgJjBHR3Ax63K8XCYwR0VwCY0RE3cx6QUQvEhgjohkDffbasQTGiGguLcaIiLpxPRI4IyQwRkQzBrd0H+N0kcAYEc318FTLTJLAGBHNpY8xIqLGzqh0RMQ+0mKMiKgzHuiviesSGCOimaHXjvWRBMaIaK7PbtdpY5bAAyV9W9IDZT7Xj3dJc4qkZ2svkPxo03IjYnow4EH3tMwUbbQYXwDeavu5MufCnZK+bvvujnTfsv3OFsqLiOnE7b2odrpoY5ZAA8+VzdllmTl/GiKisX4bfJFbGGaXNAu4H/gF4LO2/7jj+ClUs3dtoZqz4UO2N45wrhXA0FyzxwGPNq5gb+YCOyaprMmU7zXzTOZ3e63tI5qcQNI3qOrcix22T29S3mRoJTC+eDLpUOBm4I9sP1zb/0pgsFxuLwM+Y3tRawW3QNJ9L2G+22kv32vm6efvNlO0Okug7WeAbwKnd+zfZfu5sr4GmC2p178wERGTqo1R6SNKSxFJrwDeRjXNYT3Na8q8r0haWsrd2bTsiIiJ0Mao9HzgutLPeADwJdu3Sno/gO2VwHuBD0jaCzwPnO02r+HbsWqqKzBB8r1mnn7+bjNCq32MERH9oNU+xoiIfpDAGBHRYb8PjJJOl/SopM2SLpnq+rRF0mpJ2yU9PHbqmUPSQkl3SNpUHkG9cKrr1IZeHq2NybNf9zGWAaPvAqdR3Xx+L3CO7e9MacVaIOnfUz2R9HnbvzjV9WmLpPnAfNvrJR1C9WDBu2f6v1m5a+Og+qO1wIVdHq2NSbC/txiXApttP257D/BF4MwprlMrbK8Dnp7qerTN9jbb68v6bmATsGBqa9WcK3m0dprY3wPjAuDJ2vYW+uCXbH8h6VjgROCeKa5KKyTNkrQB2A6std0X32sm2t8Do7rsy1/pGUDSwVTP319ke9dU16cNtgdsLwaOBpZK6psukJlmfw+MW4CFte2jqV5yEdNY6YO7Ebje9k1TXZ+2jfRobUye/T0w3gsskvQ6SXOAs4FbprhOMYoySHENsMn25VNdn7b08mhtTJ79OjDa3gtcANxG1Yn/pZFehzbTSLoBuAs4TtIWSedNdZ1acjKwHHhr7Y3wy6a6Ui2YD9wh6UGqP9hrbd86xXXab+3Xt+tERHSzX7cYIyK6SWCMiOiQwBgR0SGBMSKiQwJjRESHBMaIiA4JjBERHf4/QyoBeR3L+pUAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Шаг по времени: 5 \n",
      "Распределение давления, МПа\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>1</th>\n",
       "      <th>2</th>\n",
       "      <th>3</th>\n",
       "      <th>4</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>20.000</td>\n",
       "      <td>19.867</td>\n",
       "      <td>19.577</td>\n",
       "      <td>19.021</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>20.133</td>\n",
       "      <td>20.000</td>\n",
       "      <td>19.804</td>\n",
       "      <td>19.604</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>20.423</td>\n",
       "      <td>20.196</td>\n",
       "      <td>20.000</td>\n",
       "      <td>19.869</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>20.979</td>\n",
       "      <td>20.396</td>\n",
       "      <td>20.131</td>\n",
       "      <td>20.000</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       1      2      3      4\n",
       "1 20.000 19.867 19.577 19.021\n",
       "2 20.133 20.000 19.804 19.604\n",
       "3 20.423 20.196 20.000 19.869\n",
       "4 20.979 20.396 20.131 20.000"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUYAAAEICAYAAAAjhV3sAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAAAe1UlEQVR4nO3dfbBdVZ3m8e9DDL4QNGBICCGKY0ca1CJoKsaipieCKKRp0W61oO1Iqz0ZLJgBBy0p2/atHIuhWlTUIh0hgiODrQVIGqN0ZHAiNYBAipfECxIZlEvSiQQhiTKE3PvMH3tf2Pfk3HvPZe/7dvJ8rF1nv6y11zo38jtrr7X3XrJNREQ874CJrkBExGSTwBgR0SKBMSKiRQJjRESLBMaIiBYJjBERLRIYIyJaJDCOA0mPSHp7ZXuppN6WNFdK+uL41y5GIsmS/iBpd7lcPtF1irH1oomuQMQUcZztzRNdiRgfaTFOUpJ+IOnfJD0lab2k11eOXSlpT6UFs1vSK9qc428l9bWk65e0tDz+YklflbSlXL4q6cVD1Kd6rp2S/pekeeWxpZJ6JX1K0uNlC/kDlbwvlvSPkn4raZuklZJeWslrSV+rpD+23Pfdyr4lkv6PpCcl3TvwHcpjP5P0d5Xtt0t6pLL9XItd0oyyDreO4p8j9jMJjJPXj4EFwGxgA3B1y/GLbc+oLE8NcZ7bqumALZVjfw8sARYCxwGLgU8PU6fbynPMBp4BPlY5djgwC5gHnAWsknR0eey/A68ry/mTMs1nKnl/B5xaCcp/B/QMHCwD8I+ALwKHAh8HrpV02DB1HcongGerOyTdKOnCEfKtL3+orpN01AsoN6aQBMZJyvZq27tsPwN8DjiuXauwpg8AX7C93fbvgM8DyzvId0C57GjZ/w+2n7H9vykC2fslCfiPwMdsP2F7F/Al4IxKvj3AWuAvJR0InAr8sHL8b4C1ttfa7re9DrgLWDaaLytpDvAR4JLqftun2b5omKz/ATgK+FOKH5YbJaUbqovlH3cSkjQN+G/A+4DDgP7y0CxgqJbhC3EE8JvK9m/KfUNZIulJ4GXAI8DfVo793vYf2pzrsDL93UWMBEDAtJZzXw58FegDfkIRLAe8GnifpL+o7JsO3FLZvlTSP5brLwIeb1P/zwFfB55o//Xas72+XN0j6TxgJ3AMcP9ozhNTR1qMk9NfA6cDbwdeQdFagSKgNGkLRdAZ8CoGX2q3ut32TOAlwHeBKyvHDpF0UJtzPQ48Dbze9sxyeUV5Sf4c2xspAuinKYJk1aPA/6jkn2n7oJZW3n8ZOAa8u03dXwe8E7h0mO/XKdP8v0VMIgmM42e6pJdIeglwIMDAdrmv2oI6mKIPbwdFsPjSGNXpGuDTkg6TNIui3++7I+SBIjD0UbQGqz4v6UBJ/x44DfiB7X7gW8BXJM2Gos9Q0jvbnPdLwE9tb2rZ/13gLyS9U9K08m+2VNKRHX/TIuB+wfbTo8iDpNdLWliWOwP4MvAYlT7Q6D4JjONnLUXL6WngJooBiKcry99U0n6H4lL0MeCXwO1jVKcvUvTV3UdxWbih3DeUt0raTXE5/5fAuZVj/wb8nqKVeDVwtu0HymOfBDYDt0vaCfwUOJoWtm+0/V/b7H+UogX9KYqBmkcpBlFG8//fHRR/131I+rGkTw2Rbw7wzxSXzw9TtN5Ps/3sEOmjCygvqo26yltnvmt7NC24iEkrLcaIiBa1RqUlHUpxmXEUxSjl+23/vk26R4BdFP1Se20vqlNuRMRYqnUpLeli4AnbF5U3yB5i+5Nt0j0CLLLd7haKiIhJpe6l9OnAVeX6VbS/TSIiYkqp22J8srxvbGD797YPaZPu/1KMWBr4J9urhjnnCmAFgA488M3T58x+wfWbrLR3omswdqY9M9E1GBvTnvjDyImmoP/HH9jjZ2rdk/nOtx3kHU/0dZT27vueucn2KXXKGw8j9jFK+inFc7Ct/n4U5Zxge0t5H9s6SQ9UniYYpAyaqwBe/Kr5PuKC80dRzNRw4O+7d8xr5kP9Iyeagl5+zVjdMTWx7vDNtc+x44k+fnHTqzpKO23uQ7NqFzgORgyMtt8+1LHyLSVzbW+VNBfYPsQ5tpSf2yVdT/GygraBMSKmFgP9dNcPYt2myxqKN6lQft7QmkDSQZIOHlgH3gFsrFluREwSxjzrvo6WqaJuYLwIOFnSQ8DJ5TaSjpC0tkwzB7hV0r3AL4Af2f5JzXIjYhLp7/B/U0Wt+xht7wBOarN/C+UroWw/TPGuv4joQsb0ddkTdHntWETU1k8CY0TEc4pXLSUwRkQMkhZjRESFgWfTxxgR8TzjXEpHRAxi6OuuuJjAGBH1FE++dJcExoioSfR12dxgCYwRUUsx+JLAGBHxnOI+xgTGiIhB+tNijIh4XlqMEREtjOjrsglHExgjorZcSkdEVBixx9MmuhqNSmCMiFqKG7xzKR0RMUi3Db50V5iPiHFniz4f0NEyHEnzJd0iqUfSJknnlfsPlbRO0kPlZ7spmo+WdE9l2Snp/PLY5yQ9Vjm2bKTv1EhglHSKpAclbZZ0YZvjknRpefw+SW9qotyImBz6UUfLCPYCF9g+BlgCnCPpWOBC4GbbC4Cby+1BbD9oe6HthcCbgT8C11eSfGXguO21rflb1b6UljQN+CbFZFi9wJ2S1tj+ZSXZqcCCcnkLcFn5GRFTXDH4Ur9XzvZWYGu5vktSDzAPOB1YWia7CvgZ8MlhTnUS8Gvbv3mhdWmixbgY2Gz7Ydt7gO9RfJGq04HvuHA7MLOchzoipriBwZdOFmCWpLsqy4p255R0FHA8cAcwpwyaA8Fz9ghVOgO4pmXfueXV6up2l+KtmgiM84BHK9u95b7RpomIKarP6mgBHre9qLKsaj2XpBnAtcD5tneOph6SDgTeBfygsvsy4LXAQooW6ZdHOk8TgbFdx0Hrays7SVMklFYM/Jr07f5D7cpFxNgaePKlk2UkkqZTBMWrbV9X7t42cIVZfm4f5hSnAhtsb3uufvY22322+4FvUVzlDquJwNgLzK9sHwlseQFpALC9auDXZNqMgxqoXkSMtX4f0NEyHEkCrgB6bF9SObQGOKtcPwu4YZjTnEnLZXRLt917gI0jfZ8mAuOdwAJJrymbsWdQfJGqNcAHy9HpJcBTA30GETG1FS+RaKTFeAKwHDix5daai4CTJT1EMch7EYCkIyQ9N8Is6WXl8etaznuxpPsl3Qe8DfjYSBWpPZRke6+kc4GbgGnAatubJJ1dHl8JrAWWAZsphtE/VLfciJgcjHi2gUcCbd9K+243KEaaW9NvoYgrA9t/BF7ZJt3y0dalkSdfyvuC1rbsW1lZN3BOE2VFxORiM+LN21NNHgmMiJo6unl7SklgjIhaTFqMERH7yItqIyIqjPKi2oiIqmL61O4KJd31bSJiAqjr3seYwBgRtRhGfKplqklgjIja0mKMiKiwlRZjRERVMfiSWQIjIiqUG7wjIqqKwZf0MUZEDJInXyIiKvLkS0REG/1pMUZEPM+GZ/sTGCMinlNcSndXYOyubxMRE6KvfF56pGU4kuZLukVSj6RNks4r9x8qaZ2kh8rPtvNCS3qknNvlHkl3VfZ3lL8qgTEiahm4XaeTZQR7gQtsHwMsAc6RdCxwIXCz7QXAzeX2UN5me6HtRZV9o8kPNBQYJZ0i6UFJmyXtU6ikpZKeqsz89Zkmyo2IyUCNTJ9qe6vtDeX6LqAHmAecDlxVJrsKePcoKzjq/LX7GCVNA75JMW1hL3CnpDW2f9mS9Oe2T6tbXkRMPk3P+SLpKOB44A5gzsB0y7a3Spo9RDYD/yrJwD/ZXlXu7zT/c5oYfFkMbLb9MICk71FE6NbAOHrTTP/MZ2ufZrLZw/SJrsKY2fWq7uyd6V/+1omuwpjo+9Httc9RjEp3/Kz0rGr/H7CqEsAAkDQDuBY43/ZOqeOge4LtLWXgWyfpAdvrO81c1URgnAc8WtnuBd7SJt1bJd0LbAE+bntTu5NJWgGsAJj2ypkNVC8ixtIob/B+vKX/bxBJ0ymC4tW2ryt3b5M0t2ztzQW2t61HMc80trdLup6i0ba+0/xVTfy8t/uLuGV7A/Bq28cBXwd+ONTJbK+yvcj2omkHH9RA9SJirPWXU6iOtAxHRdPwCqDH9iWVQ2uAs8r1s4Ab2uQ9SNLBA+vAO4CNneZv1URg7AXmV7aPpGgVPsf2Ttu7y/W1wHRJsxooOyImWIOj0icAy4ETKwO1y4CLgJMlPUQxlnERgKQjJK0t884Bbi2vSn8B/Mj2T8pjbfMPp4lL6TuBBZJeAzwGnAH8dTWBpMOBbbYtaTFFQN7RQNkRMQk0cYO37VtpfwUKcFKb9FuAZeX6w8BxQ5x3R7v8w6kdGG3vlXQucBMwDVhte5Oks8vjK4H3Ah+VtBd4GjjDduvldkRMQbbY22VPvjTySGB5eby2Zd/Kyvo3gG80UVZETD55u05EREVeVBsR0UYCY0RERV5UGxHRRtOPBE60BMaIqMWGvXlRbUTEYLmUjoioSB9jREQbTmCMiBgsgy8RERV2+hgjIlqIvoxKR0QMlj7GiIiKPCsdEdHKRT9jN0lgjIjaMiodEVHhDL5EROyr2y6luyvMR8SEsNXRMhxJ8yXdIqlH0iZJ55X7D5W0TtJD5echneYtj31O0mMtE2wNq5HAKGm1pO2SNg5xXJIulbRZ0n2S3tREuREx8exmAiOwF7jA9jHAEuAcSccCFwI3214A3Fxud5p3wFdsLyyXtW3yD9JUi/FK4JRhjp8KLCiXFcBlDZUbEZNAE9On2t5qe0O5vgvoAeYBpwNXlcmuAt49irwvSCOB0fZ64IlhkpwOfMeF24GZkuY2UXZETDy7swWYJemuyrKi3fkkHQUcD9wBzLG9tSjHW4HZw9WlJe+Ac8ur1dXtLsVbjVcf4zzg0cp2L0NEc0krBv5ofbv+MC6Vi4gXzoj+/gM6WoDHbS+qLKtazydpBnAtcL7tnaOpyxB5LwNeCywEtgJfHuk84xUY27Wh245j2V418EebdvBBY1ytiGiCO1xGImk6RWC72vZ15e5tA1eY5ef2UeTF9jbbfbb7gW8Bi0eqx3gFxl5gfmX7SGDLOJUdEWOpocEXSQKuAHpsX1I5tAY4q1w/C7hhFHlp6bZ7D9B2kLhqvALjGuCD5ej0EuCpgT6DiOgCzTQZTwCWAye23FpzEXCypIeAk8ttJB0hae0IeQEulnS/pPuAtwEfG6kijdzgLekaYClFx2ov8FlgOoDtlcBaYBmwGfgj8KEmyo2IyaGJt+vYvpX23W4AJ7VJv4Uirgyb1/by0dalkcBo+8wRjhs4p4myImJyMdDfn2elIyKeZyCvHYuIGKzbnpVOYIyI+hIYIyKqOnoOekpJYIyI+tJijIioMDij0hERrRIYIyIGy6V0RESLBMaIiIrc4B0Rsa/c4B0R0Sqj0hERgyktxoiIik5fzz2FJDBGRE3K4EtExD7SYoyIaNE/0RVoVgJjRNTThfcxNjIZVjmJ9XZJbWffkrRU0lOVSWo+00S5ETE5yJ0tw55Dmi/pFkk9kjZJOq/cf6ikdZIeKj8PGSL/KZIelLRZ0oWV/R3lr2pqlsArgVNGSPNz2wvL5QsNlRsRk0EzswTuBS6wfQywBDhH0rHAhcDNthcAN5fbg0iaBnwTOBU4FjizzEsn+Vs1EhhtrweeaOJcEbF/sr3V9oZyfRfQA8wDTgeuKpNdBby7TfbFwGbbD9veA3yvzEeH+QcZzz7Gt0q6F9gCfNz2pnaJJK0AVgC8dM4M3vjvHhvHKo6PnscOn+gqjJmneclEV2FM9L14vKZgH199BzZznlHc4D1L0l2V7VW2V+1zPuko4HjgDmDOwDz0trdKmt3mvPOARyvbvcBbyvVO8g8yXoFxA/Bq27vLSbB/CCxol7D8I60CmPmns7vsJoCILmRG80jg47YXDZdA0gzgWuB82zuljs7dLtELjh/j8jNoe6ft3eX6WmC6pFnjUXZEjINm+hiRNJ0iKF5t+7py9zZJc8vjc4HtbbL2AvMr20dSXJ12mn+QcQmMkg5XGfYlLS7L3TEeZUfE2GtoVFrAFUCP7Usqh9YAZ5XrZwE3tMl+J7BA0mskHQicUebrNP8gjVxKS7oGWErRf9ALfBaYDmB7JfBe4KOS9gJPA2fY3faiooj9WDP/NZ8ALAful3RPue9TwEXA9yV9BPgt8D4ASUcAl9teZnuvpHOBm4BpwOrKOEbb/MNpJDDaPnOE498AvtFEWRExCTUQGG3fytCTx5zUJv0WYFlley2wtk26He3yDydPvkRELZ1cJk81CYwRUV9eVBsRMVhajBERrRIYIyIq0scYEdFGAmNExGDqshfVdueT8RERNaTFGBH15VI6IqIigy8REW0kMEZEtEhgjIh4nui+UekExoioJ32MERFtJDBGRLRIYIyIGCyX0hERrbosMNZ+JFDSfEm3SOqRtEnSeW3SSNKlkjZLuk/Sm+qWGxGThItR6U6WqaKJZ6X3AhfYPgZYApwj6diWNKdSzCO9AFgBXNZAuRExWTQ3fepqSdslbazsO07SbZLul/Qvkl7eJt/Rku6pLDslnV8e+5ykxyrHlrXmb1U7MNreantDub4L6AHmtSQ7HfiOC7cDMwfmeY2Iqa+J6VNLVwKntOy7HLjQ9huB64FPtGay/aDthbYXAm8G/limHfCVgePlpFnDavTtOpKOAo4H7mg5NA94tLLdy77Bc+AcKyTdJemuPU8+3WT1ImKsNNRitL0eeKJl99HA+nJ9HfBXI5zmJODXtn/TWeX31VhglDQDuBY43/bO1sNtsrT9M9leZXuR7UUHznxpU9WLiLHSaVAs/oufNdDwKZcVHZSwEXhXuf4+YP4I6c8ArmnZd245vrFa0iEjFdhIYJQ0nSIoXm37ujZJehn8ZY4EtjRRdkRMLDGqS+nHBxo+5bKqgyI+TDF2cTdwMLBnyLpIB1IE0R9Udl8GvBZYCGwFvjxSgU2MSgu4AuixfckQydYAHyxHp5cAT9neWrfsiJgcGuxj3IftB2y/w/abKVqCvx4m+anABtvbKvm32e6z3Q98C1g8UplN3Md4ArAcuF/SPeW+TwGvKiu1ElgLLAM2U3SKfqiBciNishjD+xglzba9XdIBwKeBlcMkP5OWy2hJcysNsfdQXJoPq3ZgtH0r7fsQq2kMnFO3rIiYpBoKjJKuAZZS9EX2Ap8FZkgaiB/XAd8u0x4BXG57Wbn9MuBk4D+1nPZiSQvLWj7S5vg+8uRLRNTT4Nt1bJ85xKGvtUm7heJKdGD7j8Ar26RbPtp6JDBGRH1d9khgAmNE1DaVHvfrRAJjRNSWt+tERFR1+FTLVJLAGBH1JTBGRDxv4MmXbpLAGBG1qb+7ImMCY0TUkz7GiIh95VI6IqJVAmNExGBpMUZEtEpgjIiocB4JjIgYJPcxRkS04+6KjAmMEVFbWowREVVdeIN3E5NhzZd0i6QeSZskndcmzVJJT0m6p1w+U7fciJg81N/ZMlU00WLcC1xge4Okg4G7Ja2z/cuWdD+3fVoD5UXEJNNU0JO0GjgN2G77DeW+4ygmwJpBMWfLB9rMXY+kR4BdQB+w1/aicv+hwD8DR5X532/798PVo3aL0fZW2xvK9V1ADzCv7nkjYoowxeBLJ8vIrgROadl3OXCh7TcC1wOfGCb/22wvHAiKpQuBm20vAG4ut4fVaB+jpKOA44E72hx+q6R7gS3Ax21vGuIcK4AVAK+a9yLWLPhJk1WcFN65t3sbzr96Zs5EV2FM7GH6RFdhTLihCNDgZFjryzhSdTSwvlxfB9wE/MMoTns6xcyDAFcBPwM+OVyG2i3GAZJmANcC57dp5m4AXm37OODrwA+HOo/tVbYX2V502CunNVW9iBhL7nB5YTYC7yrX3wfMH6YW/yrp7rKBNWDOwLzS5efskQpsJDBKmk4RFK+2fd0+tbV32t5drq8Fpkua1UTZETGxBm7w7mShmC/6rsqyYvizA/Bh4BxJdwMHA3uGSHeC7TcBp5bp/+yFfqfaDWlJAq4AemxfMkSaw4Ftti1pMUVA3lG37IiYBOzRvKj28Zb+vw5O7weAdwBIeh3w50Ok21J+bpd0PbCY4hJ8m6S5trdKmgtsH6nMJlqMJwDLgRMrt+Msk3S2pLPLNO8FNpZ9jJcCZ9hddqt8xP5sDC+lJc0uPw8APk0xQt2a5qDyrhgkHUQRSDeWh9cAZ5XrZwE3jFRm7Raj7VspWtPDpfkG8I26ZUXE5NTU4IukaygGSmZJ6gU+C8yQdE6Z5Drg22XaI4DLbS8D5gDXFxewvAj4n7YHRm4vAr4v6SPAbyn6KYeVJ18ioh4DDc35YvvMIQ59rU3aLcCycv1h4LghzrkDOGk09UhgjIj6uqxjLIExImrLSyQiIlpk+tSIiKoufLtOAmNE1FLc4N1dkTGBMSLqm0KvFOtEAmNE1JYWY0REVfoYIyJajepZ6SkhgTEi6suldEREhafWfC6dSGCMiPrSYoyIaNFdcTGBMSLqU393XUsnMEZEPSY3eEdEVAnnBu+IiH10WWCsPeeLpJdI+oWkeyVtkvT5Nmkk6VJJmyXdJ+lNdcuNiEnE7myZIppoMT4DnGh7dzmN6q2Sfmz79kqaU4EF5fIW4LLyMyKmui7sY6zdYnRhd7k5vVxafxpOB75Tpr0dmFlOYxgRXUD9/R0tI55HWi1pu6SNlX3HSbpN0v2S/kXSy9vkmy/pFkk95ZXreZVjn5P0WHUW05Hq0cT0qUiaJukeivla19m+oyXJPODRynZvuS8iprwOL6M7u5S+EjilZd/lwIW23whcD3yiTb69wAW2jwGWAOdIOrZy/Cu2F5bL2pEq0UhgtN1neyFwJLBY0htakrSbXrXtX0nSCkl3Sbrrdzv6mqheRIwl01hgtL0eeKJl99HA+nJ9HfBXbfJttb2hXN8F9FCj8dVIYBxg+0ngZ+wb8XuB+ZXtI4EtQ5xjle1Fthcd9sppTVYvIsZKf4dLMV/0XZVlRQdn3wi8q1x/H4NjyT4kHQUcD1SvXM8tB35XSzpkpAKbGJU+TNLMcv2lwNuBB1qSrQE+WI5OLwGesr21btkRMTnI7mgBHh9o+JTLqg5O/2GKS+O7gYOBPUPWQ5oBXAucb3tnufsy4LXAQmAr8OWRCmxiVHoucJWkaRSB9vu2b5R0NoDtlcBaiomxNwN/BD7UQLkRMVmM4a04th8A3gEg6XXAn7dLV94Vcy1wte3rKvm3VdJ8C7hxpDJrB0bb91E0W1v3r6ysGzinblkRMQnZ0Dd29+tImm17u6QDgE8DK9ukEXAF0GP7kpZjcytXqO+huDQfVqN9jBGxn2po8EXSNcBtwNGSeiV9BDhT0q8ouui2AN8u0x4haWCE+QRgOXBim9tyLi5v9bkPeBvwsZHqkUcCI6K+hi6lbZ85xKGvtUm7haKLDtu30v7uF2wvH209Ehgjoh4DmfMlIqLK4O56JjCBMSLqMWM6+DIREhgjor4p9OacTiQwRkR9CYwREVVT612LnUhgjIh6DGQyrIiIFmkxRkRUje0jgRMhgTEi6jE49zFGRLTIky8RES3SxxgRUWFnVDoiYh9pMUZEVBn3ddfEdQmMEVFPXjsWEdFGl92u08QsgS+R9AtJ90raJOnzbdIslfRU5ZXjn6lbbkRMDgbc746WqaKJFuMzwIm2d5ezdN0q6ce2b29J93PbpzVQXkRMJs6LavdRzgC4u9ycXi5T56chImrrtsEXuYFh9nJO6buBPwG+afuTLceXUsz32ksxy9fHbW8a4lwrgBXl5tHAg7Ur2JlZwOPjVNZ4yveaesbzu73a9mF1TiDpJxR17sTjtk+pU954aCQwPncyaSZwPfCfbW+s7H850F9ebi8DvmZ7QWMFN0DSXbYXTXQ9mpbvNfV083ebKhqdV9r2k8DPgFNa9u+0vbtcXwtMl9TpL0xExLhqYlT6sLKliKSXAm+nmBi7muZwSSrXF5fl7qhbdkTEWGhiVHoucFXZz3gA8H3bN0o6G8D2SuC9wEcl7QWeBs5wk9fwzVg10RUYI/leU083f7cpodE+xoiIbtBoH2NERDdIYIyIaLHfB0ZJp0h6UNJmSRdOdH2aImm1pO2SNo6ceuqQNF/SLZJ6ykdQz5voOjWhk0drY/zs132M5YDRr4CTKW4+vxM40/YvJ7RiDZD0ZxRPJH3H9hsmuj5NkTQXmGt7g6SDKR4sePdU/zcr79o4qPpoLXBem0drYxzs7y3GxcBm2w/b3gN8Dzh9guvUCNvrgScmuh5Ns73V9oZyfRfQA8yb2FrV50IerZ0k9vfAOA94tLLdSxf8R7a/kHQUcDxwxwRXpRGSpkm6B9gOrLPdFd9rKtrfA6Pa7Muv9BQgaQbF8/fn29450fVpgu0+2wuBI4HFkrqmC2Sq2d8DYy8wv7J9JMVLLmISK/vgrgWutn3dRNenaUM9WhvjZ38PjHcCCyS9RtKBwBnAmgmuUwyjHKS4AuixfclE16cpnTxaG+Nnvw6MtvcC5wI3UXTif3+o16FNNZKuAW4DjpbUK+kjE12nhpwALAdOrLwRftlEV6oBc4FbJN1H8YO9zvaNE1yn/dZ+fbtOREQ7+3WLMSKinQTGiIgWCYwRES0SGCMiWiQwRkS0SGCMiGiRwBgR0eL/A3DC/G5YZXOGAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "pd.set_option('display.float_format', lambda x: '%.3f' % x)  # вывод трёх знаков после запятой\n",
    "for i, field in enumerate(p):\n",
    "    df = pd.DataFrame(field / 1e6)\n",
    "    df.index += 1\n",
    "    df.columns += 1\n",
    "    print(f'Шаг по времени: {i} \\n'\n",
    "         'Распределение давления, МПа')\n",
    "    display(df)\n",
    "    plt.figure()\n",
    "    plt.title(f'Шаг по времени: {i}')\n",
    "    plt.imshow(field / 1e6)\n",
    "    plt.colorbar()\n",
    "    plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}