NB6-Creating_Custom_gates.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from qiskit import(\n",
    "  QuantumCircuit,QuantumRegister,ClassicalRegister,\n",
    "  execute,\n",
    "  Aer)\n",
    "from qiskit.visualization import plot_histogram\n",
    "from qiskit.visualization import plot_bloch_multivector"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "simulator_q = Aer.get_backend('qasm_simulator')\n",
    "simulator_s= Aer.get_backend('statevector_simulator')\n",
    "simulator_u = Aer.get_backend('unitary_simulator')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Create a custom gate to use later in circuit"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAACoCAYAAAAoyHP8AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAWJ0lEQVR4nO3dfXRU9Z3H8fdMEpJAQB7CgwYCYsJTTJDEBYLUNGgR2dbWB7CobIUVMGQ9q5b1dBdRTylgKesRW49H2mrTVboK1IoVbBE0goiYimAADc8hEhTCYyBPM3P3j7s8REIySX4zdzLzeZ1zz0nu3PnNF+bmfu69v/u712VZloWIiIgBbqcLEBGR8KFQERERYxQqIiJijEJFRESMUaiIiIgxChURETFGoSIiIsYoVERExBiFioiIGKNQERERYxQqIiJijEJFRESMUaiIiIgxChURETFGoSIiIsYoVERExBiFioiIGKNQERERYxQqIiJijEJFRESMUaiIiIgxChURETFGoSIiIsYoVERExBiFioiIGKNQERERYxQqIiJiTLTTBYh/2q1Z7cjn1n7vVkc+V0TaJh2piIiIMQoVERExRqEiIiLGKFRERMQYhYqIiBijUBEREWMUKiIiYoxCRUREjNHgRwkpZceguAxqPZDYETL7QVyM01VJqDhTA5/uh+NnIK4dDE2Gnp2crkouFnFHKj6fj0WLFpGamkpcXBxDhw6lsLCQgQMHMn36dKfLM8KqqqJuwo/xrd9wYV5NDZ6Hf4rn5/OwfD4Hq2tYZTX8Zg0sWg3vfA7rdsLrm2HOCvjgS6erE6dZFryzDZ74M6wostePVVthwVvw+0KornO6Qjkn4kJl6tSpzJ07lxkzZrB69WomTpzIpEmT2Lt3L1lZWU6XZ4QrPh73XXfiffVPWJaF5fXi/cUCiIkh6meP4XKH1tde64Hn18Keby59rc4Lfy6CD3cFvy4JHX8rtnc2vA3sD31eBr99H0JwXykihdbWJcCWLl1KQUEBK1euZNasWeTm5jJ79myys7PxeDznQ+Xrr79m7NixtG/fnqFDh7JlyxaHK28+920/gIpjWBs+xPvsr7GOHiXqqSdwtQu9c0n/2A/lJ8BqZJm3PwOPN1gVSSg5UwNrihtfZs83sPNQcOqRxkVUqCxYsIBx48aRk5NTb35KSgoxMTGkp6cDkJeXx6BBg6ioqCA/P5+77roLr7dtbdFc8XG4J9yJ91fPYG3bRvT8ubg6tHe6rAZ9tBtcTSxzthaKvwpKORJiPt3f8BHKxVwu2LQnKOVIEyImVMrKyiguLmbChAmXvFZaWkpaWhqxsbGcPn2at99+myeeeIL4+HimT5+O1+tl06ZNAavN5XI1ObVYdTVRd0/E1aVLwGpr7bRzT3mjRynnTH3w0aDUoym0pv988pdYVuOpYlmwdsMWx2sN58lfERUqAL169ao3v6qqisLCwvOnvnbt2kW3bt1ITEw8v0x6ejo7duwIXrEG+N5dh++113HdMhbvG29iWf5stp1Rc/ZEkxsNgNqzJ4NQjYSa2qqTuFyNb6p8Pi81Z48HqSJpTMSEyrmQKCkpqTd/4cKFlJeXk5mZCcCZM2fo1Kn+NYqdOnWisrIyYLVZltXk1By+zZ/g/c3zRD05h6iZD8KJE1gfrA9Yba2d7rllcJMbjSg3bHrn90GpR1NoTX/9n/lNrqdudxSP3j/G8VrDefJXxIxT6d+/PxkZGcyfP5+uXbuSlJTE8uXLWbVqFcD5I5UOHTpw+vTpeu89deoUCQkJQa+5JXzbd+Cd/zRR//FT3Bl2H5F7wp14X/kTru+MDrkrvwCyU+H9L6Cm7vKd9dkpkBAX1LIkRPToBBl9YNvBhl93AR3j4fp+waxKLif0tjAB4na7WbZsGWlpaeTl5TFlyhQSExPJz88nOjqajIwMAFJTUzl69CgVFRXn31tcXMyQIUOcKt1v1r79eJ94iqgZ03DfMOr8fPdt34eTJ1t8tBJoV8TDg2MuHeR47ixueh/4UWbQy5IQck82pPS0f/722f2O8TDzJogNvQsbI5LLas5xTRiaPHkyW7duZdu2befn3XHHHSQnJ/P000/zyiuvMH/+fHbt2kVUVJRjdUbC44TP1sDmffCXf9i/D+sLo1LsjUkz+gklTPl8sLMcPt5z4ahl4nDI6qdACSURc6RyOUVFRZcMenzhhRfYvn07Xbp04bnnnmPFihWOBkqkaB8L3x104fefjIbUXgoUsbndkJYEU2+8MG9UqgIl1ERMn0pDKisrKSkpYebMmfXm9+zZkzVr1jhUlYhI2xXRoZKQkNDmBjWKiISyiD/9JSIi5ihURETEGIWKiIgYo1ARERFjFCoiImKMQkVERIyJ6EuK25JgjmwXEWkpHamIiIgxChURETFGoSIiIsYoVERExBiFioiIGKNQERERYxQqIiJijEJFRESMUaiIiIgxChURETFGoSIiIsYoVERExBiFioiIGKNQERERYxQqIiJijJ6nIhIGHn7Vmc999l5nPldCl45URETEGIWKiIgYo1CRkHLy7IWfDx6DOq9ztUjoqfFAacWF3yurnatFGqY+FXFc2TH4cBdsL4NTF20k/ns1uF2Q1AVGXAPXXw1xMc7VKc44UwMf74FP9sHhk2BZF157fAV0bg8ZfeCGAdCzk3N1is1lWRd/RSLBU1kNyz+Bz0r9Wz4+Bm6/Hv7panC5AltbWxOOHfU+CzaUwFtb/D9izU6BH2Zq58NJOlIRR+w/Cr97Hypr/H9PVR0s/Qh2fAX3jYLoqICVJw6rroOXP4AvDzfvfR/thi/KYUYu9LoiMLVJ49SnIkFXWgEvrG1eoFzss1J4eT14fWbrktBQ64EX32t+oJxz/Az8Zg18fcpsXeIfhYoE1bk90BrP5Zd59t6mT6ts/wrW7jBbm4SGt7bAviONL9PUOlJZAwXrwaMLPYJOoSJB9dYWOH626eX88bfP7Y7bUOD1wekqOzSl5fZ8A+tLzLR16AS8u91MW61lWfYFB5XVdl9ROIu4PhWfz8czzzzDiy++yMGDBxk4cCDPPfcc06dPJycnhyVLljhdYtg6VQWb9phrz+uDdTvgnmxzbTbX6SpYt9M+l38uUFJ6wpjBMCTJubraKtMh8P4XkDsEYh3a0vks+8q1wi8u7AB1bg+jB8CNA6FdGG6Bw/Cf1LipU6fyxhtvMGfOHLKysti4cSOTJk3iyJEjPProo06XF9Y+3mO+H+TTA/CjTGgfa7Zdfxw/A4v/Die+deS152vY/TXcNgzGDAl+Xf46VPIhm9+cx5H9W6irPk18px4MHHUPoyb8wpF6Kirhi0Nm26yugy0HYOQ1Ztv1h8+CpRuhaD9cfLHiibPw18/sS+gfvMm5wAuUiDr9tXTpUgoKCli5ciWzZs0iNzeX2bNnk52djcfjISsrC4Ann3ySIUOG4Ha7Wb58ucNVh4+SFna8NsbjhX1Hzbfrj1c21h+sec65sxsrt9QfqBdK9m99hzcX3sqQ7/yEqYsPMOPFY4x/6HW69x3mWE27Dl/4vzPdrhM+3mMHCjT879p3FFZtDWZFwRFRobJgwQLGjRtHTk5OvfkpKSnExMSQnp4OQGpqKosXL2b48OFOlBmWLAvKjgem7YMObLjLT9jn/xvbCLqwx1mEom3vPk9yxi0MGHk3UdHtiIpuR8/+15M6/E7Hajp4LDDtOhHslmWf8mpqONWm3Y1ftNIWhdmB1+WVlZVRXFzMI488cslrpaWlpKWlERtrn0O57777AJg3b15QanNFwEi+mNgOzPx9Zb15TV3hdbnXvz3Q71eLl3DrSzNaUV3zZdw8k9z7n290GQtY+8lB7h2VHPB6/v2V5u3jx3fszp6iN9jyzmL6po+la9LgFn2uyXX3B4++Sf/M2+rNM7GOfHWkEperYyura57Y9p15cEnTe1E1Hrh6yCjKd30UhKpax99x8hEVKgC9evWqN7+qqorCwkLGjx/vRFmRI4DB6UQou9z+jbx0+7lcsOVMXkyn7lez44OX+eDVR+jYtQ8j7niKtJwpDlYVmO/RFaB2G/3MZnzvzVm2LYiYUElMTASgpKSkXoAsXLiQ8vJyMjMznSrN7z2Atsxnwc9eg9qLxg1c7tYi5/Y+/b31yMP/No13fzetdQU2074jdid9Y1wuGJF+Fb8Nwvfb3Nu0tIvvyIjb5zDi9jmcPfkNm/78FGt/9wDJ136Pjt16+92OyXX3T5vsfoiLmVhHenTrEPS/MZ8FP//LpRdxfFuUG0q2rqeDAxeaBErE9Kn079+fjIwM5s+fzx//+EfWrl1LXl4eL730EsD5TnoJjHM3hgyEPl0D025j+iXatwFpbB/YsuxLR0Nd+yt6kDr8LizLh6fW0CCiFugdRuuH29X0d+8CMvsSVoECERQqbrebZcuWkZaWRl5eHlOmTCExMZH8/Hyio6PJyMhwusSw17+H+TZdQL/u5ttt8nNdcG82xERfPliyU2Bgr8u86KBPVi7gwOdrqK2uxPL5+Gbfp2z438dITh9LlyudS8FrArB+QGDWO3/kDIKrL7NuuoDOHeAHzl1sFzARc/oLYMCAAbz33nv15k2ePJnBgwcTHx9/fl5dXR1erxefz0ddXR3V1dXExsZGRId6II1MMX9rlWt7wxXxTS8XCH26wcNj7TEHOy4aX3FFPOQOhhsHhebdlD01Z1n/6iOcPloKbjcduyUz+IbJDB37kKN1XdUF+ibCAYOXiEe77btaOyEmCvLGwOptsHE31Pz/4NgoN2T2g9uug44OrbuBFFGh0pCioiJGjhxZb960adMoKCgAYP369QDs27ePfv36Bbu8sNK9o/3ci20HzbWZ27KLloy5qgtMz7XPnT/1hj3vyR+BO4TPAWRPmEv2hLlOl9GgMYPtm4WaMuIaZ08vtYu2b8V/awY89po97+d3hN8pr4uF8KofeJWVlZSUlFzSSf+HP/wBy7LqTQoUM+683n4uigk3pDp3auPbOre/8HMoB0qoy+hjTyZ0bg/fv85MW6118e1YwjlQIMKPVBISEvB6dRvTYLqivX2vrpc+uPzAQX+u6OndNTzPR0c6lwsmDLdvBnn09OWXa2odiY6CyTdAfDuz9UnTtE8lQZfeB/5ltH1uuSWSu8GDuXq6X7jqGAf5N0GPFj4auF00TMsJXMe/NC6ij1TEOcP62s8TX7rJfka9P9wuuGkI3JKupz6Guy4d4Kfj7EclbNjl//uu6QGTRkJicAfQy0UUKuKYq7rAI7fYHfcbSux7aTUkLgaG97f7UHrqEbERIzYG7hpuXzW4oQQ+3V9/8Ow5LmDglfa4kCFJ9s6HOEehIo6KcttHLcP6wtla+OoYHK20b5EfHwNJXaFHR3V+R7LeXeHHI+2+lsMn4dBx+55Z0W7o3skeNBmrU6EhQ6EiIaN9O0jtBalOFyIhKcpt35UhUHdmEDO0/yciIsYoVERExBid/hIJA009d0QkWHSkIiIixihURETEGIWKiIgYo1ARERFjFCoiImKMQkVERIxRqIiIiDEKFRERMUahIiIixihURETEGIWKiIgYo1ARERFjFCoiImKMQkVERIxRqIiIiDEKFRERMUahIiIixihURETEGJdlWZbTRYi0RbUeKC6DAxVw6Djs+tqeP6wvJHeDIUnQs5OzNYpzLAv2H4WSw3DwmL2uAAy6Evp0hf49YGAvcIfZrr1CRaSZaurg78WwcTdU1Ta+7IBecGsGXN09OLWJ8ywLtpbC37fbOxuN6doBvjsYRqeGT7goVESaYe838OpHUFHp/3tc2BuO718HUWGy4ZCGnamB1z6GbQeb975+iXDfKEjsGJi6gkmhIuKn7WXw0nrw+lr2/rQkmPIdiI4yW5eEhtPV8Py7cPhky96fEAv5N8OVnc3WFWzabxLxw8Fj8HIrAgVg+1fw+mZzNUno8Prgt++3PFAAKmvghXVQWW2sLEcoVESa4PHC0o/A00SgPHuvPTVm894LHbYSPtZsh9KKxpfxZ/04VQUriszV5YSICxWfz8eiRYtITU0lLi6OoUOHUlhYyMCBA5k+fbrT5UkI+mg3lJ8w196KIvC14ohHQsvJKlhTbK69LQfsvru2KuJCZerUqcydO5cZM2awevVqJk6cyKRJk9i7dy9ZWVlOlychxrJgfYnZNo+fgZ2HzLYpztm0u3WnRRuywfA6F0wRFSpLly6loKCAlStXMmvWLHJzc5k9ezbZ2dl4PB6ysrKoqanh/vvvJykpic6dOzNmzBh27tzpdOnikK+OwzenzLdbtN98m+KMf+w33+bWg1DnNd9uMERUqCxYsIBx48aRk5NTb35KSgoxMTGkp6fj8XhISUlh8+bNVFRUcPPNN3P33Xc7VLE4ranz5C118Fhg2pXgqq4LzE6H12f2lGswRUyolJWVUVxczIQJEy55rbS0lLS0NGJjY+nQoQOPP/44SUlJREVF8dBDD/H5559TXd3GL8mQFmnN1TyNOXq67e6JygWBWj+g7YZKtNMFBEtZmX3JTa9everNr6qqorCwkPHjxzf4vo0bN9KvXz/i4uICVpvL5QpY29I6N/3rEq7NnVZvXlNX8Fzu9Ydfrf97QqfO1J4N4FZJAi5pcA53zX6/3jxT68f0B2ey7d0XWl6cYf4OaYyYI5XExEQASkrq94AtXLiQ8vJyMjMzL3nP8ePHyc/PZ968eUGpUUKPp878Eeq5P05vANqW4PLWBu479ASw7UCKmBH1Pp+PYcOGUV5ezqJFi0hKSmL58uWsWrWK0tJSNm3axIgRI84vX1VVxdixY7nxxhsVKhFsQwks/8S/Zc/tgX57j7Mh3RJgzg9bXpeEhrO18F/L/Fu2OesHwMO32LdvaWsi5kjF7XazbNky0tLSyMvLY8qUKSQmJpKfn090dDQZGRnnl/V4PEycOJHU1FQFSoRL7haYdvt0DUy7Elzt2wXmfl1uF1zVRm/XEjF9KgADBgzgvffeqzdv8uTJDB48mPj4+PPzHnjgAXw+H0uWLAl2iRJiene1jyqacwNJfwzra7Y9cc6wZHtEvUnX9oZ2bXTrHDFHKpdTVFRUb9DjgQMHKCgoYN26dXTu3JmEhAQSEhIoLS11sEpxitsFN6SabfOKeHujIeFhVKq9npg0eoDZ9oIpokOlsrKSkpKSep30ffv2xbIsqqqqqKysPD8lJyc7WKk4afQA6G7wFMcPM3UL/HDSpQPkDjbX3rW9IbWnufaCLWI66kVaY98R+PUa8LXyr+W6ZPjJaNBV5OGlzgvPvNP6sSUdYuGxf7aPZtsq7S+J+OHq7vZDlFpzmiOlJ9yTrUAJRzFRMCPX7n9rqfgYu422HCigIxWRZtl5yL4N/ulmDiHIToHbs9pu56v452QVLN0IXx5u3vuu6gyTb2j7D+gChYpIs52pgVVb4ZO9UNvErVaSu8H4oTDoyuDUJs6zLPh4j31FWFNXDSbEwo0DYcyQ8HkiqEJFpIXO1sLWUth/1L6bcVWt3QGfmAB9utmPDw7UOBcJfT4LviyHksNQdgxOnAULO0j6dIX+PSC9d/iEyTkKFRERMUYd9SIiYoxCRUREjFGoiIiIMQoVERExRqEiIiLGKFRERMQYhYqIiBijUBEREWMUKiIiYoxCRUREjFGoiIiIMQoVERExRqEiIiLGKFRERMQYhYqIiBijUBEREWMUKiIiYoxCRUREjFGoiIiIMQoVERExRqEiIiLGKFRERMSY/wPRwBRGNepe8gAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 517.72x204.68 with 1 Axes>"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sub_circuit=QuantumCircuit(3,name='toggle_cx')\n",
    "sub_circuit.x(0)\n",
    "sub_circuit.cx(0,1)\n",
    "sub_circuit.cx(1,2)\n",
    "sub_circuit.s(1)\n",
    "sub_circuit.cx(0,1)\n",
    "sub_circuit.cx(1,2)\n",
    "sub_circuit.draw(output='mpl')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "#This step turns the sub_circuit to a gate\n",
    "toggle_cx=sub_circuit.to_instruction()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUIAAAEDCAYAAABAqQ1VAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAcTUlEQVR4nO3de1SUdeLH8TeIgIIKDCKE5SVIthFNLdP8CVqomKvsViqmVtp6Se13NnOzpMturmjEnq3NPf3caotdQ8oLiYmhZk5eMrM1lVXDTZNIvCABQqIMzO8Pc1a8FCbDKN/P65w5x/k+l/k8nPHD95nnGfVwOBwOREQM5unuACIi7qYiFBHjqQhFxHgqQhExnopQRIynIhQR46kIRcR4KkIRMZ6KUESMpyIUEeOpCEXEeCpCETGeilBEjKciFBHjqQhFxHgqQhExnopQRIynIhQR46kIRcR4KkIRMZ6KUESMpyIUEeOpCEXEeCpCETGeilBEjKciFBHjebk7gDSML9fBiaPuTiGNSYsQ6HSnu1PUDxWhIU4chZICd6cQuTrp1FhEjKciFBHjqQhFxHgqQhExnopQRIynIhQR46kIRcR4KsJzlJaWMnnyZEJCQmjevDl9+vRh48aN7o4lIi6mIvyBw+EgISGBzMxMUlNTycrKIjg4mAEDBrB9+3Z3xxMRF/JwOBwOd4e4GqxYsYJhw4aRnZ3N4MGDATh9+jRWq5XIyEiys7PdnPDKbMuo+zdLqmuqeSP7SVZve4vT9kp6RA7kt/ctoJVfsGtDyjUloC3cmujuFPXDiBlhTU0NqampREZG4uvrS9euXbHZbHTq1ImJEycCsHz5ciwWC/Hx8c7tvL29SUxMZM2aNVRUVLgrfoPL+Ggem/+9nFce/ZRFSWfa84VFY92cSsR1jCjC8ePHM3v2bCZNmsSqVasYMWIEo0aNYv/+/fTo0QOA3NxcrFYrHh4etbbt3LkzdrudvXv3uiO6W2Rv+Rsj+88kzNIRv2atmDAkhc++/IDDxV+7O5qISzT6IkxPTyctLY2srCxmzJhB//79SUpKonfv3tjtdmcRFhcXExgYeMH2QUFBzuUAR44cYeDAgTRv3pyuXbs2us8PK06WcrQkn8jwHs6x64JvpLlvS/YX7nRjMhHXafRFOHfuXOLj44mNja01HhERQdOmTYmOjgbOXCw5fzYIXDD2yCOPEBUVxfHjx5k6dSr33Xcf1dXVrjuAHzJc6cNmW1+n16o4VQaAX7NWtcb9fQP4vrKsvg9NrmE22/p6eW+68lFXjboICwoKyM3NZfjw4Rcsy8/Px2q14uPjA4DFYnHO+s51diwoKIgTJ06wcuVKnn32WZo1a8bEiROprq5my5Ytrj2QBtTcpwVwZmZ4rvLKEpr7tnRHJBGXa/RFCBAaGlpr/OTJk9hsNudpMYDVamX37t2cfxE9NzcXLy8voqKi2LdvHxaLheDg/149jY6OZvfu3S48ijOz1St9xMb2q9Nr+TcLICTgBv7z7b+cY4XH9/N9ZRkdw7q46AjlWhQb269e3puufNRVoy7Cs4WVl5dXazwlJYXCwkK6d+/uHEtISKCoqIicnBznWFVVFRkZGcTFxeHn50dFRQUtW9aeFbVs2ZLy8nIXHkXDu7vXRN5Z/wKFxQeoqCzjteyZ3HrTIEKD2rs7mohLNOp/obpjx4506dKF5ORkgoKCCA8PZ8mSJc57As+dEQ4dOpS+ffsybtw4UlJSCAsLY/78+eTn57No0SIA/Pz8OHHiRK3XKCsrw9/fv+EOqgEk9n+S8u+/Y9rLt1FVfYrukQN48v6F7o4l4jKN/obqvLw8Jk2axNatW7FYLDz44IO0aNGCpKQkysrKaNasmXPdkpISZs6cybJlyygvL6dbt27MmzePmJgYAE6cOEFwcDCHDh3CYrEA0KFDBxYuXEifPn3ccnx1dTk3VIvURWO6obrRF+HFjB07lh07drBz5+XfDnLPPfdwww03MG/ePBYuXEhycjL79u2jSZMmLkhaf1SEUt8aUxE26lPjS9m2bRu9evX6Wdu++uqrjBkzhsDAQCIjI1m6dOlVX4Ii8uOMK8Ly8nLy8vKYMmXKz9q+TZs2rFmzpp5TiYg7GVeE/v7+Lr8BWkSuLY369hkRkbpQEYqI8VSEImI8FaGIGE9FKCLGUxGKiPFUhCJiPBWhiBhPRSgixlMRiojxVIQiYjwVoYgYT0UoIsZTEYqI8VSEImI8FaGIGE9FKCLGUxGKiPFUhCJiPBWhiBhPRSgixlMRiojxVIQiYjwVoYgYT0UoIsZTEYqI8VSEImI8FaGIGE9FKCLGUxGKiPFUhOcpLS1l8uTJhISE0Lx5c/r06cPGjRvdHUtEXEhFeA6Hw0FCQgKZmZmkpqaSlZVFcHAwAwYMYPv27e6OJw1gwO88yD2gX3ym8XJ3gKvJ+++/j81mIzs7m8GDBwMQExOD1WolKSmJ7OxsNyd0vY++yCBr81/Zf2gHlVXfk/OC3W1ZHn+1H90j4xgd97TbMogZjJkR1tTUkJqaSmRkJL6+vnTt2hWbzUanTp2YOHEiAMuXL8disRAfH+/cztvbm8TERNasWUNFRYW74jcY/2aBDO09hUeGveTuKCINxpgZ4fjx48nMzOSZZ56hR48ebN68mVGjRnHs2DGmT58OQG5uLlarFQ8Pj1rbdu7cGbvdzt69e+nRo4c74jeY2zoNAmDHV+vdmuOVzGnkHtjAnoOfkPHRPIJbhfP64/8mfV0yq7e9RcXJEm4M78aUhJfpENoZAHt1Fa+tfIJ1/3obDw9P7o2ZTvanf+P+u55m0G0PAbBq6xss+jCZ0opj9LYmgMOBp6cXTyS+ddEcu/Zv4I1VT5F/ZPeZXxJ3TOG+mOkXvEfOZ6+u4t31L7Lm8zSOlx0iwC+ECUNSuKPzr3hiwV1cFxzB48NfB2Dt5wtZ8P7j/N9jX2BpGVZvP0OpOyNmhOnp6aSlpZGVlcWMGTPo378/SUlJ9O7dG7vd7iy34uJiAgMDL9g+KCjIuRzgueee4+abb8bT05MlS5a4PL+Hh8cVP2y29S7PWZ8e/fV8Onfoy+i4Z1gxp5w3n/iSd20vsvbzfzDn4Wwyni0kukNfnvzbACoqywBYtG4un+1dxV8e3cI/Zx2gqLSAI98ddO5z1/4NzM+cxmPDX2PZH4rpGXU3tp3vXjLD14f/TdLf72ZE7O9Y/Ptj/PHhlWRtms/az//5k/nf/OBpPvzXQp4Zs5jls8v40yM2woMjaeLZhFmjF/Hp7vdZs+0fHDyym1cyp/DUqLevuRK02dbXy3vTlY+6MqII586dS3x8PLGxsbXGIyIiaNq0KdHR0cCZiyUX++GdPxYZGcnLL79Mz549XRdaLpDz2ZuM7D+TG0Ki8PbyYUzcs3h6NuHTPSsBWPv5PxjR7wnCLB3xadqM3wx5AU+P/77FV3+eRkyX4XSLuJMmTby4s9soom64/ZKvt+KTV4npMpw7OifQxLMJN4REMazPNNZ8/o8fzelwOMja/Fcm/vJFOl7XBQ8PD1oHtKXjdV0AsLQM46nR6cx/bxrPvfUr7ot5nO43xdXDT0h+rkZ/alxQUEBubi6PPfbYBcvy8/OxWq34+PgAYLFYnLO+c50dOzszHDNmDABz5sxxVexaHA7HFe9jWwaUFNRDGDc6VvINoUEdnc89PT1pE9ieYyXfAFBU+i0hge2cy32aNqOVf2vn8+Ol33JT21tr7bNNQDsu5XDxAb74zzo27lrmHHM4amgdcP2P5iypOEbl6QrCg2+65Dq33NifMMuNfFu0j3tjpv/o/q5WsbH9cLx65e/Nq4ERRQgQGhpaa/zkyZPYbDbuvvtu55jVaiUrK+uCmWFubi5eXl5ERUU1TGgBqDWbA2gdcD2Hiw84n9fU1HDku6+dxRTcKpyj55wKn6o6SWn5MedzS6vwWqfKAEdL8muV67naBLZj0G3j+d97/npZuQP8WuPr7ce3Rfto2zryouukfziH0/ZKfnFDL17JnMrMUT8+yxTXavSnxsHBwQDk5eXVGk9JSaGwsJDu3bs7xxISEigqKiInJ8c5VlVVRUZGBnFxcfj5+TVMaDeqrqnmdFUlVfbTAJyuquR0VWW9zEovV2CLUL4t+o/z+cBbH+Ld9SkUHMujyn6a9A/nUF1t5/ZfDAEgrvtY3rW9SGHxAU5XVfJG9lPUOGqc2w/o8QAbdi3hi/98RHVNNbYd77Inf8slX3/oHVNYvyODT3avwF5dRXW1nYNHdrPjK9uP5vbw8GBo70d4feUTHDici8PhoKj0Ww4U7gLOXIh6d30Kz4xdzFOj0/nXf9bywda/X8mPSq5Qo58RduzYkS5dupCcnExQUBDh4eEsWbLEeU/guVeBhw4dSt++fRk3bhwpKSmEhYUxf/588vPzWbRokbsOoUGt/fyfpL47zvl8yKxmAPzzqQOEBrVv0Cz3xjxG6jvj+NUzAQS3Cuf/HvuCKvspnnxtIBWVpdx43S3Mm7AaP9+WACTe+RRl3x/n0b/0xNOjCffEPIal5XU09Trz0UeXjjFMSXiZP707nrLvj9PbmsAd1l/h/cPy83UI7cwfx73PmzlPk/rOOByOGq4LjmBEvyd+Mvu4+Dk082nB79/6FcUnDhPUIpTfDEkhwD+E5LdHMSXhL86r3bPuT+eZN4fS6frb6BAWXU8/PbkcHg53/KpvYHl5eUyaNImtW7disVh48MEHadGiBUlJSZSVldGsWTPnuiUlJcycOZNly5ZRXl5Ot27dmDdvHjExMRfst1+/fkybNo377ruvIQ/nZ2kMnxFerpOnyvn1s4H86REb1vZ3XHSd/32lN71uHsr9d81q4HTXvoC2cGuiu1PUj0Y/IwS46aab+Oijj2qNjR07ll/84he1ShAgICCABQsWsGDBgkvur6qqiurqampqaqiqqqKyshIfH5/Lulwv9e/E99+xN/9TukXexanT3/PqisdoE9iOTtff5lxnw86l3NppEF5NvFm97S3yCrbxRGKaG1PL1cCIIryYbdu20atXr5+17YQJE0hLO/OXZ8OGDQAcOHCA9u3b11c8+RlqHNW8mfM0sxcOx8uzKZFtb+X5cSvwatLUuc7Hu5bwp8UPU+Oo5jpLBL9/MJO2rS99dfdS0j9MZtG65IsuS354FdEd+/7s45CGZ8Sp8fnKy8tp1aoVL730Eo8++qi74zQIE0+NxbV0anyN8/f3p7q62t0xROQq0ehvnxER+SkqQhExnopQRIynIhQR46kIRcR4KkIRMZ6KUESMpyIUEeOpCEXEeCpCETGeilBEjKciFBHjqQhFxHgqQhExnopQRIynIhQR46kIRcR4KkIRMZ6KUESMpyIUEeOpCEXEeCpCETGeilBEjKciFBHjqQhFxHgqQhExnopQRIynIhQR46kIRcR4KsLzlJaWMnnyZEJCQmjevDl9+vRh48aN7o4lIi6kIjyHw+EgISGBzMxMUlNTycrKIjg4mAEDBrB9+3Z3xxMRF/FwOBwOd4e4WqxYsYJhw4aRnZ3N4MGDATh9+jRWq5XIyEiys7PdnPDn25YBJQU/vd5rK2fy6Z73OVbyDb4+/tweNYTfDHmBls2DXB9SrikBbeHWRHenqB/GzAhrampITU0lMjISX19funbtis1mo1OnTkycOBGA5cuXY7FYiI+Pd27n7e1NYmIia9asoaKiwl3xG4ynZxOeHLWQpX84zoLHdlBUWkDqO+PcHUvEpYwpwvHjxzN79mwmTZrEqlWrGDFiBKNGjWL//v306NEDgNzcXKxWKx4eHrW27dy5M3a7nb1797ojeoN6eHAyEeHd8GrSlAD/1iT0mcaO/evdHUvEpYwowvT0dNLS0sjKymLGjBn079+fpKQkevfujd1udxZhcXExgYGBF2wfFBTkXH7q1CkeeughwsPDCQgI4M4772TPnj0NejwNafu+D+kY1sXdMURcyoginDt3LvHx8cTGxtYaj4iIoGnTpkRHRwNnLpacPxsEao3Z7XYiIiLYunUrx48fJy4ujpEjR7o0v4eHxxU/bLb1l/26G3YuJXvra0wZ9nL9H5Rc82y29fXy3nTlo64afREWFBSQm5vL8OHDL1iWn5+P1WrFx8cHAIvFQnFx8QXrnR0LCgrCz8+Pp59+mvDwcJo0acKjjz7Krl27qKysdO2BNDDbjsX8eckEnn8oi8i23d0dR8SljChCgNDQ0FrjJ0+exGazOU+LAaxWK7t37+b8C+m5ubl4eXkRFRV1wf43b95M+/bt8fX1dUH6MxwOxxU/YmP71fn1PvjsTV5eOonnx63gloj+LjsuubbFxvarl/emKx911eiLMDg4GIC8vLxa4ykpKRQWFtK9+39nOwkJCRQVFZGTk+Mcq6qqIiMjg7i4OPz8/Grt47vvvmPq1KnMmTPHhUfQsDI3/oW/vT+DuRNy6Nyhj7vjiDSIRn8fYU1NDd26daOwsJDU1FTCw8NZsmQJ2dnZ5Ofns2XLFm6//XaAH2ZOsezbt4+UlBTCwsKYP38+H3zwAZs2bao1ezx58iQDBw4kJibmmijCut5HOOB3HjTx9KKpl0+t8RVzyl2UTK5Vjek+wkZfhHBmNjhp0iS2bt2KxWLhwQcfpEWLFiQlJVFWVkazZs2c65aUlDBz5kyWLVtGeXk53bp1Y968ecTExDjXsdvt/PrXv6Z169b8/e9/d8chXba6FqFIXakIG4GxY8eyY8cOdu7cednbPvTQQxw7dozly5fj5eXlgnT1T0Uo9a0xFeG18bfYBbZt20avXr0ue7uDBw+SlpaGr68vAQEBzvHdu3dzww031GdEEWkgRhZheXk5eXl5TJky5bK3bdeu3WVdjRKRq5+RRejv7091dbW7Y4jIVaLR3z4jIvJTVIQiYjwVoYgYT0UoIsZTEYqI8Yy8amyiFiHuTiCNTWN6Txn7zRIRkbN0aiwixlMRiojxVIQiYjwVoYgYT0UoIsZTEYqI8VSEImI8FaGIGE9FKCLGUxGKiPFUhCJiPBWhiBhPRSgixlMRiojxVIQiYjwVoYgYT0UoIsZTEYqI8VSEImI8FaGIGE9FKCLGUxGKiPFUhOcoLS1l8uTJhISE0Lx5c/r06cPGjRvdHUtEXExF+AOHw0FCQgKZmZmkpqaSlZVFcHAwAwYMYPv27e6OJyIupP/g/QcrVqxg2LBhZGdnM3jwYABOnz6N1WolMjKS7OxsNycUEVcxYkZYU1NDamoqkZGR+Pr60rVrV2w2G506dWLixIkALF++HIvFQnx8vHM7b29vEhMTWbNmDRUVFe6KLyIuZkQRjh8/ntmzZzNp0iRWrVrFiBEjGDVqFPv376dHjx4A5ObmYrVa8fDwqLVt586dsdvt7N271x3RRaQBNPoiTE9PJy0tjaysLGbMmEH//v1JSkqid+/e2O12ZxEWFxcTGBh4wfZBQUHO5QCjR4+mTZs2tGrVip49e/LJJ5803MGIiEt4uTuAq82dO5f4+HhiY2NrjUdERNC0aVOio6OBMxdLzp8NAheMJSUl8eabb+Lt7c3KlSu59957OXTokOsO4CIZRKRu6noJpFHPCAsKCsjNzWX48OEXLMvPz8dqteLj4wOAxWJxzvrOdXbs7Mzw5ptvxtvbG4fDQdOmTTl8+DCVlZUuPAoRcbVGPSMsKCgAIDQ0tNb4yZMnsdls3H333c4xq9VKVlbWBTPD3NxcvLy8iIqKco6NHj2apUuXcurUKaZOnYqvr69Lj0MX9kVcq1HPCIODgwHIy8urNZ6SkkJhYSHdu3d3jiUkJFBUVEROTo5zrKqqioyMDOLi4vDz83OOv/3225w4cYL33nuPXr16ufgoRMTVGvV9hDU1NXTr1o3CwkJSU1MJDw9nyZIlZGdnk5+fz5YtW7j99tuBM7Ou2NhY9u3bR0pKCmFhYcyfP58PPviATZs2OS+qnK9r164sXryYm266qSEPTUTqUaOeEXp6erJ48WKsViuPPPII48aNIzg4mKlTp+Ll5UWXLl2c63p4eJCVlcWwYcOYPn06Q4cO5ejRo6xevfqSJQhnbrr++uuvG+BoRMRVGvWM8FLGjh3Ljh072Llz52Vtd/z4cdatW8eQIUPw8vLitddeY9asWXz11VfO03ARufY06hnhpWzbtu1HZ3k/5pVXXiEsLIw2bdqQkZHBypUrVYIi17hGfdX4YsrLy8nLy2PKlCmXva3FYuHjjz92QSoRcScjT41FRM5l5KmxiMi5VIQiYjwVoYgYT0UoIsZTEYqI8VSEImI8FaGIGE9FKCLGUxGKiPFUhCJiPBWhiBhPRSgixlMRiojxVIQiYjwVoYgYT0UoIsZTEYqI8VSEImI8FaGIGE9FKCLGUxGKiPFUhCJiPBWhiBhPRWiYKVOm4OXl5e4YIlcVFaFBNmzYQEVFhbtjiFx1PBwOh8PdIcT1Tp06xZ133sl7771HWFgYdrvd3ZFErhqaERri+eef5+GHH6Z169bujiJy1VERGmDnzp18+umnjBs3zt1RRK5KKkIDbNq0id27d9OhQwfat29PdXU17du3p6yszN3RRK4K+oywjj777DOee+45Nm/eTFVVFdHR0UyfPp0RI0a4O9pl8/Ly0meEIufQfRR18NFHHzFo0CB8fX1JTEykRYsWLF26lJEjR/LNN9/w+OOPuzviFXE4HHh4eLg7hojbaEb4E+x2O1FRURQUFLBlyxZuueUWAEpLS+nZsydff/01eXl5tGvXzs1Jf75FWR9iCWjJwJjb3B1FxC30GSGwfPlyfvnLXxISEoKPjw/t2rXj/vvvZ9euXaxbt46vvvqK+++/31mCAK1atWLWrFmcPn2atLQ0N6a/MoeOFLFjz1d4euqtIOYy+tTYbrczZswY3nnnHdq2bcs999xDy5YtycvLY+nSpTzwwAN8/PHHAAwcOPCC7QcNGgSAzWZr0Nz16cPN/8LXx5s+t3Z2dxQRtzG6CKdNm8Y777zDhAkT+POf/4yfn59z2TfffENAQABvvPEGAJGRkRdsHxoair+/P/v27XNpzidf+JtL9w/wh5ev3VmtyKXMmzmxTusZW4QbNmxgwYIFxMfHs2DBggsuFlx//fXAmc8C4cyp8MW0bNnSuY6IXJuMLcKXXnoJgHnz5l31V0zr+lvtchw6UsRf3lpGXJ8exP1Pj3rfv8i1xNgiXLNmDe3bt6dr164/ut7ZmeClZn1lZWUEBgbWe75zufLUeO2mz1m76XOX7V/Eneo6iTDyUmFJSQknTpygffv2P7nu2c8GL/Y54OHDhykvL7/o54cicu0wckZ49tbJo0eP/uS6sbGxzJ07l9WrV5OYmFhrWU5OjnMdV6rvU+N/Zq7mq4OHmDl5FM18fep13yLXIiNnhIGBgdx4443s2bOHtWvXXrD8yy+/dP75rrvuomPHjqSnp/PFF184x0tLS0lOTsbb25sHHnigQXLXh0NHivh33tf8z63RKkGRHxg5IwRITk5m5MiRxMfHk5CQQEREBEePHmXz5s3cfPPNZGZmAme+l/v6668zaNAgYmJian3F7uDBg6SmptbpFPtqUVxygqCAFrpvUOQcRn/FLicnhxdffJHPPvuMyspKQkJC6NmzJ7/97W/p27dvrXW3bt160X90YeTIkW5K//PV1NTomyQi5zC6CEVEwNDPCEVEzqUiFBHjqQhFxHgqQhExnopQRIynIhQR46kIRcR4KkIRMZ6KUESMpyIUEeOpCEXEeCpCETGeilBEjKciFBHjqQhFxHgqQhExnopQRIynIhQR46kIRcR4KkIRMZ6KUESMpyIUEeOpCEXEeCpCETGeilBEjKciFBHjqQhFxHj/DzrmAPsmSQ+8AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 408.156x325.08 with 1 Axes>"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#Let's create another circuit and insert it\n",
    "qr = QuantumRegister(4)\n",
    "cr=ClassicalRegister(4)\n",
    "new_qc = QuantumCircuit(qr,cr)\n",
    "new_qc.append(toggle_cx, [qr[0],qr[1],qr[2]])\n",
    "new_qc.draw(output='mpl')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "outputstate:\n",
      " [0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+1.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j\n",
      " 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j]\n",
      "UnitaryMatrix:\n",
      " [[0.+0.j 1.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j\n",
      "  0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j]\n",
      " [0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+1.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j\n",
      "  0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j]\n",
      " [0.+0.j 0.+0.j 0.+0.j 0.+1.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j\n",
      "  0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j]\n",
      " [0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 1.+0.j 0.+0.j 0.+0.j 0.+0.j\n",
      "  0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j]\n",
      " [0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 1.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j\n",
      "  0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j]\n",
      " [0.+1.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j\n",
      "  0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j]\n",
      " [0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+1.j 0.+0.j 0.+0.j\n",
      "  0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j]\n",
      " [0.+0.j 0.+0.j 1.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j\n",
      "  0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j]\n",
      " [0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 1.+0.j\n",
      "  0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j]\n",
      " [0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j\n",
      "  0.+0.j 0.+0.j 0.+1.j 0.+0.j 0.+0.j 0.+0.j]\n",
      " [0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j\n",
      "  0.+0.j 0.+1.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j]\n",
      " [0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j\n",
      "  0.+0.j 0.+0.j 0.+0.j 0.+0.j 1.+0.j 0.+0.j]\n",
      " [0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j\n",
      "  0.+0.j 0.+0.j 0.+0.j 1.+0.j 0.+0.j 0.+0.j]\n",
      " [0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+1.j 0.+0.j\n",
      "  0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j]\n",
      " [0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j\n",
      "  0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+1.j]\n",
      " [0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j\n",
      "  1.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j 0.+0.j]]\n",
      "counts:\n",
      " {'0101': 1024}\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcoAAAEDCAYAAAC1e0SgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO3deVxU9f7H8dcMOyPIomyuISouuGe5tWiZ5qVSu5JZZrmVWVnR3i31Wll2w26Lpd5rptctrX4hpubuTVMzFTU1Rc00XAGVnWHm9wdXEpUDJsOwvJ+Phw+HM9/vmc8Zhnmf71lNk1Mm2xEREZErMju7ABERkYpMQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImLA1dkFSPkI+eEmPM/UdnYZUoVkB57i+I3rnF2GiMMpKKsJzzO1sRyv6+wyRK7KmIAxTnndySmTnfK6UjFp06uIiIgBBaWIiIgBBaWIiIgBBaWIiIgBBaWIiIgBBaWIiIgBBaWIVClhYWG0bNmSZs2aYbFYDNuGhobSvXv3cqpMKiudR3mRrHNZxI+NJ3FxIjkZOdSNqkv02GjCbwx3dmkiUgyz2cydd97J0KFD6dy5M0FBQYXP2Ww29u3bx9KlS5kyZQr79+8vfC40NJTVq1fToEEDevXqxdq1a51RvlQCGlH+j91uZ/qg6SQmJHLX+LsY9p9hWAItTOk3haOJR51dnohcQZcuXfj555+Jj4/nnnvuISgoiJSUFHbt2sW+ffuwWq00a9aMp59+ml9++YVZs2YREBBQGJJNmzZl79697Ny509mLIhWYgvJ/di/bTdL3Sdz/4f10vK8jTW9pypAZQ6gZVpOECQnOLq9c5dvymbr4Oe4dW5u7XvVh3Mz+nM047eyyRIoYN24c69ato2nTphw8eJDY2FjCw8MJDAwkKiqKyMhIfHx86NKlC9OnTycrK4sHHniAPXv2sHHjRpo2bcr27dvp0aMHKSkpzl4cqcCqRVDabDZWfbCKCR0mEBsayzvd3uHA9wd4o+MbzB8zH4BdS3ZhCbDQ7LZmhf1c3V1p17cd+9bsIycjx1nll7t5qyeyYff/8cETm5j7SsFo+u25Dzq5KpE/TJo0iddeew2bzcbf//53IiMj+cc//sGhQ4eKtMvNzWXDhg0MHz6cqKgofvjhB4KCgmjQoAG//PKLQlJKpVrso5z7xFx2Lt5Jz9ie1GtTj0ObD/H58M9JP53OLaNuASB5TzIhkSGYTKYifUOahWCz2ji5/yT12tRzQvXlb8kPU3ng9tcIDSzYNzu8zzs89HYEx1MOExLQ0LnFSbU3YMAAYmNjyc3NpX///ixevLhU/TIzM/H39y/82cfHB5vN5qgypQqp8iPKrQu3smXuFobNGUb3J7rTuFtjej7bk4YdGmKz2qjXuiD8MtMy8fLzuqy/xb/gqLnM1EwAzp88z5R+U3iuznO80+2dKrf/MiPrLCfTjtC4TvvCaWG1GuHt6cvB5EQnViYCtWrV4qOPPgLgySefLHVIXrxPcvv27WzatInQ0FDi4uIcWa5UEVV+RLkibgWRPSKJ6BJRZHqt8Fq4uLkQ1jwMKDiY59LRJACXTPoi9guCmgQx9D9D+XH+j8wYMoNXtryC2cVx6xxlcQeFdx9tTetGJd89JCPnHAAWr5pFptfw9CMz+9w11yFVx4HvDxD7YPne3ePxxx+nVq1arFy5kk8//bRUfS4NyR49ehAYGMjOnTsZPHgw48ePv2yTrbPuWiLlq7R3ianSI8q0Y2kk70mmzd1tLnsu9WgqIZEhuHoUrCtY/C2Fo8aLXZjm7e9N9vlsdi/fzR3P3YG7lzudh3TGlm/j8JbDDl2O8uTt4QMUjCwvlp6dhrenrzNKEgHAxcWFESNGADBhwoRS9blSSKakpLB//37mzZuH2Wxm5MiRjixbqoAqPaJM+z0NAN/gol/wuVm5JG1IovltzQunhUSGsGvprstGlsl7kjG7mglqHMTJAyexBFioEVij8Pmw5mGc2HfCoedalsW98RomRMDxktvV8PIjyK8+B479RESdghWM5DMHycw+R3hoq2uuQ6qOiC4RTH7TsfdtvHhkFxUVRVhYGIcPH2bNmjUl9i0uJC+YMWMGDz30EHfccQcvvvhikb66H6VcrEqPKC2BBfsXTx44WWT6qn+u4tzxc9Rt/cemyKg7o8g4k8HelXsLp+Xn5bPtq200ubkJHhYPcjNz8fTxLDIvTx/PKndE7J03jmD+mrdJTjlERvY5pi15gQ5N7tCBPOJU7dsX7DffuHFjiW1LCkmALVu2kJ+fT4sWLfD09CxmTiJVfEQZ2DCQsBZhrIhbgbe/N36hfmz/Zjt7vtsDUCQoW/RqQXincOY+MZfosdH4Bvuyfvp6Uo+mMnjaYADcvd3JSS8aitnns/GweJTfQpWD+259kfTMVEa/fz15+Tm0a3w7L94/29llSTVXv359APbt22fYrjQhCQVHwR49epQGDRoQGhp62X5KkQuqdFCazWaGzBjCgmcWsDB2Id7+3nS8ryNdh3UlYUICYS3CCtuaTCaGzRlG/Nh4vn71a3IzcqkTVYfHFj1WeFpI7fDapJ9JJyMlA0tAwWg1eU8yPZ7q4ZTlcxQXswsjo99lZPS7zi5FpNDkyZOZM2cOqamphu2aN29Ow4YNS3UxgZ49e2Kz2fjtt9/KulypQqp0UAIERQQx+pvRRabNfnQ2wU2CcfdyLzLdu6Y3MXExxMTFXHFenj6etOjZgmWTlhH9ejQ/LvgRk8lEw+sbOqp8Efmf1NTUEkMSYOXKlfTq1YvExMQSLybwyy+/lFV5UoVV6X2UxTmy7Ujh+ZNX66/v/pXje4/zcvjLrJu6jodnPuzQU0NE5OqtWbNGV9yRMlPlR5SXyknP4VTSKboO7fqn+vsE+TDqq1FlXJWIiFRU1S4oPWp4EHdaV+MQEZHS0TZDERERAwpKERERAwpKERERAwpKERERAwpKERERA9XuqFcRqTyu9uLkyVML7ikbOiKryGORa6ERpYiIiAEFpYiIiAEFpYiIiAEFpYiIiAEFpYiIiAEFpYiIiAEFpYiIiAEFpYiIiAEFpYiIiAEFpYiIiAEFpYiIiAEFpYiIiAEFpYiIiAEFpYiIiAEFpYiIiAHdj1JEpAyNCRjjlNe92nt3SulpRCkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiUsm4ubk5u4RqRedRiog4SXBwML169aJ9+/Y0atQIV1dX0tLS2LFjBxs3bmTt2rXYbLYifRo3bszSpUsZNWoUy5Ytc1Ll1YuC8hJZ57KIHxtP4uJEcjJyqBtVl+ix0YTfGO7s0kSkimjevDmvvfYa/fr1u+LocMCAAQAcPHiQDz/8kA8++ACr1Urjxo1Zs2YNYWFhjBkzRkFZTrTp9SJ2u53pg6aTmJDIXePvYth/hmEJtDCl3xSOJh51dnlSDm5/zsSuQ/91dhlSRZlMJl566SV++uknYmJiMJlMxMfH8+KLL3LXXXfRq1cvBg8ezOTJkzl48CDh4eG89957bN68mT59+hSG5Jo1a+jfv7+zF6fa0IjyIruX7Sbp+yRGzB9B89ubA9CocyMmdp5IwoQERi4Y6eQKHW/19nl8s+EjDv6+g+y8TJa9bXVaLc9OuYV2jW9j0G2vOq0GkbJiNpuZOXMmDzzwAABTp05l/PjxHDt27LK2s2bN4tlnn6VPnz68//77tG3blm+++Qaz2cyaNWvo06cPmZmZ5b0I1Va1CUqbzcaaj9awYeYG0o6lERQRRL+J/Zj/9HwiOkcQMzmGXUt2YQmw0Oy2ZoX9XN1dade3HSveX0FORg4eFg8nLoXj1fDyJ7rTKHLzsohbNMLZ5YhUGe+99x4PPPAA58+fZ8CAASxdutSwvc1mIz4+niNHjrBx40a8vLywWq2MGTNGIVnOqk1Qzn1iLjsX76RnbE/qtanHoc2H+Hz456SfTueWUbcAkLwnmZDIEEwmU5G+Ic1CsFltnNx/knpt6jmh+vJzfdM7ANiRtMapdXzw1Wh2HVrPnl83Mm/1RGrVrMP0Z3czZ9WbLP/xMzKy0mhUpy2j7n6f60JaAmDNz2NawvOs+uk/mExm+t/0DEs2TeX+Hq9yx/VDAPh287+Yu/JNzmacolOLu8Fux2x25fn7PrtiHTsPrudf377EkRM/F6xEdB7FvTc9c9ln5FLW/DwWrJnEd1tncubc7/hZghje5x06t7yH5z/tQVitCJ7963QAVmydzaeLn+WTp7cT6BtaZu+hVBzdu3fnqaeeIjc3lzvvvJP//rd0m/cbN27MkiVL8PLy4tSpU9SuXZuPP/6Ybt26XXaQjzhOtQjKrQu3smXuFkbHjyaiSwQAjbs15uiOoyQuTqRe64Lwy0zLJLhp8GX9Lf6WgudTC9bivn3rW7b/33ZO7j/JQ/9+iDZ3t3Fo/WVxN4J3H21N60Z1y6Ca8vFE3w85fHxXkU2vc1e9xYqtn/PG0CWEBFzHvFVv8eLU2/n38/uwePoyd9VbbNn7Lf984gcCfEOZnvACJ1J/LZznzoPr+fCr0UwYmkCr625ibeIXTJr/ELe2uf+KNRw+vptX/n0nL943mxua/4Vjp/fzyvTe+Flqc3uHwYb1z1j6Kj/8HM/fHviC60KjOH32GOczU3Axu/DyoLk8FteW7378nCb1OvDBV6N4ffCXlS4kD3x/gNgHnXOnjOK8MPFToOBv5uLHzmQymfj4448BGD9+/FWF5MX7JAcNGsTmzZvp3LkzgwcP5rPPPivS3tnLWRmV9o4r1eJgnhVxK4jsEVkYkhfUCq+Fi5sLYc3DgIKDea44UrhkUu1Gten3Vj/qt6vvqJLlCpZtmUHMrS9QPygSd1cPHrjtNcxmFzbtSQBgxdbPGXDL84QGhuPh5sWwPm9jNv3xEV++dSY3tforbSO64+LiSve2A4msf0Oxrxe/cQo3tfornVvejYvZhfpBkdzVZTTfbf3csE673c43Gz5ixF8mER7WCpPJRG2/uoSHtQIg0DeUlwbN4cOvR/P6Z/dw703P0q7JbWXwDklFdNttt9G0aVOOHDnC22+/Xao+l4Zknz59+P3333nllVcAGD16tCNLlktU+RFl2rE0kvckc/NjN1/2XOrRVEIiQ3D1KHgbLP6WwlHjxS5M8/b3BqDDgA4ALH9vuaPKLqIs7jPXMCECjpdBMU50Ku03QgL+OE3HbDYT7N+QU2m/AXD67DGC/BsUPu/h5kXNGrULfz5z9hhN6nYoMs9gvwYU53jKIbYfWMV/d35ZOM1ut1Hbz3jze1rGKbJzM6hTq0mxbdo0upXQwEYcO72f/jc9Yzi/iiqiSwST36xY90BMnlrw/+SUyUUel6dLR3YPPvggAJ9++ilWa8kHx10pJC/sk5w3bx7vvvsu7du3p1mzZuzZs6ewn+5H6ThVfkSZ9nsaAL7BvkWm52blkrQhqXCzK0BIZAjH9x3HbrcXaZu8Jxmzq5mgxkGOL1gKXTwaBKjtV4/jKYcKf7bZbJxIPVwYXLVq1uHkRZtac/KyOJt+qvDnwJp1imyKBTiZdqTY1w/2b8Ad1z/C139PK/z3fxPOMT12t2HdfpbaeLpbOHZ6f7Ft5qx8g1xrNs3q38gHXz1uOD+p3Dp27AjAkiVLSmxrFJIAOTk5rF69GoDrr7/eMQXLZap8UFoCC/Yvnjxwssj0Vf9cxbnj56jb+o/9dlF3RpFxJoO9K/cWTsvPy2fbV9tocnOTKn/EK0C+LZ/cvGzyrLkA5OZlk5uXfdnKQ3nw9wnh2OkDhT/37DCEBWve4eipX8iz5jJn5Rvk51u5oVkfAG5r9yAL1k4iOeUQuXnZ/GvJS9jsfxzwcHv7wazfuZDtB1aTb8tn7Y4F7DnyQ7GvH915FGt2zGPjz/FY8/PIz7fy64mf2ZG01rBuk8lEdKfHmJ7wPIeO78Jut3P67DEOJe8ECg6UWrDmHf724Be8NGgOPx1YwdLN/76Wt0oqKA8PD5o2bUpeXh67dxuvYJUUkhds27YNgFatWjmkZrlcld/0GtgwkLAWYayIW4G3vzd+oX5s/2Y7e74r2GRxcVC26NWC8E7hzH1iLtFjo/EN9mX99PWkHk1l8DTjgzeqihVbZ/HugocLf+7zshcAs146REhAw3Ktpf9NT/Pu/Ie5529+1KpZh0+e3k6eNYcXp/UkI/ssjcLaMHH4ciyeBVsL7uv+Eucyz/DEPztiNrnQ76anCfQNw821YAWnVfhNjLr7ff6x4BHOZZ6hU4u76dziHtxdr7wCdF1ISyY8vJgZy17l3fkPY7fbCKsVwYBbni+x9od7vYGXhw9jP7uHlPPHCfAJYVifd/CrEcSb/xnIqLv/WXi07sv3z+FvM6JpWu96rguNKqN3TyoCs9nMJ598gs1mIy8vr9h2bm5uLF26tMSQBPjhhx+YPn06P/74o6PKlkuYJqdMLv+hQjk7eeAkC55ZwJGfjuDt703H+zri4eNBwoQEJv46EXcv98K2mWczCy9hl5uRS52oOkS/Hk2jzo0um+8H0R/QbVg3hx/1WhYaJvTHcrzyHPVaFrJy0un7mj//eGwtLRp2vmKbJz/oxI3No7m/x8vlXF3llxFylMN9Fjm7jCKSpxas2IWOyCryuDz92aNP77zzTp588kn69ev3p86T1D5Kx6nyI0qAoIggRn9T9Cix2Y/OJrhJcJGQBPCu6U1MXAwxcTHFzi8/Lx9bvg27zU6+NZ+87DxcPVxLPLdOHOt8Zip7j2yibeMe5ORmMiX+aYL9G9C03h/7ctYnLqJD0ztwdXFn+Y+f8cvRH3n+vplOrFqkwJIlS0q1H1PKX7UIyis5su0IDTs0/FN9542Zx5a5WwA4uPEgAH/b/jcC6weWVXnyJ9js+cxY9ip/n/1XXM1uNK7bgfEPx+Pq8sdFp9ftXMg/vhiKzZ5PWGAEYx/6irq1iz86tThzVr7J3FVvXvG5N4d+S1R4tz+9HCJSsVTLoMxJz+FU0im6Du36p/oP+mgQgz4aVMZVybWqaanFx08Z77d5ZdDcMnmt+3u8rM21ItVEtQxKjxoexJ2Oc3YZIiJSCVT500NERESuhYJSRETEgIJSRETEgIJSRETEgIJSRETEQLU86lVExFH+zBVyKsIVhaR4GlGKiIgYUFCKiIgYUFCKiIgYUFCKiIgYUFCKiIgYUFCKiIgYUFCKiIgYUFCKiIgYUFCKiIgYUFCKiIgYUFCKiIgYUFCKiIgYUFCKiIgYUFCKiIgYUFCKiIgY0P0oRUSEMQFjyv01/8y9O51BI0oREREDCkoREREDCkoREREDCkoREREDCkoREREDCkoREREDOj1ERESumre3Ny1btqRmzZrk5eWRlJTEb7/9Vmz7Dh06EBgYyLJly8qxyrKhoLxE1rks4sfGk7g4kZyMHOpG1SV6bDThN4Y7uzQREafy8/PjoYce4uGHHyYqKgqzuehGyVOnThEfH89HH33ETz/9VDi9Q4cOfPfdd3h6etK1a1e2bt1a3qVfE216vYjdbmf6oOkkJiRy1/i7GPafYVgCLUzpN4WjiUedXZ6IiNMMHTqUX3/9lcmTJ9O6dWvy8/PZsWMH3333HevWrePMmTPUrl2bRx55hK1bt7Jw4UKCgoIKQ9LPz4/FixezY8cOZy/KVdOI8iK7l+0m6fskRswfQfPbmwPQqHMjJnaeSMKEBEYuGOnkCh1vWsILbNqzmFNpv+HpUYMbIvswrM/b+HoHOLs0EXECT09P5syZQ9++fQFYvXo1H3zwAUuWLCEnJ6dI28jISIYNG8aIESPo378/PXr0wMXFBR8fHxYuXMjAgQOxWq3OWIxrUm1GlDabjVUfrGJChwnEhsbyTrd3OPD9Ad7o+Abzx8wHYNeSXVgCLDS7rVlhP1d3V9r1bce+NfvIycgpbvZVhtnswosDZ7No3Bk+fXoHp88e5d35Dzu7LBFxAjc3N7788kv69u1LSkoKAwcOpHv37nz11VeXhSTA3r17iY2NpWXLlmzevBk/Pz98fHxYtWpVpQ1JqEYjyrlPzGXn4p30jO1JvTb1OLT5EJ8P/5z00+ncMuoWAJL3JBMSGYLJZCrSN6RZCDarjZP7T1KvTT0nVF9+hvZ+s/CxX43a3N1lNG/Oud+JFYmIs7zyyiv07t2bkydPcuutt/Lzzz+Xql9QUBBNmjQp/LlRo0Z4eXlx/vx5R5XqUNUiKLcu3MqWuVsYHT+aiC4RADTu1pijO46SuDiReq0Lwi8zLZPgpsGX9bf4WwqeT83EmmNl/jPz+WXNL+Rk5FCvdT36v9OfkKYh5bdA5Wjb/pWEh7ZydhkiUs6ioqJ4+eWXAbj33ntLHZIX75NctGgRDRo0oEOHDkycOJHHH3/ckSU7TLUIyhVxK4jsEVkYkhfUCq+Fi5sLYc3DgIKDeS4dTQJw0aR8az61r6tNn1f64Bvsy8r3VzJz6Exe+O8LDqu/LK7q/+6jrWndqO5V9VmfuIglm6fxj0fXXvPrS9Vz4PsDxD5Y/necMPLCxE+Bgr+Zix9XdBWx7ueeew43Nzc+/vhj1q9fX6o+F4fkhX2STZs2ZceOHQwbNoxx48Zx8uTJwvbOXsbS3r2kyu+jTDuWRvKeZNrc3eay51KPphISGYKrR8H6gsXfQmZq5mXtLkzz9vfGw+JBz9ie+IX5YXYx0214N5J/TiYvO8+xC1LO1u74griFwxk/5Bsa123n7HJEpBwFBgYyYMAAbDYb77zzTqn6XCkkrVYru3fvJj4+Hnd3dx555BEHV+4YVX5EmfZ7GgC+wb5Fpudm5ZK0IYnmtzUvnBYSGcKupbsuG1km70nG7GomqHHQZfM/tPkQAfUDcPN0c9ASlM092xomRMDx0rVdumUGU+OfZfzD8bS8rss1v7ZUTRFdIpj8ZsW6n2Dy1IL/J6dMLvK4oqsIdV88uuvSpQseHh6sXr2aX3/9tcS+xYXkBbNmzeKee+6hR48eTJw4sXB6ZfjdQDUYUVoCC/Yvnjxwssj0Vf9cxbnj56jb+o/NkVF3RpFxJoO9K/cWTsvPy2fbV9tocnMTPCweReaRmZbJoucX0efVPg5cgvL11X//ydTFsbw1fJlCUqSaat++PQCbN28usW1JIQmwadMmANq1q5xbp6r8iDKwYSBhLcJYEbcCb39v/EL92P7NdvZ8twegSFC26NWC8E7hzH1iLtFjo/EN9mX99PWkHk1l8LTBReabm5XL9Pun07ZvW9rf275cl8mRPv6/p3AxuxL7ya1Fpse/ke6kikSkvIWFFRy3kZSUZNiuNCEJcOzYMbKzswkICMDT05Ps7GyH1O0oVT4ozWYzQ2YMYcEzC1gYuxBvf2863teRrsO6kjAhgbAWYYVtTSYTw+YMI35sPF+/+jW5GbnUiarDY4seK3JaSL41n5mPzKRWeK0qNZoE+G6S3dkliIiTjRo1itjY2CueK3mxWrVq4eXlVaqLCdSpU4fs7OxKF5JQDYISICgiiNHfjC4ybfajswluEoy7l3uR6d41vYmJiyEmLqbY+c17ah52m52YycW3ERGprPLy8jh79myJ7ZYuXUrXrl3Zvn17iRcTSElJKavyyl21CMorObLtCA07NLzqfim/pbBl7hbcPN146bqXCqe/tPEl/Ov6l2GFIiIV348//ujsEhyuWgZlTnoOp5JO0XVo16vuG1AvoNIcqSUiIteuWgalRw0P4k7HObsMERGpBKr86SEiIiLXQkEpIiJiQEEpIiJiQEEpIiJiQEEpIiJioFoe9VodZQeecnYJUsXoM1W1XO1pb8lTvQAIHZFV5HFVpKCsJo7fuM7ZJYiIVEra9CoiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAQSkiImJAt9mSKmlMwBinvO7V3tNPRCo+jShFREQMKChFREQMKChFREQMKChFREQMKChFREQMKChFREQMKChF/qSaNWs6uwQRKQc6j/IiWeeyiB8bT+LiRHIycqgbVZfosdGE3xju7NLEgdq3b0/v3r1p3749DRs2xMXFhTNnzrB9+3a+//57vvnmG3Jzc4v06d69O4sWLeKBBx4gISHBSZWLSHlQUP6P3W5n+qDpnPjlBHeNv4uaITVZP209U/pN4amlT1G3VV1nlyhlrHfv3owdO5aOHTte8flbbrmFMWPGcOrUKT766CMmTpxITk4O3bt3Z/HixXh5edG7d28FpUgVp6D8n93LdpP0fRIj5o+g+e3NAWjUuRETO08kYUICIxeMdHKFUlYsFgsffvghQ4YMASAlJYU5c+awceNG9u7di9VqJTQ0lA4dOtC/f3/atm3L2LFjiYmJ4f333ycuLg4vLy+mTZvGE0884dyFERGHqxZBabPZWPPRGjbM3EDasTSCIoLoN7Ef85+eT0TnCGImx7BryS4sARaa3dassJ+ruyvt+rZjxfsryMnIwcPi4cSlkLLg4+PDsmXL6NSpE1lZWbz22mt8+OGHZGdnF2mXmJjIsmXLeOONN7jpppv45JNPaNasGVOmTMFkMjFt2jRGjhyJ3W530pKISHmpFkE594m57Fy8k56xPanXph6HNh/i8+Gfk346nVtG3QJA8p5kQiJDMJlMRfqGNAvBZrVxcv9J6rWp54TqpSzNnz+fTp06cfjwYXr37s3evXtL7LNu3TqeeeYZ4uPjcXV1JTc3l0mTJikkRaqJKh+UWxduZcvcLYyOH01ElwgAGndrzNEdR0lcnEi91gXhl5mWSXDT4Mv6W/wtBc+nZgIwa8Qs9q3dhzXbSlDjIPq+2ZfrOl5XTksj12L48OH07t2b06dP0717dw4dOlSqft27d+fLL7/E1dWVQ4cOcd111/Gvf/2Lm2++WQP28E4AAA1PSURBVGEpUg1U+aBcEbeCyB6RhSF5Qa3wWri4uRDWPAwoOJjn0tEkAJdMuv2Z2xn44UBc3V3ZvXw3M4bMYPzP4x1VPuC8O2FUJd7e3rz99tsAPP7441cVkhcO3Jk2bRovvvgiu3btolu3bsTExDBv3rwi7fW7cq4XJn4KFPweLn5c0VXGuitjzZcq7d1+qvR5lGnH0kjek0ybu9tc9lzq0VRCIkNw9ShYV7D4WwpHjRe7MM3b3xugoI+7K3a7HRdXF86fOE9edp4Dl0LKwv3334+/vz8bNmxgwYIFpepzaUiOHDmSlJQUxo0bB8CoUaMcWbKIVBBVekSZ9nsaAL7BvkWm52blkrQhiea3NS+cFhIZwq6luy4bWSbvScbsaiaocVDhtFkjZrEjfgfWHCtdh3XFzdPNocuhexxevUvXbAcOHAjAlClTStX/SiF5YTPr7NmzmTRpEt26daNOnTocO3assJ9+V86VPLXg/8kpk4s8rugqY92VseY/q0qPKC2BBfsXTx44WWT6qn+u4tzxc9Rt/ce5kVF3RpFxJoO9K/84uCM/L59tX22jyc1Nihzx+uDUB3n7yNsMnT2Uhh0aOnYh5JqZTCbat28PwIoVK0psbxSSABkZGWzcuBGADh06OKZoEakwqnRQBjYMJKxFGCviVrB53mZ+WfsLC55dwKb/bAIoEpQterUgvFM4c5+Yy5b5W9i3Zh8zHp5B6tFU+rzS57J5u7i5EHVnFKs+WHVZEEvFUqdOHWrWrMmJEyc4fvy4YduSQvKC7du3A9CiRQuH1CwiFUeV3vRqNpsZMmMIC55ZwMLYhXj7e9Pxvo50HdaVhAkJhLUIK2xrMpkYNmcY8WPj+frVr8nNyKVOVB0eW/SY4Wkh1lwrKUdSCIoIKraNOFdmZiavv/466enphu18fHz44osvSgxJgKVLl5Kdnc3333/viJJFpAKp0kEJEBQRxOhvRheZNvvR2QQ3Ccbdy73IdO+a3sTExRATF3PFeWWkZLB/3X6a92yO2dXMxs83XrYJVyqelJQUxo8v+cjk8+fPM3DgQKKjo3nyyScNT/1YvXo1q1evLssyRaSCqvJBeSVHth350/sW101bx7yn5mEymwhtFsqI+SOoEVijbAsUp1m+fDnLly93dhkiUoFUu6DMSc/hVNIpug7tetV9LQEWnkx40gFViYhIRVXtgtKjhgdxp+OcXYaIiFQSVfqoVxERkWuloBQRETGgoBQRETGgoBQRETGgoBQRETGgoBQRETFQ7U4Pkerhz9zFIHmqFwChI7KKPBaR6k0jShEREQMKShEREQMKShEREQMKShEREQMKShEREQMKShEREQMKShEREQMKShEREQMKShEREQMKShEREQMKShEREQMKShEREQMKShEREQMKShEREQMKShEREQMKShEREQO6cXM180XsF2ycuZH3Tr3n7FKqlI/7fkz66XTsdjtBjYIY+MFAPH09nV2WXKX9/93PoucWYc21EtElggFxAzC7VOzxxKIXFpG4OJHzJ89Xmr/r1KOpzBk9h7PJZzG7mGnRqwXRr0U7u6xiVexPgJSppI1J5GbmOruMKunhmQ/z/PrneeG/L+BX14/VH692dklylWw2G/Ofms+QGUN4deurZJ/P5scFPzq7rBK17duW2NWxzi7jqphdzUS/Hs3Lm14mdk0shzcfZueSnc4uq1gKymrCmmNl8bjF3D3+bmeXUiV5+XoBBV+2uZm5mEwmJ1ckV+vIT0fwDfElJDIEgBseuIEd8TucXFXJwm8MxyfIx9llXJWaITWp37Y+AK7uroS1CCP1aKqTqyqeNr1WE8smLeOGB26gRq0azi6lwsj42YXMPZf/CZxa5HHZYxeLHf87cjHKv2n3T+PXH38luGkw9/z9njKvVy5nt0HqcnfyM4r+Yq70OwSwNLfi3Sz/ivM6+/tZ/Ov4F/7sX9eftGNpZVxxgawkF9K3l+6zZ/awE9A7F5OLQ0opNbsVUr51x5Zbuve6Rrs8vK6zlTjfjJQMdibs5LEvHyu7YsuYRpTVwO+7f+fXrb9yw6AbnF1KheLVKB9bhgnrGTPWM3/8KVz62HrGjFdEvmFIAgyfM5zxe8fT8PqGbPhsg6PKlouYzOAVkV+q36Etw4RnoyuHJIDdboeLf8d2R1RcwLNBPvZcSlW3Z8N8p4ckgMkVPBrYSlWz3Qqe9UsOSWuOlRlDZnDzozcT3CTYIXWXBQVlNXBw00GO7zvO+DbjGdd6HLZ8G+NajyP7XLazS3MqswfU6JBXYju34HzDL9gi8zSb6TiwI5vnbr7W8qSUPBvl4xZU8peyz/V5mN2Lf96vjl+RzX+px1LxC/MrixIvY3IFnxtK/uy5+tnwbl66z155sLSw4uJb8nvte2NeieFuy7cxa8Qs6kTV4dbRt5ZRhY6hTa+ldOSnI3w78VsObT6EzWojtFkot4y6hbZ92zq7tBJ1faQrXR/pWvjzM7Wf4fUdrzuxoorDOzKfzJ9tWFOKX2f07ZRnOJrMPpdNTmYONUNqApAYn1i4n0scz2QC3865nPm6+KOMXQNseDU1Dpz6betzNvksx/ceJyQyhE2zN9HqL63KutxCntfZcA/NJze5+ETx7ZSHqQINZ0wuBSGYutyj2DbudfLxKMVocv6Y+XjU8OCeCRV/N4WCshT2r9/PJ/d+gquHK+36tcOjhgc74ncwc+hM0o6lVfi1oZLY7ZS4WbGqMpkLvoxSEq78h+/V2Ip7kPE2uKxzWfz7oX9jzbaCCYIbB9P/7f6OKFeK4R5kxyvCStaBK3+llSZwzC5mYibHMGPIDKy5Vhp1bkSHmA4OqLaAyVRQ1+kvzRTd5lvAo14+HvVKEThPz+fn5T9jy7fxeovXibozinsn3euAiv9XVwMb7nXyyT12hYA32UtcsQQ4+MNBNv1nE6HNQpl08yQAbhh0AzePvNkBFV870+SUyQ7cEl/55VvzeeuGt0j7PY0xy8dQN6ouUPDl+N5t75FyJIVXtrxCQL0AJ1f656WudMPV147P9VZnl+I0Kcvcyfm16B++ydVO7ZhsXCxOKkquSn46nFrgid1a9Fvao0E+AXdU3NOi0ta6kbXvkoA32an91xxc/Srm13NeionTizzAXvS99m5mpWa3kjcpVzYVaFDvPDuX7GTqfVN5tcmrPBvyLONajePz4Z/z+8+/s3/dfk4fOk27/u0KQxIKTge4/enbyc/Nr9T7o/JOm8hOcq32nwTfG/PAXPRLydLGqpCsRFxqgKX1JSt7ZnvB77YC87k+D5Nb0c+ed4v8ChuSAG4Bdrwji27KNrnbS7XPvzKq1pte8635zB45m21fbcMvzI9Wf2mFp48nJ5NOsiN+B9fHXE/SxiQAIm+NvKx/ZPeCaUkbksq17rJ0/ic3TO52LC2r72gSwLWmHUsLKxk73QAwW2zUaFW935PKqEZrK5l7XbBlFKz5WVpaca1ZcQMHwMUbarS1cn5zwWfP5GHHp13FD5waHfLIOuCCPa9gVFmjnRUXLycX5SDVOigXPb+IbV9to9PgTtzzxj14WP7YT5V6NBWvml78MPsHAGo3qn1Zf99gXzxqeHAq6ZRD60ye6vhP34mZVfQT/ifZMswc/7fek8ouI9GNjEQ3Z5dxVew5Jk58Xvk+e+d/cOP8D5XrvQ4dkVWqdtU2KJM2JrHhsw1E9ohkQNyAy66k4l+34MTjC6dQFHfdTk8fz2p/moWISFVWbYNy7ZS1AES/Hl3hLzdW2rWeq5F32sTpLz2p0T4Pn/baxHiB3Qb550wVev+QlMyaZsLF116hTq0oid1eULebf+X67OWlFvy9VPCv0WtSbYNy35p9BNQPoE7LOobtLowkixs1Zp/PxqumYzeTOHLTa/pWN9K3Vq7NJSIiZaG0g5BKtL5VdjLPZpKTnkNA/ZJP6biwb/JK+yHPnThHTnrOFfdfiohI1VA9R5T/27KRfiq9xKaNOjdiRdwK9q7eS7v+7Yo8t3fV3sI2jlTWm15TlruT+7uZoIHZmIu/wIaIiFBNR5Teft7Uuq4WJ345wb41+y57/sT+E4WPm9zchMCGgfy06CeO7jxaOD3rXBbfxX2Hi7sL1993fbnUXRbyTpvIOeyCJcqqkBQRKYVqe2WebV9tY+bQmZhdzLTs3ZJa4bVIP5XOoS2HCGkawtBZQwvbFncJu9TfUrl7/N2V6hJ2WYfMnP/BjVr9chSUIiKlUG2DEmDPyj2s/nA1R346Ql5OHj61fKjfrj43P3YzjToV3Zz669Zf+XbitxzefJh8a37hRdHb9WtXzNwrLruNSnU0oIiIM1XroBQRESmJxhUiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIGFJQiIiIG/h8/uxl7E6VJkAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 588.756x325.08 with 1 Axes>"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#Statevector part\n",
    "result3=execute(new_qc,simulator_s).result()\n",
    "outputstate = result3.get_statevector(new_qc, decimals=3)\n",
    "print ('outputstate:\\n',outputstate)\n",
    "#unitary vector part\n",
    "result1=execute(new_qc,simulator_u).result()\n",
    "unitary=result1.get_unitary(new_qc)\n",
    "print ('UnitaryMatrix:\\n',unitary)\n",
    "#Now Quasm Simulator part where we need to measure\n",
    "new_qc.measure([0,1,2,3],[0,1,2,3])\n",
    "result2=execute(new_qc,simulator_q).result()\n",
    "print (\"counts:\\n\",result2.get_counts(new_qc))\n",
    "style = {'backgroundcolor': 'lightgreen'}\n",
    "new_qc.draw(output='mpl',style=style)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "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.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}