Que pasa si el otro input.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Esta parte la puede ignorar es una manera rapida de hacer los calculos con matrices, es un pedazo del programa que uso para hacer las graficas, de la parte de N interferometros"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "from sympy import *"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "init_printing(use_unicode=True)\n",
    "theta1, theta2, theta3,theta4,gamma1,gamma2,gamma3,gamma4= symbols('theta1 theta2 theta3 thetha4 gamma1 gamma2 gamma3 gamma4 ',real=True)\n",
    "beta=Symbol('beta')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAALYAAAAzCAYAAAA6jfJvAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAJiklEQVR4Ae2c27HUOBCGD9QGwEIEe8hggQiADLhEAGQAxRO8UZABEAELGUAGLGTAZsAlA/b/NGof2SN5ZI/tI824qzySdWm1ulutlqzRhWfPnl06icDz589/RZLXpJUDRXFAehrV34ui8ouen53nZVHUr8SsHEhz4L2yuvr7/kQW+5ueW79//z6p7RHdp3PQPBfeOWgtCedcfBuCV2Uf6vmIxa4SNAU9FuF/z0T8qcc/E/rDQ1uaPC5Iu7+JzY9E2Kda2C1a74jWGwqfpGj2jL7i89/p/WuqbCzd1/+l8E0sf8k0T8tVhY/GtKt6p6qHy/lkjv4IZzHyEC0P1c+71VlsEc5i4anCPqX+qDIoJWVe6MEPGwSq+0oVGPDRxckgZPsXRqH3Wcyj2PTj2v6ktDF4/hQnj+oUW2xlYfu6zd6zNzHaLXwVOkurEIXAtRjjttDOuS+kRTvWOjmQz3ofj6kus/GfCkdZ/DjWJrVIedSo2PckoKh7oHSUF987poxYrUHg26G9Eqz2INq7hdWHfSx+F134XqQ8qlJsCQdf7r+Qq534U/JVrlkvBEp5uVM295X27uUWPqZyJcvjj6GC8IqCRQwtwGelfzBcQRkWpsBVPe+VHioc1vWWHvCYRbyvMn1+4G2VbXAo3oBvE8XvTtnXfaG+AdHgiURoj3ajs0Sk/CRJvj8shODdF70n21dekpfKY6ZijUH4Se93FZ74Om99+j8KkSn8A27oQaasM/qgWHkMUmx1FOawun6guFNkhSgnSoul/KoHJsPI26QpdKD4R/L0vNKDIr9UCGMa0PsuHxAlTfnXKAHQVXzoARpaNq/ZvwxO+rg0uAWZeOL4q8ajir2Ll8qn39cUIrcG9M4uEen077IeJxsKKA2eMZg+6OnjW7HyGKTY6iwKiwI31lnvZm2xvABlYgzBksIsFA9GXlf8kh6rp6Sk0pIH0NYPF9v+ue+TUIgwFyvEDkmfgMLy3TjtMaAXA9GKYn32DWJh+2hHuXJ4meIbuG+pTWfJaVNxDBRRBlV0QJEpKFYe2YqtjtIJGN6anpSOkpv1Jh8lMKEouoGAWbgb7KfC6J8KUXS25xgMLdybmq1fBkQ4EMJM2m6mWjKED5pRbLb8GvDp+M1s5/W5PtRB8OBZEtzs5xuEzhb9ISGinz6P4WWI5t/wJYjv6nex8rgYdGJX1KzW956CVialfFRFAQEUCkWmDv7dNwko5WYoOw2qZwJggIRgi77G6qgs7eO2UMfqhXW68T7hdctO8i4aHf8UMiihsaE/0cBkvEzgH5RcgjyGKLZNhyxmUmBl+hQGa+QGgEIsN3u0F4QQ//qh4qb4sTawTH24v3Yq4f68EU6nKOQpzjTLgDJaSe4D2qPd8wB2eZjJcKWi/Vb6WF5O0Z9i5ZGt2DBXnEBx8Om2QPl39JBPudaikMLKw18D8MHNaroEfpSPVcKlieKnjABldIJ0b/5HdWmzBUrD2mFtu7skrXIZL+DIHQQZ6PKKiH4UGT7ZLIaSx2AsL2O4hqYVK49sxfY9ZoHBVzzbgXDJeseVQKmBm3rYtIfhIVCGHRHbtWCRh/BC4N3yw3SL0wZbUTFgUDh8Hi/t3VT8V6zwgDSmeevbgGp7FzVLjA9N3LZOY4hzednlN7gYuGOhWHlkLx7puRiMG/GXomzVYUmM2a/J82WY6lEG3AxTKgTDuykt6VhoXA8FDq7oF7x91vGdymDxY/BAidSnLdyl1nZjrEJmGjPNrm3ITFT5xdQP+Igb4oyIwpSf3ctL1cPAYO3dTKh3+Aev4BPp5OPqkM4ilXTrLwOGw2bNjonyQihWHtWd7hOTGUx3Fe5lRVUfV4WBkFwzKA8hs/+eLBNK+Rjj4k1R8hA9GIL6TveJaFwMsyhz6xL+Oe2tkOZAkfK4mKa3zByNSKZk/Hys6Wzg8dNOygWYre2aEJcqj0E+dkEMx+fDJ9zafdlFo1dYLD6+M4qLxfmusPtxiDXEUjODmqoaipNHdT62id8rKFuMXYW0IqND4XysyrFjAaNxHnrFUuQhOpyPXa1iH7qirP0bxwFT7Op87HHdXWsdGwdWxT42iR9Jf1fFPhJBH1s3D06x/SJmcjnOhXdyQhdEOBdPpsB7UIothrCb0T2jMpWo2RoE/wriQOm8zt4V8R1ZL23ZnKs46o820gWOIxR5YZFoG/xJnY8VdqhpjNU6VaVLejggNSmoM+DlwE7yiKryDu0SnUl5mIusFl5nuyLq0Hppy0b6fJE85vMj9B0eREF64nij0M1qCjGGuHFjXMTRvM5W7GgvBib6Tg6slVW8yEtbsiivr1AVvF5UseeQoQYL/l7fGe6n5KucnQU/URzXBbi8CQb/0p79n3Jw5Vor1MTr3kNQXgFwxjmPvF7acqaRDBIOYC2+iPQyYbpnijdYLywyTviwV7FVxi3IxExOwnGaLipIY7bC1mk7vbvTcQqxcEd3aYv6PCmIjyzA1wuLMriaVGwxEWf/s8fBscS+6Z6/HVV9aYv6i3uCe5Fz18gPlUPJlgaMC26Vu8fFN25ulVlwysROJrJjVMSFRUvwOqnYnoH29ysE/kJPFERo1Ze2iH4GMbMSYIqyeYv/MshzysVrj0j1ygCdrWO6SkfJnaL7fjDgzCA1LSnPbne6rzj/P2VwLn5hkadxdl4nF48iwFkAhSzOEGLUDWk4N+EFOAHO0VHRbYqXc4nO0LtGWHQ6/owmcHhFmyFqv7BoEV4nFTvgO7sKB31pS9DX3CiDBou3JJgr2PfHYitjgzpG37leWBQjaEfaKF73Kra3ekx/tiGPkseAMuyeNKC6WHimSPzvOQFhmjVr2lH7WxZVacw+WNvkF8oGQX8EHKZE/SUnyvX9wTWM8pO+6SGffrcW8ZCgPJv+8cHHyqsaXvcqthjgFEZMWS9tQTvO4Lwu0WERz1e8rhFh+w+lBm7qOfoLi/oWj4zy9dIWpytbP1i/xf/oK3msFxZtiSKekH26L169jFQJvKhLW8rgyjxUlM5r0Tf4dN88nJoGK1PxUhb02C/RqYLXu3zsadRuZiwapSxU8T23FpE5TVNPDwJjceziet/6U4HHT/6urc+cZqss4/tePK97fezKOL/PpS2s9nN2StgdWmpmKJn9xfP6ICw2GiBLwjYXn8O3LO0UGuLxgn/Rbb4paJ8aRw28PojF49SCW/HVywENOrd4NFcEn4kvPA34Udm8r5GVAyVyoKu3otHpsbki+I4/g4eF1AorB2rgAF9St3T3fy26vtO4x65wAAAAAElFTkSuQmCC\n",
      "text/latex": [
       "$\\displaystyle \\left[\\begin{matrix}\\cos{\\left(\\theta_{1} \\right)} & i \\sin{\\left(\\theta_{1} \\right)}\\\\i \\sin{\\left(\\theta_{1} \\right)} & \\cos{\\left(\\theta_{1} \\right)}\\end{matrix}\\right]$"
      ],
      "text/plain": [
       "⎡ cos(θ₁)   ⅈ⋅sin(θ₁)⎤\n",
       "⎢                    ⎥\n",
       "⎣ⅈ⋅sin(θ₁)   cos(θ₁) ⎦"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "BS1=Matrix([[cos(theta1),I*sin(theta1)],[I*sin(theta1),cos(theta1)]])\n",
    "BS2=Matrix([[cos(theta2),I*sin(theta2)],[I*sin(theta2),cos(theta2)]])\n",
    "BS3=Matrix([[cos(theta3),I*sin(theta3)],[I*sin(theta3),cos(theta3)]])\n",
    "BS4=Matrix([[cos(theta4),I*sin(theta4)],[I*sin(theta4),cos(theta4)]])\n",
    "BS1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAGMAAAAzCAYAAACDtbZiAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAFSUlEQVR4Ae2c7VHcMBCGIUMBhHRAOiChghwdQFJBoAMYfsFf0kHoIAMdQAckdBA6CKGEvI/xHjpbdz5b648L3hmd9a3Vu9KubO3c+tnZ2eZahM7Pz58i2a8yS1iA0ZaeDx4A5P2VunqjnF8KfwvholTzdWecavr7jhBcqa8i5lfr2hm/VXAkad06DjZ2VQMBYX+o6gfsjJEGgsCGNx+Ssqm4P+r7vcKF8lx0rTevVf2J74nq7CjsKn5g9RXfVvy7AuW3Su8FZYdKX1q6ztN1Z4gJ7M+dnicK3xQ/UbhRHOZXisQzRnsnn8d+nrY5IAj0PkK4VxnpNT0RXOOF57YzxAh6b1PPaxiDFH/K0zA7XT1Z4fB/tsTipfjHcD8wF1jO09hYA52dYdpgojiLcE1Pdg0CZSGiIa6Ut9AuuwlDg7GN7xWKdKeMYzGCoLIJFSsMMS1eM7D1PBJ/2cqHT6Wniy3g+4fyWYzh/AD/bd4GoXB6Wic9jzzVFCvhMTKQrSDKV4oEJiDCNzsEFTSPmCOqObQVH4LK7I5QUEHRS9RlZ+RMv/Qaj7HtB0U536diisMGhDoBVAMOQWATULfsEEKJ8vIZraA8W4TUR41NDwClDvIMF2GoLwPaJhEbj1U2GBJYAI36OVA8A1JPVj/C4eABoeO/KP9YT7ML5M+Qyln5NzOZeUJljIGAZ4QVq+sljFjfxbx3xYy+0gF4GOIQJAQ0BVVlLK7KFa06tPupMENqjwCznaY4djMz7jOVgoSXMB6DPotR2zWmCorlfaRZrRCHClY9xGLhGL7wxJPVLP/sqV1oL9aU5pjPTkMItEBttS8MDYZOZcCYKrK8UIdSt09iJV+L54Xg1GCwNDf1HRrwpbryPE2xotCdRbKd0WTFFfvyTHPkdiEBbzYmqT9PYfBG+jHCDSskO5FEyvrKYiVHbZiA9fw6W2t+bsLQJNCZj+FkFEdFfVb4Wour9itjWFFVU4JXBWxJaNCn5V1ENpwHYRfwYXBXTww2z09K9zZBjV8i8cNLHOAjFK4QINLR94jn4vZ/XYWhyXAU7HVCy0ImXr2M97JDVtZzU1OVI40VKhEYhVEJUXcVRmF0h3XlSKMwKiHqrsIojO6wrhxpFEYlRN1VcBWGjouc1be7Y///GslVGIKGu4DePiesumi8X/pcPpitOqhN+ffeGU35GNsJAZedITsxUV9cpHTi7NWX5DRPu3rlu5u7g17yzsBoIwg9+dbTibNXH8LQ/Li5a9VBz2NncHnU2NkLYDVRTmB8vsY5gI+NgyLx1ImDnsfOMG87vtba3TIAc61ZvI4sOXupDuqNgECGSjglxK4BuC2caA5oh2RKFgYc5MxgN2o7e6ktt4AxL73kyTl2wNxiThe22ChPpqXUVA427xAYLsjN2eu5u+H+5nOvYhBVnUyVwhAzSN30ebZVlYdaQTj2XpHs7JU8k/Y6MKAX2TIXNbVQGAIdPY5TV6vOXu3h2FnPUeeGuqMvFIY6M4PMN6e2nL3YZew+VuCpxuH4OCQbErMVYjUj2zWmvi2/0bNKGIDUtrMXqo8wuDtpENXCeFIgGlNFlmeGnHqNaZnT1OCcvRrPtnnDThz0qoQxSGev5pg2btmJg16VMAbp7NUY0oYNpaY6cdBbaDNgQgHjPShnr4aYpjZr3UFvoTDgXoIYpGFNRbZue+HAe0arDnpVaqouz2P9BARGYSSA5910FIY3ogn9jcJIAM+76SgMb0QT+rPT1LZOC/Zqn3WXnx4Suh6bzkOgiLXqZdjbzuCDYPhnVHbxPq+/MT8NAd7oS3j/A2Xi92k3wwgvAAAAAElFTkSuQmCC\n",
      "text/latex": [
       "$\\displaystyle \\left[\\begin{matrix}0 & e^{i \\gamma_{2}}\\\\e^{i \\gamma_{1}} & 0\\end{matrix}\\right]$"
      ],
      "text/plain": [
       "⎡        ⅈ⋅γ₂⎤\n",
       "⎢  0    ℯ    ⎥\n",
       "⎢            ⎥\n",
       "⎢ ⅈ⋅γ₁       ⎥\n",
       "⎣ℯ        0  ⎦"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "M1=Matrix([[0,exp(I*gamma2)],[exp(I*gamma1),0]])\n",
    "M1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "Av=Matrix([[ 1,0],[0,beta]])\n",
    "Ah=Matrix([[ beta,0],[0,1]])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAABYAAAAzCAYAAACT1fi+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAB5ElEQVRYCe1Y7VHCQBANDgUw2kHsQLQCUwLYgdKB/Ez+YgdqCbED6YCRDrQDkBJ8L2TPkNv7gBn+3c7c7GZv92XzuNtcGJRlOcoUqapqp7gtF+LU/AtEfmH89sbCQnA76l4useoBKv6GMcOdl+7c+BngPCF6yorPIsMYVFQh1GwQf42xgO/HlxusGAD8DVbQc4wX2HOMT9j5ycBIJl8j6A8Bgc3VwutX8Wk6VPEUSWslcQVfgZuoS43xIeACMVsG9kT45bwqTmBfNR2ky459YDqBESVJvh14MhUHVSgXV4qvcfkq1rgVHHkarmtVnMDtsmKS9rjikx/RAncCt5HsH7mV9c+/s7+EgNm5bhXgMXzrzlNZIV5gJL4hYws9kUzYpOEB41F8mo5pQqyOTecOmj8W9T2utR2Jqb0EgQGwQ+hMEmK1l4pYEC0uARtWEhWJCsOAMYJbWiKxtW9gs9uN220uU6r2AgOAnewdg28Tts8cI0pCwGxAPFtkuMkzFKuOkrRBDE2JikSFYcAYx6wKOVnKgdCAaIZ3SzMBW5mNh1LsVVbDx8MgP3B4UlIlBrjpFWq2x3kMFR4YeyoBG07ORoWsihxLh68hI7jm2yMo/TwkNDhSMb+Lu3+GyNd+EBgBXOdW7h+r9njPVpljOAAAAABJRU5ErkJggg==\n",
      "text/latex": [
       "$\\displaystyle \\left[\\begin{matrix}0\\\\1\\end{matrix}\\right]$"
      ],
      "text/plain": [
       "⎡0⎤\n",
       "⎢ ⎥\n",
       "⎣1⎦"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "jones1= Matrix([0,1])\n",
    "jones0= Matrix([1,0])\n",
    "jones1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ0AAAAfCAYAAAA80BDFAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAMT0lEQVR4Ae2d23EeNRTHbY+feAqmApIOAqmA0EEIFSR0ECZPyRsTOghUwKUD6ACSDqADgp8ZZsL/J0trrT5pL5+03l18NCNLq8vR0V865+hm+/T9+/cn5gwBQ8AQMAQMgVYIvHz58o5oPff07vrwidIvz1s1YnQMAUPAEDAEDAGPwCsZmK8CGoq/VvyN/L2zkGihIWAIGAKGgCHQCIGnMjQPI1qvFL+rtPtmdCJULGoIGAKGgCHQBAF2Ob/nKJ1u8U5H1pDzwAuFf+aYtjRDwBC4QsBkxWbCHhDQPGWn80jhvTP94Kxta44LqEdbY8r4MQQ2iIDJygYH5TaxNGZDlH9feKDPPwEXHhI8lu8ufEhc24nJr9fmwdo3BPaAgMnKHkbpf89j0YZofvJyjV3OJ4pfgoTd6YDCis4PSnMOlqLbnNECwaX4X4puoRuW3BCBpcZuKboNuz5Iain+a+n6+l8r/Fz+km/8poyOGHoo/0z+pxhlGJX/Rf49YZL3NP7eU1x9eSZ+2Xou4cAM+rtzhsv4kAkjk5VxmKaWMFnJI3U0Lpqfd0WSq5vXit/HK84J1rvTFy9e/K2ED/WRdcpDcd1TWDyCUx4NkP8XZeVxWDi3nbr6HP6psjwe4Jndt/L8xuqHob5CDA2GiIcFn8vfUdpX8nSEBwe/KmzqRJM+8a6cfnzXlLiIiSZnnA8UFo8SlQf2H/m2f9D3Wx+fFPj6rDCa8z+JgSMKiddbhYsfo0H5SmFUHZOVBBSPo8nKCrgI+wMbQppYYZ72nNJPp/xyKMbk517N6ENEUPysujrlqTi7j+/lv4iKjkUvVOA71UXp/KnQGSz/jYEJL9l+VZwzQhztfktEIW/C6STGAsP3k9JqjBF0oOcuvxQ2c+ILus8VFmkrzxlahRhhyrtfrJrDhK/7RuGP8pdz6q5R1vfztuEyKF+FcaiSFWgKa+Y3K9EvGswNkxVAvUG3dVkRf8WNzNkYTqrMKqwzKHF53/EvFTrFH+VhIDBGk51oBEODECIMzin9Z/LCtw9Z9WPYYkWKkaEsvMBv7xjO15sciA4Gi90W/LR2GM2ujylxtemMqkK3Q1FIP9nqzsLU06WdYKTTprb2vRtcNBYseKqPdkWjKF+lwVGdKllRfeYRHmNR7UTPZKUaxdkEdiMrac9GjU5aIfnmOCxnkNjhzDoKgq5osaJnx8KOZ0jBYoTSY69414AwxQZJn/OdeKimUWj1sWhnj7x8vzlWyxmK2UrCt0N7YLt1tydcwHM1TP14HiUrqvtWvnh6ccwkET2TlWOAO77OnmSl18tao9MjxocmH8djCMMTvmc66iEQTODiDsPn94ya0uLdEAp7ztHeTDaPL+7xiXlNiT1XAivZ7mhQ8aDcLtLCE79p7/HEsqsUM1xmw360rMxuaaUKNifywO8dl/Nct7yS4+iAuxHuBA5W5Upj0nO/wuqbIxwMRViJd2+yldY5TxelyoMDHPTjBwcoWo7rSit96pz4NrPHZ8qDF2j2jJKrGP1QPjsp+gDfQanTNrzTDx4uENJHZ8B8He6qSP9RHuOGkcU9kP9NZdKjRpcZ/eAhRGdQovQT1YUP6KW7x099uSFjFZNK47RHuwfjGAr6tukPeARHf7oVcVTmD1/g4O5MZYq4BqKFcJO4FHitSvY4jskXc2FRWZnaiaExVd6tkxVw82O4lrzsWlayRofJLlBR3ChllG9OWfH22ilHheGJM0oRpc9E7Cl9T4s8Li5dnkIUFIIV6KDwpuxQ4Ovg7/qIHpMAvnmx9Uw+awCUjkDzV1AZvM7p2+2uFNIPjA+X953TN3yTjtJl18FTQNeGQvqCgc7dQSmrcxgQcMi5cEeQGiVo4441OvALZlknnhkv+sqfHndGRqEbe4Xsuth9wgPjzLh3fCjOU3aHg8JBXFV3yG0OlyFmK/MG5Us4gv2isjKV/7ExVf6tkhVwU5/Xlpddy8p5OvkEKMrlN5+OAegUTFo2fKtOZ5QUR/HjO4XuB4mdCa/QYmOEcGV3LEofcii+rk0K6hulCe8YGwWO76zRUR6D9qnK8fQaQxdcagzehYwkBJOHqtsZSMVRzBSjTz3eSIwcirlE90tfDqUUVXG7Hwzp6FjElaI47SEoJYcxwbh0uxp9wycu4EOZnEFlwYCxxVBeyE/BVcUO3BZxOWCyNkE4DcqX8hmnpWUFHpinjBdzrbejVVrsTFZiNK7ia8vLrmXlXBh+kGDqVrY+jXuAb5L8E01SOh2OyNJsp5CTxKDMUfIcneF4U88qOV3Vu8yRHwfKV3TihwSD1WlTHkXM+3LaR8jDy7fBulHmwU7L54HNkEPQLwsFUAbdcR5lxB/0Hsn3xsGnMz4Y8rG+g1eWL0+HdnsGWukYoLDrIR9lGBYjil45lQvGlqNJdpnH4lqNi9qmj889a8HIsnsr4V3ExdM4UV3mLgo6dfBLfu7uEUy6BUlScUy+bkJW3oonfG/MEz7dp/qxa1kR/03nhKe3trxUywqD6/tSrUNEKrUhbu6UfmB0/okzxYgTUIUoOgYst2pn9VMyFgxIKuQI7Vylrip5J97ccVw+d3IqihoFRT/d7kx0eTWXUyKTiR5bUO2CNS7d+TEpcN04qCwYB0UY6rlChR9DkzQo59IiApKhTDqucXPwhGuK6xxc1HbxH0c5zg5/DOHiSpfmg9KZNzxjH1XccbMq7zD09Rm7blx9OZOVGLBMXNiFOT8qK6reek4EWdicvMzBRWVb6pCeDckMWS/pTF//9lKuP1DIGAqOdcIgh1yOt1gp5VzpkutglZyrfBNp4t1NHIWszPk9iVO1i7F5qnhQnkuxwk4gxTNuK8UVA4sxvAyFFGcljbI72PGFMklIe7Sbc4EGjwJKLpQZ4psVfA2u1biIecYvGGP6wmJi6PebhnCh/pJuSL5MVq6Q3+KcCLKwprxU49JYh5RsSFZ+MDoHTgwhjCjfsNVHQGIXC3aXrnrUYReU7kQYqI+6glFEdVgx3rSDz3Bp79oWH6w4OU6C/yUdWDjlHDei9i/jb+IemwtFUzzTomPf0AjC0ivr232rxGy/4UGefPjr7ukCEeWFucA5dw2uLXBh4VA69gwsx2ERl7hQ67gwG5Ivk5VrwDc3JzR2yMHa8tICl2uUp8WayUrW6IiHsGLlPJd4eCJ74gXmd4Wp0kb5oHh4nQYosWPFGZSTS4eOPEaNAVzDcYGK8MeO7/TYMC1DeQbgWEd/HxQqY/Rce543cPtM8ctC+anJHHkN4cz9AzuCdExpP9T7THF+IQ3DEjvK8Kd6Am5TcY1pEK/GRTz0doSiiRGK71BoJ3ZjuMRlW8aL8qVGTFaukd7qnFhbXqpxuYZ4cqyZrJznmpTwcnzD0ZpTQghzVA7jgbHgD23GjwL+0nd2y0l9eYwMAhUMGN+r3J+IB5Q4feI4RoFz7MQ4/0VJoVjZ3bnVv74xpk/kURakk8+xI+nfyJMe+oLS5Q95li6Sf1BZ6uUcbcAD9MCy9zw5V2FiGmMW+Duoovbo88fKoG3GNozRa/KooJA5wcRzT9JJk4NPvoPBGcTV1Sj/aIqLeGKM2EXDc8kN4lKqVJvusczKl/JMVq4B3uSc0BitLS9NcbmGezDWTlb4K9P8y+qpXuVfTS1r5fK4CsM/5O/X4iMaj6A1REf5d8fKDNW/yTz4lG+BC33+Rf5OiX/lVeGi+mD/rETf0vNzfy4uwng3c2Ju32rKN8SlWoeIl1k25HzQtlnmUgi8EmF2HsXdR8OGuQ+ivT24aly0Cg27L3f/5L9PWJ0mANTiwq6u9tgzYck+MwjsaU5k2F8sqRqXGZzVykqvqbPe18iHBPeOigw9FRyhYNkgIBw52uMOBQW5mPP0aSc+Hl2svVrCtbj4/nI8yLEgfyWBIzYE5l3Mmy9XhYtocN+5C1zjvu8t7jE+Wlb8WN/InLhJbGtxmcqrx69KVtK25u50uOMI5/cpLfuehwB3PtztHLwIGyPjJwK7JM5ZmRCserhTS39nBGG7id2UmmnmjsZFHLyRZ2FE2DnhkmKwR1y6/tzCiM2J/KAfjcuqOmTueVzNOaTV7Z9zC3vuFRa5F4Au9PeIueHSnyd7HMPWPNucyM+JLeAiHmbd6Yz+u+q8gbVUQ8AQMAQMAUPAXRcc/LvqIVxm3ekMEbI8Q8AQMAQMAUNgDAGMDr/0Z84QMAQMAUPAEDgGgVk25D9bztg33PEHiAAAAABJRU5ErkJggg==\n",
      "text/latex": [
       "$\\displaystyle \\left|{i \\beta e^{i \\gamma_{2}} \\sin{\\left(\\theta_{1} \\right)} \\cos{\\left(\\theta_{2} \\right)} + i e^{i \\gamma_{1}} \\sin{\\left(\\theta_{2} \\right)} \\cos{\\left(\\theta_{1} \\right)}}\\right|^{2}$"
      ],
      "text/plain": [
       "                                                     2\n",
       "│     ⅈ⋅γ₂                      ⅈ⋅γ₁                │ \n",
       "│ⅈ⋅β⋅ℯ    ⋅sin(θ₁)⋅cos(θ₂) + ⅈ⋅ℯ    ⋅sin(θ₂)⋅cos(θ₁)│ "
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "r=BS2*M1*Av*BS1*jones0\n",
    "pd1hv=jones0.T*r\n",
    "abs(pd1hv[0])**2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY8AAAAfCAYAAAAfi0EQAAAACXBIWXMAAA7EAAAOxAGVKw4bAAALeklEQVR4Ae2dXXLWNhSGk0yuekXDCgo7oGUFhR0AXQGwAxiu4I6hO4CuAOgOSlfQwg7aHZTmutMZ+j6K5Mqy5M/+bH+Wk3NmFP1LR6/O0dGPkxx/+fLlyMgQMAQMAUPAEEgRePHixTWlPfPpN7z/UOnnp2lhixsChoAhYAgYAh6BVzIUjwMaCr9W+KPczZOQaL4hYAgYAoaAIZAg8EgG406U9krhG0q7ZcYjQsWChoAhYAgYAi0EOHX83krxkeO13jxkubhLO5P/Z44xSzMErgICpgdXYZYvzxglr5w87sm/eaIf3GGtQTzC3FujY+vTEKgIAdODiibjKrOyyxYo/5bwYc3+Fpx4MH8g1zyIkHgIEiNPD9GP9WEI1IyA6UHNs3PleCvaAskpX1px6vhW4XOQsTcPUFiYPPCz97JUu7MzGjW4FM9LtRuxbsGZEFhqrpZqd6Zh72xmKf6ntuvrP5V/V+6cOO7gxkOd3pF7Ivc+RhNm5H6R+4Kf5D2K41sKayxPxC/HvSUIzGh/E2RY/D9NwuJK6UEYuclAQKLt14qL+LohTnnaeK3wLZzC3Bp9Pn7+/PnfSvhakV7yjXC99ZfcTV8Ya+SOML2VfabK8kjOp18/yvHbiV+H+vIxGBgUHtDvyl1T2mM5mOVh/YP8WUltAgzfLDOON7M2rsbUJveDt+UXr+iUx+J/3ff9VvFPPjzI8/XZDUziX/UNi0GITy8krKvSg+kjGtaCxr0ZfRg2onlK1YKL+OjYAtI0SuS1RUo/HvRLgirIAs5OqVkEFeY08JPc/Var/ZEzZb9RXYToT/nO8Pg4hiJ8efVBYe7XIPr9kYB8vjdmICx0GLD3SptiVGiH9twDkPzZSHzR7jP5xbaV5wymfIwp5d0v34xhwtf9KP+d3PmYuklZwyIBZMHoJD2AL80188WO8P7Eeae5xUk8bk0fFseEDmrHRfwVDxYnuxDyg/tBvlvAo/Is9BiVwaQ2gsHgBIPgO1L6z+SFuPfZhWOg4gURY0FZeMGQta63fL3BntrB8HD6WeKDAYxfM8aUKfXpjKN8d2KQzzi5hhqFqW+XfoKxTbsaFFe/hsUgpKYXEtaT9MDLCHKCAdkKbUofDgjqZnHZaTwEItdMzYkjApUTx6grFuqqLXYgnCA4gfQtlBiT9Dop3sWjOLFhUXQ8iYfJbRR6faC2s1dJftxcV+UW/NELgu+H/sB2b1J9w2Jv9MZV9HO1lx6o7ie5n8f1uHrpzenDgRDbLC5DjEcHQwku104I/sNO5u4E6iH8LFTFHb/PbxknpcWnExbeMVdmuzmbqYTHJ+Y1bfmZEth9NlduCoeF/ywtPDBOfw8Glj1YMcOiCPXeelBssdIMk4H8xGwdl9P8sC5SNTgEnPcHdsNcjbDgh51x872v0hpSWRZBFkce1iHeJuKHdRZMrsFKO2/qHPk+s9dSyoMX2mwZF1cx+qF8TjaMAb7D4kzf8M44eKDHZ4zOEPk6vOWQ/k4OI4WxhG7L/aYy6RWey4x+8ODfGIYo/Uh14YP20tPcd75cn9GJm0rD9Ee/facdwyJFbaG4n+dF9WAM654fZBldCIQsNyeYqMwfvkDnXVFlijoVGs341ekDPEbjXQMTWKgSFxgbQr3GQw3wXa9b5OSHT2tZ3Fi8WVxbi7fKsDiRxyOey5OPsKFEoR0masiJgbY6f1NF7aEAzhgp/EQuu5ArnUWavwjJBDWkuDvtyGccGBEeqRtSHL5JR4E4BfB5mutDPmPhcTr3RqOshjAE4JCj8NlxalxoG9rXeMAvmHVI/BoWHVSWSxDei+vBGO7FD7qKnPOntJ2x8DzyhsgJmJsA5A8dR+cbGVSYz+edDsjvlSPVLVFV+gCTGsvamMBGdbjA1FA6HVpQYDc7WoVZwHHNwuwng5MCX03FRgVFyp4glN5HCHHTJwUVRwEQcoyGPLfQZo2H8piY71SOT34xWIHSRf1zyEh8FOiO6jaGTmGUjGKMqcUbiRGhZKV2f/Dl+BIrquJOI3xy2yhunDkgTH8oRI4MixwqC6Rp/piDpfUAHUAGz+SQo9YJQmkpuU/gVa45ZagAMgoF3aBMblPEpo8NE5sd+huiUyrWotr0AebWxgQeasQFvgbRqUp9lSspYWFg4eopLeIW1iQxLMos1lxJQfz+AjuXdJftMnf86Cyiaid+MO+tTp9yLKh8u0z/KHT4Uqu3bpTZOfn4PLDpI5TsvFAAxW+uySgj/mjvntxL4oF8Ou8YGORdYwevLF+qu2ksPA6cXqFgINlFlzAuYnHRxMVP3+6vimVxi8tG4eZUHaXFwUPowSd1iCttnBp+/BiRuVZZpWNIwimEfHD9Ta5FKhc2TFz3cuLfR6dq0wfme21MwHkyLpoPxjKXbmRtAYzmCOPxTy5DaexWS4s+wKeKy05o7OKsKnlCUPM5o1JZcAGWhdmdltQuX3kVH+pHtT6ysJ9oaqUnsfDQ3ZxmVBaMwRQasrj1CSJtbBYL8V78hzQMLEO7sHBVhDEyvMsoZ5rvTapND4KxLW0EGUwok+p0PFDkEZpNjlbUhzDe6jAB4DG4qPiculGyBbDVoROl/NtJvUjg2ojdTY5KDz2dnUuu8iHSxLsTEPnslm7KHatfjMYjhYMiLMUKu7O+BT/FFUOJUWuUV2F2fOwWOyewAtP0R78dUjubxkIDYs6CIWV8bAT6fiemiAWVD0DV6IHGGuQn/FWI3PBDmT6Z5W1kXzmqSh8qwYR5mIyL2phTN0q2AF47hPEoUaysTRkJEAsvp5L0ZIAAXm8KRgHVYed/aILP8Djt+hYf7Ow5qsP/kgQWTtHiTtR/YxxCusfmTPEUz1BkqE8bYRFI62wdC4x+6QoxHSvxPixy5edMq0oPvMyxWcnKPPInRz6y2bxhBkCUF9YB3gj2laOq9KESTIB4DlxW042s8RC47EB+l58uvggSQsS9b7pQsRsMgqbgxfFL5bgDRjjXIB4T090U8fQ6Li0DryxA+xLj5bPeHGG8XH+eN3D7XuHzXOERaVwn9OG8WSyETetUpnGiMO4roQI+u7AoVJsluUY94KMPTmqpPsNrkJnvFeYX1jAQMVGGP58TdGaoHMVt1KgPa2MCPpNxWVM3TuMZjsIYARZ9/iBh/Pj9l+LZ468fBI/lCFv4Tpw4ir4GsRhz0uBYF/rnZMQdIQsPSsJ7iNuRKY5RfCjHiYF08vn6ifSXcqSHsaBAt+WaL7GUF9NbRaiXI/qAB9oDy9ankbkKA9OYs8BfWuXSYCHcmBdOshiIEvVhUaozS7r4w9DVpAdH4gd5/0YDRO7Q66Cfr8lj4PK5JgVT9xk8aSJklHgwHL1y5Grkf9SmD4x3bUxAalZcNKbD6gZ/VZd/RRs7pb2K4xZu4zMED2H4h9ytIWX7yqiNe7S1o8yNXWX66i+dB29yc2DBOH+Ru1biWXlVY1Hi+7Kna17mkoHN60M81zPiMlk3xEvHFsS8puHT/EbBUmdAgBMYJ4HSaWCGLpomeC+hv1ppMhbaVYVdsLuX9/Ej+en1ae1Y1DpHS/M1WQZGMLglGZiMy1q6cZJOiBi5prS+T9jSKhbPICAcuTLjnplFbzHy7dMP/VVJU7HwY+S6hWsWftuZ4zkLxOd4wL5c1VjE/F6lsObG9CEz4VNx8TK/im7kTh68AYQ7zsxwLWkEAryJ8PbR+YplVxteKDi1cH/PgsgOhTen1i97KQ3BOcTpRt1Mor2xUK8f5djU4DckLNJxbwWLZgxXLLC3DFxCfYinfm9c1Mh6ujH2niu997J4/3uI8OUu8skSONEu7S/R9hJtGhb9srIE5rW1aTKQl4EacBEPo948Bv8b2thMWtgQMAQMAUPgciGg013n39D2jbDz5tFX2PIMAUPAEDAEDAEQwHjwy0FGhoAhYAgYAlcbgVG24D+GeA1KBtTalAAAAABJRU5ErkJggg==\n",
      "text/latex": [
       "$\\displaystyle \\left|{\\beta e^{i \\gamma_{2}} \\sin{\\left(\\theta_{1} \\right)} \\sin{\\left(\\theta_{2} \\right)} - e^{i \\gamma_{1}} \\cos{\\left(\\theta_{1} \\right)} \\cos{\\left(\\theta_{2} \\right)}}\\right|^{2}$"
      ],
      "text/plain": [
       "                                                 2\n",
       "│   ⅈ⋅γ₂                    ⅈ⋅γ₁                │ \n",
       "│β⋅ℯ    ⋅sin(θ₁)⋅sin(θ₂) - ℯ    ⋅cos(θ₁)⋅cos(θ₂)│ "
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd2hv=jones1.T*r\n",
    "abs(pd2hv[0])**2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ0AAAAfCAYAAAA80BDFAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAMXUlEQVR4Ae2d23EeNRTHbY+feApOBSQdBFIBoYMQKkjoIEyekrdM6CChAi4dQAeQdBA6IPiZYSb8f7KkT6tP2vt6d/HRjCytLkdHf+mco5vt048fP56YMwQMAUPAEDAE5kLgxYsXt0Trmad3x4ePlX55PlcjRscQMAQMAUPAEPAIvJKB+Tagofhrxd/K3z0LiRYaAoaAIWAIGAIzIfBEhuZBQuuV4neUds+MToKKRQ0BQ8AQMARmQYBdzh8lSqdbvNORNeQ88ELhnyWmLc0QMASuEDBZsZmwBwQ0T9npPFR490w/OGvbmuMC6uHWmDJ+DIENImCyssFBuUksddkQ5d8THujzz8GFhwSP5OOFD4lrOzH53do8WPuGwB4QMFnZwyj973ms2hDNT16uscv5XPFLkLA7HVBY0flBmZ2DpejOzmiF4FL8L0W30g1LnhGBpcZuKbozdr2V1FL8T6Xr63+n8Cv5S77xmzI6YuiB/FP5n1OUYVT+V/mPhFnek/R7T3H15an4Zeu5hAMz6O/OGS7dQyaMTFa6YepbwmSljNRoXDQ/74gkVzevFb+HV5wTrA+nz58//1sJn+qj6JSH4rqrsHoEpzwaIP8vysrjsHBuO3X12f5TZXk8wDO77+X5jdVPQ32FGBoMEQ8LvpK/pbRv5ekIDw5+UzirE036xLty+vFmVuIiJpqccd5XWD1KVB7Y3/Zt/6jvdz7eK/D1WWHMzn8vBkYUEq83Chc/Rq3ylcOoOiYrGSgeR5OVFXAR9kc2hDSxwjxtOKWf9vnlUIzJL42ayYeIoPhZdUXlqTi7jx/kv06KdkUvVOCN6qJ0/lToDJb/xsCEl2y/Kc4ZIY52v7+KOkWOocC6fq10Vz/kjQihBWju8mtE/WoV8QbdZwqrtJXnDK1CjDDl3S9WVYkWMnzdtwp/kp+KR6GFeZN8P28aLq3yVUF4kqwIZ35/gjnFHGeR+LPSpizcTFYE4nW6rcuK+KtuZM66gFJlVmHRoKTlfce/URgVv8/HQGCMejvRCIYGIcRwOKf0X8gL3z5k1Y9hi4pU37SHRwAmO9FDCNltwc/cDqMZ+5gTV5vOqCp0OxSF9JOt7iBMPV3aCUY6b2pr37vBRWPBgmfy0a5oVOWrNjiqM0lWRBcjg1wht8h248i61m4tXXRMVmrgLJe+G1nJIeg0OnmF7JvjsJJBYocz6CgIuqLF6otVGDueNgWLEWoce6n8O/nqjkzlBzvRi0ZtcOX2Co9Eu3jkpXT6zbFayVAMNqi+HdoD2627PeECnqth6sdzlKyI73SHzZyaPM/Fz2Qalcm5pzlR6cIiybvFZarROUJTk4/jMYTh8VFmdwL1MB5M4OoOw+cPNmrdzS9fwuOT79zShp/pg5VsPO5QPCi3i7TggDjtPRpQ/tqLGi6DIR8tK8I6nX8sboYcgw9mdGwFmxNl5PaOy3mpW+oUSo6jA857uRM4WpUrjUnP/QorJY5wMBRhJR7fZCstOk8XpcqDAxz00wcHKFqO62orfeqc+DanHgmwo6AP8B2UOm3DO/3g4QIhfXRCqZA63FWR/pM8AouRxd2X/11l8qNGl5n84CFENChJ+onqwgf08t3jF75cqizSql1x2qPdo3EMFX3b9Ac8gqM/cfeYlHnvCxzdB6hMFddAtBJuEpcKr5OSPY5d8sVcWFRWxAdyi/y1LuDaxlR5N05WGHz1m/FZS152LStFoyMwucxlMqKUUb4lZcXba6ccFYYnzihFJjITsTGRPS3yuOR3eQpRUAhWoIPC67Pqgq/i3/VReqdTu27CKGTwotO3210ppB8YHy7vo9M3fJOO0r2Q5ymgMzIK6QsGunQHpazoMCDgUHLhjiA3StDGjTU68AtmRSeeGS/6yp8ed0ZGoRt7hey62H3CA+PMuEc+FOcpu8NBYSuuqtvmNodLG7MT81rlSziC/aKyojZQmMg4rxufyhcXS0pvHVPl3yhZYdzV57XlZdeycg6IqROgKJfffRoGICqYtFwaV51olBRnMuOjQveDxM6EV2ipMUK4xuxYUHyxTdE40Td8Q+9CHqFurNKVljoG7QuV4ek1hi643Bh8CBlZCCYPVDcaSMVRzBSDhwZvJCYOIa7R/caXg/+kitv9oBw6xyKtlMRpD0GpOYwJxiXuavQNn7iAD2VKBpUFA8YWQwn2fXBVsSO3RVyOmJyaIJxa5Uv5jNPSssICAz4wNgqcjBeNjvJMVkCo6daWl13Lyrmw/KSJp1M+wTBwD/Ayyz/RRKXT4Ygsz3YKOUsMyhwlz9EZjjf1rJLzVb3L7PhxpHxFB57xNeGJJGlTHkXM+3LaR8jDa55YriNS22mBTZtDMV9WCqAI4nEeZcQf9B7Kx3Hwac/IlwvGhF1KjS54FfnytGi3gZvSMUBh10M+7YTFiKJXTuWCseVokpXzWFwn4wJHap9+Mm9Z4KQX5mTnropLKCgazF0WErmDX9or3T2CSVyQZBXdztGnleTrOmSlC5fIsvqxa1mhI+rDbHPC01pbXibLiu/HLDpEEOc2JM6fUgSj80+aIWac4lKIomOwSqt2Vj81Y8GAOBoKg0Nohyr1UPcoFG/uOO4oY1gCggfo9NPtzkSXV3MlJTKM8ojSahescfnOD8WES8eh+g+Srooe/WybpMFo1RYREAtl8nFNG2LccbPiOgQXlYWHYCACno6pyo82XFyV2nxQOvOGZ+wNY11pJyarvMPQ14fHdFwpZ7IS0SpHhF0Y205ZUdm550SQhc3JyxBchOycOqRhQ8qjdkg9U/Tfw2cjhkLGUHCsEwY5FOB4i11FydUuuY5WyaXK15Em3t3EUcjKnN+TOFW7GJsnigfluRQr7ARyPNO2clwxsBjDy6QQfAblSjJGs+33eGiPdksu7Bp5FFBzoUwb36zgp+A6GRe1zw4DIxD4rfUnpLfhEsosFbbJl8nKFepbnBNhbq0pL5NxEbxz6pCaDSnKDkbnyElwEUaU72ufiYCkLlV4MV31qMMuKN+JMFC3Y8EkojqsGK/bwWe4tHdtiw9WnBwnwf+SDiycck4bUfuX6Tdxj82FojmeGMja8R5VcweNICyNPN/uOyUW+w0P8uTDX7ynC0SUF+YC59xTcJ0Dl8BW37CKS18CY8oJszb5Mlk5gLq5OaGxQw7Wlpc5cJlNhxyGq1+saHRUNaxYOc8lHp7InniB+UNhrrRRPigeXqcBSupYiQfl5NKhI49RYwDXcFzWI/yp4zs/NszLUB5lNdbR3/uVyhg9157nDdy+VPwyLa/vfOfDBErvCtLixDnyasOZ+wd2SvmY0n6o96Xi/EIahiV1lOFP9QTc+uKa0iA+GZecYI/vLlx6kBhVpCpfomaycoB0q3NibXmZjMsCOuQwah2x81K+GOKYgqM1p4RgMCmH8cBY8Ic200cBf+m7uOWkvjxGBoEKBoxvlOUaDiVOn9hihvbZiXHOifJGsbK7c6t/fWNMH8ujLEgnn2NH0l/Kkx76gtLlD3nWLpJ/VFnqlRxtwAP0wLLxPLlUQWXhhd0iCrTmGLPA31EZ0aDPnymDthnbMEavyaOCQuYEbbhntqTJwSffweC04upq1H/Miku9mUZOKy6NkjN+eCyL8qU8k5UD1pucExqjteVlVlzUn8k65DBkPWL8lWn+ZXVfr/Kv+pa1cmVcheF7+XtT8RGNO/K/yt+q0fJl3tfyt5QuXufC5SG02vo2FRfVp42nbW1YXnn+D8GFcZSfQ1YWnxND+jW17Iy4TNYh4mWQDTnvYZesyPwIvBJJdh7V3UdXk1qdhF2Gu2fx3ycK86NN7oNobw9uMi4DOjkVF3Z1jWPPAW1b0f4I7GlO9O/V9JKTcVlLh5wN6buYvKXybU8Fh5C7sWWFI0d73KFgOAY7X49jMI6/+GsAbI9Roh9SYr4c7aTHo2mRTcU9n6Nx6duZOXARDe47d4FrX1y2WG5Pc+I68ZuKi+qje1bRIUN3OtxxhPP768T4/9gWdz7c7Ry9COvR2bcqwwKAMDpNpHznxKTK02L5jUZG4+IFif5yV4PxYjXIXWP+uzR7xGWjw3UtbNmcKMM8GheRW0+HDD2Pm3oWafUP59zCnvPURe4FoAv9PeJtuBzmyB7HbwmebU6U58QWcBEPg+50Ov9dddnAWqohYAgYAoaAIeBeth79u+o2XAbd6bQRsjxDwBAwBAwBQ6ALAYwOv/RnzhAwBAwBQ8AQGIPAIBvyH1Cd2DeALSxYAAAAAElFTkSuQmCC\n",
      "text/latex": [
       "$\\displaystyle \\left|{i \\beta e^{i \\gamma_{1}} \\sin{\\left(\\theta_{2} \\right)} \\cos{\\left(\\theta_{1} \\right)} + i e^{i \\gamma_{2}} \\sin{\\left(\\theta_{1} \\right)} \\cos{\\left(\\theta_{2} \\right)}}\\right|^{2}$"
      ],
      "text/plain": [
       "                                                     2\n",
       "│     ⅈ⋅γ₁                      ⅈ⋅γ₂                │ \n",
       "│ⅈ⋅β⋅ℯ    ⋅sin(θ₂)⋅cos(θ₁) + ⅈ⋅ℯ    ⋅sin(θ₁)⋅cos(θ₂)│ "
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "r=BS2*M1*Ah*BS1*jones0\n",
    "pd1hh=jones0.T*r\n",
    "abs(pd1hh[0])**2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY8AAAAfCAYAAAAfi0EQAAAACXBIWXMAAA7EAAAOxAGVKw4bAAALYUlEQVR4Ae2dWXIUORCGbYef5gnMCca+gYETYG7AcgLgBiZ4gjcCbgBzApYbwJyA5QZwAxg/T0yE5/9kSVZVS7V0dXVXtTMjZKlS+y9lpraG3fPz8x0jQ8AQMAQMAUOgjsCLFy+uiffM8w+9/0j8s/16Yvs2BAwBQ8AQMAQ8Aq9kKJ4ENBR+o/A3uaO9wDTfEDAEDAFDwBCoIfBYBuMk4b1S+FC8YzMeCSoWNAQMAUPAEKggwK7ja4XjP3Y3dechy8VZ2oH8n7mGGc8QuIoImFxcxVGfT581P9l53JN/tKc/nGFtgriEubeJiq1OQ2DCCJhcTHhwtrlpbbZA8cfqPzr7JjhwYf5ALl6IwFwHqSFP11GP1WEIzAkBk4s5jdbWtbVoCzQveWnFruOmwmf03O48QGFN5Adg5bWNVe7KG1oocKz2j1VuoRvG7onAWOMzVrk9u9c7+VjtHlquz/9U/l25M75xazceqvRE7lTuQ4oujZH7JHeOX4t7nH7PMaw+nardbPvGILCj/NmR4XIxZMLhSsmFjXtVVKeKh9p1qJZytfFG4WOcwpwa/d59/vz5P2Jc10cj+UI43vold+QTY43cFqYxs49UWi7Jefr1Wo5fJ14P+eVjMDAoXKDflbsm3hM5GsvF+mf5syS1nXPC2/KLR3WKQ/nf8B18p+/vfTrr87MqeNsn3ybTqq2GiwZAOExKLtQeFAZv+ZHvlc8nlTmrcb8qeKifC7YAnuYB87NC4u92+pGgEqLAWRlF5acwu4G/5O5XSm3+OFD0W+Vl8vyU7wyP/8ZQhJdXnxXmfA2i3tcXQSdowRLeF9/lD3FT9NVGgH8m310y5dqoOGc45WNUSe9+hJNLW+L5vN/kv5czXDxQM8FlkFyoj7zDZ94gGyzsPog3ZLFFOZRXnLOKW4rULsqdmzxcWTw0XsWNxV7bDPCD/VB+VOA+D4oeo9KZVEYwGOxg4isv8T8SVyuI1TcGKipCfVMfjsGcC2EEY1/rjVafnJGU71Z48ukvx1C9sPXlUk8wuvWqpvZtuPgR0VgPkgsVg7FAhpBRFniVY19fTWdP5WB4OBUY4yHN7Mbd8MhPnVbjoWwcM8UdR1IMO45eRyvkVVmsPFgpsQNpUpAYk8q2Wem/y30Uf070QG3Obv19/zmuyin83gbS10N9YDx1MlySEfJjtpRcqJh0h8C8iQuupIpeQbVncBmFCmc57obH4mh2MR4LuQQkx05M9EcLke0M8mEEmJzFlY2P722c2qtfXwqPU31HlTbgmT5YdcYjBoWD4j9IE/YIU9+DHunXntRwyUK+tFwIz3SOsRDpc5ScbcwYTBv3Kqpzx2O/2p3qlzrHhOb+gdUMRyIo/LAiju99xYuktCg/lCIX6xBnsOnFOoqSY7DSips8O77OQdtvXw7tQaBoe6AvKj/uYBQOaX74BAvnxkrDLgk8KIf0EP1IV30X3Mu/XPxHw3DJjjswjHB9V3fLp0sVQpq1LUx91Jvd7YTMSZ8NlwDKiL7He1S5UB3IKLLWuOhSfHEuKw755uEKPrLvDJHPwx0n/PdyyBTzF7oth0zVj7ZdZPJnkvLg+5aVbcVdOTyS8WoMNhoP5eRdr1Nu8sPTWpQakxRQK5NUaRgA4rjMdnHymagITSgHZdVlZURZ2X9TRfxOpLppI5fP/BPCzljIp1zOiFnxswOiffSNvkaFrTDPhnmaFi6x+dclmfyR9F3cOflEGALwyFF4flw3LrQHim25+Oz8FwNIH4ukdhsuRXRWHyG8R5cL1YEyd4s0hU/lsopcfBY+xbmseOYdC0PkJpK+kWf4zK8DOScbJBCPOctjjdzdJUkCTU4e1GbDI4xOT3+/a3qBHFeyCjNRcVGZiodCYqfAq6nUqCA4y+wgUOaxTpWxo28mKeUxeXmxUdlBiFcn9/RX6eIuQwmYLBBGDCJNbtJj7BAIlDv13VKY58Mhn1hFw0AcRF2/XWjxz0PPoh9pLKs5ntwuazyoj7FoIsOlCZ0Vxmkc1yEXKHpkA6NB65k7WeMhPgq8y1wuzVvKPlE9cQGoMIswsZ1sVmQWZkJTlAfDIxmgPsF9Jf4jl0GTgYEOR0/1JG4C1ZhhhY2C5UgKuiHHCr6+unaRLX8WlKfKwSjhSoIRi/TtR6AqacXHkIRdCPEI9xe5ClGXHDyOpljRIUy8g6YvGMPwukXBImF0UmOTJqTueCxAhMoGc4zHS74DeT73GBjmpmMysoAb5WTJlzV7XHw/2NFCwViywyzh3YjLRTFxDP7WdxHDkDbx40474YXgOuSibU6EtuwIH+bcMnM5lqHA1/QjCbdhNkV5mDUeGkswX5UcZG1BMr6VIMbj3wrn8gOLXFL6KJ+6kJ6I10WhXtbQEBIo7pirIUlbVFAoJQNI/pCm3pe0bPoKIaAMEsrd7bzURl6MtR1dKXmV/IDDrO/IwkV3XL0pLfWDLdQmnKRpElDiQ59njYv6wdFLxF5hlDQr8PADVgUr1IaLS6xymAudlXGlhvzH1OSCVq5sLue73I8rzMO83oQ8zBoPNX6VclCyBdkB3RP3v2zMxR0Aq/wclS6+Flbwucxr4oWdS0mZ0IyQJkzeXNO4G3EKVz47kCO5XSVEcT1WOBiXXF5WeE1l1/HFYGKQUGCOFGYHxO4ptNXHFD3qo94ShXJmjYs6B/bBqNJXDHrT72PacKGMsWgyciHMlp3Lq8BmcvIwdzw0KKuUg5ItyI49xqNEqWDGNAIbZcmupL4zQCndiAmTgPKwWl8rqU4UMMqZti4QbZIjnnTx7iYkVFzoP/cD9DlccLskimd3wPFXtnyX6ELhO2H1385TXuqskHhgxOq4jmslXYcPyggGYiG5r3sbcMF4l45PFvotRiMuuQwr4k1KLtSnZefyKuAAi0nJwxbgsTE5yBoPKRhWaV/l1xUmChVlyhlvXUGx8gsKV8GL82Ol4zgBZbUJ4lKP1Wi9H7Q1tOmOwvxwCaFKiTS8tApHd1xsg0tKfIf4lB/C1MEzxhxheFx5vlzqu6PwWS5xDx5HEqFvpWyzx0U4VXZo6ihC5F7QFTrdBZdC1kHsKcpF17lcn+8AgRFelqYqD7PFY5NysF+YBRgBlD7/IGF6+f1L39njDt8JLssRlvB7Cb4R6o2Q6kaZ/KnKORekP6Fdb4ijUfI5FkKxuGeO8ESsjvgOhgGFzk6DLaI8R+yyKLduRH20897pL8Y2R4/EJD91gWnlqXAuQ0ceY9eIOW2W2xpc1BcMPzs3xrFErbiUMg7hq20YuSnJReNc9lg+U5/djlrfzF/mKvMUPlifef5Lzw/zDSV8Wy6+xFJ8SlOUh63BQ7ivVw74V3X5r2hTJ96r9NvCVXz64CEsf8gd98mTS6sy7lFWLi7wFH/Yliak3bRPO+VWgQt9/iR3rdQnxc0Gl1IftoWvsVjVuG+FPKwQj8FyoLYs2IKmebefLgssPAoC7MRYmYXV2SiV+EK5L6G+OdBgXLTSCjtEd2flv3fk13eDc8JlDmM3pI2Dx71H5XMY98F4bEoO9uoDoYZwztn0jLOexb4bEBCeHHdx74KiG418+dRDfZOnobj4/nIUyREkv3Zmy46y+J123qebDS5p27cxrPEweUgGdigefn5vRA5yOw/OOsNZf9JNCw5AgDNgzo4XXnW1leknB7sWzuxRgqxUuHuq/PhRPCbQOnY3qmZltDQuagG/6WChgx9JuNQxmCMusT9bGlh63LdUHpbGQ/Njc3LQ95yr6QzM4sp3I8KZM8nTMTCiXMofo+yxyzRcynNmbOw3Wb6Ne3Xcp4CH2tDrzqPzf0O7pSsg65YhYAgYAoaAENCubuG/oW0CZuHOoymxxRkChoAhYAgYAiCA8eBHckaGgCFgCBgCVxuBXrbgf/1ODUoUwo02AAAAAElFTkSuQmCC\n",
      "text/latex": [
       "$\\displaystyle \\left|{\\beta e^{i \\gamma_{1}} \\cos{\\left(\\theta_{1} \\right)} \\cos{\\left(\\theta_{2} \\right)} - e^{i \\gamma_{2}} \\sin{\\left(\\theta_{1} \\right)} \\sin{\\left(\\theta_{2} \\right)}}\\right|^{2}$"
      ],
      "text/plain": [
       "                                                 2\n",
       "│   ⅈ⋅γ₁                    ⅈ⋅γ₂                │ \n",
       "│β⋅ℯ    ⋅cos(θ₁)⋅cos(θ₂) - ℯ    ⋅sin(θ₁)⋅sin(θ₂)│ "
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd2hh=jones1.T*r\n",
    "abs(pd2hh[0])**2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY8AAAAfCAYAAAAfi0EQAAAACXBIWXMAAA7EAAAOxAGVKw4bAAALWUlEQVR4Ae2dyZEUORSG6Y4+zYnFggEPGrAA8IDFAsADCE5wI8ADGAtYPICxgMUD8ACmzxMT0fN/KilbqZRyqcqsyqx+L0Ktffv1Fi1ZcHB6enrByBAwBAwBQ8AQSBF48eLFRaU98+lXvf9Q6SdHaWGLGwKGgCFgCBgCHoFXMhSPAxoKv1H4m9y1w5BoviFgCBgChoAhkCDwSAbjdpT2SuGrSjs24xGhYkFDwBAwBAyBGgKcOr7WUnzkYFdvHrJc3KVdlv8zNzBLMwTOIwImF+dx1ZczZ/EnJ4+78q8d6g93WLsgHmHu7qJj69MQmDECJhczXpx9HlqXLVD+seaPzr4ODjyY35erHkRI3AZpIE+30Y/1YQgsCQGTiyWt1t6NtWgLxJd8acWp47rCJ8zc3jxAYUvkF2D03qZqd/SBFhqcavxTtVuYhiUPRGCq9Zmq3YHTG1x8qnFv2q6v/1T+HbkT4ritGw91elvuidyHGF0GI/dJ7hQ/yXsUx5cY1pyeaNwc+6YgsKP9xZHhsloy4XCu5MLWvS6qc8VD47qqkfK08UbhY5zC3Br9Pnj+/Pk/SrikSCv5Rrje+iV3zRfGGrkjTGtln6myPJLz6ddrOX6deCnUl4/BwKDwgH5H7qLSHssxWB7WP8tfJGns3BPelF+8qlMeyv+Kn+A7xb8Pmayvz67g7ZB6uyyrsRouWgDhMCu50HhQGHzLj3yPzk9qc1Hrfl7w0DwbtoA08QH8WSOlH/T6kaAKosDZGVXKT2FOA3/J3au12h65rOy3qgvz/JTvDI+PYyjCl1efFeZ+DaLf1wTk870xE4G5MWAflDZro6LxMd5n8t0jk8INUp4znPIxqpR3P8JpFGxJ8HW/yX8v53BtKb7zLD9Pw2W1EhvJBU0IT2SCHeK9EdaftuDDIs8qby1a6LqfWzy0XsWDxWEXB/jFfiDfKfCoPIoeo9Kb1EYwGJxgYHRHSv9IXoh7n903BipWhBgLyjIWDFntesvXm5uHEazmmg5Oc3FGUr7b4clnvlxDDcLWt0s/weimXc0tbrj4FdFabyQXnlfgF5TcxqT22JBxKzDFhzSLW3fDI89SncZD1bhmqk4cUTOcOAZdrVBXbbGj4QTBCaRNQWJM0mNzvBNCUGLDougs6b7mmT36+/lzXZVT+IMVge+H/sB47mS4RCvk12wtuVDd73Ifo+Y2Dqq9qWRrketueDRZqo/xaNQSkFw7wegPG5ndCdSD2WHO4s7G59eMk9Li0wkKd8iVWffIRi7hcYrHnPbwTAnsOqurN4WD4r+cFu4Zp7/7PcvupJjhkoV9bbnItjbDRFv3+qIsHY+j+nTqMU0Ohub9gV0wVyIo/LAjrr73VVpFKovyQynysA7xNhE/rKMouQYr7bipc8H3mb2WUh5joc2acXEVkz8qw3gwNIw90BelVzu1qMwPX6DxnqIynJLAg3ZoE2Ie8WlolXr2l4f/yjCcJVcnMIxweqq74cu1GZ24qTRMf/SbPe2EwtGcDZcAyoS+x3tSueg7/DZeVh7yzYcr+Mi+26D5Orxxkv5eDpmCf6GbcshUerXtMqM/s5QHP7esbCvv3OERrVdrsNV4qCbf9TrlJj98WotSQ3kDak15qwwLQB6Pdi5PPkoXoQntoKz6nBhoq/Fvqqg9mNYZI4WfyBUZVnmMkcdn/glhZyzk0y5vJ+z4OQExPubGXCuFrTCfDfNpWnjE5l+XhPkrUrx4cvKFMATgkaPw+XFqXBgPVI1lFe39FwPIHIukcRsuRXTGzxDek8tF31FrLGx8irysfPiOjSFyU5HiyDPp8NdlOScbFFAaPMvHGrm3S4oEmp08aMyGR1idgf5R3/ICudrJKowCx1XKVGkoJE4KfDUVGxUEJ3uCUHobocyrPimoOAwNo2I05DkFWzQeynef/qpsdcpQGswCYcQgyuSYHmOHQKDcEZYbCvP5cKinpKJhIA+ir98u1PzzwCfxxVGcy26OT27XNR70x1q0keHShs6IeVrHbcgFMoGcwafwU+1knUwHBd6Hl0t8C1/eVh/VBlBhNmF0wxhqMktiRHOUB8MjWqAhwSMV/iNXQczAQoerp7SIY6AkMeywUbBcSUFX5NjBp7trl9nxp6E81U7bFVGtOT9+hKpmXJSOIQmnEPIR7i9yNVK5IBBcTXHSQZj4Dpq5YAzDV1+1ekkEYY6NTZxN39W1ABlqG8wxHi+JB/LpvGNgmLswADfayZJva/G4+HlwooWCseSEWcK7FZdVM9Ua/K14EcNQNvKrk3aUFoLbkIvv6gxX4/UwgNgXPvDcOrwcN/M1jkThLszmKA+LxmNkOcjagmh9a0GMx7+1lLMIFrmk9FE+qZCy6+ijUM96aAkJFHfN1VKkKysolJIBpH4ok84lbpu5QihtlBXK3Z28NEa+GOu6ulLxOvkFJzE9kYWH7mr3prL0D7ZQl3BSpk1AyQ9zXjQumgdXLxX2CqOkOZmGH7AqWKMuXFxhtQMvdBnoWsMdkbnJBcMdjZc75t4rW5gHvt6FPCwaDw1+TDko2YLsOh4q9b9szuoNgN1MjkoPX40dfK7yltLCyaWkTBhGKBOYNzc03kacwpXPCeSa3IEKorgeKRyMS64uO7y2tlN8MZgYJBSYI4U5AbGjDGP1OUWP/ui3RKGdReOiyYF9MKrMFYPe9vuYLlxoYyqajVwIs3V5eQxsZicPS8dDizKmHJRsQXbtMR4ligWzKiOwUZacStKTAUrpSlUwCqgOu/WtkvpEAaOcGWuDGJMc+ZSr3m5CQeWF+fM+wJzDA7cronxOB1x/Zdt3hVYK3wmrjztPdemzRkoDI3bHKa61cj0itBEMRKO473sfcMF4l65PGvNWQisuuQojpc1KLjSndXl5DDjAYlbysAd47EwOssZDCoZd2lf5qcJEoaJMueNNFRQ7v6BwFVzdH6sc1wkoq10Qj3rsRtN5MNYwplsK88MlhComyvClVbi64yESXGIiHvLj9BCmDz5jzBGGx7Xn26W/Wwqf5AoPSONKIsytVG3xuAin2glNE0WI3Bd0hUn3waVQdaPkOcpFX15O+R0gMMLr0lzlYbF47FIOjgpcgBFA6fMPEsaP378Uz153+EnwWI6whN9LEEeod0LqG2XypzrnXpD5hHG9IY9ByedaCMXiPv8lTcTuiHgwDCh0ThocEeU54pRFu6kR9dnOe6e/GNscPVQi9ekLTGufCucq9Exj7VoxZ8xye4OL5oLh5+TGOpaoE5dSxU3SNTaM3JzkopWXPZbPNGd3olYc/oVX4VPSwfrEp7/06YHfUMI35aovsZQf0xzlYW/wEO7blQP+VV3+K9rYKe1VHLdwHZ8heAjLH3LHQ+rkyqqNu7SVywtpyr/aVSaU3bXPOOXGwIU5f5K7WJqT8haDS2kO+5KutRhr3fdCHkbEY2M50FgatqCN747ibYGFJ0GAkxg7s7A7m6QT3yjvJfS3BNoYF+20wgnRvVn5+AX56WlwSbgsYe02GePG6z6g8yWs+8Z47EoODtOF0EC452z7jDOtYvEWBIQn1128u6DoJiPfPv3Q3+xpU1z8fLmK5AqSXztzZEdZ/I4n78stBpd47PsY1nqYPEQLuykenr93Ige5kwd3neGuP5qmBTdAgDtg7o4bX3V1temZg1MLd/YoQXYqvD2lPwiDgbZxulE3o9HauGgE/KaDjQ5+RcIlxWCJuFTz2dPA2uu+p/KwNh7ij93JwdB7rrY7MMsrv40IZ+4kn0yBEe3S/hRtT92m4VLmmamx32X7tu71dZ8DHhrDoDeP3v8N7Z7ugGxahoAhYAgYAkJAp7rGf0PbBkzjzaOtsOUZAoaAIWAIGAIggPHgR3JGhoAhYAgYAucbgUG24H/f4A1KwZiO3QAAAABJRU5ErkJggg==\n",
      "text/latex": [
       "$\\displaystyle \\left|{\\beta e^{i \\gamma_{2}} \\cos{\\left(\\theta_{1} \\right)} \\cos{\\left(\\theta_{2} \\right)} - e^{i \\gamma_{1}} \\sin{\\left(\\theta_{1} \\right)} \\sin{\\left(\\theta_{2} \\right)}}\\right|^{2}$"
      ],
      "text/plain": [
       "                                                 2\n",
       "│   ⅈ⋅γ₂                    ⅈ⋅γ₁                │ \n",
       "│β⋅ℯ    ⋅cos(θ₁)⋅cos(θ₂) - ℯ    ⋅sin(θ₁)⋅sin(θ₂)│ "
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "r=BS2*M1*Av*BS1*jones1\n",
    "pd1vv=jones0.T*r\n",
    "abs(pd1vv[0])**2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ0AAAAfCAYAAAA80BDFAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAMS0lEQVR4Ae2d23EXNxTGbY+f8kRMBYEOSKgA6ICQCiAdkOEJ3hjoAFIBJB1ABwl0QDoI8XMmM+T7ydL+tVpp7+vdxUczsrS6HB190jlHN9vHX758OTJnCBgChoAhYAjMhcCzZ8+uidYTT++GDx8q/fx0rkaMjiFgCBgChoAh4BF4IQPzc0BD8VeKf5C/eRISLTQEDAFDwBAwBGZC4JEMzd2I1gvFbyjtlhmdCBWLGgKGgCFgCMyCALucP3OUjrd4pyNryHngmcK/ckxbmiFgCFwgYLJiM2EPCGiestO5r/DmiX5w1rY1xwXU/a0xZfwYAhtEwGRlg4NylVjqsiHKvyU80OffgwsPCR7IVxc+JK7txOQva/Ng7RsCe0DAZGUPo/TV81i0IZqfvFxjl/O94ucgYXc6oLCi84MyOwdL0Z2d0QLBpfhfim6hG5Y8IwJLjd1SdGfseiuppfifStfX/0XhPflzvvGbMjpi6K78Y/nfYpRhVP6d/BfCJO9R/L2nuPryWPyy9VzCgRn0d+cMl+4hE0YmK90w9S1hspJHajQump83RJKrm1eK38IrzgnW5+OnT5/+o4Rv9ZF1ykNx3VRYPIJTHg2Q/zdl5XFYOLeduvhs/6myPB7gmd1LeX5j9dtQXyGGBkPEw4J78teU9rM8HeHBwXuFszrRpE+8K6cfr2clLmKiyRnnbYXFo0Tlgf113/YbfX/08V6Br88KY3b+ezEwopB4vVK4+DFqla8URtUxWUlA8TiarKyAi7Bv2BDSxArztOaUftznl0MxJr/XakYfIoLiZ9VVKU/F2X38Kv9jVLQreqYCr1UXpfOXQmew/DcGJrxke684Z4Q42n1JRCFvwukkxgLD95vSphgj6EDPXX4pnM2JL+g+UVikrTxnaBVihCnvfrFqCBO+7geFb+XPh9Rdo6zv51XDpVW+CuMwSVagKayZ36xEf5xhbpisAOoluq3LivgrbmROunBSZVZhlUGJy/uO/6TQKf4oDwOBMertRCMYGoQQYXBO6b+TF759yKofwxYrUowMZeEFfmvHcL5e70B0MFjstuBnbofRrPqYElebzqgqdDsUhfSTre4gTD1d2glGOm1qa9+7wUVjwYJn8tGuaBTlqzQ4qjNJVlSfeYTHWEx2omeyMhnFwQR2IytpzzqNTloh+eY4LGeQ2OEMOgqCrmixomfHwo6nTcFihNJjr3jXgDDFBkmfw514mEyj0OoD0c4eefl+c6yWMxSDlYRvh/bAdutuT7iA52qY+vEcJSuq+1G+eHoxZpKInsnKGODG19mTrNR6OdXo1IjxocnH8RjC8JDvgY56CAQTuLjD8Pk1o6a0eDeEwh5ytDeQzfHFPT4xrymxJ0pgJVsdDSoelNtZWrjnN+096Fl2lWKGy2DYR8vK4JZWqmBzIg/83nE5zXXLKzmODrgb4U6gsSpXGpOe+xVW3xzhYCjCSrx6k620ynm6KFUeHOCgHz84QNFyXFda6VPnyLeZPT5THrxAs2aUXMXoh/LZSdEH+A5KnbbhnX7wcIGQPjoD5utwV0X6W3mMG0YWd1v+D5VJjxpdZvSDhxCVQYnSj1QXPqCX7h5/8OXajFVMKo3THu02xjEU9G3TH/AIjv5UK+KozCdfoHF3pjJFXAPRQrhJXAq8Tkr2OHbJF3NhUVnp24m2MVXelZMVcPNjuJa87FpWskaHyS5QUdwoZZRvTlnx9topR4XhiTNKEaXPRKwpfU+LPC4uXZ5CFBSCFeig8PrsUOCr8Xd9RI9JAN+82HosnzUASkeg+SuoDF7l9O12VwrpB8aHy/vK6Ru+SUfpsuvgKaBrQyF9wUDn7qCUVTkMCDjkXLgjSI0StHFjjQ78glnWiWfGi77yp8edkVHoxl4huy52n/DAODPuFR+K85Td4aCwFVfVbXObw6WN2Yl5rfIlHMF+UVnpy3/XmCr/SskKuKnPa8vLrmXlNJ18AhTl8odPxwBUCiYtG75VpzJKiqP48ZVC94PEzoRXaLExQriyOxaltzkUX9UmBfWN0oR3jI0Cx3fW6CiPQftB5Xh6jaELLjUGn0NGEoLJXdWtDKTiKGaK0acabyRGDsVcovuTL4dSiqq43Q+GtHMs4kpRnPYQlJLDmGBcql2NvuETF/ChTM6gsmDA2GIoz+T74KpiDbdFXBpMTk0QTq3ypXzGaWlZgQfmKePFXKvtaJUWO5OVGI2L+NrysmtZORWG3ySYupWtT+Me4HmSf6RJSqfDEVma7RRykhiUOUqeozMcb+pZJaerepfZ8aOhfEUnfkjQWp025VHEvC+nfYQ8vHxrrRtlNnZaPg9s2hyCfl4ogDKojvMoI/6gd1++Ggef9oR8uWBM2KWU6IJXli9Pi3ZrBlrpGKCw6yGfdsJiRNELp3LB2HI0yS5zLK6TcYEjtU8/mbcscLrmRBEXaOFEg7mLgk4d/JKfu3sEk2pBklTskq/LkJWP4glfG/OET/epfuxaVuiE+jDbnPC01paXybLi+zGLDhHEqQ1xc6f0A6Pzb5wpZpziUoiiY7Byq3ZWPyVjwYA4GgqDQ2iHKvVQtxGKN3cc18gYloBCAnT66XZnosuruZwSGUZ5RGm1C9a4dOeHAsXF41D8B0kXRRs/2yZpMFqlRQTEQpl0XOOGGHfcrLgOwUVl4SEYiICnY6rwow0XV6U0H5TOvOEZe6fijttWeYehrw+P8bhS1GQlBiwTF3ZhbDtlRWXnnhNBFjYnL0NwEaxz6pCaDckMWS3pRF//1VIOHyhkDAXHOmGQQy7HW6yUcq50ydVYJecqX0aaeHcTRyErc35P4ljtYmweKR6U51KssBNI8YzbSnHFwGIMz6NC8BmUK8kYzbbf46E92s25sGvkUUDJhTJtfLOCn4LrZFzUPjsMjEDgt9SfkN6GSyizVNgmXyYrF6hvcU6EubWmvEzGRfDOqUNKNiQrOxidhpPgIowo31c+EwGJXazwqnTVow67oHQnwkBdrwpGEdVhxXjZDj7Dpb1rW3yw4uQ4Cf6XdGDhlHPciNo/j7+Je2zOFE3xxECWjveomjpoBGGp5fl2Pyox2294kCcf/qp7ukBEeWEucM49Bdc5cAls9Q2LuPQlMKacMGuTL5OVA6ibmxMaO+RgbXmZA5fZdMhhuPrFskZHVcOKlfNc4uGJ7JEXmD8Vpkob5YPi4XUaoMSOlXhQTi4dOvIYNQZwDccFKsIfO77TY8O0DOVRVmMd/b1dqIzRc+153sDtjuLncXl9pzsfJlB8VxAXJ86RVxvO3D+wU0rHlPZDvTuK8wtpGJbYUYY/1RNw64trTIP4ZFxSgj2+u3DpQWJUkaJ8iZrJygHSrc6JteVlMi4L6JDDqHXETnP5YohjCo7WnBKCwagcxgNjwR/ajB8F/K3v7JaT+vIYGQQqGDC+UZZrOJQ4fWKLGdpnJ8Y5J8obxcruzq3+9Y0xfSiPsiCdfI4dSX8uT3roC0qXP+RZukh+o7LUyznagAfogWXteXKugsrCC7tFFGjJMWaBv0YZ0aDP3ymDthnbMEavyKOCQuYEbbgn6aTJwSffweC04upqlH/Miku5mVpOKy61kjN+eCyz8qU8k5UD1pucExqjteVlVlzUn8k65DBkPWL8lWn+ZXVfr/Iv+pa1cnlcheEn+VtT8RGNG/Lv5K+VaPkyn0r5W0oXr3Phch9abX2biovq08bjtjYsLz//h+DCOMrPISuLz4kh/ZpadkZcJusQ8TLIhpz2sEtWZH4EXogkO4/i7qOrSa1Owi7D3bP47yOF6dEm90G0twc3GZcBnZyKC7u62rHngLataH8E9jQn+vdqesnJuKylQ06G9F1MXlP5tqeCQ8hd2bLCkaM97lAwHIOdr8cxGMdf/DUAtsco0c8xMV+OduLj0bjIpuKez9G49O3MHLiIBvedu8C1Ly5bLLenOXGZ+E3FRfXRPavokKE7He44wvn9ZWL8NbbFnQ93O40XYT06+0FlWAAQVk4TKd05ManStKr8RiOjcfGCRH+5q8F4sRrkrjH9XZo94rLR4boUtmxO5GEejYvIradDhp7HTT2LtPqHc25hz3nqIvcC0IX+HvE2XA5zZI/jtwTPNifyc2ILuIiHQXc6nf+uOm9gLdUQMAQMAUPAEHAvWxv/rroNl0F3Om2ELM8QMAQMAUPAEOhCAKPDL/2ZMwQMAUPAEDAExiAwyIb8D1vO2DcBKbq7AAAAAElFTkSuQmCC\n",
      "text/latex": [
       "$\\displaystyle \\left|{i \\beta e^{i \\gamma_{2}} \\sin{\\left(\\theta_{2} \\right)} \\cos{\\left(\\theta_{1} \\right)} + i e^{i \\gamma_{1}} \\sin{\\left(\\theta_{1} \\right)} \\cos{\\left(\\theta_{2} \\right)}}\\right|^{2}$"
      ],
      "text/plain": [
       "                                                     2\n",
       "│     ⅈ⋅γ₂                      ⅈ⋅γ₁                │ \n",
       "│ⅈ⋅β⋅ℯ    ⋅sin(θ₂)⋅cos(θ₁) + ⅈ⋅ℯ    ⋅sin(θ₁)⋅cos(θ₂)│ "
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd2vv=jones1.T*r\n",
    "abs(pd2vv[0])**2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY8AAAAfCAYAAAAfi0EQAAAACXBIWXMAAA7EAAAOxAGVKw4bAAALkElEQVR4Ae2dW3LdNhKGJZWe5slRVjD2Dpx4BXF24DgrsLMDp/xkv6WcHTizglx2kMwKEmcHyQ7i0fPUVHn+DwIQEAR4yEPyEJS6qyCAuDQaP7vRuFDS+YcPH86MDAFDwBAwBAyBHIHXr1/fU95Ln3/fx8+Uf32ZV7ZnQ8AQMAQMAUPAI/BGjuKrgIbSb5V+p/DgImRabAgYAoaAIWAIZAg8l8N4nOS9Ufq+8h6a80hQsaQhYAgYAoZABwF2Hb91cvzD+VZ3HvJcnKVdKf6zJJjlGQJ3AQGzg7vwlm/PGKWv7DyeKH5woR+cYW1BXMI82aJj69MQaAgBs4OGXsZdFuWQL1D5Q+HDnP0JOHFh/lQhXoiQeQqSIF+foh/rwxBoGQGzg5bfzp2TreoLpKd8acWu4xOlr0HG7jxAYWXywC/ey1p8Fxc0YbiWzGvxTUS35EIIrPWu1uK70LAPsllL/rl8ffuvFX+ucM0z4eTOQ50+Vnih8GOKJsIo/KzwgTgre54+7ymtsbyQvGz31iAwg/8uyLD4+zUJiztlB2HkpgMBiW7cKi6S674k5WrjrdIPCUpzavT+/NWrV/9Rxkd6GCTPhOOtvxQe+Mp4I7eFGWzsC1WXS3I+/fpWgd9O/Ci0V4zDwKFwgf65wj3lfaWAsFys/6J4URJPgOGbZcbx3aLMxUw8OR98pLh6RKcyJv+Pfd/f6/l3nx4V+fasBmbJr/aGxSjE51cS1k3ZwfwRjeOgce/GHsaNaJlareAiOXq+gDyNEn3tkPLPR/2SoCoygbNSipOg0uwG/qXwRYfr8MOVir9TW5ToT8XO8fhnHEX48uoXpTlfg+j325ukm5CZ5PCEXyjftQ9lR8TwAhh3AXRE+2oTyQbfl4qrvFXmHKZinCn13S/fVJkWCnzbd4p/UJiDh2FRwHelrFl2oPfMd/foC++MhdyPylt8cSW+i5Hk25s9LDb2IUat4yL5qhuLi6GBUeYH96XiOIH7Nkz0OJXRJB7BYbCDwQE4Uv5PlIVnH7MKx0HFCVHP9EfAaGaT+GFw7H7W+GAA5xfHmAurPp1zVOx2DIoZJ8dQkzD1fOknONu8q1HP6tewGIXU/ErCepYdSAKcBTaDTbKg6xzzzpdwFQ67sodVECgz3S0uB52HxssxU9xxJONnxzHpiIW24sUKhJUTO5ChiRJn0jlOUv3fFX5S/mIkftE5Lcb0htFT8S4eJSmfcXNcVZrwJztG3w/9ge3RpPaGxdHoTWvo39VRdqCe0t0s+rLWe5s2qOHau7OH4eEsVrpbXMY4jx5KUnyOnVD8Z73Cwxm0wwmg8NUVvy+f7JwOd79+DY9PvpNKO36pB1af8ahB6TDxX6UVJ6Tp7+mE+iepalhUYT7aDoRpqlssQKYcHVcFWqvAdKCM7N5xuSwP6yZXg0PBuX9gdcPRCBN+WBnH732VF0l1mQSZHLlYhziTTS/WmTA5BqutvGlz5vuctR0XD1b4jAG5w+RM38jOOLigJ2aMzgB9G+5yyP9BAePEWUKPFH5VnfwIzxUmP7jwj44hyT9TW+SAX76b+9TXSyeGtOmhNP3R79Bux7A4hOJC5f49r2oH6gObxLYOLrK8POgythAIXY47+aTOH75C7z5Fdao2FZgW4ubsARmT8W6BCSI0iQuCjaFB5yEGfNfrJjnF4dNaJjeUlsm1o7Sqw+REGZfZrkwxyoYRBT68qDErJXgV/6aK8g+S+mWS5i9C8oIi6dntdhQzDpwIl9SR9Izc5GNA7AL4PM05C8WMhcvp0h2NiiLhCMChROGz49y5wBs61nkgL5j1SPIaFj1U1ssQ3qvbgfrAEbhFmdIvFKoLGpVhq+g5f0rbOQvFyMjdCTtgTgLQP2wcm486qDSfzzsbUDyoR2pbo6bsASE1lq0xQYzmcEGosXQ5tqLAjitapVFcQpyY/ctgp8BXU6lTQUmP2UGgxLFP8TjTMwoOvysFvmTqrJyUlxIv5lPV4ZNfHFagfFJ/HwqyGAN6rLbR0SmNkVENGTqykZkQRlbj+6Wvh/xJE7cb4ZPbaLhp4Yg0/WEQJTIsSqiskKf3xztY2w5wBNgCTkORW3BUnYfK3Sfwqht3GcpDR6FgG9QpLYpY9LFgYrGD3Y2xKVXrUGv2gHBbY4IMLeKCXKPoUrX+UaopZWFg4egpr+Im1iwzTMpM1hxJQfz+AiuXfJXtCg/86E2i4oNTIgwZimNLnwpMqHy7TP8YdPhCxdUZ8aO28wGbIcLIrisVMPp4TEYdyQe/Jwrf8BzI53OPgUNOL0lDlTQGr6JcartrLDwO7F6h4CBZRdcwrmJxw+Lmp+f7bz0VcUvrJum4q07y0uQp7OCQLkR5/BjRuY7NKB9HEnYhlIPrrwodUr2wYOK4l53OMTbVmj3wvrfGBJxn46L3wViWso2iL0DQEuE8/lsqUB6r1dqkD/C54T5W3tTJWU3KhKKWSyblYmQAy8Tsdkviy1de1Yv6SdwnVvYvmlb5TixcdMfdjOqCMZhCYya3IUWEx26xkOzVf0jDwAp0CAvXRBijw6Mn4kI/pazW7CA429pCkDGEOrlNp+NDH6HF9GhDewjjbQ4TAJ6Ci6ovaRs1X4BYPbpQzv96uTcZHBuxyi9R7aKnt3IpNT5FnmR3CqKY1dIDhXP1i9N4rnQwhLVEYXU2NOHnuOIocWrReJVmxcdqsbcDqwhNf/TbI/HZNRYaEO8sOFLGx0Jg6HdiqljQ+ATUjB1orEF/wl+FKA0/1BnSWe5GjtWjpuyhEUx4D7NxEY8lbaPmC5C1RziPGqXGGutIgZh42ZXkOwMU8ONYMUmoDSv/UxNyhstp17fkYGXPVh351ySwcIaWdqL+o3MI+R6bKz3neIYqY2N4hEkgb7N3LHD6tSPEfKw8D2FRqr9kXlN24HWOxUpR59E/BcrRzXiHGQBRWZgHuCM4Vo+asodGMAHiJXDZzDaKzkPgsgL5TXE++aJIKBHnvvlExWowKJqSN9sv1eMMGOXcgriUzldTPOfHcXkdZGUCOpYYL5/1lgjn5frzsoHbZ0pflypPyOM4YQjn3WIhbDq7Mo0Tg3FfCVXwOYRFpdki2S3aAR99sFPL7RlZg858pjS/sIaDSIk6/PmcYDNj9Sjl0aI9bI0J+MzGZUvbuEzfcJLGCTDp8wcJ08vvv/Rc3P76QXBZjrKF78R5xtC3ICZjdhps60L/7Iw4I2TiwUi4D3ErMj3jFJ8psGMgn3K+fiL/GwXyw1gwoEcK8UsslaX0vR5oVyL6QAb4gWXn08hSg5F5vLMgX97k1mAh3Hgv7GRxEDUawqLWZpF8yYeja8kOziQP+v5PDRC9w66Dfb6ljIEr5pgUTN3nv+SJ0FGeg+MY1CPXovyjNXtgvFtjAlKL4qIxndY2+Ku6/CvaNCjvTfps6S4+Y/AQhn8oPBxTd6iOeDyB14E69w/VGWq/dhmyKSyBBeP8WeFeTWaVNY1FTe7bnq/3spQO7N4e0ne9IC6zbUOy9HxBKmueviwvFCx3AQTYgbETqO0GFugisuC+hP5apdlYaFUVVsHuXN4/nynOj09bx6LVd7S2XLN1YIKAe9KB2bhsZRsX+QuRIPeUN/QJW97EngsICEeOzDhnZtJbjTx/+qG/JmkuFn6MHLdwzMJvO7M9Z4J4nw7Y12sai1Teu5TWuzF7KLzwubh4nd/ENko7D+4AwhlnYbiWNQEB7kS4++h9xXKIh1cKdi2c3zMhskLhzqnzy17KQ3FOsbtRN7PoaCzU6zsFFjXEkYRFPu69YBHHcMcSR+vALbSH9NUfjYuYbGcbU8+58nMvex6+DxG+nEW+WAMn+MJ/Dd5r8DQshnVlDcxb42k6UNaBFnCRDJPuPEb/G9rUTVraEDAEDAFD4HYhoN1d79/QDo2wd+cxVNnKDAFDwBAwBAwBEMB58MtBRoaAIWAIGAJ3G4FJvuD/kTUNSs2RjtAAAAAASUVORK5CYII=\n",
      "text/latex": [
       "$\\displaystyle \\left|{\\beta e^{i \\gamma_{1}} \\sin{\\left(\\theta_{1} \\right)} \\sin{\\left(\\theta_{2} \\right)} - e^{i \\gamma_{2}} \\cos{\\left(\\theta_{1} \\right)} \\cos{\\left(\\theta_{2} \\right)}}\\right|^{2}$"
      ],
      "text/plain": [
       "                                                 2\n",
       "│   ⅈ⋅γ₁                    ⅈ⋅γ₂                │ \n",
       "│β⋅ℯ    ⋅sin(θ₁)⋅sin(θ₂) - ℯ    ⋅cos(θ₁)⋅cos(θ₂)│ "
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "r=BS2*M1*Ah*BS1*jones1\n",
    "pd1vh=jones0.T*r\n",
    "abs(pd1vh[0])**2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ0AAAAfCAYAAAA80BDFAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAMZElEQVR4Ae2d3XEdNRTHbY+feApOBSQdBFIBoYMQKkjowEyekjcmdBCoIEAH0AEkHUAHBD8zzIT/T5ZkrVbavXu1e+8uPmdGllYfR0d/6ejoy/bphw8fTowMAUPAEDAEDIG5EHj58uUd8Xru+d3z/lPFX53PVYnxMQQMAUPAEDAEPAKvZGC+Dmgo/Frht3L3z0Kk+YaAIWAIGAKGwEwIPJOheZTweqXwPcU9MKOToGJBQ8AQMAQMgVkQYJfze4nT6RrvdGQNOQ+8kP9nSWiLMwQMgWsETFdsJGwBAY1TdjqP5d8/0w/O2tZGXEA9XptQJo8hsEIETFdW2Cm3SaQxG6L0B8KD+fxTcOEhwRO5eOFD5LFJQn5zbBmsfkNgCwiYrmyhl/73MlZtiMYnL9fY5Xyq8BVI2J0OKByRfKfMLsFSfGcXtMJwKfmX4ltphkXPiMBSfbcU3xmbPshqKflb+fry38j/Qu6Kb9yqjI4EeiR3KfdTijKCyv0i9wE/S3uWfm8prLZcSl62nksQmMF/c2S4jHeZMDJdGYdp1xymK2Wk9sZF4/OeWHJ181rhBziFOcF6f/rixYu/FfGxPoqkNCau+/KrR3BKowLS/yKvHISFc9up68/hn8rL4wGe2X0nx2+sfhzKy8fQYIh4WPCF3B3FfS1HQ3hw8Kv8WUk8aRPvymnH97MyFzPx5IzzofzqUaLSwP6ur/uNvt/58E6eL88KY3b5dxJgj0yS9Vbh4vtoUL9yGFXGdCUDxeNounIEXIR9z4YQJ1EYpx1S/OkuvxyKMfm5UzL5EBMmflZdcfJUmN3HD3JfJlnHghfK8L3KMun8Kd8ZLP+NgQkv2X5VmDNCiHq/uw66iRxDgXX9UvGufEjbw4cXoLnLrz3KV4tINvg+l1/lrTRnaOVjhMnvfrGqyrSQ4Mu+lf+jXCsehRrmjfLtvG24DOpXBeEmXRHO/P4EY4oxziLxJ8W1LNxMVwTiIWntuiL5qhuZszGgVJhVWDQoaX7f8K/kx4nfp2MgMEY7k3gEQ4MSYjgcKf5n0sK391n1Y9jiRKpv6sOhAM0kfighuy3kmZswmrGNOXPV6YyqfLdDkU872epOwtTzpZ5gpPOq1va9GVzUFyx4mo92xaOqX7XOUZkmXRFfjAx6hd6i250j61q9tXjxMV2pgbNc/GZ0JYdg1OjkBbJvjsNKBokdzqSjIPiKF6svVmHseIYmWIxQ59hL+d/JVXdkyj+ZxC8atcmFhws8Ee/ikZfiaTfHaiVDMdmg+nqoD2zXTlvCBTyPhqnvz710RXKnO2zGVPM4lzzNPCqDc0tjotKERaI3i0ur0emhqcHH8RjK8LSXOB5BOYwHA7i6w/Dpk43aePXL5/D45Du3tOLn+mAlG487FA6T20WacUKY+p5MyH/wrIbLZMj31hVhnY4/FjdTjsEnC7pvARsTZeS2jst5qVlqFJMcRwec93In0FuVK45Bz/0KKyWOcDAUYSUe32QrLpLny6TKgwMI/umDAyZajutqK33KnPg6W48E2FHQBuQOkzp1Izvt4OECPm10SimfMtxVEf+jHAqLkYUeyv2mPPlRo0tMfvAQIhqUJP5EZZEDfvnu8TOfL50s0qJjYeqj3l4/hoK+btoDHoFoT9w9Jnn+8Bl69wHKU8U1MK34q8SlImtTtMdxTL8YC4vqiuRAb9G/wQXcUJ8q7dbpCp2vdtM/x9KXTetK0egITC5zGYxMyky+pcmKt9ducpQfnjgzKTKQGYidgex5kcYlv0uTzwSFYgU+THi7rLqQq/h3fRQ/SqrXDRj5dF4kfbvdlXzagfHh8j6SvpGbeCbdCzmeAjojI5+2YKBLd1BKioQBAYcShTuC3CjBG9rX6CAvmBVJMtNftJU/Pe6MjHzX9/LZdbH7RAb6mX6PcijMU3aHg/xBXFV2iFaHy5CwjWmD+iUcwX5RXVEdTJjoOK8bL+WKiyXFD/ap0m+VrtDvavOx9WXTunIOiCkJUCaX33wcBiBOMGm+NKwy0SgpzGDGxQnddxI7E16hpcYI5dpnx8LEF+sUjxN9Izf8LuRQ6s4qXXEp0WmfKQ9PrzF0gXJj8D4kZD6YPFLZaCAVZmImGzJ0ZCMyIZS4xvcrnw/5kyJu98PkMNoXaaEkTH0oSo0wJhiXuKvRN3JCAR/ylAwqCwaMLYYS7HfBVdl6tEZcekK2RginQf1SOv20tK6wwEAOjI08p+NFo6M00xUQ6tKx9WXTunIuLD/q4ukmn2AYuAf4Nks/0UCl0eGILE92E3IWGSZzJnmOziDe1LNKzlf1LnHkR2/yFR9kxtWUJ7KkTjkmYt6XUz9KHl7zxHwjgdpOC2yGiIn5qpKBiSAe55FH8sHvsVynH3w8/YMhTy+GFdUj8CrK5flQbwc3xWOAwq6HdCbDsBhR8JqULxhbjiZZOe+LazMuqps2PveiBSPL7q2GdxUXz+NEZRm7LCRyQl7SS3ePYBIXJFlBt3P0cSX9OoSujI2XKLLasWldkfyzjgnP79j60qwrdLBvS/McIla5DYnjpxTA6PyTJkgQp6DymejosNKqndVPzVjQIbmSo7RTJ3UVKZNkc8dx5dSdY1E8Jija6XZn4surudIksjPTfTOqXrCG8p0fgwKK/aC8YBwmwlDOZar8GBqkYXKuLSJgGfLk/ZpWh0zQrLhOwUV1V/9xlJOs/2MIF5e7Nh4Uz7jhGXvHWPer6MYov8PQl6fvYr/6nKYrXch6X8IujPlRXVHhucdE0IXV6csUXJR3zjmkY0N6HZZFnOn73ywufDIhYyg41gmdHNI43mJXUaLaJVdvlVwqfIg4ye4GjnxW5vyexKnqxdg8UzhMnkuJwk4gxzOtK8cVA4sxvAqZFGYlzWTX2/GFPJlPfdRbosCDRwE1CnmG5GYF34JrMy4Snv4Lxpi2sJgY+v2mIVwovyQN6ZfpyjXyaxwTQReOqS/NuMw8h9RsSFF/MDo9kkAoI5Nv2OqjICmlih3jVY4y7ILynQgddTdmTAIqw4rx0ISc4dLe1S05WHFynIT8SxJYuMk5rUT1X6XfhD02FwrmeOZZx77hEZSlk9fX+06RxXYjgxzpyBfv6QITpYWxwDl3C65z4MLCoXbsGURO/Souaaa5w8JsSL9MV24AX92YUN+hB8fWlzlwuUF5t9BsulI0OpIhrFg5zyUcnsieeIX5XX4+aTP5MPHwOg1QUmLFGSYnFw8fOYwaHXgM4rIe5U+J7/zYMM9DfjpgX6K9DyuFMXquPi8buH2u8FUl/67RHHkN4cz9AzuCvE+pP5T7XGF+IQ3DkhJ5+FM9AbddcU15EG7GRTJ0doTiiRFK71CoJ6UxXNK8c4ar+qVKTFdukF7rmDi2vjTjcgPxzqHZdOW8VKWUl+MbjtbcJIQyJ/kwHhgL/tDmpY9nF/OXvotbTsrLYWRQqGDA+D7K/YlkYBKnTRzHyHNEGzj/ZZJiYmV351b/+saYPpVjsiCedI4dif9WjvjQFiZd/pBn7SL5jfJSrkTUgQzwA8vO8+RSgR3j6LMgX6+I6qPNnyiBuunb0EevSaOAfMYEA889syVOhJx8B4MziKsrUf8xKy6SiT5iF43MNRrEpVaoNd5jWdQvpZmu3AC8yjGhPjq2vsyKyw3cg6H5dIW/Ms2/rN7VKf+rXfNavjKuwvAPuQet+IjHY3gN8VH6vbE8Q+UPmYaccnPgQpt/kbtTk19pTbioPNhf1vhbfHnsT8VFGG9mTExtW0v+GXFpnkMkyyQbcj5o2yxxKQReiTE7j+ruY8aKuQ+ivi1QMy5ahYbdl7t/8t8nrE4zAFpxYVfXeuyZiWSfBQS2NCYK4i8W1YzLBMladaVT1Vnna+RDintHWYaeCo5wsGQQEI4c7XGHwgS5GHn+1JMejy5WXyvjVlx8ezke5FiQv5LAERsK8z6VzedrwkU8uO/cBK5p27cW9hjvrSu+rw8yJg6JbSsuu8rq8WvSlbyuqTsd7jjC+X3Oy76nIcCdD3c7vRdhY2z8QGCXxDkrA4JVD3dq+e+MoGyH2E2pmtlob1wkwVs5Fkb4kYRLjsEWcYntuYUBGxPlTt8bl6POIVPP41rOIa1s95xb2HOvsMi9AHzhv0XMDZfuONliH84ts42J8phYAy6SYdKdzui/qy4bWIs1BAwBQ8AQMATcdUHv31UP4TLpTmeIkaUZAoaAIWAIGAJjCGB0+KU/I0PAEDAEDAFDYB8EJtmQ/wBQndg3NYEseAAAAABJRU5ErkJggg==\n",
      "text/latex": [
       "$\\displaystyle \\left|{i \\beta e^{i \\gamma_{1}} \\sin{\\left(\\theta_{1} \\right)} \\cos{\\left(\\theta_{2} \\right)} + i e^{i \\gamma_{2}} \\sin{\\left(\\theta_{2} \\right)} \\cos{\\left(\\theta_{1} \\right)}}\\right|^{2}$"
      ],
      "text/plain": [
       "                                                     2\n",
       "│     ⅈ⋅γ₁                      ⅈ⋅γ₂                │ \n",
       "│ⅈ⋅β⋅ℯ    ⋅sin(θ₁)⋅cos(θ₂) + ⅈ⋅ℯ    ⋅sin(θ₂)⋅cos(θ₁)│ "
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd2vh=jones1.T*r\n",
    "abs(pd2vh[0])**2"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Como se puede ver lo que importa es si el objeto esta en el camino en el que se envia el photon o no, es decir usar el otro input channel daria los mismos resultados (no mas que $P_{D_{1}}\\xrightarrow{}P_{D_{2}}$)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAAaCAYAAACzU/tpAAAACXBIWXMAAA7EAAAOxAGVKw4bAAALZklEQVR4Ae2c7ZXVNhCGL5wtgCwVhHRAkgoCHZCkAkgH7OEX/ONAByQVJKEDkgr46CB0kGQ7IO8jNI7sK8nytXXlu+s5x6tvafSOZkaSfffGp0+fdi3o2bNntzTuucKPLcbfxtwQWCsCm26sVTLXh68zLcLfNd27wZTvH8lYP9GYf+t5GYy9RTcENgR2u003tlVwNARk7x9rsJ+CAS92T58+fc+pYXuOg4HwvlMD61r91uB12GdN3mv2PZzHlp6uQ7XkU6vfmjKuyXNp36r3iOdm4CW2aGUEvGcOT2dLjnjH979kn9X7qowJ/J8kLtWBX8EAlWV/UnKvjMVkPbgh7/BeTH1dsk5U747qceTgCugr3+ZC+Zc+Phqo7j1Vwjh+q/j31sD3/Uppyv9Q+n5Q9kjpny19iqH4fyC+mfNFin+VcaS77ct/VfpDqm4s37e/VHgSWInPUUyYp5/XlcdF81yNbogXdP29HvR78fWkPkdlrzqr0YeaeJRgcSw9EC+PGKv4xKAGGPMHClkoL/XgIFg4v+gpIrXhhfNd2iukL9JGOAXed+AQPqiM9E4h4570C2o/zycKc07hjeaJUafOcz1gMYnUFlx/8uNNanvsyp7HLCbwpHrXAhePx5p0A8eAfhZtGpFVKZXIfoVyr4JHCRbgemw8ihyDZ/5HhcMXxRjsKVcj56r/s/pht/BRoTtp+DQGjTJOCxhHV6bwHnkKdwqJ41Ae63lFmvwToBfi0Tm6GK+aB+U7hW5nppC5cxSegq11zTiuP8tYaZjFBJ41/+uEyyzd8HixZt7owaDPIvWBzn2hMHwpOavPoHFW9hpzdXKviEcWCzBrgUeRYxBvXOXEdrtcBRVfd6gPcwYsts5QKv81ZcoLiasUjjXmICj73dfFQcEPu8lToB/Ed/Q4rnyMP0dmpwyDydwZpEeTfhzGm20cRgebVyGJCd1eN1w031m64fFiLU1eMykxqs9Q91LVDslPyn7Ncq+ERxILgG2FR6lj2BO+GGbXz4794V5hJkPtMFi043SQ2xHjKIb3m+GxFgWotXDV9TLkcRo6vbDzJ0pgFNypiALFzaizizyEGO+HQxoeo00BJrBxHXE5WDeEKdevr48hvzljFMj+2si9AItmenA2JmQxjxHnegdDzC4fY2y7kq+V3zPOSrO4ES4vqCFeUmPgrR79sYi5T+fkED2q+vLeaUR5oYFlh929vFY8SmoDP9S18an3VvmdEgV1/qJQBM+cTkJjjRODd/qhT4jrtdBZfc7t/+WdSddPWOTHxcEOT2Pf+HrhfMOmY3HGY9zoKYXGfuxWuCQxCXhrggvj1yKPeRXdmMqzeEmuZ5Wh37zjIkT3nZ75NrxTJP83Pawf5AR9qwe9Gl43u8LgT1L2aoteNZG7n1tUv1VWC48kFuDVEo+zQGCpKD94c4ZLob0QxWDhJACsM94qB1jyv1fc5StkAaIMZvwwWhjU1PWJih0o9B29KlJbxsDZdGPTZkgqpw9ekD9U3DkChfCI0WeXjoOCP/fSmzzFHSnOXa17GaiQBftCIYLsSOmoU+sqfI5g5OE3Ru4LABUMHQc8QR0/n5PFf3FwzDNK4rs1LjlM4LkJLlGwFsoU5kfRjRJ2xUt2PaucdcemD93pSGn0jXzW17kee1m+Ux5rli8cY9fCKuooJ/smchfPrfDIYQFgTfBg4DP+lJIA7HagirNb4HHGUmkz5LxEDg02CtEZeJVdKj2601cd2r3T0yM/rjuBKP5YT26HgsHHAXSnA6VZBBB8QNSJLWYcGQsdo40SfKP4LT3WTllJg0+ZEeP9Y4lB+KNP83VOWPRACU5UhzoGxkMeKWqNSw4TeG6FSwqvWfmS41F0YwKTGKSS9Zxat6zLe5pXp8eKs8mCBfS2sxNkDCgn+1Zyb4VHDgtga4VH3jFI0DBuV0ID+brdLIvAyHbFGE9OAxDfnrPzHu6IXeHIH04qvQWmNDsYdiY4BJqzQKOOQeXwTt1eufJxEnZ6oBylfaunR6pnC53TDY4IJflXIXPB0eFMen33Ovg/gVO5/D/ZizF+d1SnRH3CN47hOWkjn897Axzv2PUVuNDPHvl+WuOSwwSei3Dxc+E0Cpkj5HSYwjuJy+cuOvz/VDqKn9UbhN0JeZBvyWPoBpihj2DLRqN3XWqMEKqMNXfoerau9jZtvmAMt5zsi+Tu58A4s/XB99UKjxwWsFaEh2QJFovqwRmjZwhPmjLqMH0ZtGVRlhrLoFkyihL3SACMGcSwvhmKlGOjrtUJ5xH2QZx5QowN+Bhtd1oSP7xAL7lOUpM+eWGS2Z2mfA0WO9Q5RdWFB/CFWARjdK4KqTnZnE8eF82R670Of8UxwGwe7MeXivYoh4urqD7Abco66w2QSBxDNz5obJ6SzQpsLrqe6XAOCXdb18fWB2P7ZPHQBBbXg5uGSiJk185ii1HsxcnezjvWsCRP416U1MvUMceSMhI0tTq2KGPdcRXljCk86flKzw1VxCA9UtwcR6wteezMcv0P8WXeOBwMlCPFOb2g8MavL0kGjMe4MbI+WuIyhgl8j+KiOuBvDpM2OOzc7z9yuNC+Jq1GN4TZnPU8F6Mx2Y/KfWF92DXEYwwLsB7FQ3UW14MxxxAqXbcgBCTGkNNEaLwxOLe7SkFE9dllH5U0JoYVUOFzj+BJD+XU671UprLKbO7cxzNfexFEMeXs6LmSivZPHU/g4hTRMgjVnnF7pDxwYlcb4tqrU5igD3MAvSZ+3Na4RDGB0Ym44JxTVxq9eftEEpdY5QXzVqUbmtec9TwXlqjsJ8p9Kg9jcm+FRxQLJjcRj8X1IOkYxBi7q3cKhwYRg4mx5F6ViRmxWzNj6vLoQw/H+6HXsza1Q16OsYMczgFejafvFOdHJiyOkKjDv/6wqzTubcEkJNJWHuaHccbhU74Y4Vhcn75vxvxO8ctY5Ql5HIttfrFmrXHJYQK/RbgIp97JSu1QEPe1GZ1EaAyXSJNFstaoG6XrebjmAQRDeyjlZF8k9wMGLpF7CzxyWDDNIjxq6MFZBmSMPEb9XAOHL5P/VnrvGsIzhyNACez3AKRR1iaksTESX2pw7uCYi/H1ijKYUsg1DQvHfelEnogdPmkz+hhqTggc2RQ4uq2/9Bs6R1/UC35VCkcao4fKpA/GA1Ou7sb6i/UzzEN2SdwZQ09LXHKYMJfJuGg+OHZOXMgyRVlcUo3m5os3HNiadCO7nj2WTzRvdxpWmvWLTFin5IP1pc9/7vNtvWFg+WeR3RdLKg8pJ/vJcg87zsTH5N4KjxwWTGcyHsJ9ET1I/ndVDYDBmnulkZHV9SkSjjiksS9XRgFRPxg+5LLnmK2xylBevgRL1rG6LUPxtwgmzMHPGccPxij5Hvk6q8dlj/ErmCFZLCJ79XPy+rAUFiwT9YXuz9ID9eFuV86u4Lpb45ReiCl2VLarqskjzpzx1k6LYOKVgdPdfSbs04TDk9ep4LJ2uS3B3yKyL2Rk7XJfBAu/7hfTg5sxcDXILeXnPmeMNdvyEggIT66heNeBR69Gvn/G6T51rTbYzI6XwMTPlx0SV4P8CpdjNIbgn5A9X+8kcAn5vqpxyWPTBy/cJbDw63tRPUidGLhbtPv1q7o+jz0v7ly5q3U72ymDe8Fz2uCuFAPHLoN3PcNv1lkcxziVaJhF6GBM/Oj8ZoFNDGFHwmWIwanh0s3lCkcOlv0V1IeDsfDrY3E9SL5juMILstnU/ILmM9mhQZ/Nk/rkA4HYv/aY3XfNDmpiAt+niktNzNfSd03Zn5rca2IxRQ/Eh3vHsDmGtWjJxseGwIbAhkBjBMwxRN8xNOZtG35DYENgQ2BDoCEC7h2DvIR93w8rS31L33Ba29AbAhsCGwIbAiUIyP5zDW3v5c4Vv/gPaL18pMo1CsMAAAAASUVORK5CYII=\n",
      "text/latex": [
       "$\\displaystyle \\left[\\begin{matrix}\\beta e^{i \\gamma_{2}} \\cos{\\left(\\theta_{1} \\right)} \\cos{\\left(\\theta_{2} \\right)} - e^{i \\gamma_{1}} \\sin{\\left(\\theta_{1} \\right)} \\sin{\\left(\\theta_{2} \\right)}\\end{matrix}\\right]$"
      ],
      "text/plain": [
       "⎡   ⅈ⋅γ₂                    ⅈ⋅γ₁                ⎤\n",
       "⎣β⋅ℯ    ⋅cos(θ₁)⋅cos(θ₂) - ℯ    ⋅sin(θ₁)⋅sin(θ₂)⎦"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd1vv"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAAaCAYAAACzU/tpAAAACXBIWXMAAA7EAAAOxAGVKw4bAAALYUlEQVR4Ae2c25XUOBCGG04HwA4R7JABlwiADLhEwJABc+YJ3jiQwbARsJABuxFwyWDJAHYyYP9PI3lltyTLbavl7nGd45GsS6n0l6p0c8+1X79+rWrQq1evbqjdI4Xfa7S/tLkgMFcEFtuYq2aujlxrDcIP6u5tr8sPd+Ssz9TmDz1vvbaX6ILAgsBqtdjGMgp2hoD8/Qs19txr8HT18uXLr+walmc3GAjv4xJYl+JbQtYuz5Kyl+Td7cfyPtyGSumnFN+SOi4pcy5vlTvhue7NEku0MAJ2ZvZ3Z1O2eGz5T8mzOK/CmCD/XuJSHPgZNFBY93ul98JYDLaDa5odvkqoOznjROWOVY4tB0dAt2ydU6Vf2HhvoLIPVAjneE/xx66C5X2ud/L/0vtDL+9E7+/c+z6Gkv+R5KbPpzH5lceW7qbNf6/3b7GyoXRb/0LhXmAlOXsxoZ+2XwePi/o5G9uQLNj6Vz3Y9+TjSTx7da8ys7GHknjkYLErO5AsJ7SVvWNQBZz5I4UMlLd6mCAYOH/oySLV4cL5NvUVwot3R0wK3HcwIXxTHu8rhbS71xfUtp9nClOTwif1E6dOmdd6wGIQqS64PrftDaq768JWxiQmyKRyVwIXi8ecbIOJAfvMWjSiq1zK0f0M9V4EjxwswHXXeGRNDFb4pwq7F8U47CFHI0cq/058WC18V2h2GvYdh0YeuwWco8lT+IA0hYYUZ4v4SY8/qbjsuYZvJJiZ6EICqi/krxSalZlC+k4/h2DrWNOO4ecSZhomMUFm9f8q4TLKNoQVdsJi64Wec97H6F31sbnfFPqXkmNY+nWTulebs9N7QTySWABaDTyyJgbJxlFOaLXLUVD2cYd4uMmAwdY4SqV/JE9pPnGUwrbGTRArveMoeY79gnsQfyLZg9tx2ye2zMYYOn0Z3E/bDu3NfeKMYgIGVw0X9XeUbQiyD+KBHbF4w1bZaY0i8WpsbxSjzcpR3atN7HuW9lAIjygWwFYLj9yJYUO1EphVP6uSZxuZiQTVw2FRj90BgyBGTBSt802V54jpY6zCHNMlr9kdJWQ7Ux5Owd8VOad+lKiXygK7J6kCNfMyMEG8q4jL1rYhvPwjHxYUpZw6utmaMnR/ZfSegQU4V8Fj3adhCY8T53iHwXauhwHnVrJ3lN4agHpncNMZLqghLqlx8K4c/HDwnKezcwhuVW1+9m5EfIIkPsjDaty1T7nPSm8mGK/MP2SKkJkVmO+smcSQHT7whDhe8w3yMrX9lzuTho+fZdtl4ujuxu7act1dlF89Fac92g3uUqho266FSxQTT7YquNB+KbKYF7EN8fbHCnptPuwI9Uflo+NZedg3d1yE2L7hZetwp0j6n3poBz1B9/RgV93jZpPp/YnqXnWxqyp6t30L2rfySuERxQK8auKx9hQWi/KDN+O4FLoLUQYhkwSANc5b+QBL+mPFTbpCBiDG4JwfTguHGtsuKtuAAu9R22G1AQ8uyJ8pbiYChciI02eVzgSFfObSmzTFDSnOPYa5DFTIgH2jEEU2pPfgpNYUuIzg5MEkROYLAGV0Jw5kghp5Ll+z/zLB0c8gSe7auKQwQeYquATBmihRmO/ENtQOY42FWGOX3S4oLzmelc+4Y9GH7TRkeZLO+DrS4y7LV0pjzPKFY+hYuOGhSEr3VfQumWvhkcICzKrgQcNr/uSSAGxWoIqzWuAxzlLvzpFziewPSgyicfDKu9B7cjWjfIh6X0xs+z84fCaAZnegdwYBhBwQZUKDmYmMgY7TxgjuKn7Dyq9XQzGH7/IJae+nn+DFn9o4X+d4yWbVxI5q24mB9tBHjGrjksIEmWvhEsNrVLr0uBPbUDvYo9mdK/5CT2z1jkPKGc+xccu4fCD+jR0rziJLycZuGz9BQodSuq+l91p4pLAAtlp4pCcGKRrB3ZFQR79mNYvzduScJM6T3QDEt+esvLsrYpPZ84edSmuA6Z1VCW3iqHGmrSMhn5/ykJ3yLeNQOpOE2z2Qj9F+1tMilXMDnd0NxoaR/KuQvjDRuYu+Vr3AC7JeBNJJov1mq06C+CP3Iz2veXdk07k3YOLtO77CcOGzQZZPbVxSmCBzFi62L2e2k24iZHcYwzuKiwPK8vxb70H8XLlO2OyQO+nudRe2weoe3JgQaJe+tsY+iZDyGXPbjudLJvFFWx9uKd1n6d32gXZG24PlVQuPFBaIloWHdAkWk9rBmtYTxEwac+oIfeHVxWHnOkuvWjTKwG6RAPimBJ7ggG8V/n/FHJvYKO6cid+PDhujHNJwxoCP0za7JcnDBXrOcZKqtMkqk8RmN2VLMNihZlJUWbAGX4hB0EepAef6vPe4CASO9xr8FccB4yDdjy+7OKVwMWXFg7HQN/F2+fa978I2hso86XjuA6AvX7i7cb1re3Ci7S0e6sDkdnDdoRIJWbXjiEMUujjZWHmHKuakqd3TnHKJMm5iiTkJqroyblCG2HEUZZwpMum5peeaCuKQThTHaaeIlVmKfxdf+s2Eg4MypDi7FyZDJ6/NiQa0R7shcjxq4tKHCXL34qIy4O8mTOowYad+/5HChfolaTa2IczGjOexGPXpvlfvE9vDqiIefViAdS8eKjO5HfRNDL7RNQNCQOIM2U34zhuHc7Mp5EVUnlX2Tklt4lgBFTk3CJn0kE+51qUyhZXn+s55PP11F0Fkk8+KniOpIH/KWAIXY4gugVD1abdFSgMnVrU+rq0ymS/wcBNAq4pttzYuQUwQdCAuTM5fWh1Mv0RxSVcbnTsr21BvxoznsWAEdT9Q70Nl6NN7LTyCWNC5gXhMbgfRiUGCsbr6orDrEHGYOEvOVemYI1ZrzpmaNHjoYXvfnfVcndIhl2OsILt9QFYn033F+ZEJg8MnyvCvP9xRGncaYOIT7y7fT/fjtMOnfCFiYjE8LW/avK/4RajwgDS2xa5/oWq1cUlhgrxZuAin1s5K9TAQ87UZTALUh0ugyiRJc7SN3PHcHfMAgqPdllK6z9L7Fg3n6L0GHiks6GYWHiXsYJ0AGSePUz9Sw/5l8g+9bxxDWOGYCDAC93sA3jHWKqS2cRK/q3HO4OiLk+ucPIRSyDENA8d8zUGaiBU+787p46jZIbBlU2Dopv7C158cbVYreK83JtIQPVMiPGgPTDm66+MX4tNNQ3dR3GlDT01cUpjQl8G4qD9M7Oy40GWMkrjEKo1Nl2xMYHOyjeR4tlieqd9mN6x3xi86YZySDtYXNv21TXfjDQfLP4tsvlhSvk8p3Q/Wu884Ee/Tey08UljQncF4CPdJ7CD631XVAA5r7JFGQldXJ0s4MiH1fbnSC4j44PjQy8bE7CorD+PlS7BoGVe2Zij5JsGEPtg+M/GDMUa+QbbM7HHZEPwAE6SLSXQvPntvD1NhwTARL2x/lB2IhzldWR/guJtjl95IKFZUblVVUkYmc9qbO02CiTUGdncP6bB9J+zuvPYFl7nrbQr5JtF9piBz1/skWNhxP5kdXA+Bq0ZuKD31OWOo2pIWQUB4cgzFXQczejGy/Gmn+dS1WGMjGU+Bie0vKySOBvkVLttoHMFPXzxbbi9w8eU+1Lj0sdiDVe4UWNjxPakdxHYMnC268/VDHZ+77hdnrpzVmpXtkMat4tltcFaKg2OVwV1P9/ccDI5d7ErUzCS0NSa2dX6zwCKGsCHh0sVg33Bp+nLAka11f4D2sDUWdnxMbgfRO4YDHpDVumYHNJ/Jdh36aJnEkw8EQv/aYzTvkgxKYoLc+4pLScznwruk7vdN7yWxGGIHksPcMSwTw1ysZJFjQWBBYEGgMgJuYgjeMVSWbWl+QWBBYEFgQaAiAuaOQbOE+74fUab6lr5it5amFwQWBBYEFgRyEJD/5xja3csdKX76H+BQfKRCIs4tAAAAAElFTkSuQmCC\n",
      "text/latex": [
       "$\\displaystyle \\left[\\begin{matrix}\\beta e^{i \\gamma_{1}} \\cos{\\left(\\theta_{1} \\right)} \\cos{\\left(\\theta_{2} \\right)} - e^{i \\gamma_{2}} \\sin{\\left(\\theta_{1} \\right)} \\sin{\\left(\\theta_{2} \\right)}\\end{matrix}\\right]$"
      ],
      "text/plain": [
       "⎡   ⅈ⋅γ₁                    ⅈ⋅γ₂                ⎤\n",
       "⎣β⋅ℯ    ⋅cos(θ₁)⋅cos(θ₂) - ℯ    ⋅sin(θ₁)⋅sin(θ₂)⎦"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd2hh"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Checamos lo mismo para cuando el camino en el que se envia el foton no es el mismo en el que esta el objeto:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZYAAAAaCAYAAACU/XqBAAAACXBIWXMAAA7EAAAOxAGVKw4bAAALJklEQVR4Ae2c7bXUNhCGF84tgEAFIR0EqIDQAUkqADqAwy/4x4EOSCogpIMkFSTQQdJBktsBeR+hcWSt5LXXktfeK53jK1kfI82rGY1G8t5rnz592m01vHz58obGflPxX1vloY27ITAXgaYHcxFs7UsjcCGhfC+iXweEH2xooX6ucf+t500w/pZsCFw1BJoeXLUZXxG/shdPNZwnwZCe7V68ePEBr6U968VAc3S71vzUpL21MW8Ri1oYr51urbmqRXcpPGuNfyxd1XvMcz2wMi25QgT8biD0KEuP8rbvozTdKvQq47EpLKoAvAGiTQbSk7QmXK7JunzQgO6kh1o2V/3cFkVcJo6vvvLUnyn/0qdHRar/jSqy2N5T+ltr5Om/1Tvlv+r9QVD2WO8/2HuJ2Pf3QbTgoShtxieaDxXB4zPeU0FluKG3fNk7vX9M1RvK8zQuFc/iQe2Z303jUQqLIbxLlWmsq9CDUvyMoSOeq+vElmTAMFsLLhrHY8a0mMeiDjEEDxWzCL/Rg4FhEfpRz+igdlzYfw0NxdDj3QJGhTsjDMpHlfG+U0zfNS74WUjpv7hh9nw9VzxkVH5R3xgE6rzSA++Tg9qD5RPf5+T2QYPN41EQiwCW8kk/V2vRg/IMJih6nqvrxFZkwCBaIy6LGBbP+Pd+wgwPYhb7qcc8N9XmB9Fi5/KXYuft+HcWR8rwVlhsXZnib8hT7ILSHHn8oic0SlY8OvY0v1AcXlyNbn+g4muVO8OYqqc+Kd8pdl6GYniFr6l4Gnn6cjQtY2qsvsH4HPAogQUy53ZvU3EcWX+WHmhsjI+N2VM9b3kf2e8pqy2pE7NlYEGgVofLIoZFAHMMldp5c4w16ehGdMyYsJh3C6/yf6YsmkyOhlBuMzA7vbPw8tyO6h71Knod7aMI5Bt9J9rJoynPA0dgKUNwFF++L/qca2w3j0chLMBxFpZ50XAbill6INrvxSc6g7eKbuL9rj0sphOFZGApPFeHy1KGZQ9gTRweB7ukR3uFBzLUFoWlLd7J0A4dQ9O7/1B9jsh+PtDFSYs1PueNDQziucpYWEIvzBYxdrLHBvD67tjGtdqdCI9VYhFiLFyO1gPRCY9v2YxchrTXlm4ykJ6RteJykR5uuVwxjgHgaArhxcNAgG1XfUf5ewKtPBSGxZNLfgIX/RgIqwtNDAT3C3guyaMoXz7JIxKtLqg9Rou+6JcxETjSY9zwwJ0GMfy5jwh8G+6NyP9JD14FhoJwT8/vqsMucShwR9QZjbCi2jIO6MUe4F1fL/bawuaH0vRJ30Oe0lXBYxCLQ0CWKPdzXUUPRDuUE2S0+wgmNXY/FuqZDlINWe42aUGdPz0N9BbPqJNlpbM65dvkolPoxEEZCHhuuAQzdxGkayX5waVbBBXb5TJCjZFh8e0t/KrDwkXZt0q7MsUIIwpmiykTzgKfOw5SsTsugP5RLr5os4C/VoxAd0HvzogphgcMDB8gdEHvjJl8lAvvwS5Yd8qDD77CSx3bqagLGAkwSAU7t++U1VeCNiFcMD7njP/LmMF/L2jMVw2PLBZ74FTIEN6L6IH6Qc7YtPX0MGRJZegRcv5IaWdIFDM+jAaeM5s85M99OEOe0i4ozV2m0wHFgzJkbTLxKXRiUAbET8MlM1kXmfwq2ZqIbiesNLsfnm7h9hOFIeASPhR0hLgzECpjdzC4w1I5gXZ/uNT0PwjyXfV1w/dnFOIF/x8riGKUiwvSbpxKo4BUY1wdFmREAQXM0f3e1+XrmLDZQ73gwXVKHRaOTNMnypIKVw2PISxS+BTL0xwyB9X1QP2gf+4kQOmnenKeNAYDA9J5J3pHRgnoIoE6qQ0Tm0E2U2yE2GiN0SlV2wun0IlDMtBw2ZumzxkXmfydBIGJ/E2PCVCuapjfeRlkehp2nBXWI+0W3ijTFm0Wc7wRAr/RYNcT79Bd4YE/eEu9BVzv7KxY2BFyFueeO2/06E8PgvUvacUoul12WrVDcc6oHcKUsV1miDP+7uiNOhof9DAsr3i34PO5M8FQh2fqViWOmZPk2NR+03h4LJ57hs14sgPP4ZzFIgRN7ZFZ5CkOzOFO5aljWjYY3YYjbqj3JfQADwRZwqAocvq4Z1hUhjxQr1emfIyMMzRKUw6mv+vpBZXZZooTBozYsToFnrm5ov8aOpGVAfFxNrh4XorqxpBhYRLHLEY9QYpe2OXmDALCEAsKCjp18Y667L0iGL0gED8qg6enKL1K/7/AP4CzaDsPS+35YCC1WPzfqlLKCwDUO+/Nd2UX7p0RVV3wtQUPJRgThpSX9pvFQ2PnWLObN6VZvFlc7Ye6SvbCISxc5ZBm2Fr5yAyff4+Rs7Ap6SX0YKxumxHObRAZr9WJ9ZkyC8gjoagMCV+T7Ro6MSQDxvPmcdGcFNeN626q6/3BY2ART4XcZdzerifVeEye+rY7mTHVe3XU1gkONPR8peeaKrAwPVbalKTXpuALuzpTmBTZGFP4xOB1iq00O0UWtT3jmiLo8+iTvveCaG0aDzHEvJmhhT82CkO/+8liQeMFwir0QHya/OQMMFBYnSGZ5ShtjgydQieGZMB43jwumr/iulHbsISKjAC6IAFjYcabiRd+JuuWqxT9URt2gEsGxmgX5a5fjQGPAPefsdcM4OCUMOxE/V+G76Q9LuysYizjqmPeoWMKE9ffOh5sCnJHkzGvvA9hkapfMm81euBljo1MUuaRPz2UI5vdfamBoTJbA7iPmCNDp9CJrAycGS7FdaOaYRHwWPs/FMeLM4KGkHEfEy9i7CJNEJX8fH+gehxbILxLB+5g4CMMvMfHe3Ed6iOUxwZ4vZdpjGFz/fmxgdl9pS8z9adkc0wxhPNm8RA+PY9OfKJM7oumDECHsMg0K5K9Nj3gPgjvLtZlxmnycl9pfqiH8QgDdfgXTqYzY2UopEH6FDpxSAbOApcaunERz17BdwwEBuGmBh5exP+t96T76Bnk4h5h5FM/Au8sAksHFmo8FNxE6xtvivNIFiQUiPsXt5PTO8bykR48DfIp5yst8l/pId/4QLn455K5C9x3qku7VKAPxgA9cOS4MTbQqXZj8pgzG2Nc/2zwEF7MDR4wC0cuDGGRa1MkX+PDCK5FD3YaC/L+pZhD7tBp0823lMG0Yo5ewdN9ZUaeAjLKuxmVQRlyLfJ/TqETgzIgvs4OF/FURDeq/XdjDRAhLHE8kxe1My4Rdihv7yu7Y9gVHRZQ5iJpzI2mylkE+PpusJ7VXzrWuErhAZ8sjmCb9PJKYCEa4H7s5f3S8G6iv4IycFAnSsjAUqAWxGW2bmgszqu9WIr51s9kBPDa8B5yHsRkggcasAmgz7WG2Xj4xYIdtLsL8O87xbHHVwILjFbScK0V4A2Ma7YMTOCxhAxM6G5W1dm4eF0ophvXZ7GTaaxB3lDR0Gd4mZYt2xAQhhzDseNlF1E1+D7oq/tcuWqHRxCfi4fnEU+F4xt+CY7Lz+LxTzgcX282FqLD7ypWi2fI81bSHs/qOlFKBpbCdS4unt+iulHLY7krUO1cdSl8z7Ef7mC4a9n72uYQs15Y8HY4J0YZ2dVwv5X6XQVCtZRnpK6ODkfjoR75zQobHuIuCI+Y761g0fFwxRJHy8AEndiiDByNi+SnuG5Uu2O5YsJejV2vDHzSmTIIs/sVXT6sSP0rjtm0axCoicfWsKiB7xZoNhlIz9IacNEY3B1LMyzpOWq5DYGGQEOgITARATMsVe5YJo6lVW8INAQaAg2BM0LA3bHIyth36bBW8ncRZwRVY6Uh0BBoCDQEYgRkPzhOt7vKm0o/+w8ZN8dUojZJhQAAAABJRU5ErkJggg==\n",
      "text/latex": [
       "$\\displaystyle \\left[\\begin{matrix}- \\beta e^{i \\gamma_{1}} \\sin{\\left(\\theta_{1} \\right)} \\sin{\\left(\\theta_{2} \\right)} + e^{i \\gamma_{2}} \\cos{\\left(\\theta_{1} \\right)} \\cos{\\left(\\theta_{2} \\right)}\\end{matrix}\\right]$"
      ],
      "text/plain": [
       "⎡     ⅈ⋅γ₁                    ⅈ⋅γ₂                ⎤\n",
       "⎣- β⋅ℯ    ⋅sin(θ₁)⋅sin(θ₂) + ℯ    ⋅cos(θ₁)⋅cos(θ₂)⎦"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd1vh"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZYAAAAaCAYAAACU/XqBAAAACXBIWXMAAA7EAAAOxAGVKw4bAAALJ0lEQVR4Ae2c65EUNxDHF+oCwOcIfGSAcQTGGfCIAJwBFJ/gGwUZYEfgRwa2I+CRgcnA9mWA/z+hntJopdmZHWluZk9dNSeNHq3WX91qPWbvxufPn3dbpZcvX96S7OcKP221D03uhsBcBJodzEWw1S+NwJmU8jcxvRMw/mFDE/Vzyf2PnjeB/C3aELhuCDQ7uG4jvqL+yl88lTg/BiI927148eIDu5b2rBcDjdFFrfGpyXtrMm8Ri1oYr51vrbGqxXcpPGvJP5avyj3huRl4mRZdIQJ+NRDuKEtLeeHbKM23Cr/KeGwKiyoAb4Bp04H0IK0JlxvyLh8k0LdpUcumqp0LcWTLxPHVbc/9mdIvfXxUoPL3VJDJ9jvFH1glz/+t3sn/U+8/BHlP9P6TvZcIfXsfxIs+FOWNfOJ5XwF9fMZ7ipTHNvRrn/eL3j+myg2leR6XCmf1QfUZ303jUQqLIbxL5UnWVdhBqf6M4aM+V7eJLemAYbYWXCTHE2RabMeiBnEE9xUyCb/Rg4NhEvpZz2hSPS7s78BDIfx4N8KpcGeEQ/moPN53Cmm7xgU/EyntF3fMvl/PFQ45lT/UNg6BMq/00PfJpPpg+aNvc3L9oMLm8SiIRQBL+agfq7XYQfkOJjj6Ple3ia3ogEG0RlwWcSy+44/8gBkehEz2U495zlXnJ/Fi5fJJodvt+HcmR/LYrTDZujyF90hTuFNIHIf0VM9b3kk/hlQXnl8pDC+ujmGVqvNaic4xpjLVJvk7hW6XoZC+cpQzFU9jT1uOpyVMDdX2qeBRAgv0zK3epuI4svwsO6ANyYe+/KGHxdEWaEmbmK0DCwK6OlwWcSwCmGOo1MqbY6xJRzfiY86EybybeJX+O3nRYHI0hHGbgyH7N1+WVToyseo/msQr5H00n0TFh+KdPJpSOs6DI7CUI7hI8DqY5NuizVmTjOpvHo9CWIDjLCyHBk0yzrID1UeHeI7SlyHZKuYtZhOFdKAiFD3Wq8NlKcfSQ4EXDRw7DnYLj3mfQqqLwVKX3QnGkSMcTXz/ER5bYVS1JsKcTAfTPTaxkwzrPdcLE4vbhZGhuE1i52HBiXHafDixTvXiV4THKrEIwfZjfpQdqC5Hxb+H/NYcbzqQHp214nKWFrdcqjqO4nM0xSTODoOJ3FZJ3yp9b2JXGpMkkyeX/BAX/TgIKwtPDIP7BXYuyaMon9/bESktnLBZ8XeX/4r3SGVxWrRFu8gEcaSH3PSBOw1C+uf4+DrcG5H+qx7awIlC3+l5pzLsloaIO6LOaYQFVRc54BfvAO/6cmH/wqpj4rRJ20M7peuCxyAWY8CcW8aPdRU7mCqblwVdNhuEBbrcOaegzN+eP3bLCUGny4pnbcrXyQVXYRMHdSDoc8MlGLmzIF4ryg8u3SSo0C6XmfxwMky+8cTPxEXeA5V3eQpRRgzMJlMGnAk+dxyk7C9nyAqSR12qSxs4q1771IOUzgT+WiEK3ZHenRNTSB9wMHyA0JHe4Uc6xsXuwS5Yd0qjH3yFlzq2U1ZHOAnkS5Gd23fG6gvBG5rjWJAZ/PdIMl83PLJY7IFTIUF4L2IHY0SXLNgpev4Y3aWOlw+nwc6ZRR765z6cIY0ykOLc4TgbUDioQ19qZP9ehU0M6oD603DJDNcSjqVrWgPRrYQVZ/XD003cfqBwBFzChxM+RtY5COWxOsjuNJRnRL339mKhb9vtgBR/qie1g0CR7yrvlh7aM4on/H8tIwoxrnuq28mpOAZIMeTqsCAhIgwwx/eRL8vXMWG1+3phB9cZdZg5Mk6bGEuKrhseQ1ik8CmWpjFkDKrbwQSBcRg4kG53ond0FDLboExqwcRikMUUCyEWWmNsSsX26Cps4pAONFz2hulLwlkmfSdFYCD/0mMKlCsapne7DBI9DzvOCssRdxNvlGiTNpM5uxGI32iw6olX6C7zwB92S70JXO+svFhd4VCojhx7joX29KBY/xFXiKFjOHtllZ6jPafmCx7CFAO8zDBF9u7ojTKSCX44lle8G/l07kxw1OHdkhWJQ7BIyqb6m8bDY/Hcd9icJyvwHM5ZLELQVB+dZaEQE2O4U37qmJYFRrfgiCvqfQk7QI+QGzlZpPSOtZTmSOnoA2V7eq90nIxzNIqTD6bv9PRIebaY4oSBxdyxNoWcubGi/Ro2kdUB9eNkcPF9KWobQ46FQRwzGfUUKXphlZtzCChDrCgo+tTJO2qy94pi9EggTukTZQGcSdvtsFSfDwZSk0WvnRovXgFg3e3efDt24d45UZUFX5vwMIIxNGS81N8sHpKdY81u3BRn8maRYT/UVbRHh7BwhUOeYW2lozN8ztubkMMyA/El7OCj2uc5JJ854dwCkW5YmdieyTNCH6GiOiR8Tbdr2MSQDlifN4+LxqS4bdx0Q13vDzsGlDdFucu4vVVPqvKYNLVtdzJjivfKqK5THHjoua3nhgowMT1R3IykV6fgC6s6M5gU2xhT+onD6wxbcVaKTBp7zjXF0KfRJm3vkXhtGg91iHEzR0v/WCgM/e4niwWVF6BV2IH6afqTc8BAYWWGdJajtDk6dBU2MaQD1ufN46LxK24btR1LaMgooCMpGBMzu5l44mewvnaFoj+qwwpwSUJGuyh37UoGdgRs/5G9JoGDM8KwEbV/Gb4T97iwsoqxjIuOeYePGUxcfut4sCjIHU3GfeV9CItU+ZJpq7EDr3MsZJI6j/7pIR/d7O5LDQzl2RzAfcQcHboKm8jqwInhUtw2qjkWAY+3f68wnpxRNJSM+5h4EmMVaYqo6Jf7A5Xj2ALlXZo4e6YfIfEeH+/FZSiPUh5L9JVPk1OEY3PtednA7HvFL1OFJ6ZxTDGE82bxED69HZ36iTG5L5oyGB3CIlOtSPLa7ID7IHZ3sS0jp+nL94rzQz2cR0iU4V84mc2M1aGQB/GrsIlDOnASuNSwjbN49Aq+4yBwCOcSPLyI/0fvye2j7yAX9ygjn/pBvDMJLE1M1OxQ2CZa2+ymOI9kQsKAuH9xKzm94ywf62GnQTr5fKVF+is9pFs/MC7+uWTuAvcXlaVeimgDGeAHjhw3xg46VW9MGmNmMsblTwYP4cXYsANm4sjREBa5OkXSJR9OcC12sJMs6Ps36hx6h02bbb4lj04r5OgVPN3XlqSJ0FHezakM6pCrkf9zFTYxqAPq18nhoj4VsY1q/91YAqKEJY5n8qp2wjnCDuPtfWV3THfFhwmUsUg6c+OpfCYBvr4bLGfllw4lVyk86CeTI9gmd3klsBAPcD/28n5peDfRXkEdOGgTJXRgKVAL4jLbNiSL29WeLdX51s5kBNi1sXvI7SAmMzxQgUUAba6VZuPhJwtW0O4uwL/vFMY7vhJY4LSSjmutAG9Artk6MKGPJXRgQnOzis7GxdtCMdu4Oas7mcoS8payhj7Dy9RsyYaAMOQYjhUvq4iq5Nugre5z5aoNHsF8Lh6+j+xUOL7hl+Bs+Zk8/g3F8eVmYyE+/K5itXiGfd5K3ONZ3SZK6cBSuM7Fxfe3qG3U2rHcFah2rroUvqfYDncw3LXsfW1zqLNeWdjtcE6MMbKq4X4r9bsFlGqpnZGaOpqOxkMt8psVFjyEHQmPuN9bwaLrwzWLHK0DE2xiizpwNC7Sn+K2Ue2O5Zope7XuemPgk86UQ5jdrvjyYUXqX3HM5l2DQU08toZFDXy3wLPpQHqU1oCLZHB3LM2xpMeopTYEGgINgYbARATMsVS5Y5koSyveEGgINAQaAieEgLtjkZex79LpWsnfRZwQVK0rDYGGQEOgIRAjIP/BcbrdVZ4r/ux/KR7HVG116isAAAAASUVORK5CYII=\n",
      "text/latex": [
       "$\\displaystyle \\left[\\begin{matrix}- \\beta e^{i \\gamma_{2}} \\sin{\\left(\\theta_{1} \\right)} \\sin{\\left(\\theta_{2} \\right)} + e^{i \\gamma_{1}} \\cos{\\left(\\theta_{1} \\right)} \\cos{\\left(\\theta_{2} \\right)}\\end{matrix}\\right]$"
      ],
      "text/plain": [
       "⎡     ⅈ⋅γ₂                    ⅈ⋅γ₁                ⎤\n",
       "⎣- β⋅ℯ    ⋅sin(θ₁)⋅sin(θ₂) + ℯ    ⋅cos(θ₁)⋅cos(θ₂)⎦"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd2hv"
   ]
  }
 ],
 "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.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}