new?

Author: zingale

lectures/05-scipy/scipy-exercises-2.ipynb

@@ -67,10 +67,47 @@
   },
   {
    "cell_type": "code",
-   "execution_count": null,
+   "execution_count": 7,
    "metadata": {},
    "outputs": [],
-   "source": []
+   "source": [
+    "x = np.arange(13)\n",
+    "A = hilbert(13)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 9,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/plain": [
+       "(array([  1.20651400e-07,   9.99985193e-01,   2.00044484e+00,\n",
+       "          2.99430401e+00,   4.03835041e+00,   4.85081620e+00,\n",
+       "          6.34405283e+00,   6.55787325e+00,   8.20220729e+00,\n",
+       "          9.23611572e+00,   9.58340889e+00,   1.12467970e+01,\n",
+       "          1.19456441e+01]),\n",
+       " array([], dtype=float64),\n",
+       " 12,\n",
+       " array([  1.81383012e+00,   3.96833076e-01,   4.90294194e-02,\n",
+       "          4.34875507e-03,   2.95177714e-04,   1.56237036e-05,\n",
+       "          6.46641856e-07,   2.07632142e-08,   5.07655186e-10,\n",
+       "          9.14127833e-12,   1.14354232e-13,   8.87566335e-16,\n",
+       "          3.19134976e-18]))"
+      ]
+     },
+     "execution_count": 9,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "b = A @ x\n",
+    "import scipy.linalg\n",
+    "x = scipy.linalg.lstsq(A, b)\n",
+    "x"
+   ]
   },
   {
    "cell_type": "markdown",