{
 "metadata": {
  "name": "",
  "signature": "sha256:171668d82e68d45a998b65d0e49379e70937bfcb873c2c2a1cea475a2a8d329a"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Chapter 10 : EARTH AND ROCKFILL DAM"
     ]
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 10.1 pg : 502"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "from numpy import array,float64,round\n",
      "\n",
      "\t\t\t\t\n",
      "#Given\n",
      "K = 5.E-4;             \t\t\t\t#coefficient of permeability of soil\n",
      "Bt = 6.;               \t\t\t\t#width of top of dam\n",
      "wb = 146.;             \t\t\t\t#width of base of dam\n",
      "H = 20.;               \t\t\t\t#heigth of dam\n",
      "hw = 2.;               \t\t\t\t#heigth of water in reservior\n",
      "hs1 = 4.;              \t\t\t\t#slope on upstream side\n",
      "hs2 = 3.;              \t\t\t\t#slope on downstream side\n",
      "df = 30.;              \t\t\t\t#length of drainage filter\n",
      "\n",
      "x = wb-df-72+72*0.3;\n",
      "y = 18.;\n",
      "s = (x**2+y**2)**0.5-x;\n",
      "\n",
      "x = array([0, 10, 20, 30, 40, 50, 60, 65.6],dtype=float64);\n",
      "y = (4.849*x+5.879)**0.5;\n",
      "y = round(y*1000)/1000;\n",
      "\n",
      "print \"x                       y\";\n",
      "for i  in range(8):\n",
      "    print \"%.2f          %.2f\"%(x[i],y[i]);\n",
      "sf = K*s*10000;\n",
      "sf = round(sf*1000)/1000;\n",
      "print \"Seepage flow per unit length of dam = %.2fD-6 cumecs/metre length of dam.\"%(sf);\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "x                       y\n",
        "0.00          2.42\n",
        "10.00          7.37\n",
        "20.00          10.14\n",
        "30.00          12.30\n",
        "40.00          14.14\n",
        "50.00          15.76\n",
        "60.00          17.23\n",
        "65.60          18.00\n",
        "Seepage flow per unit length of dam = 12.12D-6 cumecs/metre length of dam.\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 10.2 pg : 502"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "\t\t\t\t\n",
      "#Given\n",
      "K = 3.E-3;            \t\t\t\t#coefficient of permeability\n",
      "nd = 25.;              \t\t\t\t#number of potential drops\n",
      "nf = 4.;              \t\t\t\t#number of flow channels\n",
      "lf = 40.;             \t\t\t\t#filter length\n",
      "H = 52.;              \t\t\t\t#heigth of dam\n",
      "fb = 2.;              \t\t\t\t#free board\n",
      "\n",
      "# Calculations\n",
      "q = K*(H-fb)*nf/(nd*100);\n",
      "\n",
      "# Results\n",
      "print \"Discharge per meter length of dam = %.5f cumec/metre length.\"%(q);\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Discharge per meter length of dam = 0.00024 cumec/metre length.\n"
       ]
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 10.3 pg : 503"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "\t\t\t\t\n",
      "#Given\n",
      "x = 4.;\n",
      "#Given scale\n",
      "An = 14.4;       \t\t\t\t#area of N recmath.tangle\n",
      "At = 6.4;        \t\t\t\t#area of T recmath.tangle\n",
      "Au = 4.9;        \t\t\t\t#area of U recmath.tangle\n",
      "L = 12.6;        \t\t\t\t#length of arc;\n",
      "gamma_m = 19.;    \t\t\t\t#unit weigth of soil\n",
      "gamma_w = 9.81;  \t\t\t\t#unit weigth of water\n",
      "fi = 26.;         \t\t\t\t#effective angle(degree)\n",
      "co = 19.5;       \t\t\t\t#cohesion value\n",
      "\n",
      "# Calculations\n",
      "#consider 1m length of dam\n",
      "SumN = An*x**2*gamma_m;\n",
      "SumT = At*x**2*gamma_m;\n",
      "SumU = Au*x**2*gamma_w;\n",
      "Le = x*L;\n",
      "F = ((Le*co)+(SumN-SumU)*math.tan(math.radians(fi)))/SumT;\n",
      "F = round(F*100)/100;\n",
      "\n",
      "# Results\n",
      "print \"Factor of safety for slope = %.2f.\"%(F);"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Factor of safety for slope = 1.41.\n"
       ]
      }
     ],
     "prompt_number": 5
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 10.4 pg : 503"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "\n",
      "#check section for:\n",
      "#Stability of d/s slope against steady seepage\n",
      "#Sloughing of u/s slope against sudden drawdown\n",
      "#Stability of the foundation against shear\n",
      "#Seepage through body of dam\n",
      "\n",
      "#Given\n",
      "#Dimensions\n",
      "H = 20.;            \t\t\t\t#Heigth of dam\n",
      "Bt = 6.;            \t\t\t\t#top width of dam\n",
      "s1 = 4.;            \t\t\t\t#u/s slope\n",
      "s2 = 3.;            \t\t\t\t#d/s slope\n",
      "fb = 2.;            \t\t\t\t#free board\n",
      "#Properties of materials of dam\n",
      "gamma_d = 17.27;       \t\t\t\t#dry density\n",
      "wc = 0.15;             \t\t\t\t#optimum water content\n",
      "gamma_s = 21.19;       \t\t\t\t#saturated density\n",
      "gamma_w = 9.81;        \t\t\t\t#unit weigth of water\n",
      "wavg = 19.62;          \t\t\t\t#average unit weigth under seepage\n",
      "theta = 26.;            \t\t\t\t#average angle of internal friction(degree)\n",
      "co = 19.13;            \t\t\t\t#average cohesion\n",
      "K = 5.E-4;              \t\t\t\t#coefficient of permeability\n",
      "#properties of foundation materials\n",
      "gamma_f = 17.27;      \t\t\t\t#average unit weigth\n",
      "cof = 47.87;          \t\t\t\t#average cohesion\n",
      "fi = 8.;               \t\t\t\t#average angle internal friction\n",
      "t = 6.;                \t\t\t\t#thickness of clay\n",
      "FOSp = 1.5;           \t\t\t\t#permissible factor of safety of slope\n",
      "PS = 8.E-6;            \t\t\t\t#permissible seepage\n",
      "#(a) Stability of d/s slope against steady seepage\n",
      "An = 302.4;            \t\t\t\t#area of N diagram\n",
      "At = 91.2;             \t\t\t\t#area  of T diagram\n",
      "Au = 98.4;             \t\t\t\t#area of U diagram\n",
      "Le = 60.;               \t\t\t\t#length of arc\n",
      "SumN = An*gamma_s;\n",
      "SumT = At*gamma_s;\n",
      "SumU = Au*gamma_w;\n",
      "F = (Le*co)+(SumN-SumU)*math.tan(math.radians(theta))/SumT;\n",
      "F = round(F*100)/100;\n",
      "print \"Parta:\"\n",
      "print \"Factor of safety for slope = %.2f.\"%(F);\n",
      "print \"Safe\";\n",
      "\n",
      "#(b) Sloughing of u/s slope against sudden drawdown\n",
      "h1 = 15.;\n",
      "b = 80.;\n",
      "P = gamma_s*H**2*math.tan(math.radians(45-(theta/2)))**2/2+gamma_w*h1**2/2;\n",
      "sav = P/b;\n",
      "smax = 2*sav;\n",
      "Ne = (gamma_s-gamma_w)*b*H/2;\n",
      "R = Ne*math.tan(math.radians(theta))+co*b;\n",
      "fs = R/P;\n",
      "fs = round(fs*100)/100;\n",
      "print \"Partb:\"\n",
      "print \"Factor of safety w.r.t average shear = %.2f.\"%(fs);\n",
      "print \"Safe\";\n",
      "sr = 0.6*H*(gamma_s-gamma_w)*math.tan(math.radians(theta))+co;\n",
      "FS = sr/smax;\n",
      "FS = round(FS*100)/100;\n",
      "print \"Factor of safety w.r.t maximum shear = %.2f.\"%(FS);\n",
      "print \"Safe\";\n",
      "\n",
      "#(c) Stability of the foundation against shear\n",
      "h1 = 26.;\n",
      "h2 = 6.;\n",
      "gamma_m = (wavg*(h1-h2)+gamma_f*h2)/h1;\n",
      "l = (gamma_m*h1*math.tan(math.radians(fi))+cof)/(gamma_m*h1);\n",
      "fi1 = math.tan(math.radians(l));\n",
      "P = (h1**2-h2**2)/2*gamma_m*math.tan(math.radians(45-(fi1/2)))**2;\n",
      "sav = P/b;\n",
      "smax = 2*sav;\n",
      "s1 = cof+gamma_f*h2*math.tan(math.radians(fi));\n",
      "s2 = cof+gamma_m*h1*math.tan(math.radians(fi));\n",
      "as1 = (s1+s2)/2;\n",
      "fs = as1/sav;\n",
      "fs = round(fs*100)/100;\n",
      "print \"Partc:\"\n",
      "print \"Factor of safety w.r.t overall shear = %.2f.\"%(fs);\n",
      "print \"Safe\";\n",
      "\n",
      "gamma_av = (wavg*0.6*H+gamma_f*h2)/(0.6*H)+h2;\n",
      "s = cof+gamma_av*0.6*H*math.tan(math.radians(fi));\n",
      "fs = s/smax;\n",
      "fs = round(fs*100)/100;\n",
      "print \"Factor of safety w.r.t overall shear = %.2f.\"%(fs);\n",
      "print \"Unsafe\";\n",
      "\n",
      "#(d) Seepage through body of dam\n",
      "s = 2.;     \t\t\t\t#measured\n",
      "q = K*s*100000/100;\n",
      "print \"Partd:\"\n",
      "print \" Seepage through body of dam = %.2fD-5 cumecs/m length of dam\"%(q);"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Parta:\n",
        "Factor of safety for slope = 1149.17.\n",
        "Safe\n",
        "Partb:\n",
        "Factor of safety w.r.t average shear = 2.16.\n",
        "Safe\n",
        "Factor of safety w.r.t maximum shear = 1.24.\n",
        "Safe\n",
        "Partc:\n",
        "Factor of safety w.r.t overall shear = 1.18.\n",
        "Safe\n",
        "Factor of safety w.r.t overall shear = 0.69.\n",
        "Unsafe\n",
        "Partd:\n",
        " Seepage through body of dam = 1.00D-5 cumecs/m length of dam\n"
       ]
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 10.5 pg : 507"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "%matplotlib inline\n",
      "import math \n",
      "from matplotlib.pylab import plot,show\n",
      "from numpy import zeros\n",
      "#design upstream impervious blanket\n",
      "\t\t\t\t\n",
      "#Given\n",
      "Zb = 1.2;       \t\t\t\t#thickness of blanket\n",
      "Zf = 8;         \t\t\t\t#dismath.tance of blanket from foundation\n",
      "kb = 0.06;      \t\t\t\t#coefficient of permeability of blanket material\n",
      "kf = 72;        \t\t\t\t#coefficient of permeability of foundation soil\n",
      "Hw = 10;         \t\t\t\t#heigth of water in reservior\n",
      "Xd = 40;\n",
      "\n",
      "a = (kb/(kf*Zb*Zf))**0.5;\n",
      "Xo = 1.414/a;\n",
      "\n",
      "#we vary value of x\n",
      "x = [0, 25, 50, 75, 100, 125, 151.8, 300]\n",
      "Xr = zeros(8)\n",
      "ho = zeros(8)\n",
      "r = zeros(8)\n",
      "\n",
      "for i in range(8):\n",
      "    e = math.exp(2*a*x[i]);\n",
      "    Xr[i] = (e-1)/(a*(e+1));\n",
      "    ho[i] = Xr[i]*Hw/(Xr[i]+Xd);\n",
      "    r[i] = Xr[i]*100/(Xr[i]+Xd);\n",
      "\n",
      "print \"x                     Xr              ho         reduction qpercent\";\n",
      "for i in range(8):\n",
      "    print \"%.2f        %.2f        %.2f        %.2f\"%(x[i],Xr[i],ho[i],r[i]);\n",
      "\n",
      "#graph is plotted between r and x.\n",
      "#after around 130m length there is only slight increase in head dissipated(ho)\n",
      "plot(x,r)\n",
      "show()\n",
      "L = 130;\n",
      "print \"Thickness of blanket = %.2f m\"%(Zb);\n",
      "print \"Length of blanket = %i m.\"%(L);\n",
      "\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The history saving thread hit an unexpected error (OperationalError('disk I/O error',)).History will not be written to the database.\n",
        "x                     Xr              ho         reduction qpercent"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "\n",
        "0.00        0.00        0.00        0.00\n",
        "25.00        24.56        3.80        38.04\n",
        "50.00        46.67        5.38        53.85\n",
        "75.00        64.78        6.18        61.83\n",
        "100.00        78.50        6.62        66.24\n",
        "125.00        88.28        6.88        68.82\n",
        "151.80        95.35        7.04        70.45\n",
        "300.00        106.53        7.27        72.70\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAXUAAAEACAYAAABMEua6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAF9RJREFUeJzt3XuUVOWdr/EHuQgCcokKRG6CUTAxhmPiqPFSoPGok6hz\nczlr5iwy8cTRiUbNJONlYuxZKzlRZrlMVsw4epg4aOItMZKYZHKUhvICSlBxUBBQEkQTuhEhGpDI\npev88e62m6Karq6qXbv2ruezVq3atbsu78vWb+/+7bfeFyRJkiRJkiRJkiRJkiRJklLhOmAl8CJw\nL3AgMBp4DFgLPAqMTKx1kqSyTQZ+TQhygAeA2cAc4J+ifdcAN9W9ZZKkPhsNrAFGAQOAR4BPAauB\nMdFzxkaPJUkpcAnwB2ATcE+0b2u3n/creixJalBTgVXABwhn6g8Df8u+Ib6lzu2SJJUwoJeffxxY\nArwVPf4xcBLQRii7tAHjCGfx+5g6dWph3bp1tWmpJDWPdcCRlbzwgF5+vho4ERhCKLOcSThzf4Rw\nwZTofn7JVq1bR6FQyOztxhtvTLwN9s3+2b/s3QhVkor0dqb+38DdwLNAB/A8cCcwHHgQuBhYD1xY\naQMkSbXTW6hDGL44p2jfFsJZuySpgfRWftF+5HK5pJsQmyz3Dexf2mW9f9XoF/P7F6L6kCSpTP36\n9YMK89kzdUnKEENdkjLEUJekDDHUJSlDDHVJyhBDXZIyxFCXpAwx1CUpQwx1ScoQQ12SMsRQl6QM\nMdQlKUMMdUnKkHLmU5ck9VFHB+zYAe++G27dt3t7XA1DXVJT2bWr8rDty3N37oQhQ+Cgg7ruO2/7\ne/yBD1TXP+dTl5S4QiGEYLmBWU3YdnTA0KG9h2tfgrjUawcPhn4VJmw186kb6pJ61NEBf/xjPOFa\n/HjgwNqEaW/PHTiw8rCtl7hD/Wjg/m6PpwA3AN8HHgAm0bX49O+LXmuoSzHYvTu+M9nuj997L5xx\nxn1WO2QIDLAY/L56nqkfAPwWOAG4AthMWJT6GmAUcG3R8w11NY3OEkI9wnbPnvqc1Q4eDAc4Rq7u\n6hnqZxHO0k8FVgOnA+3AWCAPTCt6vqGuhrRrF6xbB2+8UX3ZoPt2//7xn9UedBAMGtT4JQRVrppQ\n7+sfPBcB90XbYwiBTnQ/ppIGSHEpFKC9HdauhTVr9r5t2ADjx8PEiT1fNBsxAsaNKz94O+u1UpL6\n8ptgEKH0cgzwJrCVUHLptAUYXfQaz9QVux074JVX9g3utWtDnfboo/e+HXUUTJ0KBx6YdMul0up1\npn4O8Bwh0KGr7NIGjAM2lXpRS0vL+9u5XI5cLldBM9XsOjpCqaQ4uNesgU2bYMqUENZHHw2zZsFl\nl4Xtasf8SvWQz+fJ5/M1ea++/Ca4H/gvYF70eA7wFnAz4QLpSLxQqiq9886+Z9tr1oQz8REj9j3r\nPvpomDTJkRPKlnpcKB0KvAYcAfwh2jcaeBCYiEMaVYHt22HxYli4EJ55JoT3O+90nXEXl0yGD0+6\nxVJ9+OUjpcJ778HSpSHEFy6E55+HGTNCueSTn4Tp0+Hwwx1CJxnqaki7d4fg7gzxp5+GadNCiHcG\n+bBhSbdSajyGuhpCRwe8+GJXiD/5ZBgyOGsWzJwJp50Go0b1/j5SszPUlYhCIVzI7AzxRYtg9Oiu\nM/FcDg47LOlWSuljqKtuXnutK8QXLgzfoDzjjHAmPnMmTJiQdAul9DPUFZuNG8MZeGeIb9/edSY+\na1YYH+7X1aXaMtRVU8uXw/e+B62t0NYGp5/eFeLHHGOIS3Ez1FW1PXvgZz+DW28NE11deimccw4c\nd1wosUiqn3pO6KWM2bYN7roLvv3t8JX6q6+Gv/gLJ6aS0spQb1IbNsB3vhMCfeZMuPtuOOkkSytS\n2vndvSazdClcdFH4JueePbBsGfzwh3DyyQa6lAXW1JvA7t3w8MOhXt7WBl/8Inzuc3DwwUm3TFIp\n1tRV0ttvw9y5ocwyYQJ8+ctw/vle+JSyzFDPoF//Olz4vOeeMILlRz+Cj3886VZJqgdr6hlRKIS5\nVv78z+GEE8LSaitWwA9+YKBLzcQz9ZTbuTNc6Lz11jAX+VVXhTP0oUOTbpmkJHihNKW2bIE77oDv\nfjcsInH11XDuuc5FLmVBNRdKjYCUWbMG/uEf4MgjwwyJP/95+Dr/pz9toEuy/JIaixbBLbeEceWX\nXgqrVsHYsUm3SlKjMdQb3K5d8I//GOZluf76UD8fMiTpVklqVOWG+khgLvBhoAD8HfAK8AAwiZ4X\nnlYV2tvhwgvDl4Sefx5Gjky6RZIaXblV2G8DvwCmAx8FVgPXAo8BRwGt0WPVyLJl8IlPhGlvf/IT\nA11Secq5ujoCWA5MKdq/GjgdaAfGAnlgWtFzHP1SgbvugmuugTvvhAsuSLo1kuot7mkCjgDeBO4C\njgOeA64CxhACneh+TCUNUJedO8PQxNZWePxxmD496RZJSptyQn0A8D+Ay4FlwLfYt9RSiG77aGlp\neX87l8uRy+UqaGb2tbXBX/5lmNN86VIYMSLpFkmql3w+Tz6fr8l7lXN6PxZ4mnDGDnAKcB2hHDMT\naAPGAYuw/FKRZ56Bv/or+Pzn4atfdby51Ozi/vJRG/A64YIowJnASuARYHa0bzYwv5IGNLu5c+G8\n8+Df/g2+9jUDXVJ1yv1NcBxhSOMgYB1hSGN/4EFgIj0PafRMvQc7d8KVV0I+D/Pnh6/6SxK48HTq\nbNwY6ueHHQbz5rlYhaS9OfdLiixZEsafn3suPPSQgS6ptpwmoI7uuANuuCGMQ//TP026NZKyyFCv\ng/fegyuugMWLw+1DH0q6RZKyyvJLzH772/BV/y1bwtBFA11SnAz1GD31VFha7vzzw+yKw4cn3SJJ\nWWf5JQaFAtx+O/zLv4TRLWefnXSLJDULQ73G/vjHsDLRsmVhpMvUqUm3SFIzsfxSQ6+/DqedBtu3\nw9NPG+iS6s9Qr5EnnoA/+ZPwpaL774dhw5JukaRmZPmlSoUC3HYbfP3rcM89cNZZSbdIUjMz1Kuw\nYwdcdhksXx7KLVOKlxGRpDqz/FKhDRvg1FPDF4uWLDHQJTUGQ70CmzfDySfDRRfBvffC0KFJt0iS\nAmdprMBnPwujRsGttybdEklZFPcapeqmtRUWLYKVK5NuiSTty/JLH+zYAX//92GVIocsSmpEll/6\n4PrrYd06eOCBpFsiKctc+agOVqyAM88M92PHJt0aSVnmykcx27MHLrkEvvENA11SYyv3Qul64B1g\nD7ALOAEYDTwATKLnhacz4fbbYdAguPjipFsiSftX7un9b4DjgS3d9s0BNkf31wCjgGuLXpf68svr\nr8OMGfDkkzB9etKtkdQM6lV+Kf6A84B50fY84IJKGtDorrgCLr/cQJeUDuWWXwrAAkL55Q7g/wJj\ngPbo5+3R40x5+GFYvdrRLpLSo9xQ/ySwETgUeAxYXfTzQnTbR0tLy/vbuVyOXC7X1zYm4u23w1n6\nvffCgQcm3RpJWZbP58nn8zV5r0pqNjcC24DPAzmgDRgHLAKmFT03tTX1L3wBdu2CO+9MuiWSmk3c\nNfWDgM4lk4cCZwEvAj8FZkf7ZwPzK2lAI1qyJJRebr456ZZIUt+UU34ZAzzc7fk/AB4FngUeBC6m\na0hj6u3cGcakf+tbYdIuSUoTv1Fa5BvfCAtePPII9Iv7X0eSSnCagBpZuzbMk/7cczBpUtKtkdSs\nnCagBgqFMAPjV79qoEtKL0M98p//Cdu2hWGMkpRWll+ATZvg2GPhl78MUwJIUpKsqVfpb/4GDj8c\n5sxJuiWS5HJ2VfnlL8NolxdfTLolklS9pg717dvhssvg3/8dhg5NujWSVL2mLr985SuwcSN8//tJ\nt0SSulh+qcDzz8Pdd1t2kZQtTTmkcffuMBXATTfBYYcl3RpJqp2mDPXbboPhw+Gzn026JZJUW01X\nU3/tNTj++DAT41FHJd0aSdqX0wSUqVAI86RffbWBLimbmupC6Q9/COvXw49/nHRLJCkeTVN+2boV\nPvxheOghOOmkpFsjST1zmoAyXHIJDBwI3/1u0i2RpP1znHovnngCfvELWLky6ZZIUrwyf6H0vffC\nWfp3vgMjRiTdGkmKV+ZD/ZvfhOnT4c/+LOmWSFL8yq3Z9CcsNP0G8BlgNPAAMImuRad/X+J1idbU\nX34ZTjsNli+H8eMTa4Yk9Uk9xqlfCawCOhP6WuAx4CigNXrcUDo6QtmlpcVAl9Q8ygn18cC5wFy6\nfnOcB8yLtucBF9S+adWZOzfM8XLppUm3RJLqp5zRL7cCXwEO7rZvDNAebbdHjxvGxo3wz/8MCxdC\n//5Jt0aS6qe3UP80sAlYDuR6eE6BrrLMPlpaWt7fzuVy5HI9vU3tXHllKL0ce2zsHyVJVcvn8+Tz\n+Zq8V2+F+P8D/C9gNzCYcLb+Y+AThJBvA8YBi4BpJV5f9wuljzwCX/oSrFgBQ4bU9aMlqSbq9Y3S\n04EvE0a/zAHeAm4mXCQdSemLpXUN9W3bwlQAd90Fs2bV7WMlqabqOUtjZ0LfBHwKWAvMih4n7oYb\nYOZMA11S88rM3C/LlsFnPgMvvQSHHFKXj5SkWDifOnDVVfCv/2qgS2pumThTf/NN+NCHYPNmGNAU\nU5RJyrKmP1NfuDBMB2CgS2p2mQj11lY444ykWyFJyTPUJSlDUh/qv/kNbN8exqdLUrNLfah3nqX3\ni/uSrySlQGZCXZKU8iGNHR0wblz44tHEibF9jCTVVdMOaXzpJTj4YANdkjqlOtQtvUjS3gx1ScqQ\n1NbUd+0K87ysW+d8L5KypSlr6r/6FUydaqBLUnepDXVLL5K0r1SH+plnJt0KSWosqaypb98OY8ZA\nezsMHVrzt5ekRDVdTf3JJ+H44w10SSqWylC3ni5JpfUW6oOBpcALwCrgm9H+0cBjhIWnHwVGxtXA\nUhYssJ4uSaWUU7M5CHgXGAA8BXwZOA/YDMwBrgFGAdeWeG3Na+qbN4ehjJs3w8CBNX1rSWoIcdfU\n343uBwH9ga2EUJ8X7Z8HXFDJh1di0SI49VQDXZJKKSfUDyCUX9qBRcBKYEz0mOh+TCytK8HSiyT1\nrJylmjuAjwEjgP8HzCz6eSG6ldTS0vL+di6XI5fL9bWNe2lthcsvr+otJKmh5PN58vl8Td6rrzWb\nG4AdwP8GckAbMI5wBj+txPNrWlN/7TU44QRoa3OlI0nZFWdN/RC6RrYMAT4FLAd+CsyO9s8G5lfy\n4X3l0nWStH+9lV/GES6EHhDd7gFaCcH+IHAxsB64ML4mdlmwwPHpkrQ/qZkmoFAIS9c98wxMnlyT\nt5SkhtQU0wSsXBmmBTDQJalnqQl1hzJKUu9SE+rO9yJJvUtFTb1z6bpXX4VDD61BqySpgWW+pr5s\nGUyZYqBLUm9SEeqWXiSpPIa6JGVIw9fUt2+HsWNh40YYNqxGrZKkBpbpmvpTT8GMGQa6JJWj4UPd\n0osklc9Ql6QMaeia+ltvhaGMLl0nqZlktqa+aBGccoqBLknlauhQt/QiSX3T0KHuJF6S1DcNG+ob\nNsDbb8NHPpJ0SyQpPRo21FtbYdYsOKBhWyhJjadhI9N6uiT1XUMOaSwU4IMfhCVL4IgjYmiVJDWw\nuIc0TgAWASuBl4AvRvtHA48Ba4FHgZGVNKCUVatgyBADXZL6qpxQ3wVcDXwYOBH4AjAduJYQ6kcB\nrdHjmrD0IkmVKSfU24AXou1twMvA4cB5wLxo/zzgglo1yqGMklSZvtZsJgOPAx8BNgCjur3Plm6P\nO/W5pr57d1i67pVXXOlIUnOqpqY+oA/PHQY8BFwJ/KHoZ4Xoto+Wlpb3t3O5HLlcbr8f8uyzMGmS\ngS6peeTzefL5fE3eq9zfBAOBnwH/BXwr2rcayBHKM+MIF1OnFb2uz2fqX/86bN0Kt9zSp5dJUmbE\nPfqlH/AfwCq6Ah3gp8DsaHs2ML+SBhTzIqkkVa6c3wSnAE8AK+gqsVwH/Ap4EJgIrAcuBH5f9No+\nnam/+y6MGePSdZKaW9w19afo+Yy+pmNUFi+G444z0CWpUg01TYBDGSWpOg0V6tbTJak6DTP3y5Yt\nMHlyWLpu0KB4GyVJjSwTy9l1Ll1noEtS5Rom1C29SFL1DHVJypCGCPXXXw819Y9+NOmWSFK6NUSo\nu3SdJNVGQ8SopRdJqo3EhzR2Ll331FMwdWrMrZGkFEj1kMaXX4bBg2HKlKRbIknpl3iod5Ze+sX9\nN4MkNYGGCXVJUvUSral3Ll23Zk2YcleSlOKa+nPPwcSJBrok1UqioW7pRZJqK9FQX7DAUJekWkqs\npr5jBxx2GPzudzB8eMytkKQUibum/j2gHXix277RwGPAWuBRYGRfP3jx4jDXi4EuSbVTTqjfBZxd\ntO9aQqgfBbRGj/vE0osk1V45of4ksLVo33nAvGh7HnBBXz+4tdX1SCWp1iq9UDqGUJIhuu/ToMSt\nW8PY9BNPrPDTJUkl1WL0SyG6lS2fh5NPduk6Saq1ARW+rh0YC7QB44BNPT2xpaXl/e1cLkcul2PB\nAksvktQpn8+Tz+dr8l7lDpmZDDwCHBs9ngO8BdxMuEg6ktIXS0sOaZw2De67D2bM6GtzJSn7qhnS\nWM6L7gNOBw4hnKF/DfgJ8CAwEVgPXAj8vsRr9wn1N96Aj30MNm1ypSNJKqWaUC+n/PLXPeyvqIDS\n2gozZxrokhSHukerQxklKT51nSagUIDx4+Hxx+HII2P+ZElKqdRMvbt6NQwY4FqkkhSXuoZ6Z+nF\npeskKR51D3Xne5Gk+NStpr5nT1i67uWXYezYmD9VklIsFTX1556Dww830CUpTnULdYcySlL86hrq\n1tMlKV51qanv2AGHHhqWrjv44Jg/UZJSruFr6kuWhKXrDHRJilddQt3SiyTVR11C3fVIJak+Yq+p\nb91aYMIE2LwZDjww5k+TpAxo6Jp659J1BrokxS/2ULeeLkn1E3uoW0+XpPqJvaY+enSBTZugf/+Y\nP0mSMqKha+ozZxroklQv1Yb62cBq4BXgmlJPsPQiSfVTTaj3B24jBPsxhAWqpxc/KcuTeOXz+aSb\nEJss9w3sX9plvX/VqCbUTwBeBdYDu4D7gfOLn5TltUiz/B9WlvsG9i/tst6/alQT6ocDr3d7/Ea0\nby8uXSdJ9VNNqBdq1gpJUk1Ucx59ItBCqKkDXAd0ADd3e86rwNQqPkOSmtE6oO7F6wHRB08GBgEv\nUOJCqSQpPc4B1hDOyK9LuC2SJEmSetPrl5JSaD2wAlgO/CraNxp4DFgLPAqMTKRllfke0A682G3f\n/vpzHeF4rgbOqlMbq1Gqfy2EUVrLo9s53X6Wpv5NABYBK4GXgC9G+7Ny/HrqXwvZOH6DgaWEkvUq\n4JvR/oY9fv0J5ZjJwECyU2v/DeEfvbs5wD9F29cAN9W1RdU5FZjB3qHXU3+OIRzHgYTj+ip1XLS8\nQqX6dyPwpRLPTVv/xgIfi7aHEUqg08nO8eupf1k5fgAHRfcDgGeAU6jR8Yuj42V9KSmlikcLnQfM\ni7bnARfUtzlVeRLYWrSvp/6cD9xHOJ7rCcf3hPibWJVS/YPSI77S1r82wv/kANuAlwnfEcnK8eup\nf5CN4wfwbnQ/iHAivJUaHb84Qr2sLyWlUAFYADwLfD7aN4bwJz7R/ZgE2lVLPfXng4Tj2CnNx/QK\n4L+B/6Drz9s0928y4S+SpWTz+E0m9O+Z6HFWjt8BhF9c7XSVmmpy/OII9ax+KemThP+4zgG+QPjz\nvrsC2ep7b/1JY19vB44g/Gm/EbhlP89NQ/+GAQ8BVwJ/KPpZFo7fMOBHhP5tI1vHr4PQj/HAacDM\nop9XfPziCPXfEi50dJrA3r9l0mpjdP8m8DDhz592Qv0PYBywKYF21VJP/Sk+puOjfWmzia7/WebS\n9SdsGvs3kBDo9wDzo31ZOn6d/fs+Xf3L0vHr9Dbwc+B4Gvj4ZfFLSQcBw6PtocBiwhXoOXSN7rmW\ndF0ohXCMii+UlupP54WaQYQzpXXEv8BKLUxm7/6N67Z9NXBvtJ22/vUD7gZuLdqflePXU/+ycvwO\noat0NAR4AjiDBj9+WftS0hGEf9QXCEOsOvs0mlBnT+OQxvuA3wE7CddA/o799+d6wvFcDfzPura0\nMsX9+xwhKFYQarLz2fsaSJr6dwrhz/cX6BredzbZOX6l+ncO2Tl+xwLPE/q3AvhKtD8rx0+SJEmS\nJEmSJEmSJEmSJEmSJElSI/r//MLloWR9quQAAAAASUVORK5CYII=\n",
       "text": [
        "<matplotlib.figure.Figure at 0x10e9d5350>"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Thickness of blanket = 1.20 m\n",
        "Length of blanket = 130 m.\n"
       ]
      }
     ],
     "prompt_number": 1
    }
   ],
   "metadata": {}
  }
 ]
}