From 2d7c25e5aab25959f5304131c5fc95fd7d5abe7c Mon Sep 17 00:00:00 2001
From: Satoki TSUJINO <tsujino.satoki@gmail.com>
Date: Tue, 2 Apr 2024 17:56:43 +0900
Subject: [PATCH] temporary commit.

---
 Lilly1966/OSE_solve.ipynb | 242 ++++++++++++++++++++++++++++++++++----
 1 file changed, 217 insertions(+), 25 deletions(-)

diff --git a/Lilly1966/OSE_solve.ipynb b/Lilly1966/OSE_solve.ipynb
index 36b858d..9adcb3e 100644
--- a/Lilly1966/OSE_solve.ipynb
+++ b/Lilly1966/OSE_solve.ipynb
@@ -74,7 +74,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": 1,
    "id": "7e7d6392",
    "metadata": {},
    "outputs": [
@@ -619,26 +619,218 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": 10,
    "id": "d6476cfa",
    "metadata": {},
    "outputs": [
     {
      "data": {
-      "text/plain": [
-       "Figure(PyObject <Figure size 700x500 with 1 Axes>)"
+      "image/svg+xml": [
+       "<?xml version=\"1.0\" encoding=\"utf-8\"?>\n",
+       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"600\" height=\"400\" viewBox=\"0 0 2400 1600\">\n",
+       "<defs>\n",
+       "  <clipPath id=\"clip6900\">\n",
+       "    <rect x=\"0\" y=\"0\" width=\"2400\" height=\"1600\"/>\n",
+       "  </clipPath>\n",
+       "</defs>\n",
+       "<path clip-path=\"url(#clip6900)\" d=\"\n",
+       "M0 1600 L2400 1600 L2400 0 L0 0  Z\n",
+       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
+       "<defs>\n",
+       "  <clipPath id=\"clip6901\">\n",
+       "    <rect x=\"480\" y=\"0\" width=\"1681\" height=\"1600\"/>\n",
+       "  </clipPath>\n",
+       "</defs>\n",
+       "<path clip-path=\"url(#clip6900)\" d=\"\n",
+       "M208.77 1487.47 L2352.76 1487.47 L2352.76 47.2441 L208.77 47.2441  Z\n",
+       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
+       "<defs>\n",
+       "  <clipPath id=\"clip6902\">\n",
+       "    <rect x=\"208\" y=\"47\" width=\"2145\" height=\"1441\"/>\n",
+       "  </clipPath>\n",
+       "</defs>\n",
+       "<polyline clip-path=\"url(#clip6902)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
+       "  208.77,1487.47 208.77,47.2441 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6902)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
+       "  637.568,1487.47 637.568,47.2441 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6902)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
+       "  1066.36,1487.47 1066.36,47.2441 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6902)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
+       "  1495.16,1487.47 1495.16,47.2441 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6902)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
+       "  1923.96,1487.47 1923.96,47.2441 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6902)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
+       "  2352.76,1487.47 2352.76,47.2441 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6902)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
+       "  208.77,1487.47 2352.76,1487.47 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6902)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
+       "  208.77,1160.15 2352.76,1160.15 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6902)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
+       "  208.77,832.824 2352.76,832.824 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6902)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
+       "  208.77,505.499 2352.76,505.499 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6902)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
+       "  208.77,178.174 2352.76,178.174 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6900)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
+       "  208.77,1487.47 2352.76,1487.47 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6900)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
+       "  208.77,1487.47 208.77,47.2441 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6900)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
+       "  208.77,1487.47 208.77,1470.19 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6900)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
+       "  637.568,1487.47 637.568,1470.19 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6900)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
+       "  1066.36,1487.47 1066.36,1470.19 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6900)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
+       "  1495.16,1487.47 1495.16,1470.19 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6900)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
+       "  1923.96,1487.47 1923.96,1470.19 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6900)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
+       "  2352.76,1487.47 2352.76,1470.19 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6900)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
+       "  208.77,1487.47 234.498,1487.47 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6900)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
+       "  208.77,1160.15 234.498,1160.15 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6900)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
+       "  208.77,832.824 234.498,832.824 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6900)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
+       "  208.77,505.499 234.498,505.499 \n",
+       "  \"/>\n",
+       "<polyline clip-path=\"url(#clip6900)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
+       "  208.77,178.174 234.498,178.174 \n",
+       "  \"/>\n",
+       "<g clip-path=\"url(#clip6900)\">\n",
+       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 208.77, 1541.47)\" x=\"208.77\" y=\"1541.47\">0.0</text>\n",
+       "</g>\n",
+       "<g clip-path=\"url(#clip6900)\">\n",
+       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 637.568, 1541.47)\" x=\"637.568\" y=\"1541.47\">0.2</text>\n",
+       "</g>\n",
+       "<g clip-path=\"url(#clip6900)\">\n",
+       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 1066.36, 1541.47)\" x=\"1066.36\" y=\"1541.47\">0.4</text>\n",
+       "</g>\n",
+       "<g clip-path=\"url(#clip6900)\">\n",
+       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 1495.16, 1541.47)\" x=\"1495.16\" y=\"1541.47\">0.6</text>\n",
+       "</g>\n",
+       "<g clip-path=\"url(#clip6900)\">\n",
+       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 1923.96, 1541.47)\" x=\"1923.96\" y=\"1541.47\">0.8</text>\n",
+       "</g>\n",
+       "<g clip-path=\"url(#clip6900)\">\n",
+       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 2352.76, 1541.47)\" x=\"2352.76\" y=\"1541.47\">1.0</text>\n",
+       "</g>\n",
+       "<g clip-path=\"url(#clip6900)\">\n",
+       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 184.77, 1504.97)\" x=\"184.77\" y=\"1504.97\">-1.00</text>\n",
+       "</g>\n",
+       "<g clip-path=\"url(#clip6900)\">\n",
+       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 184.77, 1177.65)\" x=\"184.77\" y=\"1177.65\">-0.75</text>\n",
+       "</g>\n",
+       "<g clip-path=\"url(#clip6900)\">\n",
+       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 184.77, 850.324)\" x=\"184.77\" y=\"850.324\">-0.50</text>\n",
+       "</g>\n",
+       "<g clip-path=\"url(#clip6900)\">\n",
+       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 184.77, 522.999)\" x=\"184.77\" y=\"522.999\">-0.25</text>\n",
+       "</g>\n",
+       "<g clip-path=\"url(#clip6900)\">\n",
+       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 184.77, 195.674)\" x=\"184.77\" y=\"195.674\">0.00</text>\n",
+       "</g>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"718.024\" cy=\"173.278\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"616.255\" cy=\"417.543\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"2276.95\" cy=\"224.244\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"957.25\" cy=\"341.185\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1641.93\" cy=\"1473.49\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"2216.3\" cy=\"260.989\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"998.826\" cy=\"490.868\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1642.41\" cy=\"1390.35\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"2155.63\" cy=\"297.73\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1642.96\" cy=\"1309.65\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1226.95\" cy=\"451.466\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"2094.95\" cy=\"334.466\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1643.59\" cy=\"1231.36\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"2034.27\" cy=\"371.199\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1308.46\" cy=\"553.452\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1973.56\" cy=\"407.928\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1912.85\" cy=\"444.655\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1467.75\" cy=\"527.969\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1852.12\" cy=\"481.389\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1644.32\" cy=\"1155.45\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1645.17\" cy=\"1081.92\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1791.28\" cy=\"518.189\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1646.16\" cy=\"1010.73\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1728.59\" cy=\"554.801\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1572.57\" cy=\"603.955\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1672.82\" cy=\"580.93\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1647.31\" cy=\"941.85\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1648.65\" cy=\"875.254\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1661.63\" cy=\"628.224\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1650.22\" cy=\"810.906\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1654.91\" cy=\"688.574\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1652.12\" cy=\"748.761\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"803.092\" cy=\"244.816\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"664.851\" cy=\"439.197\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"2276.93\" cy=\"224.219\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1101.42\" cy=\"359.154\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1642.37\" cy=\"1431.59\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"2216.22\" cy=\"260.924\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1032.03\" cy=\"525.278\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"2155.47\" cy=\"297.612\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1642.88\" cy=\"1349.68\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"2094.68\" cy=\"334.286\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1643.47\" cy=\"1270.2\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1349.47\" cy=\"448.499\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"2033.83\" cy=\"370.948\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1644.14\" cy=\"1193.11\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1972.92\" cy=\"407.601\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1328.42\" cy=\"589.683\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1911.93\" cy=\"444.251\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1644.92\" cy=\"1118.41\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1850.8\" cy=\"480.866\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1573.89\" cy=\"518.443\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1789.78\" cy=\"516.913\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1645.83\" cy=\"1046.06\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1735.97\" cy=\"551.999\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1646.88\" cy=\"976.046\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1648.11\" cy=\"908.325\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1697.03\" cy=\"596.494\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1586.1\" cy=\"639.762\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1649.54\" cy=\"842.869\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1663.3\" cy=\"655.044\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1651.29\" cy=\"779.627\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<circle clip-path=\"url(#clip6902)\" cx=\"1654\" cy=\"718.226\" r=\"14\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<path clip-path=\"url(#clip6900)\" d=\"\n",
+       "M1989.93 251.724 L2280.76 251.724 L2280.76 130.764 L1989.93 130.764  Z\n",
+       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
+       "<polyline clip-path=\"url(#clip6900)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
+       "  1989.93,251.724 2280.76,251.724 2280.76,130.764 1989.93,130.764 1989.93,251.724 \n",
+       "  \"/>\n",
+       "<circle clip-path=\"url(#clip6900)\" cx=\"2097.93\" cy=\"191.244\" r=\"21\" fill=\"#009af9\" fill-rule=\"evenodd\" fill-opacity=\"1\" stroke=\"#000000\" stroke-opacity=\"1\" stroke-width=\"3.2\"/>\n",
+       "<g clip-path=\"url(#clip6900)\">\n",
+       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:start;\" transform=\"rotate(0, 2181.93, 208.744)\" x=\"2181.93\" y=\"208.744\">y1</text>\n",
+       "</g>\n",
+       "</svg>\n"
       ]
      },
