{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Random variables from multivariate normal distribution" ] }, { "cell_type": "code", "execution_count": 19, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 22, "metadata": {}, "outputs": [], "source": [ "# example: 2d normal distribution \n", "mu = np.array([0., 0.]) # mean\n", "cov = np.array([[1, 0.8], [0.8 , 1]]) # cov. matrix\n", "\n", "# Cholesky decomposition, A = lower triangular matrix\n", "A = np.linalg.cholesky(cov)\n", "\n", "# number of random vectors\n", "n = 1000\n", "\n", "# random numbers following standard normal distr.\n", "# Ndim rows, N columns (here Ndim = 2)\n", "z = np.random.normal(size=(2, n)) \n", "\n", "# matrix with Ndim rows (here Ndim = 2): \n", "# [[mu1, mu1, ...], [mu2, mu2, ...]] \n", "mean = np.repeat(mu, n).reshape(2,n)\n", "\n", "# matrix of random vectors\n", "# vector i: (rand[1][i], rand[2][i]) \n", "rand = mean + np.dot(A, z)" ] }, { "cell_type": "code", "execution_count": 21, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQAAAAD8CAYAAACYVXqwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy86wFpkAAAACXBIWXMAAAsTAAALEwEAmpwYAAAsdElEQVR4nO2df3Bd1XXvv+v+kHBoXHeEE2KMhxB+BD9cbCKUKEyNCMSGAkGJ25kkbZUCQXYHaJwQ3LpM3vCGPjuBkCjBppGIzUNT+tLOuPwMlB8OCia+iSPHJiYQCKHBJcAE9HBNJ7GwdPf7Y3l173Pu+XV1z/2he9Zn5ox0zz3n3H1l77XXXj/JGANFUbJJrtkDUBSleagAUJQMowJAUTKMCgBFyTAqABQlw6gAUJQMk5oAIKI8Ee0hogfSeqaiKPUlTQ3gcwCeTfF5iqLUmVQEABEtBHARgG+n8TxFURpDIaXnDAFYB+CdYRcQ0SCAQQA4+uijP/D+978/pY9WFMXP7t273zDGzI+7rmYBQEQXA/iNMWY3EfWFXWeMGQEwAgDd3d1mfHy81o9WFCUEInopyXVpbAHOBvAxIvoVgO8A+AgR/WMKz1UUpc7ULACMMeuNMQuNMScA+CSA7xlj/rzmkSmKUnc0DkBRMkxaRkAAgDFmDMBYms9UFKV+qAagKBlGBYCiZBgVAIqSYVQAKEqGUQGgKBlGBYCiZBgVAIqSYVQAKEqGUQGgKBlGBYCiZBgVAIqSYVQAKEqGUQGgKBlGBYCiZBgVAIqSYVQAKEqGUQGgKBlGBYCiZJiaBQARHUVEu4joKSL6GRH9rzQGpihK/UmjJuAkgI8YY/6LiIoAniSih4wxP0zh2Yqi1JE0yoIbY8x/HXlZPHKYWp+rKPWiVAI2buSfrfi8RpJKVWAiygPYDeAkAJuNMT9K47mKkjalEnDeecDbbwMdHcD27UBvb+s8r9GkYgQ0xkwbY5YCWAigh4hO919DRINENE5E46+//noaH6soVTM2xpN1epp/jo211vMaTapeAGPMAQCPA7gg4L0RY0y3MaZ7/vzYnoWKUhf6+nilzuf5Z19faz2v0aTRHHQ+gMPGmANENAfARwF8peaRKUod6O1lNX1sjCdrrep62s9rNGnYAN4D4M4jdoAcgH8xxjyQwnMVpS709qY7UdN+XiOpWQAYY34KYFkKY1GUtqFUmh1aQaq9ARVFmV2eAQ0FVpSUqdYz0Mw4AtUAFCVlxDMgGkCUZ6DZ2oIKAEVJmWo8A0HaggoARZnlJPUMVKMt1AMVAErbMVss8EDz4whUACizjqgJ3uw99UxoZhyBCgBlVhE3wZu9p55tqBtQmVXEudjSis2fzSm+1aAagDKriDOapbGnjtMyZpONIQ4VAEpL459sSSZ4rXvqqG3EbLQxRKECQGkZ/JM9bLLV22gWpWW0m41BBYDSEriTPZ8HLr+czzdjskVpGc3226eNCgClJXBX1ulpYHgYKBaBwpH/oY2cbFF7/Gb77dNGBYDSEsjKeugQYAwf09PAlVcCixZVP9lmaqhLssefzfn/flQAKE0hzLg3OgrccQcwNWVX/6STWJ7Z1QWsXTszQ1277fFjMcY0/PjABz5glOyyc6cxc+YYk8/zz507K99fs8aYjo7wa6KeWSwak8uxHpHPG7NhQ3pjmy0AGDcJ5qIGAikNJy6Yp7eX1X6xByTJqffbEHK5ymCgkRFg5Ur+GYZoIjfeaDWHdg4K0i2A0nCSWNKrtbb7rx8aAiYm7PZhZARYvZqvfeQR/rlkSaXbUV6vX8/XtJvf308aVYGPBzAK4N3gjkAjxphv1PpcpX1JGsxTjbU97vpt27yvt2wB9u3zCowgu0G72wTS0ACmAFxrjPkJEb0TwG4ietQY80wKz1baFJlEotrXI6LPZdUqu/IDwIIFwO7ddmJv2xa8Ldm/n7cSQHv4/f2kURX4VQCvHvn9LSJ6FsBxAFQAKKFUq1oHRQnGRQ3u28cTe9UqYHCQnyOvlywBHn7YXr9qFbBjB78uFIB/+zfgS19id2SxyO7IgYH2Wv0BpOsFAHACgP0A5ga8NwhgHMD4okWL6moBVVqfDRvY0p7EUu+3zA8PV75escJr+e/vl2gCPhYv5uv8z92wwVr6xftQKHjvzeWq8yS0Ami0F4CIfg/ANgBrjTEHAwSNtgZT/pu+Pl5pifhnlGo9OsoBQkHq+uQkcPXVwKOPAuUyP6+jA3jlFe8znnmGjYDiAQgKFHK9Dy65XPup/kJa3YGL4Ml/lzHmX9N4ptL+GOP9GUSpxIFBck0+71XXiThoSN4nYoMeAOzaVfm8bdtY/fdvFwAbRFQs8nvyeZs3t6Hqf4Q0vAAEYAuAZ40xX6t9SEq7EbTajo3xSishv2HW9dFROxmJOElocNC68A4cAL76VSsAymVgzx7gH/6BXw8NAc8+a5+3alWlZX90FLjzTisQbr2VnwG06b7fIQ0N4GwAfwFgHxHtPXLu74wxD6bwbGWWE2bsi/LzuyG9/tVf3pdJed55lRrEHXfwxBVBsXy5N7R4/35vkhHgFQgTE1aAtDtpeAGeBEApjEVpM0ol4IYbeJ9eLnv96GF+e1dgEPF9AP9uDHD77bxay71vv83n5X2AJ7t8ztiYVzu46ip+nc+zZX/ZMl7t29nVF4VGAip1QSayTP5cjidXVxeH1bpJQC5i8DPGhvMS8e/T015B4moRhYLdTriT2L3GfYYgwT+FQhu7+iJQAaDUBVmdZfKffz7vv6Oy9EolYOtWu2IXi8A3v8kr9GuvAQ8+aCd4Vxd/hhvyK5/rjxWQa/xZgoBV/QH2AGRp8gMqAJQ64d/j33BDfFjt2Bir78Jll/HPLVv4HgnIWbYsXJDE1e5z4/8Br/EvS6q/oAJAqQthe/yoBJ+uLq96Pncu79lFKBw+zKv0xIQ3DuCGG/jwCxPZfkxOem0P7nVSgyCrqABQ6kJYoE1Uws7EBG8XZNuwd69XIOTzVmhI9aBymYOAduzwagIHDth7y2UWLkHj6+qyWoAYF4O6DbVLCTA/KgCUmgiaHDNNoe3rAzo7K+PzJyfZEHjRRRzfv2cPcNZZwJNP2oDdQ4esJgAAX/+6fS4RCxd3zDI+gAWEMcHbEk0HVpQQwiZH2F6/VALOPdda3S+6CDj2WGt5dzWEri6etENDPOG3bAHuvRe4557gsRgDPPYYC4zPfMYbzutqDoB3fC7+6/zXajqwojiETQ4xAE5Osiov6vfoKJ8DeD8vk/nb3waeeMK7P3cFy8qVfH0QRDYGQFyEAGsSk5M8qTdtqiztXSh4BYBEGfond7uVAfejAkAJJG7fWypVRtTJ5OjtBa65hkN0p6bYYr9kCbvygpiaYuHghgm7gsWf2CNIMM/cuazyi7HwtdcqKwK536mryxs9mMuxh0Gucb9vu5UBryBJymDahxYFbW3iCmMOD9vCmx0dnELrXrNzp031lXTaNWuM6ez0ptm6x5o14Z8/PBx879Kltnio/71isXJMQUVDczljenqqK0A6G4AWBVVmSlTRzlKJXXOHD7PKPTVVGUAzOlqpXgNeH79LPs92AMFfmHNwkAOC5DnC3r3cQETUfpfDh73uPTel2C0a2tkJnHlmdQVI2wndAigAvCp/0L5X3t+/P9w1J/hV/fnz+VyhYC3uQi4H3HZbcJUfV3X31/QTolKJ3e8WFGHoRhAmDQhqO5dgEjUh7UO3AK1FkMrvVstx3+/sZHU5l+PKOVJlx73eX41HjkLBmOXLvep3f3/lZ7hq+Lp19vqgI5fzbjfk6Oiwz3CrDxF5txvu38CtDpT079SqIOEWQDUAJVDlX7/ernAf/zjwu9/Z6/3tuvzuwJUrgz9nagr4wQ/493yeV//77+fafCtX2iSgQ4dYZd+3D7jpJu8zFi4EPv1pYN486yrs67P1/+bPB15/nWMIZPx+jcbdbghJCpC2o0tQBYAS6eoaGfH63okqM+bciXHoEJ8rFoNdd2IbILK/T04C993nrRC0dSuwdGnl/S+/DHzjG8Djj1da691KPzt28HkREGlY8tvRJagCIKO4LjEJuPG7zYDKvfeZZ3r3611dbBcQA50xwEMPcZDPvfcG79Hd3H3BtSsALBwWLAge++QkawZ33+097xdEf/VXfL6zkwWANPuYKW3pEkyyT0j7UBtAc5G9rLsXd/e07n54eNi7tx4eDr4/n+f9tVTlXbPG7pcLBf6Zy7ELbuFCey2RMaedVrmHl/EMD7ObLmifH1Tl1x3XbK7qWytQN6AShpurD3iLbMh+/ktf4p9LlrCrbcUK/jk4GHy/MWzll358AwO8Wl5yiY3SEy3hlVdswY+jjuJAoY4Ofr9YBPr7OZwX4M/70Y9YpT/mGO/32LLF27dPVuhTT/VeR9Qe6no9SKsq8FYAFwP4jTHm9DSeqdQPN1TXrdbT1xduEJTGGm4EoEzsXI7VbInbB2zSzv33272+VOMxR0p4dXfzPaJKy3ZD7AH+7LwDB7zfY3wc+KM/4mtFzQeA55/3XnfJJW2irteBtGwA/wfAJnCPQKWFCaqS49/7RxXrFCNbPs9CYdkyrz/9r//axvtHYQy35tq3jw837184dMha2kdHK9937QaS8y/PFvJ5YN26+PGE0XZ+fx+pCABjzBNEdEIaz1LqR5LUVn9Gnkwqf5YfwK7AwUF+7ugo8P3vx09+d3JOT3OzjiCjIOBV3cPyCAQ3IEkSgYh49Z8p7Z4KDCA9IyC4LdjTEe9ra7AGERbUUks7rrBgnZ07o2P8aznWrbNjKRaTXSvXS+5BWNBOksCfav5erQZaLRDIGDMCYAQAuru7EwRwKjMhatVK4sd2Q379zTMWLeIsv717baDNxz8evOrnjpiX/e49IHzFF3p6gCuuYAPkxo3c4ScsHViYN8+rri9axFuGoKCdpCt7O/r9/WgcQJsRFa0W58d2J4ZY9AH+fcsW24KLCPje9ziHP6xARy4HnHIKd+XxT/aoyQ9wrIEb1OO//j3vAV591b7O53m74o79wguDU5Xj/kYuben396ECoM2IW7WiQl7diWEMB+KceCKH1rrttYxhYXDXXeHjmJrihpxxzJsHzJnjndD+sQTd88YbrBUQAWefzcFHEko8Pc1BSLkcv7d4sff+alb2JCHCs5ok+4S4A8D/BfAqgMMAXgZwRdT1GghUX5Lsb8PuCwqkqffhb8e9bl1lAlKh4A0ekt+THP5AJ2NsS3F/MFG7gIQ2ADJx+lgd6O7uNuPj4w3/XKUSf0hwVxdw883ACy8EX5/PW/9/vcjlvC7Gri6OKfj+972aSBBuSLJLPs/1Bdavz4Z1n4h2G2O6467TLUDGcA1l+/YBV19t9/ZSGuuDH/QKAInQu/xyjvC75x7bkdcVCGmtJeUyRx0WCsBJJ/FY/C29gpCApGuuAW65JbwwaDtm9c0UFQBtQpKAFb+RT6zkgjTReOIJnkynnsp+9HnzvKm/t97K1+dyds+dy3mbedaKMbzHj1vxhUIB+OxnbabiwYMsRMRo6Rb8zIJ1PykqANqAJCqtv1NvnHuuXOaVt7+fX4+N2Zx7eYaLhAQ3krlzgbfesrt9tzTZwIC3yo+/5NjQEH8Xt25AFlEB0AbEqbRBnXqLRWvNz+WAL3yBV82RETu5p6fZ/3/nnZWTPklKb9q4NQGNsZPfzWUQolx4pZLtLbhjB7scsyoEVAC0OElU+ziV1s3eE9X+ne/k6/zq/ZYt3on8zDPWveYyZw7w29/O7DvJdiHIxedCBLz//ZU2ANE0ZPKff35lb0Ag3IWnNgCHJK6CtA91AyYjrAZdkJsvyvXnPsfvcnPdYBs2NMYFuHCh163nHuLiy+e99QZ7erzX5fMzr803m2r7zRQkdAOqAGhhgopZzvQ/r8THv+td3om0YoX3Gr+AaPQhdfr9/nl/odH+/pnFOrjftZb7W52kAkC3AC1MXx+7ryQy7447+PxM1dctWypj6letstl8ANsCbr45PZdetZTLwI9/zAZH2ZuXSsCDD9prikVO8fXv66sJ2W37CL+EqABoYXp72X0l7izJh5+JC+ummyonf38/T7Jzz7UJPR0dwHXXVScE/AbBd7wj2D4gwUZxGONt0HHDDfa7E3GiEMCJQhLPcNVVLDykMIhO7mSoAGhxxJ0l+e3LlvG5pKvdyAiv/D/+sfd8Ps8JM2vXerP5Dh9mw+B111WW5PZDxCW9PvEJb16AW0LcxT/5/YLj2GOBN9/kyd7RYRN8XO9FZyf/DdzCJNLVB7CFQVQAJENrArY44rPO53kSrF3L5926/WGMjHDBjV27vBONCLj2Wo6Y27XLe49k1u3dGz+2RYtYg1i+nFVy1zofx/LlXM7L5UMf4nLf0hJsYsJ6L6SEmHtetkF+b4J0I1biUQ1gFjAxYYN3JDc/iQYQ1E6rWOR22RMTwTn25bJVpwGeePk88N73Ar/4hffal17iY9cu3k74e/dFcdRRnGn4xBPe8/69uWsDeeopPhfW3ltMhNKNWLWAeFQDmAWInz+f5//4W7faqr2lUvh9q1Z5X/f3c0LN4CA/U1piu0jDT1G5P/pRNgy++GL0GF95pbpIwHe8g1V5dwwPPeT9PmIDcZuLinp/2WX2fC4HvO99NjQ5aw0+a0EFQAsjJa8B2y33ssvsnndykg1kfiEg9/3yl3ZSFgq2OKY8U7YTcdx8c3zQzptvxlftcbnnHk5EuugiO5EPH66cuAMDtmR4oWCNngMDrEWIUDz9dG9Z8izH91dFEl9h2ofGAcQTFQQU19QjrDnGmjXcNFMCad7znur99NXk4QcF+PjPL13qfe1vNjo8zGMm8jb8lGvc2n+dnTZWIutAG4PMboLCVQEb437++Tb7zi2J7W/aIRBxZd233+bX09OVVXiSMNP4AGOCtwgvvuhV5ScmvM1J3HTl6WmvhtDb6639NzXlTQhS4lEB0KK4+/6gRJdVq+xkLJet5Vvu8xvkaimPnRaXXGKzC4WDB61w6Ozk7yFZi7LVyeXCVfuov5MSjwqAFkVWenGJ+Ve1iQm7ouZyXDHHtResXs0GNiKeWBdeCHz3u95ndHZyBd5GIDaIu+/mn66AImKNZmiI7RKPPebNWrz4Ym5JHtXHIOzvpMSQZJ8QdwC4AMBzAF4A8Ldx16sNoJJqY9NdG0FHR3ANfPeZa9ZU7snXras8n/Zx2ml2X75zJ8fwu81BxR4hY5XcB8kJ6OiYWd5DtX/LdssLQKOSgQDkAfwSwIkAOgA8BWBx1D0qALzUkuCzYQNPMH8DC9eI1t9vzPz5lZOzWAw36s3U2Oc/1qyxYw1r7lEsepuMEPHP/n5vx+EkjTmq/Vu2a2ZgUgGQRiBQD4AXjDEvAgARfQfApQASFIVWgHCDX1ywjwTNlEre6jddXcA550S75Yiscc3PH/wBcNZZwCOP1Pa9ikVbiWdsLHw8Ytzr67PjmZ7mLYu8dl2AUVSb65/12gBpCIDjAPyH8/plAB/0X0REg+D2YFi0aFEKH9v6JM1Q8xf0OHCAJ/D0dLLkFn/1m6DEH2HePK6kE1XE8803gUcfTfAFjyDRen6uuMKOOyo8N5eznYkl6s8f5XfZZckmZrX1/rJeH1Bbg9WJakpP+xtyup1yJydt6K+/k6+k8UrjzGOP5Z+vvBI+LrfFNhEwfz7wm99UXhcmHFzmzweOPhr41a+C35871/4+MRHeEuwLX6jsTJzPWy3FX9Mvimq7+WSh+08kSfYJUQeAXgAPO6/XA1gfdU8WbAAzbSzpr8ojAS7uuVyODXhBTTnzeWP+7M+852YS8BN3FAo2qCjqGtcgGXR9Luf927gGuXY0zjUKNNAIWADwIoD3whoB/0fUPVkQADMxRonRTiL5CgU2hIkgSWqky+dZQEjnm+Hhyuvnzq1NAMybF28o9E/u4WHvdwnq2KOkQ1IBUPMWwBgzRURXA3gY7BHYaoz5Wa3Pne34VUvAFrCIasrZ0eHtwAuwMcy/xzaG985BlXinp4Hnnwcefpiffe65lar3wYO1fT93K1EscgMPfw3/zk7vntqtByDpvUNDGVS7W4kkUiLtIwsagEucNuD3fxcKVvXv6OBz/tj+QoFX+Sg33vBw/f38gDHLl/NnuecWL+bxuSq8uPlcDWD5co3frwfQXIDWIczNJ7jhrFL4Q7L93IIY/f3ACSewJf/DH+YU2DPOCP5MY9iYKAbCpFST0y/84AcciejG+j/zDHsjrr/epi1LGq9QLnM9gG99i7WUqNRmpT6oAGgAcfHqbjjrpk22r70gjS9OOYUt7gcO8MRZvTq6cs/0NHsGOjqSjZMIOO20xF/rv5HtRWdn8HtustLAQOX3AzSHv1moAGgASeLVe3u5zNfgoLcIhjS+2L49WZkuF1mR165l4ZPk+qBefETRmkFnJ7v8liyxrkgXN1mptxfYvLlyPEGCUeoaqGZQP7QkWIMIK0MdFCy0bJldVctlrp4DsFHQH51HxBO3qwv4/d/n9Fpj7Pnbb7eVcuIIK/oRdS8RayZRBUSJrAGwVOLfb7uNtw2vvcZCQ5p6Cllo4d0KqABoIqUST/zDh9mSLmGokuknFv5772WL/vbtbAe45x4+n8uxJX3PHuD1121ATy4HnH02782lpVa9GncaY2v1uUjQj1QO7utL3sR0bAzYv99rN0laB1GpDhUAKVNNg4rRUVugQ/6T9/baopfynruPXreOhcHkJKvRRx1VGfZbLgNPPmnVdlnBxcCYRBuohTPOAN79bmDpUm/vwY0bkzUxlfblsk3I57kpikQFqjaQHmoDSBG3kk1Qwc6ke1ophumSz9uJdM01fG5qilf5IMST4L5Oa/KLR6K/P7iw6NNPcy7B0JBXEMYZQ11vydQU/w1uvJF/StUfNRamiwqAFBkd5U66Qf9Rg4TDwAAb0IjsRBLhMDDAHXhlL//5z9v4/1tusZNZKubEGerSIpfjpiE9Pbz3DwtEMsZqNYL0ODjlFC4zvm+f9z6/gBgYYMOoFAbVqj91IEmwQNpHOwYC+YNc/AUsw3IDpFiHBPzk8xz+u3MnB9cUi96Q2bDAHqL6d/adN4+De6RIZ5J73HoAbt0COfxNQMPi/zUvoDqg3YEbS1AnX5eoaMCgttxS4dYVGmvWVLbJTnI0ouU3wEVHgir9yHcPEhpud2IlPZIKADUCpoQ/r9yfvirq77Zt7M5zjVh9fZVx/WIAdNNj77ijutr7RNzY49Chyg489eCNN3haSxHPz3/ea9E3AfYHf/OSZlNtl+FZTxIpkfbRjhqAMdHqq1u/XjQA93p/plxnp/eaIPU57pB03P7++qz4QZpFLmezEINqFnZ2cvx/T0+l+t9s2qk8GFQDaDxBwT5i/Dt0yK6Ak5McnffUU17X1o4d1mjmBsaI8W/LlvgOPS6f/Swb2u67b+bfqVgM1jqCgoukAvENN3gt+gBX9V20qLVX1iyWB1MBkCJB6qP8p3InS7nM7brlnPxnC+r4K8/s6gq2uIeRy3EfwJGR6u5zIeLWXfffXyl4gp556aUcp+Cv7iNbolafTFksD6YCICVkpZ+c5Mm3eTPH9Xd1WVdeocDBMePj3u67Yf/Z3MCYXC756i/2hKC4fj9hZboAXv2DYvvD7uvp8Wot9Sy1VY+9ehbLg6kASImxMZ78EnBz9dV8fu1afp3PA7feygkzMqnzeQ5yCYqD94fDxq3iImQuuYRrAu7alWzcxxzDYcSCqPa5HI9Xnu1SKHAdv699zdYuzOcrC3+G5T/USj3zBOo15lZFBUBK+C3509Ns8Xfz+ScmKguAulVyAO9/bpl4Uas0AJx8MvDv/86f893vAscdl3zc7uQHrDkP4HbdDzxgJ7l06N20ibWb/n5OArr/fr5n7VoWcPWeQFncq9cLjQSskrBwXklzLRZtn7tVq8Ij2Pbv55Bef9iwaBISDitRdVH84hdWSzh8OLxKbxwywaXl9n332ckP8DjKZSu0entta7FymQ2dbuRfvdB+gCmSxFUQdgD4UwA/A1AG0J30vtnqBkziJvK7AoNe+4Ni3MhAf2mtoCMooCaNTj5z5vDnb9jArrogF5//e/ur/Yr7st5oZGA0aJAb8GkAnwAwXONzZgVJVE//HtLt3rNxY2VQjKy6+/fzNQ89FD8OIs79dwtznnEGv45b/QsFHn+hYBOGiKwFH+BV/MknvfctXw5ccEGlcUwSl4aH+TtNTTVGJc/aXr1uJJEScQeAMWRYA4hbjYIaeUrc/9KlNqCmUAgOrjnmmPjVu1gMXrX9R08PBxVJnoF7vwTv+LUJCemt9u+Sxt9bV/mZgUbmAiQRAOC2YOMAxhctWlT3P0C9CFPpo/7z+xOB+vvDG2UGHf5GH2FHkkhBUePXrKmc6D09lc/I570Re3HRjmlV+G2nqLxmkFQAxG4BiOgxAEHe4OuNMfdWoWm0RWswv+qZZFvgDzA59thwn76/aEcux0a+5cuBl15idf3ll6270cV9pqQIyznxKJTLNs/AH1uwYAFHDoo7E+CtwpIl/Huc+00alN55Z+2uObX0N4ZYL4Ax5nxjzOkBR+LJ38709dk+dlK0w4+/KOjAQHCJrmIRuPZa73vlMvv0n3iCBcD+/fw5H/tYdJkvN9Aon+cc/qOOsve89hrwyU967znlFOAzn+EyY3Kd7OmB6PLmcaXPq0Ut/Y1B4wBSwPXXh+HXHDZv5mAhMcJdcgkb4UZHo4N+ZAX3NwAN6hIkmoQxPCHPOouNe9PTXFdQtATRNm65hX8WCtZY6E6+qFDZtMNosxiV1wzImJlr40T0cQC3ApgP4ACAvcaYlXH3dXd3m/Hx8Rl/biuxcSP78qenebW68UaO6U+CdPcFbPrwOefY5JtisbK0lyC+ern2U58C/vmfK4N24qIIgwSHeAV6euxEdluchU3KzKXStjBEtNsY0x17YRJDQdrHbPUCBFGLscp/rz9td/lyW2WIyJiTT/a+v3Sp9RzMmcPGQjHsdXSw8W7Fimij4AknBMcQSEWjNI1xatVvHNB04MYQp6pGrYo33QT87nf8e5Ba/8YbbJADeFr+8pfe9596ildricL7p3+y6rwx1njn7yXg8utfszYhWxERAdPTdh+fRnlurfPfmqgASIGoph9h/+lHRmx9f4AnX18fsHu3DdR57jnv8/y7NWNsEpBfjZ+a4rz8E0+MziVwJ75sG8pl7z7eX5VoJuW51arfmmguQB3x/6cfHbV5BNu2ea89+WTOvjOG9/4f/nBwS/AgursrDZDGAI89BmzdWtmLb+lSzlVwm5EC1rV45ZV2crsejFrKc6tVvzVRDaAOuEU8ZPWUNl3G8OS75hqvan7qqcDzz9vMwUOHvM+cOxd4661gIXDmmdZ/L626nnvOGhCvvNJeK6nH7hivusoaD43hyj1h4czi6692EqtVvzVRAZAyfrV/aIhbd91+u13RJyeBgwc5nfaVV4ArruD9+sMP2644CxbwT5mY7uSXjjkiTAYG+HAt9dJyLJ+vnPRA5bZFXJKdneETu9ZJrPH7rYcKgCpI4ubyq/0TE7yiuis3Edf3m5riSb5nDwuA7dt5m7B1K6fiuvt6//2bN/Oz3bHIz1LJG5uwb599rvj23f374CB/fpKJrZO4zUjiKkj7mI1uwKTuMLkul2MX3fBw5Tl/0g6RfaabNxB1+PsOCMPDxpx0knXtyWe6rj6p3KvuuPYFCd2AagSMwC3+kTTUVer/i3Ft7Vo+v3078Pd/zyG9ixd77zHGPlOMZRKKK5Z5/z1BjIwAq1cDL7zAz5T6/P6egOUyGwiD+hcq2UIFQAj+Xn5i0EtixZ6YsBF40tW3t9dW/ZVed0IuZ58pAqS729sd9+KLvZ+xbJl3rBs38rbC5cQTuXxXsWj7D/b0WLdhkCBL2sBUaQ/UBhBC0F4+zgAmob3PPON1re3axe+5e/WxMW9dQHlmqcRag9tHYGoK2LvXTtxczpbl8lcOdrnuOt7by3OI2OC4b1+wJV+DdbKHCoAQ/MktXV3xk7+vz6bautx7L1f6uewya5H3G9OiKgZ1dHB9wR07vBO3VOJgHzd9t78f+O1v+frBQX6m1BWcno4WZBqskz1UACDYuu+v3rt2bfTKODYW3rfPGJ6kw8PBufIjI9YNVyxa1b9Q8AoN11IP2D4EohV0dHgbcwDJs/RKJZtqDGiwTmZIYilM+2glL8BMKvpIAU//c9zimMUiJ/d0doYXAJX7CgWvR6C/Pz5pxh1TnFXfTcIJ+r7uOelKrB6C2Q00GSgZM6noI+q3qzXIvt7f229khI1ze/ZUxtjL5/v9/fffD1x4YfhWwx9l2NHBWwE32MevzcjvGzcGezPcPn7+SECljUkiJdI+ZpsGINdFraJxz+7oCF5Zd+60Kb/uUSwGX+t+rpTwrqY+YZwGoPX32gOoBpCMmYS3JjWWudcBlSvryAgnBX3wg2zgc331ko4bZaSbmPAWH0latjzo+2qcfjbJvAAAvJ183ddCUHx/R4dNvvH3xBOiDHB/8zdcD0AoFtndJwE8QTH5cQa9pAa/oHBeDfHNJioAEO//Dlp5h4as5d7fE8/dh7srK8B78K4u4Ktf9Y7hpJOA44/nVN1586zbEUjecVcz7pSqSbJPCDsA3Azg5wB+CuBuAPOS3NdKNgBj4q38QXvksHuimofI+aDmH8Uiny8WjVm3zuYOSMOOOOpVbkvLeM1O0IjGIABWACgc+f0rAL6S5L5WEwAz7fkXdE+YYAhL8iHi5CBXKORyXtdhoVDZjy9oLNJtqL9fm3NknaQCoKYtgDHGrTb3QwB/UsvzmkUS1Tmo51/QPbIPn5zkvbzYB+S81AAUzjoL+PKXuRqwuAON8ZbxKpetQS9ouzI25g0dvucejjx8/HFtzqFEk2Yy0OUAQltbEtEgEY0T0fjr/qb0LYCbrJOEsNoAkswjXXfWrrV5ANu3c6iuyxVX8HubNnHkXy7HDTy++EX72jUIBk3Krq7KSkHanENJQiqtwYjoegBTAO4Ke45pk9ZgQLTRUOr9Sfy9mw3Y2wvcfbd1/0m8PhBclKO/P1zDcC39Y2OVhT9rnbAi4IaGKguPKG1Ekn1C1AHgLwGUALwj6T2tZgOoFrexZi5XaQD019lPYsSrhqQ2gJka8HTvP/tBI2wARHQBgHUAzjHG/LZWYTQbKJU4tNfdn8s+X9RzdyV2U3fTIok9opbUXt37Z4da4wA2AegE8ChxEbofGmPW1DyqFmZszBbqBFj1npiw2XSFAgsAydCLKrIJVLYHi5toUbaHtFJ70+7zp7QutXoBTkprILMF18oPcARfV5ddbfN53s8vWxa/dy6VuO23CJQ77oi23Luruj9VOGycWsJbiUIjAauktxf45jd5G7BgAeffB8X8i3EvitFRrzbhGgyDGBuz+f/T0+H1BWScWsJbiUMFQJWMjHAjjXKZ1ft169JTmfP56Hu7uipTh6PUe53EShxaFLQKSiWO/5+aqiz4Ke2zqjG2DQywECHiyX/bbdH3Tkx46/5JuTDdoyszRTWAKhgb8/bry+Xs5JvJatvby3v+pGp6Xx8LjCQ2AEVJggqAKpAJODnJK/amTbVPvmoEhxrnlLQh448hbQDd3d1mfHy84Z+bBq4bDtDJqLQmRLTbGNMdd51qACHE+dv9gTbXXMO1+1etSt5nT1GajQqAAJJE0bmuv8lJW93nkUc4NkAKgGpzDaWVUS9AAEFRdH7cTDk/hw/H9xBUlFZANYAAkvj1XYPcgQOV9f2CSoArSquhAiCApNZ214L/vvfZFF+1ASizBfUCKEobktQLoDYAaEtsJbtkfgugLbGVLJN5DcBv8R8dVW1AyQ6Z1wBci38+zzn5U1OqDSjZILMawMgIsHIlsG+fzeS7/HKe/OrDV7JCJgXAyAiwejVH7a1ezUJg/XrOrNMy2EqWqEkAENGNRPRTItpLRI8Q0YK0BlZPtm0Lfj3TvH5Fma3UqgHcbIz5Q2PMUgAPAPiftQ+p/qxaFf662gYhijKbqbUo6EHn5dEAZkXDD6nX52/OoShZo2YvABH9bwADAP4TwLk1j6hBDA7qxFeU2C0AET1GRE8HHJcCgDHmemPM8eC2YFdHPKelewMqShZJLReAiBYBeNAYc3rctZoLoCj1pSG5AER0svPyUgA/r+V5iqI0llptAF8molMBlAG8BKCt24IpSrtRqxdgVfxViqK0KpmMBFQUhVEBoCgZRgWAomQYFQCKkmFUAChKhlEBoCgZRgWAomQYFQCKkmFUAChKhlEBoCgZRgWAomQYFQCKkmFUAChKhlEBoCgZRgWAomQYFQCKkmFUAChKhlEBoCgZRgWAomSYVAQAEV1LRIaIjknjeYqiNIaaBQARHQ9gBYD9tQ9HUZRGkoYG8HUA6zBL+gIqimKpqSz4kfZgvzbGPEVEcdcOApBufJNE9HQtn92iHAPgjWYPok6063dr1+91apKLYluDEdFjAI4NeOt6AH8HYIUx5j+J6FcAuo0xsX9MIhpP0rZottGu3wto3++W9e8VqwEYY84P+YAlAN4LQFb/hQB+QkQ9xpjXqhyvoihNYMZbAGPMPgDvktfVaACKorQGzYoDGGnS59abdv1eQPt+t0x/r9TagyuKMvvQSEBFyTAqABQlwzRdALRbGDER3UxEPyeinxLR3UQ0r9ljqgUiuoCIniOiF4job5s9nrQgouOJ6HEieoaIfkZEn2v2mNKEiPJEtIeIHoi6rqkCoE3DiB8FcLox5g8BPA9gfZPHM2OIKA9gM4ALASwG8CkiWtzcUaXGFIBrjTGLAXwIwFVt9N0A4HMAno27qNkaQNuFERtjHjHGTB15+UNwfMRspQfAC8aYF40xbwP4DoBLmzymVDDGvGqM+cmR398CT5bjmjuqdCCihQAuAvDtuGubJgDcMOJmjaEBXA7goWYPogaOA/AfzuuX0SaTxIWITgCwDMCPmjyUtBgCL6zluAtrygWII0kYcT0/v15EfS9jzL1HrrkerGbe1cixKdVBRL8HYBuAtcaYg80eT60Q0cUAfmOM2U1EfXHX11UAtGsYcdj3EojoLwFcDOA8M7sDLX4N4Hjn9cIj59oCIiqCJ/9dxph/bfZ4UuJsAB8joj8GcBSAuUT0j8aYPw+6uCUCgdopjJiILgDwNQDnGGNeb/Z4aoGICmBD5nngif9jAJ82xvysqQNLAeKV504A/88Ys7bJw6kLRzSALxpjLg67ptlGwHZkE4B3AniUiPYS0beaPaCZcsSYeTWAh8FGsn9ph8l/hLMB/AWAjxz5d9p7ZNXMFC2hASiK0hxUA1CUDKMCQFEyjAoARckwKgAUJcOoAFCUDKMCQFEyjAoARckw/x9ry4canl1irwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots()\n", "ax.set_aspect(1)\n", "ax.set_xlim(-4., 4.)\n", "ax.set_ylim(-4. , 4.)\n", "plt.plot(rand[0], rand[1], 'b.');" ] }, { "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 }