{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## TP - Diffusion thermique\n",
    "Modélisation de l'expérience vue en TP\n",
    "On suppose que le problème est à deux dimensions et on écrit l'équation de la diffusion thermique en régime permanent sous la forme :\n",
    "$$\\frac{\\partial^2T}{\\partial x^2}+\\frac{\\partial^2T}{\\partial y^2}=0$$\n",
    "On fixe la température à 50°C en x=0, la température de l'air extérieur est de 20°C."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "\n",
    "N=30\n",
    "Te=20\n",
    "T=20*np.ones((N,N),float)\n",
    "\n",
    "T[0,:]=50\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "On écrit ensuite :\n",
    "$$\\frac{\\partial^2T}{\\partial x^2}+\\frac{\\partial^2T}{\\partial y^2}=(T_{i+1,j}+T_{i-1,j}+T_{i,j+1}+T_{i,j-1}-4T_{i,j})/(\\Delta x)^2$$\n",
    "Ce qui permet de calculer $T_{i,j}$ (en prenant $\\Delta x=1$):\n",
    "$$T_{i,j}=(T_{i+1,j}+T_{i-1,j}+T_{i,j+1}+T_{i,j-1})/4 $$\n",
    "Il faut respecter les conditions aux limites : continuité du flux pour toutes les surfaces en contact avec l'air, par exemple en y=0 :\n",
    "$$\\lambda \\frac{\\partial T}{\\partial y}= h(T(x,0)-T_e)$$\n",
    "Ce qui donne :\n",
    "$$T_{i,0}=(a*T_{i,1}+T_e)/(1+a) $$\n",
    "avec \n",
    "$$a=\\frac{\\lambda}{h\\Delta y}$$"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "a=380/18\n",
    "nbiter=800\n",
    "k=0\n",
    "Ta=T.copy()\n",
    "\n",
    "for k in range(nbiter):\n",
    "    \n",
    "    for i in range (1,N-1):\n",
    "        for j in range (1,N-1):\n",
    "            T[i,j]=(Ta[i-1,j]+Ta[i+1,j]+Ta[i,j-1]+Ta[i,j+1])/4\n",
    "            T[i,0]=(a*Ta[i,1]+Te)/(1+a)\n",
    "            T[i,-1]=(a*Ta[i,-2]+Te)/(1+a)\n",
    "            T[-1,j]=(a*Ta[-2,j]+Te)/(1+a)\n",
    "    k+=1\n",
    "    Ta=T.copy()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "On affiche ensuite la répartition de température."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.colorbar.Colorbar at 0x1046f9470>"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAASoAAAD8CAYAAADAKumpAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAG0ZJREFUeJzt3W2sZVWd5/Hv7966RWEBAwxKaqCcaiNvCBPLpMKYYDqo\no412p8E3RCbSdNpM+aLH6ISkx+aNdMwkxPj4wpC+KOly2icSpTXE0dC0hDHpgIVN8yBOY+xypFJW\nyYMDFPVw7zn/eXF24bkP+7/23efcc/ap+/uQk3vP2Wevve/l3FVrrf3baykiMDPrsrlpn4CZWYkr\nKjPrPFdUZtZ5rqjMrPNcUZlZ57miMrPO2zbtEzCzs5ukQ8DLQA9Yjoh9ki4GvgnsAQ4BN0bEi3Vl\nuEVlZpPwjojYGxH7qucfBx6IiCuAB6rntVxRmdk0XA8cqL4/ANyQvVmTTKZfcvF87Nm9MLHjmW01\nh361xHMv9DRKGX/wjp3x/Au9Ru999PFTTwEnh15ajIjF4fdI+lfg/zHo+v11RCxK+m1EXFhtF/Di\nmefrGWmMStJ1wBeAeeBLEXFH9v49uxd45Ae7RzmkmSWu/oNfjVzG8y/0eOQHb2z03vldz5wc6s7V\neXtEHJb0BuB+ST8b3hgRISltMbXu+kmaB74IvBe4ErhJ0pVtyzOzbgig3/C/RuVFHK6+HgPuBa4G\njkraBVB9PZaVMcoY1dXAzyPiFxFxGvgGg36nmc2wIFiKXqNHiaSdks4/8z3wHuBJ4LvALdXbbgG+\nk5UzStfvMmC4nfks8B/XOdH9wH6AN17mNITZLGjaWmrgUuDewTAU24CvRcT3Jf0YuEfSh4BfAjdm\nhWx6zVENrC0C7HvLDs8pY9ZxQdAb00W2iPgF8JZ1Xn8eeFfTckapqA4DwyPjl1evmdmM69OtNsUo\nFdWPgSsk/R6DCuoDwH/Odni+v42vvvxv1922e+H59GAX6FTttgXVN1OzbfOF/xnZAN58ywvA0wiu\nzU/hmDC4Fj1pbTssveSjUCqzR/2HYSnq/49n216Kc9Jj/mpp/b+j5/tH0v2aCKB3tlRUEbEs6b8C\nP2Dwt3B3RDw1tjMzs6k5m1pURMT3gO+N6VzMrAMCWOrYFOW+DGdmKwRx9nT9zOwsFfmY3TS4ojKz\nFQbJ9G5xRWVmqyi9kjkNrqjMbIXBYPoWrqhO9Lfz+Kvrz57wm+3np/teOP9q7bbXzdVnrHZoqXbb\nduVJnwUt127LMlhzaXYrb1SXsl1tZOezmfpJTqit0r/0vSSplp1PnoXK/0xOR31S7WTUT2v0ar8+\nK/Xb3uvSYx45vf6MKCf6/5Lu18QgR7WFKyozmw39rdyiMrPuc4vKzDovUNqFngZXVGa2hrt+ZtZp\ngdILBNPgisrMVhgEPrdw1+9Ub55/Pb7+9BQvL+9I9z1/28n6bfP12143d7p22zlz9dEFgIUkvrA9\niS5k+5XiCZk89tCxex4KssHaUWIN2djKUtJKOJ1EELL9AE71swjC9tptL/fqP/Olv4fnTu1c/1x6\n42kJeTDdzDotQvQ2IQM3CldUZrZG3y0qM+uywWB6t6qGbp2NmU3dlh9MN7PZ0HOOysy6bMsn05f6\n8/z6+AXrbnt1uf4yLsAFC/URhJ3b6iMI585n29rHE9pumy/MZNA2gjCtGRLaajuTQSm6kF2tymIG\nbbcBnOjVxxNO9Oo/18eTz/xLS3k84aVT629f6o8nnrAZM1+Mwi0qM1thcFOyKyoz67BAxVbkpLmi\nMrMVIvIu9DR062zMrANEv+GjUWnSvKR/knRf9fx2SYclPVY93lcqwy0qM1shGHuL6qPA08DwlbTP\nRcSnmxbgFpWZrdFjrtGjRNLlwB8CXxrlfCbaour153jx1XPX3Xa6cNf38YX6S7k7F+ojCDuSCMKO\n+foZEADOmavfvm0uiyfURwWy/SCfXWFO7WZIaLvfqNpOvpbtV/rjWE4uzy8lrYRsv1P9/M/kZK9+\n+8kkunB8qf4zfWKpfj+A46fW37fXH73tEWicE+d9HvgLYPXqLR+R9CfAQeDWiHgxK2Skn0rSIUlP\nVP3Mg6OUZWbdMFgua1ujB3CJpINDj/1nypH0R8CxiHh01SHuBN4E7AWOAJ8pndM4WlTviIjnxlCO\nmXXChhYgfS4i9tVsuwb442qwfAdwgaS/jYgPvnYk6S7gvtJBPEZlZisEg2R6k0daTsRfRsTlEbEH\n+ADwDxHxQUm7ht72fuDJ0jmN2qIK4O8l9YC/jojF1W+omoL7ARZe/29GPJyZTcImz/D5KUl7GdQf\nh4APl3YYtaJ6e0QclvQG4H5JP4uIh4bfUFVeiwDnvvnfzdZ8uWZbUITGfq9fRDwIPFh9f/NG9x+p\nooqIw9XXY5LuBa4GHsr3MrMuGwymd+sWmtbVpqSdks4/8z3wHhr0Nc2s6wZzpjd5TMooLapLgXsl\nnSnnaxHx/WyHfl+cePWcdbf1evkPfWqh/lRPLNRnTnZsS1aLmc8zTeckOattLbNS2+YK07ykU7m0\nm+Zlfko5qmzytaxrkWV4SrdtLCc5oiwrtZycz6kkJwWwlGQATy7X73tqqd02gNOn1v/M9/ujjy0N\nBtPPkonzIuIXwFvGeC5m1hGe5sXMOm3MyfSxcEVlZmt4cQcz67QIWBrDPYPj5IrKzFYYdP1cUZlZ\nx21yMn3DJltR9UXvxPqXck/2CpedF+ovAZ/eXr/txHx9dGH7tjyesC2JLywkMYP5ZNso8YSs3LTM\nGZvmJZuqZJR4QlZu1tVZLkxBdHo5iT0ksZul0/V/fstL+THjdE25jieY2dbgrp+ZzYCm86FPiisq\nM1thcNWvW/f6uaIysxUc+DSzmeCun5l1mq/6hdDJ9fu+UYgnLCXbs5kX5ufrL+mf3pZf7s/2zaIC\n+bY8KpDtq3T2hPYRhLYfyVFCD9kfQmSr0BQS073k8ny2b7qtMLNHb7ndvv2lpNy6+EFFdfuOqYLx\nVT8z67QIpdPeTIMrKjNbY2t3/cys8zxGZWYzwRWVmXWac1RmNhO2do6qD3Mn1/8F9AuXnSO57NxP\nLg/3kwhCrxBP0HwSB0hiBHNJBCHbBoUIQrJvtl8XZRGEbIGCbL/Svvm2+s9QKTqTff5It9WXq6X8\nmHOna7a3m2BjhYh8FoppcIvKzNZw18/MOs1jVGY2E0pd7ElzRWVma3RtML1bI2ZmNnURgzGqJo8m\nJM1L+idJ91XPL5Z0v6Rnqq8XlcpwRWVmq4hef67Ro6GPAk8PPf848EBEXAE8UD1PuaIyszUi1OhR\nIuly4A+BLw29fD1woPr+AHBDqZziGJWku4E/Ao5FxFXVaxcD3wT2AIeAGyPixWJZAfOn1v/hVFg9\no9+rzwnFtmxbkp2ZL9TTSY6qn2wjyzuVclTZvtl+o+So2u47woBr9iHPzibL0xW3Z/msLCtVyFFl\n25VlpZJtc6Uc1VJNmWOI023wXr9LJB0cer4YEYtDzz8P/AVw/tBrl0bEker7XwOXlg7SpEX1N8B1\nq17bcNPNzGZEDMapmjyA5yJi39DjtUpK0pkGzqO1h4oIGkxtVqyoIuIh4IVVL2+46WZms6OPGj0K\nrgH+WNIh4BvAOyX9LXBU0i6A6uuxUkFtx6gaN90k7Zd0UNLB3vHjLQ9nZpMSYxpMj4i/jIjLI2IP\n8AHgHyLig8B3gVuqt90CfKd0TiMPppeabhGxeKZZOL9z56iHM7MJ2EDXr407gHdLegb4T9XzVNvA\n51FJuyLiSNOmm5nNjnEn0yPiQeDB6vvngXdtZP+2LaoNN93MbDYMWkvjiSeMS5N4wteBaxlchnwW\n+ASDpto9kj4E/BK4sdHR+jBXE0+gX7hsn07zklx2TqML+TEjiSCkc98nEYMoxBPyDEK2LZsCJj/k\nZkm7BtmHPN2vcNCW8YRsehQV4gmtIwi9pMya+MFr+9bFF8YwzQvM4OwJEXFTzaYNNd3MbHaMMP60\nKXxTspmtEKg4keWkuaIyszU61qByRWVmq4TnozKzWdCxJpUrKjNbY0u3qBQwf7pmW2n2hOzycRIz\n6CeXliO5dAyFeMJ8smMWQSiMUUZ2+3vL6ELH/nEcaBlBUOkPKLs8n3zGlEUFSvGEZN80gpDOnpAe\nkrnlmjLHNXtC4e9x0tyiMrOVgpGm8NkMrqjMbA3nqMys+1xRmVm3TfY+viZcUZnZWm5RmVmnRXlu\n+kmbbEUVMFcTTyjV4EouO2cRhCy6kEYMgJhPog1ZHGAu2ViYPSHdNzvmKJ+rKSzukB4yiyeUZgfI\nIghZxCU7ZhIxKG1PIwg1EYMmx6yNL4ytJbSVKyozmw3u+plZ57miMrNOc+DTzGaBA59m1n1b+qqf\nmc2EcdzcPE6uqMxspUaLrE9Wd6Z5Kfxi+smZts1YlXNUybaWU7mkOSnYnKxUcb9Naua3/LC3zVhB\nIWeV5ahaTgFT2p5O85Ltl2SsoD5HNZ6WkDyYbmYzYCu3qMxsRoxpfcBxcUVlZit1MEfVrcW7zKwT\nFM0exXKkHZIekfTPkp6S9FfV67dLOizpserxvqwct6jMbK3xjVGdAt4ZEa9IWgB+JOl/Vds+FxGf\nblKIKyoz2zQREcAr1dOF6rHhanDy07ws1ZzjCKuLZDECJT9htrINAMkl4ixmEFmHurQKTRJ72Lx4\nwibZhHhCabWiPILQblvp58hXmkm2jRJPWK77O8r3a2oDMYdLJB0cer4YEYsrypLmgUeBNwNfjIiH\nJb0X+IikPwEOArdGxIt1BymOUUm6W9IxSU8Ovbah/qWZzZBgcAtNkwc8FxH7hh6La4qL6EXEXuBy\n4GpJVwF3Am8C9gJHgM9kp9RkMP1vgOvWef1zEbG3enyvQTlmNiui4WMjRUb8FvghcF1EHK0qsD5w\nF3B1tm+xooqIh4AXNnZKZjbLxnjV7/WSLqy+Pxd4N/AzSbuG3vZ+4Mn19j9jlDGqRv1LSfuB/QAL\n5100wuHMbGLGd9VvF3CgGqeaA+6JiPsk/U9Je6sjHQI+nBXStqK6E/hkdZBPMuhf/tl6b6z6rIsA\nr3vD7o4F881sXWP6S42Ix4G3rvP6zRspp1VFFRFHz3wv6S7gvjblmFn3NO3WTVKrikrSrog4Uj0t\n9i9f2y+y1TPy30w6Q8K2JCqQXZIuRQXS2RPabStFBULtVqEZJYLQNvYw0oe57So0pWO2XMEm3Vaa\nPaHlvrURgwbHrIsvjK2CmbWJ8yR9HbiWQV7iWeATwLUb6V+a2WyZuRZVRNy0zstf3oRzMbOumLWK\nysy2mLNljMrMznKuqMys69J7H6fA81GZWee5RWVma23prl+aoyrsmmZg6nfuz4+wCk2WZWmZo0oz\nVtB+FZqzKEe1WavQtF+hpnDMdBWadlmptjmqsVQwHkw3s5ngisrMOs8VlZl1mejeVT9XVGa2kseo\nzGwmuKIys87b0hVV1E9tUWpq9pNpJ1pHFwr98GylmTSekEUMRogntI4ulLTdd4QPc+sIwijxhGzf\nlp8hKKxC03JbFmtI9538KjQT4RaVma3lisrMOi181c/MZoFbVGbWdR6jMrPuc0VlZp3WYhXkzTbR\nikpJPCFK19eTVWqyfbPoQnZJGoC55Jit4wmFn3MTZkgYKbowgtbdh5HiCclsBS1XqCl9TtJ4QrpC\nTf0JZWVm+46jy6YxlTNOblGZ2RquqMys+1xRmVnndayi8pzpZrZS/G5Z99KjRNIOSY9I+mdJT0n6\nq+r1iyXdL+mZ6utFWTmuqMxsrWj4KDsFvDMi3gLsBa6T9Dbg48ADEXEF8ED1vJYrKjNbQ/1mj5IY\neKV6ulA9ArgeOFC9fgC4IStnCrMn1GxK4geDN2QRhCRGkCzgoGRGBmi/SEMaB0giD8V9HU8ol9ly\n9oTWi0JQiES0jC7M0OwJl0g6OPR8MSIWV5QlzQOPAm8GvhgRD0u6NCKOVG/5NXBpdpBiRSVpN/CV\nqqCoTuQLki4GvgnsAQ4BN0bEi01+MjPrsI0FPp+LiH1pcRE9YK+kC4F7JV21antIedXYpOu3DNwa\nEVcCbwP+XNKVbLCPaWYzZHxjVL8rMuK3wA+B64CjknYBVF+PZfsWK6qIOBIRP6m+fxl4GriMDfYx\nzWw2nEmmj+mq3+urlhSSzgXeDfwM+C5wS/W2W4DvZOVsaIxK0h7grcDDQKM+pqT9wH6A7edeuJHD\nmdmUlGY13YBdwIFqnGoOuCci7pP0j8A9kj4E/BK4MSukcUUl6TzgW8DHIuIl6XcjtFkfsxpYWwQ4\n76LdHYuRmdkaY7wpOSIeZ9C4Wf3688C7mpbTKJ4gaYFBJfXViPh29fKG+phmNjvG1fUbl2JFpUHT\n6cvA0xHx2aFNG+pjmtkM2YTB9FE06fpdA9wMPCHpseq124A72EAfE0ARzC3VTfNS2DnZnmaastVr\n5ttnmrLpWtKVZkqz2WzCKjRTilGl2q5CU/xXvHVWqt30MFDKSiXljpKjqskjqviH1MzMzZ4QET+i\n/rPeuI9pZjNk1ioqM9tivAqNmXWdZ/g0s9kwprGucXFFZWZruEVlZt221VehIUZYPSONJ7SLIJRW\nvsnjCdl18Gy/9JCEWoYJRsggtJ0GZqR/dVvuW7r8ng4Ct40ulOIJrad5SVahqVmt6TV15zuuaV48\nmG5mXeeKysy6LfBgupl1nwfTzaz7XFGZWZc58Glm3RcxzonzxmLiFVXdXeGlqMBc0hZNZ09oOesC\nlOIJyY5JxKB0zKzNvRkr1EzNKDMkJPJ4QruZDKYRT8j2K+07Ft2qp9yiMrO13PUzs24LYKt3/cxs\nBnSrnnJFZWZruetnZp235a/6mVnHefYEUN1d4dlsBJBfm08u+WeXnYvxhPlkW7ZoRPKzlG72bL1I\nQ9tZF6al5b1km7a4QxZdGCUq0HZBidLiDps4e8Ig8NmtmsotKjNby7MnmFnXuUVlZt3WwTGqRku6\nm9lWMrjXr8mjRNJuST+U9FNJT0n6aPX67ZIOS3qserwvK8ctKjNba3xdv2Xg1oj4iaTzgUcl3V9t\n+1xEfLpJIa6ozGylMS5AGhFHgCPV9y9Lehq4bKPluOtnZmtFNHtsgKQ9wFuBh6uXPiLpcUl3S7oo\n27fYopK0G/gKcCmDIbbFiPiCpNuB/wL8pnrrbRHxvbSsCOZ661fVxRVhsulasrxTOldL++yWkkBP\nPgVM6ZBtg1Qtp4fZRK1vw0ingCn9P0v2bb1CTfuVb9pmpcrHrFvNaUxdtubFXCLp4NDzxYhYXP0m\nSecB3wI+FhEvSboT+GR1pE8CnwH+rO4gTbp+Y+ljmtnsUL9x3++5iNiXliUtMKikvhoR3waIiKND\n2+8C7svKKFZU4+pjmtmMCMYW+JQk4MvA0xHx2aHXd1V1C8D7gSezcjY0RtWmjylpv6SDkg6eXjq+\nkcOZ2RSIQNHs0cA1wM3AO1dFET4l6QlJjwPvAP5bVkjjq35t+5hVf3UR4ILzL+tYjMzM1jWmsa6I\n+BHrj66m49mrNaqoxtHHNLMZ0rFbaIpdv6yPOfS2Yh/TzGbEmTGqJo8JadKiOtPHfELSY9VrtwE3\nSdrL4Mc6BHy4yQFrp3nZrOhpdpl3vpQVSLZlKYJsapnCdCxZ7CGTRSK6OAFM60BhaUWYNMfSrtxi\nVCCLGbQst3jMur+jMdnAVb+JaHLVbyx9TDObFRsPc24230JjZisFrqjMbAZ0q+fnisrM1vLEeWbW\nfa6ozKzTIqBm8oBpmfgqNHVxAS0X9s0SX8nsCcwlMyAU/tVIZ2XIrvmnMyCUrq/nm2t369bnqmwT\nZlaAwq83W2kmixEUVoRpG20YZfaE2tjNuBpCblGZWee5ojKzTkt6PtPiisrMVgmIbo0luKIys5WC\nLT6YbmazwWNUZtZ5W72iUt3iDoUZZ9rOAJBFECKJLgzekJSbrRmRXOtW4ZitYw8jaLv4Q+vFG0pa\nxgiAdBA4jy7UbyrFWNIoQcuFH0qD2XV/R+Phm5LNrOsCmLVpXsxsC3KLysy6bavfQmNm3RcQzlGZ\nWec5mW5mnecxKjPrtIgtftXv+Al45Il1N+nq/5DumuWsshhQ9g9DMT6UZpqybVnIapRpXpJ9S5mw\ntofcLC27FsUcVds81Ag5KlquQpNmvkqD2TV/R8SJfL+m3KIys24Loteb9kms4IrKzFbyNC9mNhM6\nFk8oLuluZltLANGPRo8SSbsl/VDSTyU9Jemj1esXS7pf0jPV14uyclxRmdlKUU2c1+RRtgzcGhFX\nAm8D/lzSlcDHgQci4grggep5LXf9zGyNcQ2mR8QR4Ej1/cuSngYuA64Hrq3edgB4EPjvdeUoJngZ\nUtJvgF8OvXQJ8NzETqDM55Pr2vlA985p2ufz7yPi9aMUIOn7DH6OJnYAJ4eeL0bEYk25e4CHgKuA\n/xsRF1avC3jxzPP1TLRFtfoXKOlgROyb5DlkfD65rp0PdO+cunY+bUTEdeMuU9J5wLeAj0XESxrK\nGkZESHlCzmNUZrapJC0wqKS+GhHfrl4+KmlXtX0XcCwrwxWVmW2aqlv3ZeDpiPjs0KbvArdU398C\nfCcrZ9qD6ev2ZafI55Pr2vlA986pa+czbdcANwNPSHqseu024A7gHkkfYjBufWNWyEQH083M2nDX\nz8w6zxWVmXXeVCoqSddJ+j+Sfi4pTaRO6HwOSXpC0mOSDk7pHO6WdEzSk0Ovbeg2gwmcz+2SDle/\np8ckvW+C5zOWWzEmcD5T+x2dzSY+RiVpHvgX4N3As8CPgZsi4qcTPZGV53QI2BcRUwvqSfp94BXg\nKxFxVfXap4AXIuKOqkK/KCJq07sTOJ/bgVci4tOTOIdV57ML2BURP5F0PvAocAPwp0zhd5Scz41M\n6Xd0NptGi+pq4OcR8YuIOA18g0GcfkuLiIeAF1a9fD2D2wuovt4w5fOZmog4EhE/qb5/GRi+FWPi\nv6PkfGwTTKOiugz41dDzZ5n+/+AA/l7So5L2T/lchl1a3SsF8Gvg0mmeTOUjkh6vuoYT64oOq27F\neCvwMB34Ha06H+jA7+hs48H0gbdHxF7gvQzu7v79aZ/QajHoo087S3In8CZgL4MbTT8z6RNYfSvG\n8LZp/I7WOZ+p/47ORtOoqA4Du4eeX169NjURcbj6egy4l0H3tAs2dJvBZouIoxHRi8Gib3cx4d/T\nOG7F2Ozzmfbv6Gw1jYrqx8AVkn5P0nbgAwzi9FMhaWc1GIqkncB7gCfzvSZmQ7cZbLYzFULl/Uzw\n9zSuWzE2+3ym+Ts6m00lmV5dsv08MA/cHRH/Y+In8btzeRODVhQMbin62jTOR9LXGczPcwlwFPgE\n8HfAPcAbqW4ziIiJDHDXnM+1DLo0ARwCPjw0PrTZ5/N24H8DTwBnZmy7jcG40MR/R8n53MSUfkdn\nM99CY2ad58F0M+s8V1Rm1nmuqMys81xRmVnnuaIys85zRWVmneeKysw67/8D3fuJpzT9q6UAAAAA\nSUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x10ba12828>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.imshow(T)\n",
    "plt.colorbar()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "celltoolbar": "Raw Cell Format",
  "colab": {
   "name": "python4tp.ipynb",
   "provenance": [],
   "toc_visible": true
  },
  "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.6.1"
  },
  "toc": {
   "base_numbering": "0",
   "nav_menu": {
    "height": "369px",
    "width": "618.333px"
   },
   "number_sections": true,
   "sideBar": true,
   "skip_h1_title": false,
   "title_cell": "Table des Matières",
   "title_sidebar": "Sommaire",
   "toc_cell": true,
   "toc_position": {
    "height": "calc(100% - 180px)",
    "left": "10px",
    "top": "150px",
    "width": "165px"
   },
   "toc_section_display": true,
   "toc_window_display": true
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