-     "metadata": {},
-     "output_type": "display_data"
-    },
-    {
-     "data": {
-      "text/plain": [
-       "Figure(PyObject <Figure size 700x500 with 1 Axes>)"
-      ]
-     },
-     "execution_count": 4,
+     "execution_count": 10,
      "metadata": {},
      "output_type": "execute_result"
     }
@@ -649,10 +841,10 @@
     "#  Plot  #\n",
     "##########\n",
     "\n",
-    "using PyPlot\n",
+    "using Plots\n",
     "\n",
-    "rc(\"font\", family=\"IPAPGothic\")\n",
-    "fig = figure(\"pyplot_majorminor\",figsize=(7,5))\n",
+    "#rc(\"font\", family=\"IPAPGothic\")\n",
+    "#fig = figure(\"pyplot_majorminor\",figsize=(7,5))\n",
     "\n",
     "x = real.(Dcal[1:Ntot])\n",
     "y = imag.(Dcal[1:Ntot])\n",
@@ -660,21 +852,21 @@
     "#ax = gca()\n",
     "\n",
     "#cp = scatter(x,y)\n",
-    "scatter(x,y)\n",
-    "xlim(0.0,1.0)\n",
-    "ylim(-1.0,0.1)\n",
+    "scatter(x,y,xlim=(0.0,1.0),ylim=(-1.0,0.1))\n",
+    "#xlim(0.0,1.0)\n",
+    "#ylim(-1.0,0.1)\n",
     "\n",
     "#ax = gca()\n",
     "\n",
-    "xlabel(\"Real(λ)\")\n",
-    "ylabel(\"Imag(λ)\")\n",
+    "#xlabel(\"Real(λ)\")\n",
+    "#ylabel(\"Imag(λ)\")\n",
     "#if i == 1 #div(ne,2)-1\n",
     "#   plt.colorbar(cp)\n",
     "#end\n",
-    "grid(\"on\")\n",
+    "#grid(\"on\")\n",
     "\n",
-    "fig.canvas.draw()\n",
-    "gcf() # Needed for IJulia to plot inline\n",
+    "#fig.canvas.draw()\n",
+    "#gcf() # Needed for IJulia to plot inline\n",
     "\n",
     "#ttext = \"OSE spectra: Re = $(Re), α = $(α)\"\n",
     "#title(ttext)\n",