example of interpolation for Julia

interpolation_test.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "# Julia 内挿の例"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import Interpolations"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import Winston;"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import Color;"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### tとfの列を与えて内挿"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "t=sin((1:10)/20*pi);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "f=t.*t;"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAAEsCAIAAADfNCTgAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3de1xU1d4/8DVcZABFLg6CclWGQoQjJf0sjpkKx+Q8vSyLUoGnrI6VgQakmYShmdIFtcyjpSkdLEU8PnU8gjcIL0M+KAeJmzPcGa4Od5gBZJh5/pgOPw7mba/N7M3web/6Y/aevdb6rpevPuzL7L0FWq2WAAAAU0ZcFwAAMLohRgEAqCBGAQCoIEYBAKggRgEAqCBGAQCoIEYBAKggRgEAqCBGAQCoIEYBAKggRgEAqCBGAQCoIEYBAKiYcF3Af6ipqbl27RrXVQDAWOHl5eXl5UXZCb/2RteuXXufW6pUqv3799MM1NbWxqytQqH44YcfGA/9xRdfMG5bVVX1008/MW7+1ltvMW5bVFR07tw5xs1pZp2Tk5Odnc24Oc2ss7Ky8vPzGTenmfWpU6fKysoYN6eZdWpqan19PbO2arV6z549jIdOSkrq6Ohg1ra2tnbdunUP2mr9+vXMhvsPWj4JDQ29zy3b29tXr17NeKAZM2bU19cza1tVVbVx40bGQ9//HG93/fr1Tz/9lHFzNzc3xm2zsrK++eYbxs1pZv3jjz8eO3aMcXOaWSclJZ09e5Zxc5pZ79q1Kycnh3Fzmllv3rxZKpUya9vX17dy5UrGQ0dFRTU1NTFrW1RU9MgjjzxoK5p/o0H82hsFABh96JOYRff/l0GtVldXVzMeKCwsrL29nVnbvr6+2tpaxkNXVFQwbqtSqRobGxk3X7JkCeO2XV1dCoWCcXOaWbe1tbW2tjJuTjPr5ubmjo4Oxs1pZt3U1KRUKhk3p5l1fX19b28vs7YajaayspLx0HK5vL+/n1nburq6v/zlLw/aipW9UX5dYrp/xsbGLi4ujJuPHz/eyIjhnvi4ceOmTp3KeGh3d3fGbc3Nzc3NzRk3Hz9+PE1bmuY0s7a2tmbcltDN2s7OjmZomlnb29vTDE0za0dHR8ZtBQKBm5sb4+ZOTk6M25qYmNDMmgYO6gHAQPSoezgZV68xGh0dvWXLFn2OeCcrVqwQCoVcV6Fvr776KtclcACzHiNyW3OvWl/lZGg9xahEItm0adPu3buVSqV+Rry7uXPnmpqacl2Fvi1YsIDrEjiAWY8Fao16Xea6XOPcH2/8qP/R9RSjZ8+eLS0tpTnnAgBwJ/uu7StSFIksRO+efbdvoE/Po+spRjdv3nzkyJE5c+boZzgAGDvaetvis+IXui9MCUkpbyvf+ctOPReAS0wAMLrFZcZ19HXsenrXfLf5zz383MeXPq7vYngXFjP8+sGTTCbT/UbEwsLCz89Pt/LQoUNj8DwmANxdr7o3vyl//LjxX+d+/dbst2bazySEfP6nz73/6h2bGXtoyaHbm0RERLS3txNCSkpKWltbCSEWFhb0lfArRj09PXNycriuAgBGgZ1XdibnJztOcDQWGD/z0DO5Dbm69c94PvO3/L+t9l/tP8V/WJOvvvpq2JqwsDD6SvgVowAA96O+q37bpW3dt7pLmksIIX9K/tOwDdamr5W8KhEIBHooBjEKAKPPxoyNao16jtOcq3VXDy89PMli0u3b9Kh7LExZOGa/J73G6LFjx/Q5HAAYpNyG3ORfk9//4/thvmG+e30vVl/865//ymE9uFIPAKOJVqtdm77W3tJ+fcD6hyc9/PZjb3+T+01+E/PHwtJDjALAaHK44LBELvk06FMrMytCyIfzPrQ1t406HcVhSTg3CgCjhrJf+f7598ePG1/XWfeJ5BPdSi+R189VP58oObHUayknVSFGAWDU2H5pe11XHSHk/Yz3h3217ty6YHGw0ISDRw4hRgFg1Ni6YOvWBVu5rmI4nBsFAKDCr73R3t7eixcvEkIsLS0fffRRrssBAMNUVFTU0tJCCGHl0Z38itH+/v4bN24QQmxtbRGjADBCamtrq6urCSF9fSw8VY9fMTphwoRVq1ZxXQUAGLhFixbpPugOfynh3CgAABXEKAAAFcQoAPCIVqvluoQHhhgFAL4oay3bcpEXLw9+IIhRAOCLtafXJlxOqGyv5LqQB4MYBQBeSC9LTytN61P3rTu7jutaHgxiFAC4p9ao151d52HrsWnepr+X/P1cxTmuK3oAiFEA4N7unN1FiqIvF3+5ce5GD1uPqNNRao2a66LuF2IUADjW2tO69eLWoGlBiz0WjzMelxCYUKQoOvCvA1zXdb/4dRdTS0tLVFQUIcTR0XH9+vVclwMA+rAxY2NXX9fu4N26xee9ng+aFhSbGRviHWJnbjcSIx44cKCoqIgQ0tDQQN8bv2LU2tp648aNhBATE34VBgAj5NemX3U7no8feHxwZd9An6pftfXi1p2Ldo7EoC+++KLubvrIyEj63viVVsbGxiKRiOsqAEB/chtyX3vktd/9qn+gv1fdOxJPYraystJ9YGWPjV8xCgBjzcpZK1fOWsl1FVRwiQkAgApiFACACmIUAIAKYhQAgApiFACACmIUAIAKYhQAgApiFACACr9+fq9QKMLCwgghTk5OCQkJXJcDAIYpMTExLy+PECKXy+l741eMikSiw4cPc10FABi4mJgY3QfdfhslHNQDAFBBjAIAUEGMAgBQQYwCAFBhLUZlMtmyZcucnJyCg4NlMtntG+Tk5Lz88stTpkxZvny5RCJha1wAAG6xE6NqtTo0NNTV1VUikQQEBAQGBqpUqqEbaDSa5cuXe3t7X7582dPTc/ny5f39/awMDQCc6xvo47oELrETo1lZWVVVVdu3b3d1dY2Nje3q6kpLSxu6QXFxsVwuj4yMnDZt2rvvvltXV1dWVsbK0ADALY1W8+fv/9zY3ch1IZxhJ0alUqm/v7+R0W+9eXp6FhcXD93A29vbw8Pj7bffTktLe+ONN3x8fLy8vFgZGgC49bf8v2VUZmz6eRPXhXCGnRgtLy+3sbEZXHRzc+vs7By6gUAgCA8PT0pKWrNmzZEjR5566ilWxgUAbnXd6tqYsdHM2Oxg3sG8xjyuy+EGO3cxCYXClpaWwUWVSiUWi4ducPLkyc2bN0ulUrFYXFpaOnPmzJCQkICAgGH9dHd3//DDD4QQMzMzd3d33Uo/Pz+BQMBKnQDAro8vftykbDoXfi70RGhkWuSllZf4/H9rQUGB7qpMbW1td3c3IaStrY2+W3Zi1N3dPTU1dXCxoqLC19d36AbZ2dlPPPGELlvFYvGTTz556dKl22O0o6MjJSWFEGJpaTl79mzdSl9fX7xvGYCHKtoqdl3Z9cqsVxa4L9gyf8uqk6tSi1Nf9H6R67ruKDs7W6lUEkIKCwt1e35D9/+Y07Lh5s2bZmZm6enpGo3mwIEDtra23d3dWq02NTVVJpNptdrk5GRzc/OSkhKtViuXyy0tLaVS6e39hIaGslIPAOjBkiNLJmybUN9Vr9VqBzQDs7+Z7bzDWXlLyXVdD4CVzGHn3KhIJEpOTl66dKmzs3NMTExKSoqlpSUhJDIyMjMzkxASGhr65ptvzps3TywWu7m5vfbaa56enqwMDQCcyKjM+En605r/t0ZoImzrbevo64h7Mk7eKU/8JZHr0vRNoNVq2eqrt7e3tLTUy8vrTsfgGo0mPz/f0dHRwcHhdzcICwvDE54A+E+tUft97Vd4s/D2ryxNLW9E3HCyctJ/VQywkjlsnnMUCoU+Pj532cDIyMjPz4/FEQGAE7WdtWG+d3zEXHVH9WiJUVbg0g0APDA3a7f3At7jugq+wKNJAACoIEYBAKggRgEAqCBGAQCoIEYBAKjw60p9W1tbXFwcIWTy5MkRERFclwMAhunw4cNSqZQQ0tTURN8bv2J0woQJr776KiFk3LhxXNcCAAYrKChI90wPXZhS4leMmpiYDD7YCQBghEyePFn3gZU9NpwbBQCgghgFAKCCGAUAoIIYBQCgghgFAKCCGAUAoIIYBQCgghgFAKCCGAUAoMKvu5gUCkVYWBghxMnJKSEhgetyAMAwJSYm5uXlEULkcjl9b/yKUZFIhFfaAcBIi4mJ0X3Q7bdRwkE9AAAVxCgAABXEKAAAFcQogCFr7G6Utci4rsLAIUYBDFlsZmzUmSiuqzBwiFEAg5XbkJt0PSmtNC2tNI3rWgwZYhTAMGm12rXpa+0t7R+yeyj6THS/pp/rigwWYhTAMB0pPCKRS7Yv3P5V8FfSFunu/93NdUUGCzEKYIB61D0bMzY+4vjIf//hvwOnBf5Z/OfNFzY3djdyXZdh4leM9vf3y2QymUxWXV3NdS0Ao1jC5YSajppdT+8yEhgRQnY9vatP3fdh1odc18UX9fX1uqjp6+uj741fN4MqlcqjR48SQuzt7d98802uywEYlWo6aj6TfOY80Tlbnp0tz9at9LD1+PZf3741+61ZDrO4LY8PLl68KJPJCCEdHR30vQm0Wi19L2wJCwvDPfUAlJYdX5ZSlPK7X81znZf1SpZ+y+E1VjKHX3ujAEDv6AtHj75wlOsqxhB+nRsFABh1EKMAAFQQowAAVFiLUZlMtmzZMicnp+DgYN0lsGHa2tpWrVrl6Ojo4uKyY8cOtsYFAOAWOzGqVqtDQ0NdXV0lEklAQEBgYKBKpRq2TVRUVF9fX0FBwc6dOzdt2sTKs/sBADjHzpX6rKysqqqq7du3GxkZxcbGfv7552lpaS+88MLgBuXl5d99911tbe2kSZOef/55c3PzW7dusTI0AAC32IlRqVTq7+9vZPTbvq2np2dxcfHQDW7cuDF9+vTs7Ozjx4/PnDkzPDzczc2NlaEBALjFzkF9eXm5jY3N4KKbm1tnZ+fQDcrKyhQKxZ49e4KCgvLz8728vAoKClgZGgCAW+zsjQqFwpaWlsFFlUolFouHbtDb29vT05OamioSiV5//fXHH3/86NGjPj4+w/qRy+XBwcGEEAsLC39/f93KmJgYExPcJgAAtPbs2dPd3U0IKSws1EXW0OBijJ14cnd3T01NHVysqKjw9fUduoGLi8ukSZMmTZqkW3R2djY3N7+9H1tb2/DwcEKImZmZi4uLbqWxsTErRQLAGLdw4cL+/n5CyKxZs5RKJSHk0KFDLPSrZcPNmzfNzMzS09M1Gs2BAwdsbW27u7u1Wm1qaqpMJtNqtc3NzVZWVt9//71Wq5VKpdbW1pWVlbf3Exoayko9AAD3g5XMYefcqEgkSk5OXrp0qbOzc0xMTEpKiqWlJSEkMjIyMzOTEGJnZ5eUlPTWW2+5uLjMmDEjMjISl5gAwDCwds4xJCTkmWeeKS0t9fLyGjyV2dDQMLjBc889t3jxYplM5uLiYm1tzda4AADcYvPSjVAovP2q0bANhp0zBQAY7XBPPQAAFcQoAAAVxCgAABXEKAAAFcQoAAAVxCgAABV+3auuUCjCwsIIIU5OTgkJCVyXAwCGKTExMS8vjxDCyoOP+RWjIpEIL1gGgJEWExOj+6Dbb6OEg3oAACqIUQAAKohRAAAqiFEAACqIUYARV3ATr8wxZIhRgJF1uuz04sOLlf1KrguBkYIYBRhB/Zr+6DPRdV11CZfxO2iDhRgFGEF/vfrXkuYSH3ufxOzEqvYqrsuBEYEYBRgprT2tH134KHBaYHpYurGR8fpz67muCEYEv2J0YGCgvr6+vr5eoVBwXQsArdjM2I6+jl1P75o6Yep7Ae+lFqdeqL7AdVFACCGtra26qNG9KJSSQKvV0vfCluDg4JkzZxJCHBwcoqOjuS4HgLkiRdGsfbMiHovYuWgnIaRX3eu1x8taaH1t1TVjAd4ZzrFDhw6VlJQQQnJzczMyMih749c99ba2tp9++inXVQCwYE36mgHNgMN4h29yv9Gt8Z/in1qc+u2/vl316Cpua4OVK1fqPrByTz2/YhTAMPzPjf/JrMwkhGw4v2HYV3E/x70086WJZhO5qAtGBGIUgH0BzgHla8rv9K2pkak+i4GRhhgFYJ+9pb29pT3XVYCe8OtKPQDAqIMYBQCgghgFAKCCGAUAoIIYBQCgghgFAKDCrx88dXZ27tmzhxBia2u7fPlyrssBAMN06tSpqqoqQkhLSwt9b/yKUTMzs0ceeYQQYmFhwXUtAGCwpk+fbmtrSwihv6Ge8DBGH3/8ca6rAAAD9/DDD+s+sLLHhnOjAABUEKMAAFQQowAAVBCjAABUWItRmUy2bNkyJyen4OBgmUx2p806OzvFYnFSUhJb4wIAcIudGFWr1aGhoa6urhKJJCAgIDAwUKVS/e6Wq1evrq6u5tWbSwAAaLATo1lZWVVVVdu3b3d1dY2Nje3q6kpLS7t9s8OHD7e2tj722GOsDAoAwAfsxKhUKvX39zcy+q03T0/P4uLiYdtUVFTEx8cfPHhQIBCwMigAAB+wE6Pl5eU2NjaDi25ubp2dnUM3UKvVK1as2LFjh4ODAysjAgDwBDt3MQmFwqG3pqpUKrFYPHSDbdu22djYTJs2rbCwUKVS1dXVVVdXu7q6DutHJpO5u7sTQiwsLPz8/HQrDx06ZGqKd9cAAK2IiIj29nZCSElJSWtrK2HpLiZ2YtTd3T01NXVwsaKiwtfXd+gGXV1djY2N4eHhhJCysrL9+/crFIovvvhiWD+enp45OTmslAQAMMxXX301bA0rL1hm56D+2Weflcvlp0+f1mq13377bWNjoy5Gjx8/XlpaSgj57LPP8v7N398/Pj7+9gwFABiN2NkbFYlEycnJS5cutbW17e7uPn78uKWlJSEkMjIyPj5+2AE+AIAhYe0JTyEhIc8880xpaamXl5eJyW/dNjQ03L5lZmYmW4MCAHCOzQflCYVCHx8fFjsEAOA/3FMPAEAFMQqGpqajhusSYGxBjIJBqeuqm3NgTmN3I9eFwBiCGAWD8t659xq6G+J+juO6EBhDEKNgOLLl2T8U/OA60fVg3sG8xjyuy4Gxgl8xqtFoOjs7Ozs7lUol17XAKKPRaqLORE21mnrl9SuTLSdHpEXgeYxwJz09PbqoGRgYoO+NX28GbW1t3bBhAyFkypQpH3zwAdflwGiSdD0ppy7nh+d/cBjvsGX+lr+c/Mux4mMveb/EdV3AR0lJSQUFBYSQxkYWTqMLePUXOyws7PDhw1xXAaNP162uh3Y/5GbtJnlVIhAINFrNnANzGrobpBFSC1MWnj0BhoqVzOHXQT0AMx9f/Lihu2HJw0syKjPOV5zPrMz8L8//qu2s/Tz7c65LA8PHr4N6AAbK28p3XdlFCNlwfsOwrz6RfPKq36tOVk5c1AVjBWIURr3x48ZLXpPc6VszEzN9FgNjEGIURr3JlpMnW07mugoYu3BuFACACmIUAIAKYhQAgApiFACACmIUAIAKv67Ud3d3f/fdd4QQa2vrJUuWcF0OABimzMxMuVxOCGlra6PvjV8xamxsPGXKFELI+PHjua4FAAyWra2t7qEkpqam9L3xK0bNzc2DgoK4rgIADNysWbN0H3SHv5RwbhQAgApiFACACmIUAIAKYhQAgApiFACACmIUAIAKYhQAgApiFACACmIUAIAKv+5iUigUYWFhhBAnJ6eEhASuywEAw5SYmJiXl0cI0d1ZT4lfMSoSifCCZQAYaTExMboPuv02SjioBwCgghgFAKCCGAUAoIIYhQew79o+Vb+K6yoA+AUxCverqr0q+kz0J5JPuC4EgF9Yi1GZTLZs2TInJ6fg4GCZTHb7BtevX583b561tbVYLN67dy9b44LerDu3rkfd85nks6r2Kq5rAeARdmJUrVaHhoa6urpKJJKAgIDAwECV6j8O/W7durVw4cLFixcXFhbGxsZGRESUlZWxMjTox4XqC8eLj6/2X20kMNpwfgPX5QDwCDsxmpWVVVVVtX37dldX19jY2K6urrS0tKEb5OXlmZqabtiwwcnJ6ZVXXnn88cdTUlJYGRr0QKPVxJyJcbd2T/xT4vqA9SlFKReqL3BdFABfsBOjUqnU39/fyOi33jw9PYuLi4du4OHhkZGRofusVCqLiormzJnDytCgB9/kfpPbkPvZnz4TmgjXBaxzs3Z75/Q7A9oBrusC4AV2YrS8vNzGxmZw0c3NrbOzc+gGdnZ23t7ehJCioqK5c+fOnj17wYIFv9uVWq1Wq9UDAwOaf2OlQmCss68zPit+vtv8572eJ4SYm5h/EvjJ9cbrh/IOcV0awIPRarW6VBkYGNBFjVarpe+WnZtBhUJhS0vL4KJKpRKLxcO20Wq1H3zwwd69e+Pi4tauXSsQCG7vRyaTPfTQQ4QQCwsLPz8/3cpvv/2WlZegAjPxWfFNyiY7C7ug5P//0lYjgVFsZmyId8hEs4kc1gbwQNauXdve3k4IKSkpaW1tJYQIhUL6btmJUXd399TU1MHFiooKX1/foRtoNJo33nijpqamuLjYwcHhTv14enrinnpeUagUv9T+8qjjo4SQtp62wfV+Dn6EkNSi1NcfeZ2z4gAe0JdffjlsDSv31LMTo88++2xkZOTp06cXLVp08ODBxsZGXYweP378D3/4g1gsvnz58j/+8Y+rV68KhULdXwMzMzNzc3NWRoeRI7IQ/fLaL1xXAcBr7JwbFYlEycnJS5cudXZ2jomJSUlJsbS0JIRERkZmZmYSQi5dunTz5k1XV1ebf/vggw9YGRoAgFsCVs6w6vT29paWlnp5eZmYMNzJDQsLw0E9AOgNK5nD5vNGhUKhj48Pix0CAPAf7qkHAKCCGAUAoIIYBQCgghgFAKCCGAUAoMKvN4Mqlcq///3vhBArK6ugoKB7bg8AwEB2dnZDQwMhpKOjg743fsUoIUT3mKjBh0UBALDOyMhIFzK/+3CPB8WvGLW0tHzuuee4rgIADNzggzp1h7+UsNMHAEAFMQoAQAUxCgBABTEKAEAFMQoAQAUxypmuW11r0tdwXQUA0EKMcmbbpW27c3ZnVWVxXQgAUEGMcqO8rXznLzsJIXhTMcBohxjlxrtn3x1nPG734t35TfkH/nWA63IAgDl+3cWkUCh0L+pzcnJKSEjgupyRklmZ+eONHxMCEyIei/iH9B8bMzaGzAixNbflui6AsSIxMTEvL48QIpfL6Xtj811M9MbCu5gGtAN++/yU/crit4vNjM2KFEWz9s2KfCxyx6IdXJcGMOawkjk4qNe3fdf2Fdws2LFoh5mxGSHEW+T9mt9ru3N2FymKuC4NAJjg10G9wWvrbfvw5w8FRBB9Jjr6TLRuZU9/j1qjXn9u/akVp7gtDwAYQIzq1bX6a8/PeP5O39Z31U+ZMEWf9QAAPcSoXgVNCwqahsdRAxgUnBsFAKCCGAUAoIIYBQCgghgFAKCCGAUAoIIYBQCgwq8fPDG4p76xu9HOws7UyHSESwMAw8HuPfX8ilGRSPSg97dGnYl6bOpjUXOiRqgkADA8MTExug+6/TZKo/ugXiKXpBSmbLmwRaFScF0LAIxRozhGNVrN2vS1DuMdlLeUcZlxXJcDAGPUKI7RQ9cP5Tbk7nx6Z8RjEfv/tT+3IZfrigBgLBqtMdp1qysuM+4J5ydenPFi/FPx9pb275x+h1fPTgWAMUJPMSqTyZYtW+bk5BQcHCyTyeg7/OjCR03Kpl1P7xIIBFZmVvFPxV+uuXzixgn6ngEAHog+YlStVoeGhrq6ukokkoCAgMDAQJVKRdNhWWvZl//75RynOR29Hecrzp+vOO9m7SayEK07u65X3ctW2QAA90MfMZqVlVVVVbV9+3ZXV9fY2Niurq60tDSaDt89+27fQF+2PDsoOUj339OHn1aoFJXtlYm/JLJVNgDA/dDH70alUqm/v7+R0W+R7enpWVxczLg3jVYTNy8ubt7vX5o3NzFn3DMAAAP6iNHy8nIbG5vBRTc3t87OTsa9GQmMHnV8lI26AABYoI8YFQqFLS0tg4sqlUosFv/ulm1tbdu2bdM18fb21q0MCgoa3JNly9dff/3yyy8LhUJ2u+W5Xbt2vfPOO1xXoW+Y9RjR0dFx4sSJlStX3mWbrKysvr4+Qkh5eXl7ezshpKmpiX5ofcSou7t7amrq4GJFRYWvr+/vbtnf319TU0MIsbCwEIlEupUj8TOm69evr1ixYqzF6LVr17gugQOY9RjR19dXUFBw922am5t7enoIIbW1ta2trbpW9EPrI0afffbZyMjI06dPL1q06ODBg42NjXeKUXt7+3379umhJAAYg1544YVha0bNPfUikSg5OXnp0qXOzs4xMTEpKSmWlpaUffb19Z0+fZpx85qamlu3bjFr29HRkZWVxXjon376iXHbmzdvXrlyhXHz2tpaxm1ramquX7/OuDnNrGUyWUlJCePmNLP+9ddfKysrGTenmfXVq1cbGhoYN6eZ9aVLl3R7agxoNJp//vOfjIc+f/48419D9vT06I5l9U9PP78PCQlpbW1NT09vbm4ODAyk77C3t/fkyZOMm1dVVTGO0fb29nPnzjEeeuj5jQfV0NBw6dIlxs2rq6sZt62srLx69Srj5jSzLikpKSwsZNycZtZ5eXllZWWMm9PMOjs7myYKaWb9888/Nzc3M2urVqtPnGB+F0xaWlp3dzeztkqlkuZvHg39PShPKBT6+PjobTgAAP0Q8Oo+9IiIiPvcSVSr1TKZbMaMGcwGysjI+OMf/2hmZsagbU9Pj1wu9/T0ZDb0r7/+eqdTw/fU1dXV3Nzs7u7OrPmZM2cWLVrErG1ra6tSqXR2dmbWnGbWN2/e1Gg0Dg4OzJrTzLqurk4oFNrZ2TFrTjPr6upqa2vriRMnMmtOM+uysjJHR0dmZ940Gk1RURHjHaYbN25MmzZt3LhxDNoqlcqcnJz58+c/UKv+/v5Dhw4xGG4ofsXowMAAKxfOAADuh6mpqakp7bsz+BWjAACjzmh9UB4AAE8gRgEAqCBGAQCoIEYBAKggRgEAqBh+jN79/SXXr1+fN2+etbW1WCzeu3cvJxWOhPt5a0tnZ6dYLE5KStJvaSPlnlNua2tbtWqVo6Oji4vLjh079F/hSLjnrHNycl5++eUpU6YsX75cIscx2K4AAAQxSURBVJHov0K9iY6O3rJlCydDG3iM3v39Jbdu3Vq4cOHixYsLCwtjY2MjIiJo7vzjj/t8a8vq1aurq6sN4xdv9zPlqKgo3UOAdu7cuWnTJrlczkmpLLrnrDUazfLly729vS9fvuzp6bl8+fL+/n6uqh05Eolk06ZNu3fvViqV3FSgNWjnzp2bNGnSwMCAbtHa2jo1NXXw2ytXrkyePHlwMSAgYOvWrfoucQTcfdY6ycnJixcvDggIOHjwoN4LZN89p6z7A1lbW6tbPHXqVFlZmb6rZNs9Z11QUGBqaqpSqbRabWdnp5GRUXFxMQeFjrBNmzYtW7bM2dl5/fr1nBRg4Hujd39/iYeHR0ZGhu6zUqksKiqaM2cOB1Wy7Z5vbamoqIiPjz948KBAIOCiQPbdc8o3btyYPn16dnb2Sy+99NFHH82YMWP69OlcVMqme87a29vbw8Pj7bffTktLe+ONN3x8fLy8vLiodGRt3rz5yJEjHP7Pa+Axevf3l9jZ2emesV9UVDR37tzZs2cvWLCAgyrZdvdZq9XqFStW7Nixg/Fd6jx0zxfVlJWVKRSKPXv2BAUF5efne3l53fMRv/x3z1kLBILw8PCkpKQ1a9YcOXLkqaee0neJY4OBx+g931+i1WpjY2Pnzp0bHh5+5swZw9g7u/ust23bZmNjM23atMLCQpVKVVdXR/NQNZ645z90b29vT09Pamrq66+/fvz48VmzZh09elTvZbLsnrM+efLk5s2bpVJpWVmZTCbbu3evYV9l4oqBx6i7u3t5efng4rD3l2g0mlWrVl27dq24uDgqKor1Nz5x5e6z7urqamxsDA8PDw8Pl8lk+/fvN4DL1nefMiHExcVl0qRJkyZN0i06Ozubm4/6l8jec9bZ2dlPPPGELlvFYvGTTz5J87xauCNOzsjqzc2bN83MzNLT0zUazYEDB2xtbbu7u7VabWpqqkwmu3Dhgr29fXV1ddu/6U7Gj3Z3n/XQLefPn28Yl5juNGXtv2fd3NxsZWX1/fffa7VaqVRqbW1dWVnJZcVsuOc/dHJysrm5eUlJiVarlcvllpaWUqmU46JHTEhICFeXmAw8RrVa7bFjx8zNzadOnTpx4sRz587pVjo4OOzbt2/r1q3D/qhER0dzWy1b7jLroZsZTIxq7zBl7ZBZnzhxwsrKytnZ2djYOC4ujrtK2XT3f2iNRhMVFWVvb+/h4WFsbLxmzRpOix1ZHMbomHhQXm9vb2lpqZeXl4mJ/p72z7kxOOt7Trm3t1cmk7m4uFhbW+u5tpFzz1lrNJr8/HxHR0dDuqjIK2MiRgEARo6BXFQBAOAKYhQAgApiFACACmIUAIAKYhQAgApiFACACmIUAIAKYhQAgApiFACACmIUAIAKYhQAgApiFACACmIUAIAKYhQAgApiFACACmIUAIDK/wH1s0TfawrP9AAAAABJRU5ErkJggg==",
      "text/plain": [
       "Winston.FramedPlot(...)"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(Winston.plot(t, f,\"g^\"))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 線形内挿"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "itp = Interpolations.interpolate((t,),f, Interpolations.Gridded(Interpolations.Linear()));"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "tnew=linspace(0,1,20);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAAEsCAIAAADfNCTgAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deUATd/4+8Hc4JCQeBImCy63BIoJHxdrSalXUldautbJVkZ9u17WH4ln91UX5ovVbbetRV7t2PZAWd7WCtrUKVjywigf1KHIo4RAEFEEQOcIVMt8/0mVZL5jMkJmE5/VXEuYzn2dKeZxkMjMShmEIAAAMZSF0AAAA04YaBQDgBDUKAMAJahQAgBPUKAAAJ6hRAABOUKMAAJygRgEAOEGNAgBwghoFAOAENQoAwAlqFACAE9QoAAAnVkIH+C+3b9++fPmy0CkAoLPw9vb29vbmuBJx7Y36+/tv2bJly5Yte/fubeeQDRs2nDp1itUsRUVFBw8eZJtty5YtbIds3bq1ubmZ1ZAdO3bExcWxGlJZWRkdHc1qCBm0OXv27KmqqmI15MCBA7t27WI1RKvVfvnll6yGkEGbExsbe+fOHVZDEhMTN27cyHYiA7IZMCQhIUGtVrMacvXq1cjISLYTGZBt27ZtbP8Qfv7556tXr7IaUlRUtGzZsnYufPjwYX3VvPzyy6xmeTJGTNzd3dkOmTp16u7du1kNyczMXLt2LduJQkJC2A6ZNWtWU1MTqyHvvffexx9/zGpISUnJkiVLWA1hDNqcRYsWlZaWshoSGRk5b948VkMaGhr+9Kc/sRrCGLQ5a9asuXnzJqshO3bs+OMf/8h2IgOyGTBk69atFy5cYDXk6NGjY8eOZTuRAdlmz57d2NjIasjevXvj4+NZDcnIyBg6dCirIYxBnfM4ce2NAgCYHEsD9uo7TlJS0rRp01gNycnJ8fPzc3d3b/+Qbt26DRgwoGvXrqwmGjJkiEKhYDVk8ODB9vb2rIYUFRW5u7uz+rBGKpX6+vp2796d1UQGbI6Pj49SqbSwYPFPb1lZWY8ePZ5//vn2D7GwsBg8eLCdnR2rbAZsjpeXl5OTk5UVi8MDDx8+ZBiG7dtAA7IZMMTDw8PZ2blLly7tH6LRaMrLy8eNG9fR2fR/CBKJpP1DnJ2d3d3dbWxs2j+kqalJrVZPmjSJVTYDOudx4jrExLbaiKhbt262trashkilUkdHR7YTeXh4sB3Cqtz1unbtKpPJWA2xsrJydnZmO5EBm+Pi4sJ2iFwub2pqYjVEIpEY8N/NgM1xcnJiO8TW1taA/0UNyGbAEKVSyXaIjY1Nt27d2I4yzh8C2/0PIrKysjLgt2PAkMfhTT0AmIk6bZ0g8xp1b3TJkiV2dnYRERE8rnP06NG9evXicYXCeuGFF6RSqdApeDNw4MDGxkahU/DG09OT7YcnYubo6BgUFCR0Ct5cqbjyi90vgkxtpL3R5OTkiIiIrVu31tbW8rvmgQMHmlONqlQqA947i5arq2u/fv2ETsGb3r17+/j4CJ2CN3Z2dkOHDhU6BT+0Ou2yU8uuWF75/ub3xp/dSDV6/Pjx7OxsAz6NAgBo01eXv8ooy1DKlB8e/7ChucHIsxupRlevXr1v374RI0YYZzoA6Dwe1D+ITIr8i9ULSdZ/Gf5z7uYLm40cAIeYAMC0RZyOeNjw8H/LBw1Y/Mne7ySfnP3kbs1dYwYQ1xee1Gq1/usUMplsyJAh+hf37NljbW0taC4AEKMDGQfyHuR9dfmr94e9r0zuTkQWEosmXdNfT/51zx/2PL78/PnzKysriejGjRsVFRVExPb7hU8krhr18vJKSUkROgUAmIA6bd3yxOUFDwtsLG0m9Z+U6mZpGfQ8SSSTiv71Teo3H/h/4N/H/5Eh27Zte+SVmTNnck8irhoFAGinz5M/L3hYQEQNzQ3jY8Y/8tOFCQuT30lmdeqUwVCjAGB6iqqKPkv+7I3+bzRoGxLzEr9585ve8t6PLFOnrZNZ8/CevU1GrdEDBw4YczoAMFfLE5frGN3fJv5Nq9P6fOmTlJ+0c9JOocLgSD0AmJgLRRf2p+9fHrDcrYdbX0XfBS8siLoWdfmOYFd8R40CgCnRMbpFxxb9rvvvlgX8dpHmVaNW9Zb3XnhsIcMwgkTCZ6MAYEq+Tv06pTilu013v+1+LS9WN1afLzx/IPPA2z5vGz8SahQATEl+Zf7c5+c+8Uc5FTlGDqOHGgUAU7L61dVCR3iUuGpUp9Ppb5pmaWkpl8uFjgMA5qmurk5/QXG299p7InHVaEVFxUcffUREffr0WblypdBxAMA8RUdHp6WlEVFJSQn3tYmrRh0cHP7+978LnQIAzNz777+vf8DLyaD4whMAACeoUQAATlCjACBGl4ovCR2hvVCjACA6P2T98Ma+Nx42PBQ6SLugRgFAXBqaGz48/mFpbenHZz4WOku7oEYBQFy2XNySU5EzwnnE1pSt6nK10HHahhoFABG5V3vvf8/+75vPvXl4+mGZtWzJT0uETtQ21CgAiMiKEyu6VtVvHLBIWa1bNXLV0eyjCTkJQodqg7i+ft/Y2JiamkpEUqm0f//+QscBAKO6evfq16lfp5/p67FuFA0eHHYlZceVHUt+WhLoGWhtwed9LW/duqU/77y+vp772sRVo3V1dcePHyciBwcH1ChAp8IwzIKEBQzDSEhCRKW1pXsubOrv0P9w1uG///L3hS8s5HGua9eu5ebmElF1dTX3tUmEutDpE82cOXPv3r1CpwAAAexL3zfj4AwiCswjt0oql9H3z/32I4VUoQ5TO8gceJ+Ul84R194oAHRab/R/o+L/V7R+JarVY7m1eC/5hhoFAFGQW8vF3JXPgCP1AACcoEYBADhBjQIAcIIaBQDgBDUKAMAJahQAgBNxfeFJo9F8//33RNS9e/cxY8YIHQcAzNPFixf1N7PTnxLKkbhqlGEYrVZLPN31FADgiZqbm/VVw8tpnOKqUblcPnXqVKFTAICZCwgI0D/Qv/3lCJ+NAgBwghoFAOAENQoAwAlvNapWq6dNm+bs7BwUFKRWP+H2KSkpKbNmzerTp8/06dOTk5P5mhcAQFj81KhWqw0JCXFzc0tOTg4ICAgMDNRoNK0X0Ol006dP9/HxOXfunJeX1/Tp05uamniZGgDE6fq960JHMBJ+ajQpKSk/P3/dunVubm7h4eHV1dXx8fGtF8jMzCwsLAwLC/P09Pzwww+Li4tzcnJ4mRoARKhMUzb669EpxSlCBzEGfmo0KyvL39/fwuK3tXl5eWVmZrZewMfHp1+/fvPmzYuPj3/33Xd9fX29vb15mRoARGjVqVUVdRXLjixgNBoy97ee/NRobm6uQqFoeeru7v7IuQESiSQ0NDQ6OnrBggX79u179dVXeZkXAEQo9V7qrqu7BvUeFL3ikkQupz//WehEHYufr99LpdLy8vKWpxqNRqVStV7gxx9/XL16dVZWlkqlys7OHjhwYHBwcMs3YFsUFxdPnjyZiORy+bBhw/QvLliwwNLSkpecAGAEi44tsre1Pz37dN0qJ6IGrU4rkvN8duzYUVtbS0Tp6en6ytKfEsoRP1vn4eERGxvb8jQvL8/Pz6/1AufPn3/ppZf03apSqUaOHHn27NnHa7RHjx76s5hsbGw8PT31L7Z8VgAA4hebGZuUn/SP1/+hkCoqly5acfIzHx9mptCp9EaMGKE/uO3p6am/J+i+fft4WC/Dh9LSUhsbm4SEBJ1Ot2vXLnt7+5qaGoZhYmNj1Wo1wzAxMTG2trY3btxgGKawsFAul2dlZT2+npCQEF7yAIAg6prq3L9wH/zVYK1Oq38l9FCodK301oNbguZ6Kl46h5+9UaVSGRMTM2XKFH2BxsXFyeVyIgoLC4uMjFSpVCEhIVevXh01alT37t1v3bo1b948Ly8vXqYGAPHYeGFjfmX+O0Pe2X11t/4VD4VHvbZ+WeKy2ODYZ481XXzep76+vj47O9vb29vK6sntrNPpUlNTnZycHB0dn7gA7lMPYLqKq4uf2/ZcTWPNE3+aNDtplNsoI0dqk+juUy+VSn19fZ+xgIWFxZAhQ3icEQDEQylT3l58+2k/lVpJjRnGmERyAA0ATF4Xyy5dLLsInUIAOAgOAMAJahQAgBPUKAAAJ6hRAABOxHWIqby8fMGCBUTk5OS0YsUKoeMAgHnasWNHeno6Ed25c4f72sRVowqFIjIykoie9s1TAADuZsyY0djYSETz5s3jvjZxtZWFhYW9vb3QKQDAzHXt2lX/gJfLHuGzUQAATlCjAACcoEYBADhBjQIAcCKuQ0wAIEZNTbRxIxHRmDE0fLjQaUQHNQoAbWlsJP33uD/7DDX6OLypBwDgBHujANAWuZzq6oiIrK2FjiJG4qrRsrKymTNnEpGzs/P69euFjgMA/yY1q4sub9y48dq1a0RUWFjIfW3iqlGlUombiABAR1u6dKn+gX6/jSN8NgoAwAlqFACAE9QoAJCmSSN0BBOGGgXo7DLLMt868JbQKUwYahSgs1t4bOGxnGNHs48KHcRUoUYBOrXvb35/Iu+EzFq25Kcljc2NQscxSahRgM6rsblxeeLyAcoBB/94UF2u3pqyVehEJgk1CtB5bbqwKbsie/OEzb/v9/vXvV5fnbS6pKZE6FCmBzUK0Endq7237ty6Kd5TxvcdT0Rf/P6LxubGiNMRQucyPeI6iwkngwIYzUcnPqpprHl/2Pt5D/KISEKSYJ/gqGtR7/u/P8RxiNDpOha/J4NKGIbhvha+zJw5EyeDAhjBL3d+GbFrhI7RPf6jkW4jz8w+Y/xIguClc8S1NwoAxlHXVLf9te1P++mD+gcKqcKYeUwaahSgMxrpNnKk20ihU5gJHGIC6EzS02nYMBo2jC5cEDqK+cDeKEBnUlNDV64QET18KHQU88Hb3qharZ42bZqzs3NQUJBarX58gQcPHsydO9fJycnV1XXTpk18zQsALCgU9Npr9Npr1KuX0FHMBz81qtVqQ0JC3NzckpOTAwICAgMDNZpHLxizePHihoaGtLS0zZs3R0RE8PI9AwBgp39/OnKEjhyhoUOFjmI++HlTn5SUlJ+fv27dOgsLi/Dw8A0bNsTHx0+dOrVlgdzc3K+//rqoqMjBweGtt96ytbVtbMTZuwBgDvip0aysLH9/fwuL3/Ztvby8MjMzWy9w8+bNvn37nj9/Pi4ubuDAgaGhoe7u7rxMDQAgLH7e1Ofm5ioU//mWmbu7e1VVVesFcnJyysrKvvzyy3HjxqWmpnp7e6elpfEyNQCAsPjZG5VKpeXl5S1PNRqNSqVqvUB9fX1dXV1sbKxSqZwzZ86LL764f/9+X1/fR9ajVqs9PDyISCaTDRny2+loe/bsscZtXQGAs/nz51dWVhLRjRs3KioqiEgmk3FfLT816uHhERsb2/I0Ly/Pz8+v9QKurq4ODg4ODg76py4uLra2to+vx8vLKyUlhZdIAACP2LZt2yOviOjOoJMnTy4sLDx27BjDMLt37y4pKdHXaFxcXHZ2NhGNHz++trZ23759RKRWqxMTE3lJDwAgOH5qVKlUxsTETJkyxcXFZenSpd9++61cLieisLCwU6dOEVHPnj2jo6Pff/99V1fXAQMGhIWF4RATAJgH3s5iCg4OnjRpUnZ2tre3t5XVb6u9e/duywJvvvnmxIkT1Wq1q6urnZ0dX/MCAAiLz5NBpVLp40eNHlngkc9MAQBMHS5NAgDACWoUAIAT1CgAACeoUQAATsR1vdHm5ubS0lIisrKysre3FzoOgIilp9Px40REc+dS165CpzExDx8+bGhoICKtVst9beKq0crKSv0NQR0dHZcvXy50HAARu3CBli4lInr7bdQoW3FxcRkZGURUVlbGfW3iqtGePXviis4A0NH+/Oc/6x+I6GRQADC22bPp4UN6+JD69BE6Smcnrr1RAGgva2vClc/EAXujAACcoEYBROq+5v71e9eFTgFtQ40CiFTE6YiwhDChU0DbUKMAYpRWmrbjyo6fC34+kHFA6CzQBtQogBgtOrZIYavw6+23PHF5nbZO6DjwLKhRANE5dOPQqVunPh798fbXtt9+eHvD+Q1CJ4JnEdcXnhoaGi5dukREMpns2ZcuBTBXjc2NH534yEfpM2foHCsLq2Cf4HVn180aNMu1h6vQ0cyHWq1+8OABEWk0Gu5rE9feaENDQ0pKSkpKCm6/DJ3WhvMbsiuyN/9+s5WFFRFtGL9BIpH89eRfhc5lVtRqtb5q6up4+MBEXHuj3bt3DwvDoUnovO7W3F1/br2DzOFk3smTeSf1L7r2cP1X2r/mDZ/3ovOLwsYzG6+//rr+gf7tL0fiqlGATm7FiRXVjdXVjdWfJn/6yI8WJiy8OOeihURc7yCBUKMAohI9OTp6crTQKYAd1CiA0AoLqbSUrKxo0CCho4Ah8AYBQGiff07DhtHo0ULnAAOhRgEAOMGbegChvfsujRtHXboInQMMhBoFEJqPD/n4CB0CDIc39QAAnIhrb7S2tnb//v1E1KNHj4kTJwodBwDM09mzZ4uLi4mosrKS+9rEVaMSiaRr165EJJPJhM4CAGZLKpXqq8bCgod35OKqUZlM1nKSFgBAB/H399c/0L/95QifjQIAcIIaBQDgBDUKAMAJbzWqVqunTZvm7OwcFBSkVquftlhVVZVKpYqOjuZrXgAAYfFTo1qtNiQkxM3NLTk5OSAgIDAw8GnXlP7ggw8KCgoYhuFlXgAAwfFTo0lJSfn5+evWrXNzcwsPD6+uro6Pj398sb1791ZUVAwfPpyXSQEAxICfGs3KyvL392/5BpaXl1dmZuYjy+Tl5UVGRkZFRUkkEl4mBRCdhASKiqKEBKFzgFHx873R3NxchULR8tTd3b2qqqr1AlqtdsaMGZs2bXJ0dORlRgAx+uwzSkqikSMJ5+B1JvzUqFQqLS8vb3mq0WhUKlXrBT755BOFQuHp6Zmenq7RaIqLiwsKCtzc3B5ZT2Fh4YQJE4hILpe3fD922bJlVlbiOk0AAEzR1q1ba2pqiCgtLU1fWfr7g3LETz15eHjExsa2PM3Ly/Pz82u9QHV1dUlJSWhoKBHl5OTs3LmzrKxsy5Ytj6ynZ8+es2bNIiIbGxsXFxf9i5aWlryEBOhw335L9fUklQqdA55swoQJTU1NROTv719bW0tEu3fv5mG9DB9KS0ttbGwSEhJ0Ot2uXbvs7e1ramoYhomNjVWr1Y8sPHr06KioqCeuJyQkhJc8AADtwUvn8LM3qlQqY2JipkyZoi/QuLg4uVxORGFhYZGRkY+8wQcAMCe8feYYHBw8adKk7Oxsb2/vlo8y7969+/iSp06d4mtSAADB8XnoRiqV+vr68rhCAADxwzn1AM/C4Iw7aAtqFOBZNl/cfKn4ktApQNRQowBPVVpbuubMmkXHFmGfFJ4BNQrwVCtPraxqqLpYdHFv2l6hs4B4oUYBnuxaybWoa1FLX1oa4BKwPHF5VUNV22OgUxLXSZb379+fO3cuEfXp0ycyMlLoONCpLTq2qKes58qRK3MqcobvHP5Z8mdrx6wVOhTwY9u2bdevXyci/f1BORJXjfbs2fOLL74gnm7XB2Cwg5e/uab++W9/+KqHTY/nnZ6f6Tdzw/kNswfP7mffT+howIO//OUvzc3NRDRnzhzuaxNXW0kkEplMJpPJpDgrGYRTp60b8+qfqtbRrJjr+lc+G/eZjZXNipMrhA0GfLGxsdFXDS97bOKqUQAx+Dz5cx2jI6LdV3fbf2pv/6m99zZvTZMmLjMuKT9J6HQgOuJ6Uw8gOK1OW1lfeXzWy9YN2uq+DsE+fVr/NK007VX3VwWKBiKFGgX4L1YWVpsmbKIJQucA04E39QAAnKBGAQA4QY0CAHCCGgUA4AQ1CgDAibiO1JeVlc2cOZOInJ2d169fL3QcADBPGzduvHbtGhEVFhZyX5u4alSpVO7di0vpAEDHWrp0qf6Bfr+NI3HVKADPmpvp9m0iIgcH6tZN6DRgnvDZKJi10lLy9CRPT/rXv4SOAmYLNQoAwAne1INZs7Oj6GgiohdfFDgJmC/UKJg1W1uaNUvoEGDm8KYeAIAT1CgAACeoUQAATlCjAACciOsQE04GBQAjwMmgAACc8HsyKN7UAwBwIq69UYCnqqqirCwiIm9v6tpV6DQA/4G9URCjgzcONjPN//XSxYs0fDgNH07XrwsUCuDJUKMgOneq78z+fvbOKzuFDgLQLrzVqFqtnjZtmrOzc1BQkFqtfnyBX3/9ddSoUXZ2diqVavv27XzNC+ZnxckVNY01EacjKusr//Pq0KF05AgdOUIDBggXDeAJ+KlRrVYbEhLi5uaWnJwcEBAQGBio0WhaL9DY2Dh27NiJEyemp6eHh4fPnz8/JyeHl6nBzFwqvhSTGjN1wNSKuorVZ1b/5wcODvTaa/Taa2RnJ1w6gCfgp0aTkpLy8/PXrVvn5uYWHh5eXV0dHx/feoFr165ZW1t/9NFHzs7Os2fPfvHFF7/99ltepgZzwjDMwoSFfbr12fOHPXOfn7stZVtGWYbQoQDawE+NZmVl+fv7W1j8tjYvL6/MzMzWC/Tr1+/kyZP6x7W1tRkZGSNGjOBlajAnMddjLhVfWh+4vmuXrmvHrO1h02PRsUVChwJoAz81mpubq1AoWp66u7tXVVW1XqBnz54+Pj5ElJGR8corrwwbNmzMmDG8TA1mo6axZsXJFSOcR4T4hhCRva39qlGrTuSdOJp9VOhoAM/Cz/dGpVJpeXl5y1ONRqNSqR5ZhmGYlStXbt++fdWqVQsXLpRIJI+vR61We3h4EJFMJhsyZIj+xT179lhbW/OSE8Rs/bn1d6rvqOxV7x19T/9KU3OThcRiyU9LxnmO62LZRdh4YAbmz59fWVlJRDdu3KioqCAimUzGfbX81KiHh0dsbGzL07y8PD8/v9YL6HS6d9999/bt25mZmY6Ojk9bj5eXV0pKCi+RwLSUacq+u/mdp8KzsKqwsOo/pzm727lrddrYzFj9LioAF9u2bXvkFRHdGXTy5MlhYWHHjh2bMGFCVFRUSUmJvkbj4uIGDRqkUqnOnTt3+PDhX375RSqV6v81sLGxsbW15WV2MDH5+XTjBhHR+PFkaal/TSlTZnyAo0lgkvj5bFSpVMbExEyZMsXFxWXp0qXffvutXC4norCwsFOnThHR2bNnS0tL3dzcFP+2cuVKXqYG03PwIAUFUVAQ/fe34gBMFG/n1AcHB0+aNCk7O9vb29vK6rfV3r17V/8gPDw8PDycr7kAAMSDz0uTSKVSX19fHlcI5ikkhF5+mYhILhc6CgAPcIUnMDpHR3r6YUYAk4NLkwAAcIIaBQDgBDUKAMCJuD4b1Wq1BQUFRNSlSxcnJyeh4wCAeSorK9Nfha6pqYn72sRVo9XV1Tt27CCi3r17L1iwQOg40Jb9+yklheztCd8CBpNy7NixmzdvElHrs9gNJmEYhvta+DJz5kzcGdSUzJ5NX39Nbm6Uny90FABD8NI5+GwUOLC1pW7dqFs3oXMACAk1Chxs305VVZSWJnQOACGhRgEAOEGNAgBwghrt7HBteQCOUKOdWn5lfvCB4B+yfhA6CIAJQ412assSl9Vp65b+tLShuUHoLACmCjXaeZ0pOBOXGfe61+u5D3LzJgwnPz/CKQ8A7InrLKb6+vqkpCQiksvl/v7+QscxZzpGt/SnpR52HrHBsdMPTm/c8SPdaSY3N6FzARhDWlqa/vyl2tpa7msTV402NTXl5OQQkUKhQI12qB1Xdly5eyXuj3FSK+nG8RsPex5pdnUb+tJLQucCMIa7d+/evn2biBoaePg4S1w12q1btzlz5gidwvxV1ldGnI4Y7T76Le+3iMhT4Xkv8kP/5M8u/L+xw4XOBmAE48eP1z/Qv/3lSFw1Csax5syaMk1ZD2mPj058pH+lurFax+gW/7T43J/OSSQSYeMBmBbUqNmpqiKGIamUbGye+PNblbf+mfZPhVRxJv/MmfwzLa8rpIobZTcOqw//of8fjJUVwBygRs2OgwM1NdGKFfTJJ0/8uYedx70P7xk5FIAZwxeeAAA4wd6o2dmwgZqbaTiOFQEYCWrU7OAr9ADGhTf1AACcoEbFqqiIxo2jcePo9GmhowDAs4jrTX1NTc2ePXuISKFQTJ48Weg4/NA0aWytbFl/GbO2lk6cICJ6552OSAXQmZ04caKwsJCIHjx4wH1t4qpRS0tLV1dXIuratavQWXjz8c8fD1AOCPULZTfM1pb0p2YqlR2RCqAzUyqV+j0ba2tr7msTV43a2tqOHTtW6BR8ynuQt/nCZgeZw5vPvdm1C5t/G1xdKTm5w3IBdGqDBg3SP9C//eUIn412rA+Pf2hpYXmn+s76c+uFzgIAHQI12oFO3Tr13c3vNqnm7yl54eGWT/PyrgidCAD4hxrtKM1M86JjizwVnu8wg2d9dXHrYe2271YIHQoA+GekGlWr1dOmTXN2dg4KClKr1caZVFhfXf4qrTRt04RN1ha/fYadmJeYmJcobCoA4J0xalSr1YaEhLi5uSUnJwcEBAQGBmo0GiPMK6CKuor/Of0/Q52G+vbyvfXywIKbF9XpZyrcei07vqyZaRY6HQDwyRhH6pOSkvLz89etW2dhYREeHr5hw4b4+PipU6caYWqhrD6zuryuvLyuvO/f+rZ+/U5d6c4rO98b9p5QwQCAd8ao0aysLH9/fwuL3/Z8vby8MjMzjTAv/zQa0mhIIqGePZ+xlI7RDeo96B+v/+OJP+0pe9ZYADA5xqjR3NxchULR8tTd3b2qqsoI8/Lv009pzRoiIoZ5xlIWEot3huDUI4DOwhg1KpVK9Tfh09NoNCqV6olLFhQUjBw5kv77zqARERFWVuI6TQAATNHnn39eXV1Nre4MystxGmPUk4eHR2xsbMvTvLw8Pz+/Jy7Zu3fvsLAwIrK2tnZ0dNS/aGlpaYSQ7fLaa6RUEm5VBGCa3n77ba1WS0T379+vq6sjos2bN3NfrTFqdPLkyWFhYceOHZswYUJUVFRJScnTalQqlY4aNcoIkQw0fDguhwxguvSX7CAiT09P/dmubzQAAAigSURBVIOdO3dyX60xalSpVMbExEyZMsXe3r6mpiYuLk4ulxthXgAAIzDS1++Dg4MrKioSEhLu378fGBhonEmfZdUqGjUK16ADAO6Md+hGKpX6+voabbo2ZGTQzz8TH5caBIBOrrMeAff1pQcPqG/ftpcEAHimzlqjq1cLnQAAzASu8AQAwInp743GxFBlJfXvT+PHCx0FADoj06/RNWsoJ4dmzECNAoAgxFWjZWVloaGhROTs7Lxu3Tqh4wCAedq0adO1a9eISH9/UI7EVaNKpfKbb74hIha3I/7lF2puJhubDowFAOZlyZIlDMMQkX6/jSNx1SixKlA9O7uOCQIA5ox11TwdjtQDAHCCGgUA4MTkazQ+Pj47O1voFLw5c+aM/pNv83D16tWzZ88KnYI3WVlZCQkJQqfgTXFxcesrWJq6hw8f7tmzR5CpRffZKFu5ubk9n3lLD9NSWFjo4OAgdArelJSUVFZWCp2CNxUVFXl5eUKn4E11dXVWVpbQKXjT0NCQlpYmyNQmvzcKACAscdVoUVER2yG3b9++f/8+qyHl5eXJyclsJ/rhhx/YDvnxxx+ZZ9616XFFRUUlJSWshtTW1p44cYLVEDJocxITE9neceHu3btsf6c6ne7IkSOshpBBm3Pu3LmKigpWQ8rKym7fvs12IgOyGTDk8uXLxcXFrIZUVlYWFBSwnciwPwSdTsdqSHp6em5uLqshdXV1Bvx2DOicx4mrRg34pebn59+7d4/VkNLS0qSkJLYTGfAp0sGDB5ub2d2VvqCggO3vtaamxoAP7AzYnPj4+NraWlZDioqK2P6frdVqDx06xGoIGbQ5p0+fLisrYzXk3r17+fn5bCcyIJsBQy5evMj2m+QVFRW3bt1iO5EB2Q4dOsT2DyE1NVWtVrMaUltba8DmGNA5jxNXjQIAmBwJ23edHcrX17e+vp6I5HL58Pbd9ejy5cu9e/d2cXFp/yy1tbUlJSV9WV5s9Pr160+7hdQzhvj6+rL6lu/169dlMlm/fv3aP6SxsTEvL++5555jm43t5ty8edPT07NLly7tH6JWqxsbGwcOHNj+ITqdLiMjg+0Vvg3YnNzcXEdHR1b3syksLCwtLX3++ec7OpsBQwoKCnr06GHH5myUsrKy7Ozsl156qaOzpaWlDRw4kNUfwp07d6ytrZVKZfuH1NbWpqSkjB49uj0LX79+veV2xdy/6iOuGm1ubm5oaBA6BQB0FtbW1tbW1hxXIq4aBQAwOfhsFACAE9QoAAAnqFEAAE5QowAAnKBGAQA4MeEaVavV06ZNc3Z2DgoKYnvCgxi0mf/XX38dNWqUnZ2dSqXavn278ROy0s5fR1VVlUqlio6ONmI0Q7S5OQ8ePJg7d66Tk5Orq+umTZuMn5CVNjcnJSVl1qxZffr0mT59ugGnSovBkiVL1qxZI8jUplqjWq02JCTEzc0tOTk5ICAgMDCQ7enewmozf2Nj49ixYydOnJienh4eHj5//vycnByh0rap/b+ODz74oKCgQORfs2vP5ixevFh/SaHNmzdHRETwckufDtLm5uh0uunTp/v4+Jw7d87Ly2v69OlNTU1CpTVAcnJyRETE1q1b2Z6szBvGNCUmJjo4ODQ3N+uf2tnZxcbGChuJlTbzX7x4sXfv3i1PAwIC1q5da9SIbLTz1xETEzNx4sSAgICoqCjjBmSnzc3R/5NWVFSkf3r06NGcnBxjp2y3NjcnLS3N2tpao9EwDFNVVWVhYZGZmSlAUENFRERMmzbNxcVl+fLlggQw1b3RrKwsf39/C4vf8nt5eWVmZgobiZU28/fr1+/kyZP6x7W1tRkZGSNGjDB2ynZrz68jLy8vMjIyKiqKx3vgdJA2N+fmzZt9+/Y9f/7822+//fHHHw8YMIDtucXG1Obm+Pj49OvXb968efHx8e+++66vr6+3t7cQSQ20evXqffv2CfgHYqo1mpubq1AoWp66u7tXVVUJmIetNvP37NnTx8eHiDIyMl555ZVhw4aNGTPG2Cnbrc3N0Wq1M2bM2LRpk6Ojo9HTsdbm5uTk5JSVlX355Zfjxo1LTU319vYW6oLB7dHm5kgkktDQ0Ojo6AULFuzbt+/VV181dkQTZ6o1KpVKW64sQEQajUalUgmYh6325GcYJjw8/JVXXgkNDf3pp5/EvBPX5uZ88sknCoXC09MzPT1do9EUFxfzcoGyDtLm5tTX19fV1cXGxs6ZMycuLm7w4MH79+83esz2anNzfvzxx9WrV2dlZeXk5KjV6u3bt5voUSahmGqNenh4tL6qa15eHturzgirzfw6nW7u3LmXL1/OzMxcvHhxyzsycWpzc6qrq0tKSkJDQ0NDQ9Vq9c6dO8V8dLvNzXF1dXVwcGi53YuLi4utra1RI7LR5uacP3/+pZde0nerSqUaOXKkOd1ByxgE+USWu9LSUhsbm4SEBJ1Ot2vXLnt7+5qaGqFDsfC0/LGxsWq1mmGYM2fO9OrVq6Cg4MG/6Y8AiFObm9Pa6NGjRX6Iqc3NuX//fvfu3f/5z38yDJOVlWVnZ3fr1i1hMz9Dm5sTExNja2t748YNhmEKCwvlcnlWVpbAodkLDg4W6hCTqdYowzAHDhywtbX93e9+16NHj8TERKHjsPbE/I6Ojl999RXDMGvXrn3kH7wlS5YImrcNz96c1sRfo0w7NufQoUPdu3d3cXGxtLRctWqVcEnb5dmbo9PpFi9e3KtXr379+llaWi5YsEDQsAYSsEZN+0J59fX12dnZ3t7eVlYmeYtTU8//iM62OfX19Wq12tXVldXFkoXS5ubodLrU1FQnJyeTOAwoKqZdowAAghP1gQsAAPFDjQIAcIIaBQDgBDUKAMAJahQAgBPUKAAAJ6hRAABOUKMAAJygRgEAOEGNAgBwghoFAOAENQoAwAlqFACAE9QoAAAnqFEAAE5QowAAnPwfJQtZKqpergsAAAAASUVORK5CYII=",
      "text/plain": [
       "Winston.FramedPlot(...)"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(Winston.plot(t, f,\"g^\",tnew,itp[tnew],\"r.\"))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 複素列の内挿"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "10-element Array{Complex{Float64},1}:\n",
       " 0.987789+0.155797im\n",
       " 0.952633+0.304122im\n",
       " 0.898704+0.438555im\n",
       " 0.832171+0.554519im\n",
       " 0.760245+0.649637im\n",
       "  0.69021+0.723609im\n",
       "  0.62863+0.777705im\n",
       " 0.580823+0.81403im \n",
       " 0.550621+0.834755im\n",
       " 0.540302+0.841471im"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "g=exp(t*im)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "itc = Interpolations.interpolate((t,),g, Interpolations.Gridded(Interpolations.Linear()));"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### 実部、虚部、絶対値、偏角の順に表示"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAAEsCAIAAADfNCTgAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deUBU9d4G8O8MIMMqICOgrMpMF1nSFE3xlprm1fLeMjUIyez22ntD3NPKJDSvaa9LppYZLr1QmkPdupYoaC6JmuFCIMKwiSAoIyDbsDjL+8d0ecmNOXMOc87A8/lrZji/c54T8XhmzpnzE+n1egIAAFOJ+Q4AAGDZUKMAAKygRgEAWEGNAgCwghoFAGAFNQoAwApqFACAFdQoAAArqFEAAFZQowAArKBGAQBYQY0CALCCGgUAYMWa7wB/cO3atczMTL5TAEBPERQUFBQUxHIlwjoaHT16NNMh69ev/+mnnxgNKS8v/+abb5huaPPmzUyHbNmyRavVMhqyY8eOlJQURkNu3769Z88eRkPIpN3ZvXt3fX09oyH79+9PTExkNESj0Wzbto3REDJpdxQKRUVFBaMh6enpGzZsYLohE7KZMCQ1NVWpVDIacuHChYSEBKYbMiHb1q1bmf4hnDx58sKFC4yGlJeXv/nmm4yGENGECROYDrkPvZD4+/szHTJt2rSdO3cyGpKbm7t69WqmG4qOjmY6ZNasWXfu3GE05L//+7/ff/99RkNu3LixaNEiRkP0Ju3OggULqqqqGA1JSEiIjY1lNKS1tXX27NmMhuhN2p1Vq1bl5eUxGrJjx44ZM2Yw3ZAJ2UwYsmXLljNnzjAa8uOPPz711FNMN2RCtldeeaWtrY3RkOTk5IMHDzIacvny5ccee4zREL1JnXMvYR2NAgBYHCsTjuq7zvfff5+Tk5OWlvbbb7+NGjXKmCGFhYVhYWH+/v7Gb8XJyWnQoEGOjo6Msg0ZMsTV1ZXRkMGDB7u5uTEaUl5e7u/vz+jDGolEEhoa6uzszGhDJuxOcHCwVCoVixn806tSqXr37j106FDjh4jF4sGDB7u4uDDKZsLuyOVyLy8va2sGpwfq6ur0ej3Tj55MyGbCkICAAG9v7169ehk/RK1WV1dXM31Xa/IfgkgkMn6It7e3v7+/ra2t8UPu3LmjVCqnTJlizMK7du368ssv09LSbt++/fe//934rdyXsE4xeXp6LlmyhIhsbGyMHOLk5GRnZ8doKxKJxNPTk2m2gIAApkMYlbuBo6Ojvb09oyHW1tbe3t5MN2TC7vj4+DAd4uDgcOfOHUZDRCKRCf/dTNgdLy8vpkPs7OyY/utLJmUzYYhUKmU6xNbW1snJieko8/whMD3+ICJra2vjfztTp06dPHkyES1YsIDphu6zafar4JCVlZUJBQcAwEj72x1Gb0ceRFg1aoKxY8f27duX7xScGTFihEQi4TsFZ0JCQtra2vhOwZkBAwYw/fBEyDw9PQ1HZN2Ds7Pz9OnTedm0xddoSEgI3xG4JJPJ+I7AJV9fX74jcMnDw8PDw4PvFJxxcXF57LHH+E7BGYlEMnLkSF42jTP1AACsWPzRKB04QPX1FBBAxp3ZBwDgluXX6KJFVFhIL72EGgUAXuBNPQAAK8I6GtVoNNeuXSOiXr16GXvl05Ej1NZG3ej8KQB0NZVK1dzcTERMr2u+L2HVaENDw/bt24nIw8Nj/vz5Ro3x8+vaTADQ7aSmpubl5RFRdXU1+7UJq0ZdXV3XrFnDdwoA6OZefvllw4OZM2eyX5uwatR8Pv6Yfv2VfH3pn//kOwoAWLaeeorp+HFKTqYDB/jOAQAWz6w1umjRolWrVplziw/k6UkDBhDze20AANzFTG/qMzIyDh8+vGXLlkWLFplni5345BO+EwBAN2Gmo9G0tLSCggITbk0GACBwZjoaXblyJRHNmDHDPJsDADCbnnqKyTQffkju7sT8/rgA0I0J64InpVJpuLe2vb39kCFDDC/u3r3b+Jvhd62mJmJ+sW5Nc42zrbO1WFj/qQF6oLlz596+fZuIrly5UlNTQ0RMJ5u4L2H9bcvl8nPnzvGd4sFGjKDYWGIypQwRvXP0nUHSQfNGzOuiUABgpK1bt971Ci6/N7vJk4nh3cKzbmYlXkjsLekdHRbdx65PF+UCAB7hs9GutfDQwt6S3vWt9fHH4vnOAgBdwqxHo/v37zfn5niXkpty7Oqx3U9srM45t//f27NDXgn1Dec7FABwDEejXaVV2/r20bcf9Xj05Vv9F8/b98tnus3JcXyHAgDuoUa7yvrT6wtrCj/6y0di0e//kX+5/st3ed/xmwoAOIdTTF2ioqFi7am1I/qPcLJ1+k3WZv/FJq1e23p9/ZK0JZNkk2ytbPkOCACcQY12ibePvt3Y1vjL9V+G7Rj2hx/U0kdnP1oWsYynXADAPWHVaHNz89GjR4nI0dFxxIgRfMcxkV6vf/nRl2PCYu77U4deDmbOAwB3ycrKunXrFhE1NjayX5uwalSr1RrmYnJ1deU7i+lEItFTAU+xXUt19e+3oZo6lYKD2acCgHYqlaqsrIy65VxMjo6Os2fP5juFMNy6RfHxRESBgahRAG6NHz/e8MDw9pclnKkHAGBFWEej8P8eeYQ0GiIiMf6pAxA01KiAWVnxnQAAOocjHQAAVlCj3c6TT9KoUfTZZ3znAOgpUKPdzpkzdOYMlZYaufgXWV/o9fouTQTQveGz0W5nwgS6c4dkMmOWzVXl/te//0ssEj/oywIA0CmRoI5EJk6cKJVKicjb23vt2rV8x+n+JiZPTCtK6+fUL39uvmMvR77jAJjJhg0bLl68SERlZWUnTpxguTZhHY1KpdLk5GS+U/QUB5QH0orSXhn8yheXvvjg1Af/HPdPvhMBmMnixYsNDziZRASfjfZQbdq2JWlLgtyDdkzZEfNozIbTGwprCvkOBWCRUKM91EdnP1JWKzf9ZZON2ObDCR++e5JKpz9N90z4BQCdQo32RFVNVWt+XvO3R/42ceBEIvJw8HittM9TP5Wo/oVPVAAYE9Zno2Ae7xx9p6617pb61gzFDMMrb9lpdE50Xl0wWa+1EuHbUwAMoEZ7nKLaohOlJwa4DqhsrKxsrDS8OP1VRyJHItLk//v5Pz3Pa0AAC4Ma7XEGug4siCvgOwVA94HPRgEAWBHW0Wh9ff327duJyM3NbcaMGXzHAYDuKTU1tbS0lIiqq6vZr01YNdqrV6+QkBAicnDAhEUA0FX8/PycnJyI6PDhw+zXJqwalUgko0eP5jsFAHRzgwYNMjwwvP1lSVg1Chbm1CkqLCRHR5o2je8oALxBjQILiYn0xRfk54cahZ4MZ+oBAFhBjQIL69dTSQmdOsV3DgA+4U09sODuTu7ufIcA4BmORgEAWEGNAgCwIqw39SqVynAzakwiAgBdp+MkIuzXJqwaxSQiAGAG3E4iIqwahR7h9m1SqYiIAgNJJOI7DQBb+GwUzG7nTpLLSS6nxka+owBwgLMaVSqVkZGR3t7ekydPViqV9y5w+PDhxx57zMPDY+rUqYWFmD0NHqaqqYrvCADG4qZGNRpNdHS0n59fRkZGRETE+PHj1Wp1xwWam5unTZsWFRV18uTJgICA559/vqWlhZNNg+V55hlKTqbkZLKzu+/PtXrtxOSJv1b8auZcACbScyE9Pd3d3V2r1Rqeuri4KBSKjgvs378/LCzM8FitVjs6Op45c+be9URHR3OSByzaJ79+Qgk0etdonU7Hdxbo5jjpHG6ORvPz88PDw8Xi39cml8tzc3M7LlBXV+fi4tL+VCQS5eTkcLJp6GZqW2rjj8V7OHicunZqf+5+vuMAdI6bGi0qKnJ1dW1/6u/vX19f33GBadOmnT9//quvvqqsrFy9enVjY2NtbS0nm4Zu5r1j79U21x6OORzeL3xJ2pKmO018JwLoBDcXPEkkko734ler1TKZrOMCLi4un3322cqVK+fNm/fss89GREQEBQXdu56SkpJhw4YRkYODg+EBEa1bt87aGhdm9QhXbl3Znrn9H+H/eNTj0c2TNkfsjFh/ev17T77Hdy7oJlauXGk4wsvJyTFUllarZb9abuopICBAoVC0Py0uLg4LC+u4QG1t7ejRo/Pz84lIq9V6eHgYJgu5S//+/ePj44nI2tq6/fAWHdpzLDq8yLGXo6E3R3qPjAyJXHdq3SuDX/Hr7cd3NOgOYmNjDb1ZV1fX1tZGRAkJCexXy01DPffcc3FxcYcOHZo4ceKuXbtu3LhhqNGUlJRHH31UJpNdv359+PDhSqWyf//+mzZtCg4O9vO7zx9G+1xM0AMdUB44VHgoJizm0o1LhlcmDJywL2ff20fe/uqFr/jNBt2D+39uSObh4WF4IJFI2K+WmxqVSqVJSUlTp051c3NrbGxMSUkxzEkXFxeXkJAgk8lCQkJmz5799NNPG869Hj58WISvr0AHbdq2xYcXE1HSb0lJvyV1/NG+nH2xw2MjfCJ4igbQCc7eL0+fPn3KlCkFBQVBQUHtb8MrKyvbF9i2bZtarS4vL5fL5VxtFLqNpjtNnzzzyYN+6mbnRpmZZPge9CefUHCw+ZIBdIbLjx0lEkloaOhDFrC3t0eHwn25SlzHDxj/sCUupNHJk0REdXXmiQRgJJy9AQvh5kZjxhAR9e7NcxKAP0KNgoUYNoyOHeM7BMB94A5PAACsCOtoVKPRXL16lYh69erVr18/vuMAQPdUVVVluH2S4epRloRVow0NDYmJiUTk4eERFxfHdxwA6J7S0tLy8vKIqKamhv3ahFWjrq6uq1ev5jsFAHRz7XOHcDKJCD4bBQBgBTUKAMAKahQAgBXUKAAAK8I6xQTAMbWadu8mIhozBt/Ehy6CGoVura6O5s4lItq+HTUKXQRv6gEAWBHW0ahKpTJcxuXt7b127Vq+44Dl8/CgqioiIicnvqOAgGzYsOHixYtEVFZWxn5twqpRqVSanJzMdwroRsRikkr5DgGCs9hw71pcfg8AIASoUQAAVlCjAACsoEYBAFhBjQIAsIIaBQBgBTUK8DBbz23NrMjkOwUImrCuG21ubk5PTyciR0fHkSNH8h0HerCaGtLpaqj5vePvBUuDT7xyQiQS8Z0JOHPp0iWVSkVEjY2N7NcmrKNRrVZbUVFRUVFRZfjmCQBf5HKSSnNmP1PTXPPztZ/35+7nOxBwqaamxlA1d+7cYb82YR2NOjo6zpo1i+8UAL+7rLocNzzubPnZJWlLpsin2NvY850IuDFu3DjDA8PbX5aEdTQKIBTr1m17+U8Hhzi+N+a9j/7y0fX66+tPr+c7EwgUahTgPv4V4TZ3QN6kVz/oY9dnlM+oGcEz1p5aW1pXyncuECLUKMDd2rRty9KXDZIOmjN0juGV9U+vF4vE7xx9h99gIEzC+mwUQAg2nd1UUFPgKnF9ZOsj7S9q9dq92Xtjw2NH+YziMRsIEGoU4A90ep2qSdV+HHqXrJtZqFG4C2oU4A/EIvH6p3E2CRjAZ6MAAKwI62gUk4gAgBlgEhEAAFYwiQgAgICgRgEAWOGsRpVKZWRkpLe39+TJk5VK5b0LnDt3btasWf369YuKisrIyOBquwAA/OKmRjUaTXR0tJ+fX0ZGRkRExPjx49VqdccFdDpdVFRUcHDwqVOn5HJ5VFQUJzdWARCWbdtowQLaupXvHGBeei6kp6e7u7trtVrDUxcXF4VC0XGB7OxsGxsbtVqt1+vr6+vFYnFubu6964mOjuYkDwA/xozRE+mfeILvHGAsTjqHmzP1+fn54eHhYvHvx7ZyuTw3N7fjAsHBwYGBgbGxsdOmTUtOTg4NDQ0KCuJk0wAC4uhIvXuTkxPfOcCsuKnRoqIiV1fX9qf+/v719fUdFxCJRDExMcuXLz958mRRUdH8+fM52S6AsBw4wHcC4AE3NSqRSKqrq9ufqtVqmUzWcYEDBw6sXLkyPz9fJpMVFBSEhIRMnz49IiLirvVUVFRERkYSkYODw9ChQw0vvv7661ZWVpzkBICebM+ePYbTNpcvXzZU1rVr19ivlpsaDQgIUCgU7U+Li4vDwsI6LnD69OlRo0YZulUmkz3xxBM///zzvTXq6Og4duxYIpJIJO1FjDlwAIATISEhbW1tROTq6mp4x/yvf/2Lg/Wy/3hVr9dXVVXZ2tqmpqbqdLrExEQ3N7fGxka9Xq9QKJRKpV6vT0pKsrOzu3Llil6vLysrc3BwyM/Pv3c9OMUEAOYkoFNMUqk0KSlp6tSphgJNSUlxcHAgori4uISEBJlMFh0dfeHChSeffNLZ2bmkpCQ2NlYul3OyaQAAfon0ej1X62ppaSkoKAgKCrK2vn8763S6rKwsLy8vT0/P+y4wc+ZMfKceAMyGk87h8tYkEokkNDT0IQuIxeIhQ4ZwuEUAAN7hO/UAAKygRgEAWEGNAgCwghoFAGBFWHe/xyQiAGAGmEQEoHvZt4+OHCEHB9q8me8oPQUmEQHoXk6fpp07KSmJ7xxgImEdjQL0RC4u1L8/dbhHWjudXicW4VhH6PAbAuDbqlVUXk7Z2Xe9fLTk6OZf8DbfAuBoFECINDrNvNR5FQ0VM8NmSu2lfMeBh8HRKIAQfZr5aa4qt761/t2f3uU7C3QCNQogODXNNSuPrxw/YPyCxxckXkg8X3me70TwMKhRAMFZcWxFXWvdpombEsYk9HXoOz91Pod3YgPOCeuz0erq6oULFxKRl5fX0qVL+Y4DwINcVe6O8ztiw2ND+oYQ0coxK1//4fWUKynTB03nO1r3kZiYePnyZSKqrKxkvzYu7zfKXlRU1Mcff0xE1tbWrve7/gOg25uYPDGzIlMZp+xj14eIdHrdiMQRNxtv5s3Ns7ex5ztdN1FfX9/a2kpEcXFx+/btY7k2YR2NWllZSaU4KQk913d536UVpRGR+4fud/1o45mN7z6B003ccHZ2Njx40D3mGRFWjQL0cDebbi6LWHbfH+lJj6vxhQk1CiAgrw99ne8IwBj+ZQMAYAU1CgDACmoUAIAVfDYKYJmysujgQSKiuDhydOQ7TY8mrBrVaDQlJSVE1KtXr/79+/MdB0DAzp2jd94hInr5ZdQoUzdv3lSr1UTU1tbGfm3CqtGGhoZdu3YRkYeHx9y5c/mOAwDdU3p6en5+PhHV1tayX5uwvsU0c+ZMTCICYBSdju7cISKyteU7igXjpHOEdTQKAMYSi1GgAoEz9QAArKBGAQBYQY0CALCCGgUAYAU1CgDACmoUAIAVYV3wpFKpZs6cSUTe3t5r167lOw4AdE8bNmy4ePEiEZWVlbFfm7BqVCqV4vJ7AOhqixcvNjwwHLexhDf1AACsoEYBAFjhrEaVSmVkZKS3t/fkyZOVSuVdP/3444+H/NH8+fO52jQAAI+4+WxUo9FER0ePGzdu3bp1ycnJ48ePz8vLs7f//8lgn3rqKW9v7/anS5YsGThwICebBgDgFzdHo8ePH7969eoHH3zg5+e3fPnyhoaGg4Ybyv5HcHDw1P+4devW0KFD582bx8mmAYCBmzdp61baupWuXeM7SvfBzdFofn5+eHi4WPx7Kcvl8tzc3PsuWVhYuGLFCsOlBgBgbiUlFBdHRBQYSL6+fKfpJrg5Gi0qKnJ1dW1/6u/vX19ff98lly9f/sYbb/Tr14+T7QIAJ3R6Hd8RLBg3R6MSiaS6urr9qVqtlslk9y5WWFj4ww8/fPrppw9aT21t7Zo1awwrDA4ONrw4YcKE9uNcAGBl2DBSqYiIevduf628vnzNz2s+eeYT3lKZy/Hjx1tbW4moqKjo9u3bRHTz5k32q+WmRgMCAhQKRfvT4uLisLCwexfbsmXLX//6Vzc3twet586dO9euXSMie3t7qVRqeFFQ9+cHsGzW1uTuftdry44s25u9d9bgWSP6j+AllNncunWrubmZiMrLy2tqaojI0Kps6blQVVVla2ubmpqq0+kSExPd3NwaGxv1er1CoVAqle2LyWSyffv2PWQ90dHRnOQBACNlXMsQJYjEK8WPJz6u0+n4jmNunHQON2+WpVJpUlLS1KlTfXx8Fi9e/PXXXzs4OBBRXFzcTz/9ZFimtLS0oKBg9OjRnGwRANjT6XULDi3o79z/k2c+OVt+NjkbX8U2BWffqZ8+ffqUKVMKCgqCgoKsrX9fbWVlZfsCfn5+erw9BxCS3Zd2/1rx694X9r4Y/GJSVtLS9KV/e+RvzrbOfOeyMFyeupFIJKGhoe0dCgBC1tDWsOKnFaN8Rr0Y/KJIJNo8aXNVU9WHGR/yncvy4Aw4QA+1+uTqm003P/rLRyKRiIiGeg2dGTZz/en1hTWFfEezMDhyBOiJCmsKN5/dTEQTkya2v9imbWvVtr515K2UGSn8RbM8qFGAnijvVt6Cxxc86Ke3W267SFzMmceioUYBeqJn5c8+K3+W7xTdhLBqFJOIAIAZYBIRAABWMIkIAICAoEYBAFhBjQIAsCKsz0YBQIg0GjKctAgPp//cwRLaoUYBoDOtrTR7NhHRhx+iRu+FN/UAAKzgaBQAOmNvTyUlREQPvud6TyasGq2vr9+6dSsR9enTJyoqiu84AEBERCIR+fvzHYJLP/74Y0lJCRF1nP3IZMKqUVtb22HDhhFRxznuAQC4FRgY2KdPHyJqv688G4Kr0ccff5zvFADQzT3yyCOGB5wcseEUEwAAK6hRAABWUKMAAKygRgEAWEGNAgCwghoFAGAFNQoAXUKv1/MdwUyEdd0oJhEB6DZe/ferWydvdbBx4DvIfWASEQAQum+vfLvn0h4fZ59VY1fxneU+MIkIAAhai6blzfQ3RSRaf3p9aV0p33G6HGoUADi28czG4triXX/bJRaJl6Yv5TtOl0ONAgCXbjbdXJex7oWgF14Z/Mqy0cv2X95/ovQE36G6FmoUALi0NH1pq6Z17fi1RPT21xU/fudw5q2ZWr2W71xdCDUKAJw5X3k++bfkJaOWBLoFEpF16uHJl5r6XynfdXEX39G6EGoUALih1+vjDsbp9LqPf/nYbZ2b2zq3Y9blyr5W153o3Z/erWut4ztgVxHWBU8AYLkKawtDPUJDPULbX9m38fcHzxGdrzg/LmAcP8m6mEhQ3zSYNGmSXC4nIi8vr7feeovvOADQPX3++ec5OTlElJ2dzf4G+MI6GnVxcYmPjycia2thBQOA7iQyMrKtrY2I5s6dy35twmorKysrwwQpAABdx8nJyfDAysqK/dpwigkAgBXUKAAAK5zVqFKpjIyM9Pb2njx5slKpvHeB2traOXPmeHl5+fr6bty48d4FAAAsETc1qtFooqOj/fz8MjIyIiIixo8fr1ar71pm4cKFra2t2dnZmzZtio+P5+T+VAAAvOPmFNPx48evXr36wQcfiMXi5cuXr1+//uDBg9OmTWtfoKio6IsvvigvL3d3d3/hhRfs7OwMp8kAACwdNzWan58fHh4uFv9+bCuXy3NzczsukJeXN3DgwNOnT6ekpISEhMTExPj7+3OyaQAAfnFTo0VFRa6uru1P/f396+vrOy5QWFioUqm2bds2c+bMQ4cOrVmz5ty5c6GhoXetp6Wl5dixY0RkY2Pj5eVleHHAgAEikYiTnADQk5WWlmo0GiJSqVTNzc1E1NjYyH613NSoRCKprq5uf6pWq2UyWccFWlpampubFQqFVCp97bXXRo4cuW/fvntrtKqqKiEhgYgcHBzCw8MNL65YsQJX4wMAewqFoqGhgYhycnJu3bpFRE1NTexXy009BQQEKBSK9qfFxcVhYWEdF/D19XV3d3d3dzc89fHxsbOzu3c9vr6+mEQEALrIkiVL7npFQJOIPPfcc2VlZYcOHdLr9Tt37rxx44ahRlNSUgoKCojo6aefbmpq2rt3LxEplcr09HRO0gMA8I6bGpVKpUlJSVOnTvXx8Vm8ePHXX3/t4OBARHFxcYav/ffp02fPnj3/+Mc/fH19Bw0aFBcXh1NMAPAQOVU5fEcwFmefOU6fPn3KlCkFBQVBQUHtH2VWVla2L/D8889PmjRJqVT6+vq6uLhwtV0A6H7SitJe/f7V/Lh8Yc7PfBcuvwwqkUhCQ0MfcjpIIpGEhYWhQwHgITQ6zcLDC683XF97ai3fWYyC79QDgLAUTok4tjT37F7HDac3XL19le84nUONAoCA1DTXXL16qW8TPWbrZyW2soj5mXE9JgAIyIpjK2igZnj4DDf5o8sidCuOrTh+9fgY/zF853oYYU0iMnHiRKlUSkTe3t5r11rGxyIAwJVcVe6j2x+NDY/96C8fEVGLpmXQtkHOts7nXz9vJeLg/srtNmzYcPHiRSIqKys7ceIEy7UJ62hUKpXi8nuAHmv+oflandbT0XPH+R2GV4b1G6bIVey8sHPO0Dkcbmjx4sWGB5xcwC6sGgWAHuu7vO+OFB8horePvn3Xj1YcWzEjeIaLRKAX+aBGAUAQRvmMKppX9KCf9rLqZc4wjKBGAUAQ+jr07evQl+8UpsAFTwAArKBGAQBYQY0CALCCGgUAS9WiaWnT8j+rm7BOMTU3Nx86dIiInJycIiIi+I4DAIL2P6f/x8HGYdHIRUwHXrhwoaqqiogMN8NnSVg1qtVqDXf2N8yXAgDwINcbrq87tc5abD0zbCbTU/z19fWGqtFqteyTCKtGHR0dcVd8ADDGW0fe0ul1zZrmd396d8eUHYzGjhkzxvDA8PaXJXw2CgCW50z5mS9/+/LNiDfjhsftvLgzsyKTxzCoUQCwMDq9bsGhBf2d+y+NWLriyRXu9u6L0xbr9fprddd4yYMaBQALk/Rb0rnr59aOX+tg49DbtvfqcatPlp6c/f3syJRIXm5ZhxoFAEvS0Nbw9pG3B0kHPe79eHFtcXFt8Vj/sQNcB3yR9cWZ8jNfZn9p/kjCOsUEAPBwa35eU9lYWdlYGfhx4L0/fevIW88HPW/mifBQowBgSUZ6j9w/fX/HV2403pifOn+efvgLTb77z6WsHbL2/bHvmzMSahQALMlfH/nrXa/M/n62jZVNfO1jbhs/jRCLHEat//uQv/u7+JstkrBqVKVSGa4bxSQiAGCMC5UX/jfrf5dFLHM7JiYisUhsLbZemr70riPWu/XdwYUAAAgzSURBVHScRIR9BmHVKCYRAQDj6fX6uNQ4G7FNeP/wn17qZTX1GyLR04XJilzFydKTT/g98aCBmEQEAICIaG/O3tNlp4lo6tdT7/rR/EPzM+dkcjsR3oOgRgHAUg3vPzxzzgO/v9SmbbOztjNDDNQoAFiqQLf7XPNkfhZ/+f3BgwcLCgr4TsGZEydOGD757h4uXLjw888/852CM/n5+ampqXyn4Mz169cVCgXfKThTV1e3e/duXjZt8UejRUVFffr04TsFZ8rKytzd3flOwZkbN27cvn2b7xScqampKS4u5jsFZxoaGvLz8/lOwZnW1tbs7GxeNm3xR6MAAPwSVo2Wl5czHXLt2jXD7VeNV11dnZGRwXRD33//PdMhBw4cYHqjhPLy8hs3bjAa0tTUdOTIEUZDyKTdSU9PV6vVjIZUVlYy/Z3qdLoffviB0RAyaXdOnTpVU1PDaIhKpbp2jfE9hEzIZsKQzMzM69evMxpy+/bt0tJSphsy7Q9Bp9MxGpKTk1NU9MA56++rubnZhN+OCZ1zL2HVaGFh4ZYtW7Zs2fLVV18ZOeTq1as3b95ktJWqqqrjx48zzWbCp0jffPMN03trl5aWMv29NjY2mvCBnQm7c/DgwaamJkZDysvLmf6frdFovv32W0ZDyKTdOXbsmEqlYjTk5s2bV69eZbohE7KZMOTs2bNMrySvqakpKSlhuiETsn377bdM/xCysrKUSiWjIU1NTcbvzg8//GCoGk7OrAjrs1ErK6vhw4cTkb29Pd9ZAKDbksvlUqmUiKysOLiwVMTL7fkeJCIiIjg4mNGQzMxMDw8PHx8f44c0NTXduHFj4MCBjDb022+/hYWFMR0SGhoqEokYDbG3tw8MZHAZR1tbW3Fx8Z/+9Cem2ZjuTl5e3oABA3r16mX8EKVS2dbWFhISYvwQnU53+fLl0NBQRtlM2J2ioiJPT08HBwa3AiorK6uqqho6dGhXZzNhSGlpae/evV1cXIwfolKpCgoKRo0a1dXZsrOzQ0JCGP0hVFRU2NjYGJrOSE1NTefOnRs7diyjbBcvXvz1118ZDbmXsGpUq9W2trbynQIAegobGxsbGxuWKxFWjQIAWBxhnWICALA4qFEAAFZQowAArKBGAQBYQY0CALBiwTWqVCojIyO9vb0nT57M9AsPQtBp/kuXLj355JMuLi4ymezTTz81f0JGjPx11NfXy2SyPXv2mDGaKTrdndra2jlz5nh5efn6+m7cuNH8CRnpdHfOnTs3a9asfv36RUVFmfBVaSFYtGjRqlWreNm0pdaoRqOJjo728/PLyMiIiIgYP348069786vT/G1tbU899dSkSZNycnKWL18+d+7cwsJCvtJ2yvhfxxtvvFFaWirwy+yM2Z2FCxcabim0adOm+Ph4Tqb06SKd7o5Op4uKigoODj516pRcLo+Kirpz5w5faU2QkZERHx+/ZcsWpl9W5ozeMqWnp7u7u2u1WsNTFxcXhULBbyRGOs1/9uxZDw+P9qcRERGrV682a0QmjPx1JCUlTZo0KSIiYteuXeYNyEynu2P4J628vNzw9McffywsLDR3SqN1ujvZ2dk2NjZqtVqv19fX14vF4tzcXB6Cmio+Pj4yMtLHx2fp0qW8BLDUo9H8/Pzw8HCx+Pf8crk8NzeX30iMdJo/MDDw6NGjhsdNTU2XL19+/PHHzZ3SaMb8OoqLixMSEnbt2sXoS4G86HR38vLyBg4cePr06RdffPH9998fNGgQ0+8Wm1OnuxMcHBwYGBgbG3vw4MHXX389NDQ0KCiIj6QmWrly5d69e3n8A7HUGi0qKnJ1dW1/6u/vX19fz2MepjrN36dPH8PtBS5fvvznP/952LBh48aNM3dKo3W6OxqN5qWXXtq4caOnp6fZ0zHW6e4UFhaqVKpt27ZNmDAhKysrKCiIrxsGG6PT3RGJRDExMXv27Jk3b97evXvHjBlj7ogWzlJrVCKRVFdXtz9Vq9UymYzHPEwZk1+v1y9fvvzPf/5zTEzM4cOHhXwQ1+nurFmzxtXVdcCAATk5OWq1+vr16ybc6dJsOt2dlpaW5uZmhULx2muvpaSkDB48eN++fWaPaaxOd+fAgQMrV67Mz88vLCxUKpWffvqphZ5l4oul1mhAQEDHu7oWFxczvesMvzrNr9Pp5syZk5mZmZubu3DhwvZ3ZMLU6e40NDTcuHEjJiYmJiZGqVR+/vnnQj673enu+Pr6uru7t0/34uPjY2dnjhkoTdPp7pw+fXrUqFGGbpXJZE888UR3mkHLHHj5RJa9qqoqW1vb1NRUnU6XmJjo5ubW2NjIdygGHpRfoVAolUq9Xn/ixIm+ffuWlpbW/ofhDIAwdbo7HY0dO1bgp5g63Z1bt245Ozt/+eWXer0+Pz/fxcWlpKSE38wP0enuJCUl2dnZXblyRa/Xl5WVOTg45Ofn8xyauenTp/N1islSa1Sv1+/fv9/Ozq5///69e/dOT0/nOw5j983v6em5fft2vV6/evXqu/7BW7RoEa95O/Hw3elI+DWqN2J3vv32W2dnZx8fHysrqxUrVvCX1CgP3x2dTrdw4cK+ffsGBgZaWVnNmzeP17Am4rFGLftGeS0tLQUFBUFBQdbWwrqNv5EsPf9detrutLS0KJVKX19fRjdL5kunu6PT6bKysry8vCziNKCgWHaNAgDwTtAnLgAAhA81CgDACmoUAIAV1CgAACuoUQAAVlCjAACsoEYBAFhBjQIAsIIaBQBgBTUKAMAKahQAgBXUKAAAK6hRAABWUKMAAKygRgEAWEGNAgCw8n+eDVcLM00JhwAAAABJRU5ErkJggg==",
      "text/plain": [
       "Winston.FramedPlot(...)"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(Winston.plot(t, real(g),\"g^\",tnew,real(itc[tnew]),\"r.\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAAEsCAIAAADfNCTgAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deVgTd/4H8E84JAHLHQXl1lBRoNqKWlm1KrbVrbtW6xYLVP3V6lrBu2pLtXg8oq1HXXVtPaoWfKpCbbe2oFURq3iviBxKuAUUuVSOcGd+f6TLsnjAZIbMBN6vZ/8gyXy/856lvJ1kMjMShmEIAAC0ZSB0AAAA/YYaBQDgBDUKAMAJahQAgBPUKAAAJ6hRAABOUKMAAJygRgEAOEGNAgBwghoFAOAENQoAwAlqFACAE9QoAAAnRkIH+B937969fv260CkAoKvw8PDw8PDgOIm49kb/9Kc/sR2yadOmuLg4VkMKCgp++OEHtivatm0b2yHbt29vampiNWT37t3R0dGshjx69OjAgQOshpBWm7N///6KigpWQ44ePbp3715WQxobG3fu3MlqCGm1OVFRUffu3WM15NSpU5s3b2a7Ii2yaTEkNjZWqVSyGnLjxo2wsDC2K9Ii244dO9j+Ifz+++83btxgNaSgoODjjz9mNYSIxo0bx3bIUzBi4uLiwnbIO++8s2/fPlZD0tLS1q1bx3ZFAQEBbIdMnz69oaGB1ZC///3va9euZTWkqKho8eLFrIYwWm3OwoULi4uLWQ0JCwubN28eqyF1dXUzZ85kNYTRanPWrFlz584dVkN27979t7/9je2KtMimxZDt27dfunSJ1ZBff/117NixbFekRbYZM2bU19ezGhIZGRkTE8NqSGpq6ssvv8xqCKNV5zxJXHujAAB6x1CLvfqOEx8f7+/vz2pIZmamt7e3i4tL+4e88MIL/fv37969O6sVDRo0yMrKitWQgQMHWltbsxpSUFDg4uLC6sMaqVTq5eVlbm7OakVabM6AAQPkcrmBAYt/ektKSiwsLF555ZX2DzEwMBg4cKClpSWrbFpsjru7u729vZERi8MDjx8/ZhiG7UdPWmTTYoirq6uDg0O3bt3aP0SlUpWVlbF9V6v1H4JEImn/EAcHBxcXFxMTk/YPaWhoUCqVEydOZJVNi855krgOMbGtNiJ64YUXZDIZqyFSqdTOzo7tilxdXdkOYVXuGt27dzc1NWU1xMjIyMHBge2KtNgcR0dHtkPMzMwaGhpYDZFIJFr8/6bF5tjb27MdIpPJtPhPVItsWgyRy+Vsh5iYmLzwwgtsR+nmD4Ht/gcRGRkZafHb0WLIk/CmHgA6iZrGGkHWK669US2MHj26R48eQqfgzdChQ6VSqdApeOPp6VlfXy90Ct64ubmx/fBEzOzs7CZMmCB0Ct78u/zf1yyvCbJqva9RT09PoSPwSaFQCB2BT05OTkJH4FPPnj179uwpdAreWFpavvzyy0Kn4KSoqsjY0NhGZtOobvw47uNUw9Sf7vw0qd8kHccQV43W1taeO3eOiMzMzAYPHix0HAAQtU/OfGJqbLpzws6vr3+dWpIqN5Uv/W3peMV4E8M2jk2lpKSUlZURUXV1NfcY4vpsVHOsTalU5ubmCp0FAETt2r1r3yV9t/vfu8/lnfs8/vM5xq/GG3845PesrZe2tjm2sLBQUzV1dXXck0gYhuE+C18CAwMjIyOFTgEAYscwzJ/2/ynnYU6DuqGyrrKJabr38P/kW3erDSSW67qnh6Tbd2/XNzF46Rxx7Y0CALTH9ynfX8y/uH7s+mW+y+qa6t5/6X25qZyIDCQGDeqGT898qsswqFEA0DPVDdXLTy/36eXz/kvv/5b1GxGdyz13c9rolIRjKeejJ7pP/C7pu2v3dHfUXlyHmAAAiKi2sXbX9V2Lhi166qsbL2wsrCg8OvXoceXx09mniSjrYdagaL+WyyyIXZDwfwmsTp3SGmoUAERn86XNq+NXv9Hnjf7y/q1eynuct+nipr7WfZMfJN96cGtq/6lEdP7u+eLq4gOTDrT8SLSmscbUmN05gdpBjQKAuBRWFm64sKFB3bD45OITgSdavbrs1LKaxpqM8ow5v8xp9VJcTtz+v+7XVcz/Qo0CgLh8cvqTRnXj4lcXb7m05RflL2+5v9Xy1fCx4eFjw5860EAizMEe3taqVCr9/f0dHBwmTJjw1MvHXr16dfr06b169Zo2bVpCQgJf6wWAzuRyweXUk5HfGL39hekkzx6eC08srGv6n692ulm5Pet/LpYugmTmp0YbGxsDAgKcnZ0TEhJ8fX39/PxUKlXLBdRq9bRp0wYMGHDhwgV3d/dp06axvfAPAHR6DMMsPLEw/ILJ+yu+NwyZ/9WbX2U9zNp+ZbvQudrAT43Gx8fn5uaGh4c7OzuHhoZWVlbGxMS0XCAtLS0/Pz8kJMTNzW3p0qWFhYWZmZm8rBoAOo3vbn13pfBK82Glsa5jJ7pPXHNuzf2q+8IGez5+PhtNT0/38fFpvqCvu7t7WlpaywUGDBjQt2/fefPmvfPOO5GRkV5eXtxvIwUAnUlVfdWnZz6VkOSz8SbmI31qTQxzIsY9qHpQWV+56uyqPRP3CB3wmfip0aysrJYXxHZxcWl17zOJRBIUFBQaGvr7779nZWUtWLCAl/UCQKdxMutkf3n//vL+hUSF/3myZ/eePbv3vPv4bnF1cQ8zkV4Sk58alUqlmsulaKhUqlYXfDt+/Pjq1avT09MVCkVGRoanp+fUqVN9fX1bzaNUKjXX1jY1NR00aJDmyf379xsbG/OSEwBEa4rHlCkeUzp0FcHBwY8ePSKi27dvl5eXExHbm008FT816urqGhUV1fwwOzvb29u75QIXL14cPny4plsVCsXIkSPPnz//ZI26u7tfvXqVl0gAIKx5MfM2v75ZaiSiy5Dv2LGj1TOBgYHcp+XnENOkSZPy8/NPnDjBMMy+ffuKioo0NRodHZ2RkUFEAwYMuHz58p07d4iooKDg0qVLkydP5mXVACBCsZmx/7z2z62X275mXSfAz96oXC6PiIiYPHmytbV1VVVVdHS0mZkZEYWEhISFhSkUioCAgBs3bowaNcrc3DwnJ2fevHnu7u68rBoAxEZzAhIRrT+/fvpL03u90EvoRB2Lt7OYpk6dOnHixIyMDA8Pj+ab1t6//8fXFCQSyZYtWzZt2pSUlGRvb6/FjTkBQF9sv7L9TumdnRN21i5ZkHVmVK9pK+iDD4QO1YH4PHdKc8P059z428DAYNCgQehQgE6sRFWy9ve1f1b8+SOfj2ZmdB9xJrM09gehQ3UsXG8UAPi0Mm5ldX315jc2E5F5H497loZnKm6K6i4bvEONAgBvkh4k7b2xN2RoyIs2LxKR4YWLp+L3+fveP5R8SOhoHQhXeAIA3iw8sbCJaTqccvinOz9pntHsh644veJtj7fNjM0ETddRUKMAwI/cR7nuNu7uNk//Es6N+zdGOI3QcSTdQI0CAD9cLF2+eesboVMIAJ+NAgBwghoFAOBEXG/qGYbRXO/ZwMBAKhXRqbgAXZpaTZcvExG5uFCvznBKUl1dXVNTExGp1Wrus4mrRsvKyhYuXEhEvXr1CgsLEzoOABARUU0NaS4k9MUX9PHHQqfhwZ49e27dukUtzrTkQlw1amtru3v3bqFTAHQ5FXUV5ibmQqfQneDgYM0PIrrCEwDor8SixMlHnnvFNamUTp6kkyfpb3/TVSh9Iq69UQDQMYZhFsQuOH/3/LHbxyZ7PKNMDQ3p9dd1m0ufYG8UoEs7mnb0/N3z5ibmS39bWttYK3QcvYQaBei6ahprlp9aPshu0PFpx3Mf5W65tEXoRHoJNQrQdW28sDHvcd5Xb3410nnklP5Twi+E36u8J3Qo/YMaBeiiCioKNl3cNM1z2kjnkUS05Y0takb9yZlPhM6lf1CjAF3UslPL6prqZr8yO/thdvbD7IamBn9P/4ikiCuFV4SOpmdwpB6gK3nwgI4eJaLrPr0PpxxmiBl9cHSrRRbELrj0wSWJRCJEPr2EGgXoSnJyaP58IrKI/MeRqUeetVRFfYWFiYUOY+k31ChAV6SwUSj6vyl0ik4CNQrQlQweTKWlRETmXejUz44mrhp9/Pjxtm3biMjGxoaXc10B4H8YGZGNjdAhhPfzzz/n5OQQUanmHxVuxFWjUql0+PDhRCSTyYTOAgCdloeHh729PRH9/vvv3GcTV42amJj4+PgInQIAOjmFQqH5gZc9NnxvFACAE9QoAAAnqFEAAE7E9dkoALRXUxM1NBAR4a5lQsPeKIB++vZbkslIJqPCQqGjdHWoUQDxalA3CB0B2oY39QAi9dOdn1JLUkNHhD795SFDKDyciMgCJ78LDDUKIEZ1TXUfn/q4qKpo5sCZvV542q3hX3qJXnpJ57ngKfCmHkCMtl3ellmeWddYt+L0CqGzQBtQowCi86D6wfrz6yf1m7TMd1nkrcgLdy8InQieBzUKIDqfnP6kprHmi3FffDriU0cLx4UnFqoZtdCh4JnE9dloSUmJ5sJODg4OGzZsEDoOgAASixIPJh1c5rtMYa0gonVj1r3/4/uHkg8FeQcJHa3z2Lx5c2JiIhHl5+dzn03CMAz3WfgSGBgYGRkpdAoAwTAMM+rAKGWZUhmiNDcx1zwzYv+IrIdZ6cHpmmeAR7x0jrj2RgG6uCOpR87fPU9EFhtaf43pi4Qv1o1ZJ0QoaANqFEBoBQVUXExGRuTt3ahu3OD39I+zpEZShmFwpzkRQo0CCO2LL2j7drKyovLyQG/c9EH/4Eg9AAAn2BsFENrs2eTnR926CZ0DtIQaBRCapyd5egodArTH25t6pVLp7+/v4OAwYcIEpVL55AIPHz6cPXu2vb29k5PTli1b+FovAICw+KnRxsbGgIAAZ2fnhIQEX19fPz8/lUrVaplFixbV1dUlJydv3bp11apVvHzrFQBAcPy8qY+Pj8/NzQ0PDzcwMAgNDd20aVNMTMw777zTvEBWVtbBgwcLCgpsbW2nTJkik8nq6+t5WTUAgLD4qdH09HQfHx8Dgz/2bd3d3dPS0loucOfOnT59+ly8eDE6OtrT0zMoKMjFxYWXVQMACIufN/VZWVlWVlbND11cXCoqKloukJmZWVJSsnPnznHjxiUlJXl4eCQnJ/OyagAAYfGzNyqVSsvKypofqlQqhULRcoHa2tqampqoqCi5XD5r1qxXX3318OHDXl5ereapqKj45ptvNBM2zzBs2LDm/VwAAK1dv35d83FiXl6eZlevtLSU+7T81Kirq2tUVFTzw+zsbG9v75YLODk52dra2traah46OjrKZLIn56mqqjp79iwRmZmZVVdXa54cOnQoLyEBoItLSUnRHP1OTU3V7Pm1et+sJYYPxcXFJiYmsbGxarV679691tbWVVVVDMNERUUplUqGYUpLS83NzQ8dOsQwTHp6uqWlZU5OzpPzBAQE8JIHQBgffcQMHcrMnSt0DmgvXjqHnzfLcrk8IiJi8uTJjo6OS5YsOXLkiJmZGRGFhITExcURkY2NzYEDB+bOnevk5NS/f/+QkBAcYoJOKC2Nrlyh1FShc4BO8XYW09SpUydOnJiRkeHh4WFk9Me09+/fb17g7bffHj9+vFKpdHJysrS05Gu9ACIyZAiZmND/fqIFnR6fJ4NKpdInjxq1WsAb/4VBJ7Zxo9AJQAA4Ag4AwAlqFOB5KusrGTHdaAdECDUK8DzLTy0/knpE6BQgaqhRgGdKLk7ec2PP8tPLVQ2tL7UD0Aw1CvBMi04skhnJ8h/nf3nxS6GzgHihRgGe7tjtY2dyzmwct/Fdz3c3XtiY9zhP6EQgUqhRgKeob6pfcXrFAPmAD1/+8MtxXxpIDD45/YnQoUCkUKMATxERsdQ8JWPXS58aGRg5mDssHb70cMphzR3kAVqRiOrLHG+88UavXr2IqHfv3uvWrRM6DnRR96vum/Tsba1iKmYFlX4RRkSqBtWoA6PcrNyuzLpiIMHOh97btm3bzZs3iSgnJyc+Pp7jbOK6pZ1cLt+9ezcRSSQSobNA1/XpmU83EUNEEbcigv8R0fx8eU35waSDMwfOFC4a8CM4OFitVhPRjBkzuM8mrholImNjY6EjQJemZtRvub+Vuq23YX2jonePowrHlq/amtoKFQx4ZGhoaGhoSDztsYmuRgGEZSAxmOIxhTymCB0E9AY+5QEA4AQ1CgDACWoUAIAT1CgAACeoUQAATlCjAACcoEYBADhBjQIAcIIahU6trIyGDqWhQ+mnn4SOAp2W6M5iamhoICKJRNJ8l2YA7dXX09WrREQPHggdBUSkqalJc049L9dmEldVlZSUzJ49m3CFJ+CLVEpvvklE5OQkdBQQkR07dmiu8FRYWMh9NnFdKC8wMDAyMlLoFADQVfDSOfhsFACAE9QoAAAnqFEAAE5QowAAnKBGAQA4QY0CAHCCGgUA4AQ1CgDAibjOYgJ4pnv36LffiIj+/GeSy4VOA/Bf2BsFMfrxzo9NTNP/PJWSQjNn0syZlJEhUCiAp0ONgujcr7o//cfp+27sEzoIQLuI6019RUXFP//5TyKytrb29/cXOg4I49Mzn1bWV648u/Jdz3ctTCz+eHbkSMrJISKytxcwG3QOMTExubm5RFRWVsZ9NnHtjZqYmAwcOHDgwIH9+vUTOgsI49q9a98lffd2v7fLVGVrzq357wtSKbm4kIsLmZgIlw46CTc3N03VyGQy7rOJa2/UxMRk+PDhQqcAwTAMs/jkYrmp/MCkAx+f+nj7le2zX5n9os2LQueCzqZ5R83U1JT7bOLaG4Uu7lDyoQt3L2wct9HcxHzdmHVm3cyWnFwidCiANqBGQSxqGmtC40JfsX8lyDuIiOSm8pUjV/6a8WtsZqzQ0QCeR1xv6qErCz8ffvfxXRdLl7m/ztU809DUICHJkpNL/Nz8jA2MhY0H8CyoURCF+1X39yXus5JaJT9ITn6Q3Py8pdSyqKro4M2Ds16eJWA8gOdAjYIo2He3L1zMw11xAHQPn40CAHDCW40qlUp/f38HB4cJEyYolcpnLVZRUaFQKA4cOMDXegEAhMVPjTY2NgYEBDg7OyckJPj6+vr5+alUqqcu+dFHH+Xl5YnqdqSgazdv0o4dtGMHNTQIHQWAB/zUaHx8fG5ubnh4uLOzc2hoaGVlZUxMzJOLRUZGlpeXDxkyhJeVgr46c4ZCQigkhGprhY4CwAN+ajQ9Pd3Hx8fA4I/Z3N3d09LSWi2TnZ0dFhb27bffSiQSXlYKACAG/NRoVlaWlZVV80MXF5eKioqWCzQ2Nr733ntbtmyxs7PjZY2gx+bNo5ISKimh7t2FjgLAA36+8CSVSlteKEWlUikUipYLrF+/3srKys3NLSUlRaVSFRYW5uXlOTs7t5pHqVS6uroSkamp6aBBgzRP7t+/39gYX73uRKRSkkqFDgFdUXBw8KNHj4jo9u3b5eXlxNM59fzUqKura1RUVPPD7Oxsb2/vlgtUVlYWFRUFBQURUWZm5p49e0pKSrZt29ZqHnd396tXr/ISCQCglR07drR6JjAwkPu0/LypnzRpUn5+/okTJxiG2bdvX1FRkaZGo6OjMzIyiOjLL79M/A8fH5+wsLAnOxQAQB/xszcql8sjIiImT55sbW1dVVUVHR1tZmZGRCEhIWFhYa3e4AMAdCa8nQw6derUiRMnZmRkeHh4GBn9Me39+/efXDIuLo6vlQIACI7Pc+qlUqmXlxePEwIAiB/OqQcA4AQ1CgDACWoUAIATXG8UOHjwgCoqyNiYXFyEjgIgGOyNAgfLl5O7O732mtA5AISEGgUA4ARv6oGDv/+dxo3DFUagixNXjdbV1V27do2IZDKZp6en0HGgLcOG0bBhQocAYC0jI0NzjZKamhrus4mrRmtray9evEhENjY2qFEA6CC3b9/Oyckhourqau6ziatGLSwsFixYIHQKAOjk/vKXv2h+0Lz95QiHmAAAOEGNAgBwghoFAOAENQoAwAlqtKtrUONm8QCcoEa7NIZhJh+ZXKoqFToIgB5DjXZph5IP/aL8ZfW51UIHAdBjqNGuq7qhesXpFcYGxl9f/7py/Fjq04dmzxY6FID+QY12XeHnw+9V3jv27jFLqWXRnWuUnU1Pu3cWADwfarSLyq/I33p5a6B34Fvub30+6vPDzpX5fxlFY8YInQtA/6BGu6jFJxcT0box64hors/cqHe9XhubXzf/I6FzAegfcZ1TX1pa+uGHHxJR7969w8LChI7TaSXkJ/yQ9sPaMWudLJyIyFBi+NWbX439buy2y9uW+S4TOh1Ah9uxY0dSUhIRFRYWcp9NwjAM91n4EhAQsHfvXiIyMDAwMTEROk7n1MQ0Dd49+GbRzSdfMjcxTw9Ot+tup/tUALpUX1/f1NRERLNmzTp06BDH2cS1NyqRSGQymdApOrl7lff8Pf39Pf2f+mre4zzUKHR63bp10/wgkUi4zyauGgUdcDR3XO67XOgUAJ0HDjEBAHCCGgUA4AQ1CgDACWoUAIAT1CgAACc4Ut/p1NcTw5CRERkaCh0FoEvA3min0707SaW0cqXQOQC6CtQoAAAneFPf6axbR01NNHy40DkAugrUaKezDNcWAdApcdVoSUlJYGAgETk4OGzYsEHoOADQOW3evDkxMZGI8vPzuc8mrhqVy+WRkZFCpwCATm7JkiWaHzT7bRzhEBMAACeoUQAATlCjAACcoEYBADhBjQIAcIIaBQDghLcaVSqV/v7+Dg4OEyZMUCqVTy5w8+bNUaNGWVpaKhSKXbt28bVeAABh8VOjjY2NAQEBzs7OCQkJvr6+fn5+KpWq5QL19fVjx44dP358SkpKaGhocHBwZmYmL6vutBobqbCQCguppkboKADwPPzUaHx8fG5ubnh4uLOzc2hoaGVlZUxMTMsFEhMTjY2NV6xY4eDgMGPGjFdfffXIkSO8rLrTysoiBwdycKCffhI6CgA8Dz81mp6e7uPjY2Dwx2zu7u5paWktF+jbt++ZM2c0P1dXV6empg4bNoyXVQMACIufk0GzsrKsrKyaH7q4uFRUVLRcwMbGxsbGhohSU1ODgoIGDx48ZswYXlbdadnZ0e7dRERDhggdBQCeh58alUqlZWVlzQ9VKpVCoWi1DMMwn3322a5du1auXLlgwQKJRPLkPPfv3w8KCiIiMzOzl19+WfPkBx98YNgFL+RuYUEffih0CIBOJTIyUnPYJjU1tby8nIhycnK4T8tPjbq6ukZFRTU/zM7O9vb2brmAWq2eM2fO3bt309LS7OzsnjWPTCYbMmQIEUml0hdffFHz5FMLFwCArb59+9bX1xORiYmJ5h1zcXExD/MyfCguLjYxMYmNjVWr1Xv37rW2tq6qqmIYJioqSqlUMgxz7ty5Hj165OXlPfwPlUr15DwBAQG85AEAaA9eOoefQ0xyuTwiImLy5MmOjo5Lliw5cuSImZkZEYWEhMTFxRHR+fPni4uLnZ2drf7js88+42XVAADCkjAMw9dctbW1GRkZHh4eRkZaflYQGBiI640CgM7w0jl8XrZZKpV6eXnxOCEAgPjhnHoAAE5Qox0uLicu+2G20CkAoKOgRjtWo7pxfuz8pb8tFToIAHQU1GjH+vr616klqT/e+TEuJ07oLADQIVCjHai8pvzz+M8n2o18S+Xw9d6/N9VUC50IAPiHGu1Aq86uqqir2GH01+NfFBxdnxH98wahEwEA/1CjHSWtJG33v3fPHTzXycJJ88z2K9vLasqePwoA9A5qtKMsPrm4e7fun7/2OQ0fTj/+mH/wH2mm1WvPrRU6FwDwjM+zmLh744035HI5ETk4OGzYoMdvgX9O//mvh//qauk6uNfg5ifjcuIq6iqS5iZ52HoImA0ANm/enJiYSET5+fnnzp3jOJu4arRznAyqZtTjD43PLH/6XVJe7/P6rj/jVlQAoiC6k0FBw0BicDLwpNApAEBH8NkoAAAnqFEAAE5QowAAnKBGAQA4QY0CAHCCGgUA4AQ1CgDACb43ykZ2NmVnk0RCY8cKHQUAxAI1ysbBg7RmDRGRmE79AgBh4U09AAAn4tobra+v11wvQCaT9evXT+g4T/jgAxo3jiQSoXMAACfZ2dmPHz8mopqaGu6ziatGa2pq4uLiiMjW1laMNerkRE5OQocAAK5u3bqVlZVFRNXVPNyTQlw1amFhsWTJEqFTAEAnN2nSJM0Pmre/HOGzUQAATlCjAACcoEYBADhBjQIAcIIaBQDgBDUKAMAJahQAgBPUKAAAJ+L6+r3uXLhA+flkaUnjxwsdBQD0W1et0S1b6McfycsLNQoAHOFNPQAAJxJGTJfOfPPNN93c3IioV69en332WQeu6d49qqwkqZScnTtwLQAgSrt27UpOTiai27dvnz17luNs4npTb21tvWHDBiIyNDTs2DX16tWx8wOAiM2YMaOhoYGI5syZw302cdWogYGBubm50CkAoJOTyWQymYx42mPDZ6MAAJygRgEAOEGNAgBwghoFAOAENQoAwAlqFACAEx3VqFKp9Pf3d3BwmDBhglKp1M1KAQB0QBc12tjYGBAQ4OzsnJCQ4Ovr6+fnp1KpdLBeAAAd0EWNxsfH5+bmhoeHOzs7h4aGVlZWxsTE8DZ7RARt306//cbbhAAAbOjiLKb09HQfHx8Dgz8q293dPS0tjbfZ16yhzEx67z16/XXe5gQAaDdd7I1mZWVZWVk1P3RxcamoqNDBegEAdEAXe6NSqbSsrKz5oUqlUigUT11SqVS6uroSkamp6aBBgzRP7t+/39jY+JmzX7tGTU1kYsJnYgDojIKDgx89ekREt2/fLi8vJyJTU1Pu0+qiRl1dXaOiopofZmdne3t7P3VJd3f3q1evspvd0pJLNgDoOnbs2NHqmcDAQO7T6uJN/aRJk/Lz80+cOMEwzL59+4qKip5VowAAekcXe6NyuTwiImLy5MnW1tZVVVXR0dFmZmY6WC8AgA7o6Ov3U6dOLS8vj42NLS0t9fPz43HmmJiYjIwMHicU1rlz5xITE4VOwZsbN26cP39e6BS8SU9Pj42NFToFbwoLC1t+2qbvHhHIXq0AAAcMSURBVD9+vH//fkFWrbvLNkulUi8vL96nzcrKsrGx4X1aoeTn59va2gqdgjdFRUWaT/Q7h/Ly8uzsbKFT8KaysjI9PV3oFLypq6vT3BdE93BOPQAAJ+Kq0YKCArZD7t69W1paympIWVlZQkIC2xX961//Yjvk+PHjbO8YWFBQUFRUxGpIdXX16dOnWQ0hrTbn1KlTbM/ivX//PtvfqVqt/uWXX1gNIa0258KFC5qvvLRfSUnJ3bt32a5Ii2xaDLl+/XphYSGrIY8ePcrLy2O7Iu3+ENRqNashKSkpWVlZrIbU1NRo8dvRonOeJK4a1eKXmpub++DBA1ZDiouL4+Pj2a5Ii0+Rfvjhh6amJlZD8vLy2P5eq6qqtPjATovNiYmJqa6uZjWkoKCA7X/ZjY2Nx44dYzWEtNqcs2fPlpSUsBry4MGD3NxctivSIpsWQy5fvpyfn89qSHl5eU5ODtsVaZHt2LFjbP8QkpKS2F7DqLq6WovN0aJzniSuGgUA0Dviuk+9r6/vgAEDWA25fv16z549HR0d2z+kurq6qKioT58+rFZ069Yttl93vXXrlpeXl0QiYTXE1NS0b9++7R9SX1+fnZ3dr18/ttnYbs6dO3fc3Ny6devW/iFKpbK+vt7T07P9Q9RqdWpqKtujkVpsTlZWlp2dHavv3uXn5xcXF7/yyisdnU2LIXl5eRYWFpZszkYpKSnJyMgYPnx4R2dLTk729PRk9Ydw7949Y2NjuVze/iHV1dVXr14dPXo0q2yJiYnXrl1jNeRJ4qrRpqamuro6oVMAQFdhbGz8vHPN20dcNQoAoHfw2SgAACeoUQAATlCjAACcoEYBADhBjQIAcKLHNarvN21uM//NmzdHjRplaWmpUCh27dql+4SstPPXUVFRoVAoDhw4oMNo2mhzcx4+fDh79mx7e3snJ6ctW7boPiErbW7O1atXp0+f3qtXr2nTpmlxqrQYLF68eM2aNYKsWl9rVN9v2txm/vr6+rFjx44fPz4lJSU0NDQ4ODgzM1OotG1q/6/jo48+ysvLE/nX7NqzOYsWLdJcUmjr1q2rVq1ieyKmLrW5OWq1etq0aQMGDLhw4YK7u/u0adMaGhqESquFhISEVatWbd++ne3Jyrxh9NOpU6dsbW2bmpo0Dy0tLaOiooSNxEqb+S9fvtyzZ8/mh76+vuvWrdNpRDba+euIiIgYP368r6/vt99+q9uA7LS5OZp/0goKCjQPf/3118zMTF2nbLc2Nyc5OdnY2FilUjEMU1FRYWBgkJaWJkBQba1atcrf39/R0XHZsmWCBNDXvdGOvWlzx2szf9++fc+cOaP5ubq6OjU1ddiwYbpO2W7t+XVkZ2eHhYV9++23rE4KFESbm3Pnzp0+ffpcvHjx3XffXbt2bf/+/dmeW6xLbW7OgAED+vbtO2/evJiYmDlz5nh5eXl4eAiRVEurV6/+/vvvBfwD0dca1febNreZ38bGRnN5gdTU1BEjRgwePHjMmDG6TtlubW5OY2Pje++9t2XLFjs7O52nY63NzcnMzCwpKdm5c+e4ceOSkpI8PDyEumBwe7S5ORKJJCgo6MCBA/Pnz//+++9fe+01XUfUc/pao+2/abM4tSc/wzChoaEjRowICgo6efKkmHfi2tyc9evXW1lZubm5paSkqFSqwsJCXi5Q1kHa3Jza2tqampqoqKhZs2ZFR0cPHDjw8OHDOo/ZXm1uzvHjx1evXp2enp6ZmalUKnft2qWnR5mEoq816urq2vKqrs+5abM4tZlfrVbPnj37+vXraWlpixYtan5HJk5tbk5lZWVRUVFQUFBQUJBSqdyzZ4+Yj263uTlOTk62trbNt3txdHSUyWQ6jchGm5tz8eLF4cOHa7pVoVCMHDmyM91BSxcE+USWu+LiYhMTk9jYWLVavXfvXs09R4UOxcKz8kdFRSmVSoZhzp0716NHj7y8vIf/oTkCIE5tbk5Lo0ePFvkhpjY3p7S01Nzc/NChQwzDpKenW1pa5uTkCJv5OdrcnIiICJlMdvv2bYZh8vPzzczM0tPTBQ7N3tSpU4U6xKSvNcowzNGjR2UyWe/evS0sLE6dOiV0HNaemt/Ozu7rr79mGGbdunWt/sFbvHixoHnb8PzNaUn8Ncq0Y3OOHTtmbm7u6OhoaGi4cuVK4ZK2y/M3R61WL1q0qEePHn379jU0NJw/f76gYbUkYI3q94XyamtrMzIyPDw8jIx0d4tTHul7/la62ubU1tYqlUonJydWF0sWSpubo1ark5KS7O3t9eIwoKjod40CAAhO1AcuAADEDzUKAMAJahQAgBPUKAAAJ6hRAABOUKMAAJygRgEAOEGNAgBwghoFAOAENQoAwAlqFACAE9QoAAAnqFEAAE5QowAAnKBGAQA4QY0CAHDy//d77uLUcAXdAAAAAElFTkSuQmCC",
      "text/plain": [
       "Winston.FramedPlot(...)"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(Winston.plot(t, imag(g),\"g^\",tnew,imag(itc[tnew]),\"r.\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAAEsCAIAAADfNCTgAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3de1xUdf4/8DfDcBtARWDkmoxTCuiqqUQrrevli1lqW7nb6paIW9/twVaul0gjKbAMy5D99muJlJRdm1raviFZXssLWancRLnYAgo6DoLcB2YY5nK+f5y+fPsplzmemTkcfD0f/jGX8/mc92dmeHnmXD7jxDAMAQDA7ZIIXQAAgLghRgEAeEGMAgDwghgFAOAFMQoAwAtiFACAF8QoAAAviFEAAF4QowAAvCBGAQB4QYwCAPCCGAUA4AUxCgDAC2IUAIAXqQ37Sk9Pnzhx4tKlS299qrm5OTs7u7CwMCoqav369a6urkSk0Wi2b99eXl4eEhLy0ksvRUREEBG7GNtKJpNlZGTYsEIAAJuzzdaoRqNRqVRJSUl1dXX9LrBy5cqysrK4uLiCgoK4uDgiMplMc+bM6e7uXrNmDcMw8+fPN5lMRJSTk+Pl5aVUKpVKpUKhsEl5AAD2Y5ut0cTExKqqKnYb81YnT54sKipqaGiQSqXR0dGBgYFpaWmtra16vX7nzp1E9NBDDwUHBx8+fHjx4sU1NTUHDx709va2SWEAAPZmm61RlUpVUlISFRXV77OVlZUxMTFSqZSIAgICfHx8iouLfX19c3Jy2AU6Ojq0Wm1YWJhWqzWZTHv27ElISPjqq6/MZrNNygMAsB9b7hsdSHV1tb+/f99dpVLZ1tYWFhYWFhZGRLW1tU899dTChQsnT55cWlra1tZWW1sbGRmZkpLyyiuvnD592t3d/aYOc3Nzr169an0BZrO5vb3d19fX+iY9PT1SqZSNfit1dXV5eXlZvzzDMDqdztPT0/omnZ2dLi4uHh4e9quKiLq7uzlVZTAYJBKJi4uL9U0aGxvHjRvHqSquAzEajQzDDPQNySarIO4DMRqNFovFzc3N+ibd3d0ymczJycn6JlwH0tPTYzAYRo8ebX0TnU7n7u4ukXDYDuNaldls7u3t5fRpb21tHTVqFKc/26CgoD/84Q/WL98vR8RoQEDAz1PPYDD07fTcunVrRkbGxo0bN2zYQEQRERH19fUhISFEtGLFiuDg4EOHDj366KM3dbhx48Zdu3YRkVQqZUPW2dlZJpMNVEBlZWVaWtrevXutr/njjz+ePn16ZGSk9U02b978xhtvWL+8VqvNzMzcuHGj9U2Sk5MjIiIee+wx+1VlNBrffPPN1157zfomhw4d8vLyeuCBB6xv8vDDDx84cMD65Yn7QE6cOGE0GmNjY61vkpycnJqayikauA7kzJkzjY2NixYtsr7Jm2++uXbt2kE+3rfi+lrl5eUVFxe//fbb1jfJyMiIi4vjtGnCtaqamppTp05x+rQ//fTTCQkJs2bNGmgBnU7Hfsft6ekxGo1E9J//+Z/iiNGoqKjc3Fz2ttlsVqvVU6dOZRjmhRdeqK6uLisrCw4OZp9tbGxkx0ZEfn5+ISEh/f5f1NXV9dJLLxGRXC5nj++HhYWtXbt2oAIMBoNMJpsyZYr1NQcGBiqVSk5NxowZw2l5dgOZUxMfH5/AwEC7VtXb2+vj48OpSUVFxahRozg1cXV15bQ8cR/I5cuXe3t7ub68U6ZM4RSjXAfS1NTk4uLCqcnYsWMjIyM5bcdxfa3Onj07evRoTk38/PzCw8M5bYlzrYqIampqODXx9PRUKBSDNHnvvfdqamqI6Mcff7x+/ToRdXV1cSqpX3aM0fz8/BkzZoSGhs6aNauioqKwsHDmzJk7duyQy+VyufzEiRP5+fklJSVeXl56vZ6IXF1d1Wr10qVLy8vLg4KC8vLytFrtggULbu3Zx8entLTU+kqcnJw4fSciIolEchtNOC3v5OR0G02GYVW38VpxXZ5uayAOeNNvY3kHvLzD87VyzB/I4As8//zzNz1yzz33cFpF/xjbmTdv3rvvvtt318fH55NPPmFv79mzRyaThYaGyuXysrIyhmG2bNlyUyU7d+5kGCYhIWH8+PH33HOPp6fnxx9/3O+KVq1axakwtVqdlJR0e4MaVj788MOTJ08KXYUNcH0Hh62RMZDTp09nZmYKXYUNpKam1tbWcmpik3fQllujx44d+/nd1tbWvtvx8fFPPPFEXV1deHg4+z9McnJycnLyrZ1kZmbq9fq6ujqlUsnp+AAAgCAcsW+UJZPJrDxi4+Hhwe7xBAAY/pwYhhG6Bs7i4+P7zjm1htls7urq4nQ+x/DU3d0tlUo5nS4zPLW2to4dO1boKmxgZAykt7fXaDRyOstteOrs7PT09HR2dra+Cdcw6ZfjtkYF5OzsPAIylIhGwAedNQKihzUyBuLq6joydqCNGjVKkPVihicAAF5EuTV67dq1Z555hojuu+++P/3pT0M3yMmhvDxyd6f/PX0VAO5ku3fv/v7774mI0/WQAxFljAYHB2dlZZH1Z6JVVNAXXxCX60AAYASLj49np5p7+umn+fcmyhglIk6XzVJwME2fTlwuzgWAEUwikbAbYbdxJUg/vfHvQgTWrqXSUvr+e6HrAIAR6M6IUQAAu0GMAgDwghgFAOAFMQoAwIsoj9QbDIaLFy8S0ahRo4KCgoQuBwBEpqGhoaOjg4gMBgP/3kQZo52dnfv27SOiyMjIRx55ROhyAEBkSkpKLly4QERsmPIkyhj19/fftGmT0FUAgFgtXrx48eLFRMR+r+UJ+0YBAHhBjAIA8IIYBQDgBTEKAMALYhQAgBfEKAAAL6I84UmtVsfHxxNRdHR0QkKC0OUAgMhkZ2efOnWKiK5cucK/tzviJ+0AAPplkzDBl3oAAF4QowAAvCBGAQB4QYwCAPCCGAUA4AUxCgDAiyjPG+3u7j5+/DgRjRs3LjIyUuhyAEBkLl682NDQQERdXV38exNljBqNRvakWWdnZ6FrAQDxaWlpYTOkt7eXf2+ijNExY8asWrVK6CoAQKxiYmJiYmKIiP1eyxP2jQIA8GL3GE1PT9+/f3+/TzU3N2/btm3ZsmXbtm3r27TWaDTr1q2LjY1dvXp1VVWVvcsDAODJjjGq0WhUKlVSUlJdXV2/C6xcubKsrCwuLq6goCAuLo6ITCbTnDlzuru716xZwzDM/PnzTSaT/SoEAODPjvtGExMTq6qqXF1d+3325MmTRUVFDQ0NUqk0Ojo6MDAwLS2ttbVVr9fv3LmTiB566KHg4ODDhw+zvzwFADA82TFGVSoVEc2fP7/fZysrK2NiYqRSKREFBAT4+PgUFxfPmjWrb7aVjo4OrVYbFhZmvwoBAPgT7Eh9dXW1v79/312lUtnW1hYWFsbmZm1t7VNPPbVw4cLJkycLVSEAgDUEi9GAgICrV6/23TUYDAqFgr29devWjIyMjRs3btiwod+2586du/vuu4nI399/0qRJRKRUKpOTk+1fNQCI2LZt29gfpq+urm5sbCQimUzGv1vBYjQqKio3N5e9bTab1Wr11KlTGYZ54YUXqqury8rKgoODB2o7ffp0TNsMAFxt2rTppkfY39HgydHnjebn57MbobNmzaqoqCgsLLRYLDt27JDL5XK5/OTJk/n5+R999NHYsWP1er1erzebzQ6uEACAE0fH6OrVq7/77jsi8vb2zsrKmjt3blhY2DvvvPPpp58S0bfffqtWq+Vyuex/7d6928EVAgBwIvBvMel0urq6uvDwcImEQ6Djt5gAwCZsEiYCX1Mvk8kwRRMAiBquqQcA4AUxCgDAiygnymtvb//kk0+IaPz48bNnzxa6HAAQmTNnzly6dImI2tvb+fcmyhiVSqVjx44lIi8vL6FrAQDx8fT0ZDPEJlO/izJGvby8HnzwQaGrAACxmjJlypQpU4iI/V7LE/aNAgDwghgFAOAFMQoAwAtiFACAF8QoAAAviFEAAF5EecKTWq1mZwmMjo5OSEgQuhwAEJns7OxTp04R0ZUrV/j3JvAMT7cHMzwBgE3YJEzwpR4AgBfEKAAAL4hRAABeEKMAALwgRgEAeEGMAgDwIsrzRpubmzMyMoho0qRJDz/8sNDlAIDIHD58uLKykohu3LjBvzdRxqiXl9f8+fOJyMfHR+haAEB8IiMjAwICiOj06dP8exNljLq7u0+bNk3oKgBArEJDQ0NDQ4nIw8ODf2/YNwoAwAtiFACAF8QoAAAviFEAAF4QowAAvCBGAQB4EeUJT5i2GQD4wLTNmLYZAGwD0zYDAAjPll/q09PTJ06cuHTp0lufam5uzs7OLiwsjIqKWr9+vaur60Ct2MXY2zKZjL12HgBg2LLN1qhGo1GpVElJSXV1df0usHLlyrKysri4uIKCgri4uEFa5eTkeHl5KZVKpVKpUChsUh4AgP3YZms0MTGxqqrq59uYP3fy5MmioqKGhgapVBodHR0YGJiWlqZQKPptVVNTc/DgQW9vb5sUBgBgb7aJUZVKRUTsrEu3qqysjImJkUqlRBQQEODj41NcXKxQKG5tpdVqTSbTnj17qqqqlixZsmjRImdnZ5tUCABgJ444xFRdXe3v7993V6lUtrW19btkTU1NW1tbbW1tZGRkSkrKzJkze3p6HFAhAMBtc8R5owEBAVevXu27azAYBtrpGRERUV9fHxISQkQrVqwIDg4+dOjQo48+etNilZWV999/PxHJ5fKIiAgiuuuuu5577jl7DQAARoQPPvjg0qVLRHTx4sXGxkYiMplM/Lt1RIxGRUXl5uayt81ms1qtnjp1ar9LNjY2Go1G9rafn19ISEi/swEqlcpXXnmFiNzc3EaNGkVEA+2WBQDo88QTT/T29hKRVqtlv+m+8cYb/Lu1Y4zm5+fPmDEjNDR01qxZFRUVhYWFM2fO3LFjh1wul8vl/TZRq9VLly4tLy8PCgrKy8vTarULFiy4dTE3N7cpU6bYr3IAGJH6fi9j3Lhx7A13d3f+3doxRlevXp2Zmbl8+XJvb++srKy5c+f6+voaDIajR48O1CQmJmb58uWzZ892dXXVaDS7du1iD0wBAAxbjrsYVKfT1dXVhYeHSyRDHNfS6/V1dXVKpXKgr+q4GBQAbMImYeK4bT2ZTBYZGWnNkh4eHuyBIwCA4Q/X1AMA8IIYBQDgRZQHcMxm840bN4jI3d0dl40CAFd9JzyJ5rxRm7tx48bbb79NRNOnT3/yySeFLgcARObAgQNFRUVE1NTUxL83TNsMAHcuTNsMACA8xCgAAC+IUQAAXhCjA3vkEVq8mPbuFboOABjWRHmk3kEOHCCzme69V+g6AGBYQ4wOLCqKTCYKCRG6DgAY1kQZo2q1Oj4+noiio6MTEhLstZoffrBXzwAgqOzs7FOnThHRlStX+PeG80YB4M6F80YBAISHGAUA4AUxCgDAC2IUAIAXxCgAAC+IUQAAXkR53mhPT8+FCxeIaMyYMaGhoUKXAwAio1ar29raiIidvJknUcZoV1fXkSNHiCg8PBwxCgBclZeXV1RUEJFWq+Xfmyhj1M/Pb8OGDUJXAQBitWjRokWLFhER+72WJ+wbBQDgBTEKAMALYhQAgBfEKAAAL4hRAABeEKMAALyI8oQnB03bDAAjFKZtxrTNAGAbmLYZAEB4do/R9PT0/fv39/tUc3Pztm3bli1btm3btt7eXitbAQAMK3aMUY1Go1KpkpKS6urq+l1g5cqVZWVlcXFxBQUFcXFxVrYCABhW7HiIKTExsaqqytXVtd9nT548WVRU1NDQIJVKo6OjAwMD09LSFArF4K0AAIYbO8aoSqUiovnz5/f7bGVlZUxMjFQqJaKAgAAfH5/i4mKFQjF4KwCA4UawQ0zV1dX+/v59d5VKJTv9HwCAuAh23mhAQMDVq1f77hoMBoVCYWXbpqamt956i4h8fHzYVmPGjImKirJHnQAwYhQVFbGba/X19S0tLUTU2NjIv1vBYjQqKio3N5e9bTab1Wr11KlTrWxrMBjq6+uJSKfTeXp6EpHFYrFTnQAwYrS3tzc1NRGRWq1mb4hy9vv8/PwZM2aEhobOmjWroqKisLBw5syZO3bskMvlcrncyk5CQ0MzMzPtWicAjDz/8R//cdMj7PWQPDk6RlevXp2Zmbl8+XJvb++srKy5c+f6+voaDIajR486uBIAAJuwe4weO3bs53dbW1v7bsfHxz/xxBN1dXXh4eESiWSQVgAAw5bAU5PIZLLIyEhhawAA4APX1AMA8IIYBQDgBTEKAMALpm0GgDsOpm3GtM0AYBuYthkAQHiIUQAAXhCjAAC8IEYBAHhBjAIA8IIYBQDgRZTnjba3t+/du5eIwsLCfvWrXwldDgCIzPfff19bW0tENvnRDVHGqFQqDQoKIqKxY8cKXQsAiM+YMWPYDHFxceHfmyhj1MvLa8GCBUJXAQBiFRkZyc4tx36v5Qn7RgEAeEGMAgDwghgFAOAFMQoAwAtiFACAF8QoAAAvojzhCdM2AwAfmLYZ0zYDgG1g2mYAAOEhRgEAeEGMAgDwghgFAOAFMQoAwAtiFACAF1GeN9rc3Lx9+3YiioiIWLJkidDlAIDIHDx4sLy8nIiampr49ybKGPX29n744YeJaPTo0ULXAgDiM3Xq1LvuuouIiouL+fcmyhh1c3ObPHmy0FUAgFgFBwcHBwcTkbu7O//ebBmj6enpEydOXLp06a1PNTc3Z2dnFxYWRkVFrV+/3tXVdaAH2UfYVjKZLCMjw4YVAgDYnG0OMWk0GpVKlZSUVFdX1+8CK1euLCsri4uLKygoiIuLG+TBnJwcLy8vpVKpVCoVCoVNygMAsB/bbI0mJiZWVVWxm5O3OnnyZFFRUUNDg1QqjY6ODgwMTEtLu3Llyq0PKhSKmpqagwcPent726QwRzMayWAgIvLyEroUAHAQ22yNqlSqkpKSqKiofp+trKyMiYmRSqVEFBAQ4OPjU1xc3O+DWq3WZDLt2bMnISHhq6++MpvNNinPcbKyyNubvL2ppUXoUgDAQRxx3mh1dbW/v3/fXaVS2dbW1u+DNTU1bW1ttbW1kZGRKSkpM2fO7OnpcUCFAAC3zRFH6gMCAq5evdp312AwKBSKtra2Wx+MiIior68PCQkhohUrVgQHBx86dOjRRx+9tU+LxUJETk5O7N2+GwJ74AHavp2IyNNT6FIAoB99U4PacI5QR8RoVFRUbm4ue9tsNqvV6qlTpzo7O9/6YGNjo9FoZB/08/MLCQnx8PC4tcNz585NmjSJiPz9/dkbSqVy8+bNDhjLEO69l+69V+giAKB/b7311sWLF4mourq6sbGRiPpNGK7sGKP5+fkzZswIDQ2dNWtWRUVFYWHhzJkzd+zYIZfL5XK5h4fHrQ9+9913S5cuLS8vDwoKysvL02q1CxYsuLXn6dOnY9pmAOBq48aNNz3C/o4GT3aM0dWrV2dmZi5fvtzb2zsrK2vu3Lm+vr4Gg+Ho0aNE1O+DMTExy5cvnz17tqurq0aj2bVrF3sMCgBg2HLcj4jodLq6urrw8HCJRDL4g3q9vq6uTqlUDnQGFX5EBABswiZh4rhtPZlMFhkZac2DHh4eERERjqoLAIAXTJQHAMALYhQAgBfEKAAAL6I8Dt7Q0LBhwwYiuvfee5966imhywEAkfnnP//JziSn0Wj49ybKGJXL5UlJSUTk5uYmdC0AID5Lly6NjY0lorVr1/LvTZQx6uzs7OvrK3QVACBWnp6enp6eROTs7My/N+wbBQDgBTEKAMALYhQAgBfEKAAAL4hRAABeRHmk3mQysWd7eXh4+Pj4CF0OAIhMe3u7Tqcjor4JjvkQZYy2tLS89957RDRt2rTf//73QpcDACLz9ddfl5SUEFFzczP/3hw3UZ4NYaI8ALAJm4QJ9o0CAPCCGAUA4AUxCgDAC2IUAIAXxCgAAC+IUQAAXkR53qharWZ/XTo6OjohIUHocgBAZLKzs0+dOkVEV65c4d8bzhsFgDsXzhsFABAeYhQAgBfEKAAAL4hRAABeEKMAALwgRgEAeBHleaN6vZ6dK3Ds2LFhYWFClwMAIlNfX9/S0kJEer2ef2+ijFGdTseeOjtx4kTEKABw9e9//7uqqoqIuru7+fcmyhj19fVds2aN0FUAgFjFxsbGxsYSEfu9lifsGxXesL2QbLgVNkzqEaoMx6zX5mvh06E1bX++jFBvzTCN0fT09P379wtdhSN09XatzFspdBX9aOxufP7g80JX8f95+ounm3U2+OUcPhq7G5878Jzj1/vd1e8yTmeIbi295t4V/73Cwlhuo63JYlrx3ytMFtPgiz134DmNVsPefvXEq+cbz9/GungadjGq0WhUKlVSUlJdXZ3QtTjCm9++qbqg2ndxn9CF3GzT15veL3z/B/UPQhfyk2OXj+05t+e1E68JW8bLX7+cVZTl4JfFwljWHFyz5eSWpu4mu67lhQMv2HYtfz3919yK3H+U/eM22r5f9H5uRe6ukl2DLHPm2pmsoqyXv3mZiC63X37n+3fWHV53m7XyMOxiNDExMT093dXVVehCHOFS26WM0xlO5PTikRcNZoPQ5fyfQk0h+9Ffe2jt7W1K2JaZMa87vM6JnHYW77zQdEGoMoo0RX8v+zsR/eXgXxz5suwu3V3SUNJp6Nx8bLNd11J6vdSGa7nedX3rt1udyCnpmyRtr5ZT2xZ9S8qJFCdyevX4q209bf0uwzDM2kNriWhv2d4z1868eORFg8lw7PKxvIt5Nqiei2EXoyqVqqSkJCoqSuhCHCHxaKKLxGXv43tr22ozfrD7VzYrsZ/OcZ7jMhdnnr12du/5vUJXRDuLd55vPJ/9SPYY9zHrDgmwuUFEDMP85dBfxnmOe3/J+33/zThAh6Fj87HNc8bPWf/L9WzS2Wktrxx7pW8t566f49/nK8deMZgMn/z2k8buxq0FWzm1ffX4q1qD9p+//Wd7T3vqidR+l/nowken1af/tvhvQd5Bv//X7z+v+jxT8fxvTRO3/2udozdKmGFp3rx577777kDPrlq1yoG12Ms3l76hFEr7No05c2bzq7+asd5To9UIXRTDMIzqvIpSKOdcjsVi+WX2LwPfCew0dApYT6u+1e9tvwV/X8AwTGZhJqVQXlWe48tgX5Y9pXssFsvsD2c77GV58ciLklRJSUNJe0+7fLt8zp459ljLhsMbnFOdbbiWIk2RJFWy8ehGxmTKSI79zUqX+tITVrY933g+/C/OO1OWMl988eev/uyyxaXyRuVNy2gN2uD04Ohd0RaL5R9l/6AUmvBfEyw+Yxiid6Mp7ds0K9dlkzARZYxOmzZNqVQqlcr7779/1apVq1at2rJliyPL489kMU19f+qE/5qgN+oZhYIhUk13jt8XL3RdjM6ouyvjrpkfzDRbzAzDFGuKJamSl79+WcCS1hxc45zqfL7xPHPT6+ZA/b4sm77eZO/11rTWuL3u9qf9f2Lv7ireRSmUW55r27VUt1S7ve727P5n2bs7i3dSCn1a8eltd2ixWH61+1fjto/r6OlgtFqGiCHKeWqKlc0X7l345kIPtlWLtsn3LV/2P9GfS/omySnF6bT6NMMwfzv7N0qhgHcC2Bj9apHS602va53Xbu05LS2NTYzZs2ezGfKLX/zitofZR5QxOgK2Rtmtqt2lu1v1reaw8QzRmfmTJKmSM+ozwhaWfCyZUuirf3/Vqm9l//3u09+5ve5W21orSD0VTRUuW1zi8uL66vms8rOftuIdSKiXZcnHSzy3el5svsiutFnXHPm3yPEZ43VGnV3XEvFeBJ+1fHzhY0qhjB8yWvWtbTeusoH44kI6Wnt0yLafV31OKVTwdCzbqrXrxusnX6cU2v/j/r5lLrdddn/DfVnuslZ9a317vd/bfoq/KiiFcl7/rfaff6/45hOXLS6r8lZZU6pNwmSYzn4/f/78xx577IUXXuj3WbHPft/W0zbx/03sO3fnfjV5GOm6F1X50y9DfvndH79zcnISpLD6jvqI9yL0pv+7PM6zl6ZfJyKaMHvxP57+0vElLfpo0eHaw7c+7u3q/eMLPwZ6BTqghltflj6PhT/2+e8/t9N6D9ceXvTRon6fSp2b+uqvX7XJWg7VHHpI9VC/T22ZtyV5TjLXDnVGXfh74Vc7r7J3JQzN0hARXRlNfhOmlD5bKpUMeNWPwWyI/FvkpbZLAV0U3ElEVBz001MTfSdeSLjg6uxKRL/71+8+q/xs8DIkTpLTz5yOChriKItNwkSUVzGJXWNX49b5A+5xb+1p9fXwdWQ9fVp0LX9d9NefP+L349XH33yDiI5MHq836T2kHo6sR9urfTzi8ccjHu/32RvdNxwToy26FtVd68JOnCOic/EPGmXuP3/Wfi+LEzl9sOSDfp/ydPW01VokThLbruV61/XNcwY81n9DN9i71tjVuDFm40DPNnU3hYwK6TH1xE6IjZ0QO2QleqMNrpe3xjDdGh2c2LdGxaS0lGbMICLat49+8xuhqxHOzp307LNERNeuUVDQUEuDaGBrFOzv7rvpyBEiounThS5FUBIJOTsTEQm0v8X2Ojro88+JiObNI+vn9ykqotpakslo6VK7VSY+iFEYlLc3xQ797Wnke+YZeuYZoYuwqWvX6I9/JCLKzeUQox98QNnZFBSEGP05xCiA+DU3U08PubmRv7/QpdyJRBmjWq32wIEDRBQUFDT9Dv+yeYdTqejTT8nZ+afvpyPAt9/S008TEf3rXzRtmrWtnnySjhyh+++nH6y+2P+ee4idtoJT8r7+Oq1fTy4uHJoMS+fPn1er1UTU2dnJvzdRxqjFYunq6iKinp4eoWsBQVVV0Rdf/LTXcmTo7qbqaiIiW8zKPhgXFxo/nnOrgAAKCLBDNY7W09PDZojFYoO5EUQZo6NHj37iiSeErgKGgaAgmj6dpKL8GPcvKIhWrCAi8vPj0GrdOvrd72jcODsVNfLcd9999913HxGx32t5GkGfP7gD/fnP9Oc/C12ETU2dSh9/zLnVov7P0gfHQIwC2E1FBen1NHYsTZggdClgR8NuojyAkeO3v6WoKHr5ZeJ1gPMAAAk+SURBVKHrAPvC1igMD93dtGEDEdHy5TR3rsDFAHCBGIXhoaeHPviAiGjy5JETo1lZpNXi4tERDzEKw4NE8tOBZk+bTbohvF//WugKwBFEGaNqtTo+Pp6IoqOjExIShC4HbMHHh65fd9C6cnLIaKRp0+i++xy0RhhmsrOzT506RURXrlzh3xtmeII7j6cn6XT04ou0fbvQpYDAMMMTcLRnD12+TAoFrV4tdCkAIwdi9E7y97/TyZP061/bPUYzMuibbygggLKz7bui21NTQxYLeXsLXQeMEIhRcSospG3biIjeeYcUCmtbSSTk5EQS+58sXFpKX33FoTAHC3TEnPlw50CMipNG89OcRpsH/LWGfhw7ZqdybqZQ0MyZFBLioNUBCAoxKk6jRtGUKUREHg79cSRrpaZSaqrQRQA4CGJUnObNowsXhC4CAIhEGqNtbW27d+8mogkTJswdMVe8AICjfPvtt9XV1UTU2trKvzdRxqirq+vdd99NRHK5XOhaAEB8xo0bx54y7+rqyr83Ucaop6fnnDlzhK4CAMRq4sSJEydOJCL2ey1Pd8REeRqN5tVXXxW6ChvIyclhr2ATu2dGyq9sjoyBFBYW7ty5U+gqbGDr1q2XL192/HrviBhlGMZsNgtdhQ1YLBab/HSM4Ewmk9Al2MbIGIjFYhkZfyBms1mQq9vviBgFALAfUcZocXExp+Vv3Lhx9uxZTk327t1bUVHBqcmmTZs4La/Vardu3cqpSXFx8cWLFzk14VqV0WhMTk7m1OTAgQMFBQWcmnB9B4n7QI4fP3748GFOTV5++WWu2zJcB3L69Ol9+/ZxarJlyxadTsepCdfXqrq6uqioiFOT7du3t7S0cGpyG1V9+OGHnJqcPXu2oaGBU5Pb+CjeSpSHmPQcf3vWbDZz/SC2t7dz/fXm6xzneTObzTdu3ODURK/XGwwGTk24VsUwTGNjI6cmWq3WycmJUxOu7yBxH0hXV1dvby+nJo2NjQzDcBoL14HodDquP4ze1NTEdU8O19eqt7eX6x9Ic3Mz1x0aXKsyGAxtbW2cmuh0Oq5V3cZH8VaijFGGYTi96z09PRaLhVOT3t7enp4eTk1MJhOn5fV6PdcmJpOJ6yee6yp6e3u5NjEYDAaDgVMTru8gcR+IwWC4vddKwmXOAa4DcUxVt/Gmc93UMBqNer3erh/Fnp4eo9HIqYnZbOb6Z2uTfaminG80NTWV0xZTb29vU1NTCJdLvLu6utzc3FxcXKxv0traOnbsWOuXt1gsnZ2dY8aMsb5Jc3Ozm5ubN5epibhWRURtbW0+Pj7WL6/X652cnNzd3a1vcvnyZQXHiUu4DoTdyvDgcrHsbbxWXAdiMBjMZrNMJrO+SXt7++jRozltI3MdSFdXl06n43QWdkdHh7e3N6dw51qV0Wjs6enh9GnXaDS+vr5ubm7WN/H19X399detX75fooxRAIDhQ5SHmAAAhg/EKAAAL4hRAABeEKMAALwgRgEAeBn5Mdrc3Lxt27Zly5Zt27aN6/nYAhq8bI1Gs27dutjY2NWrV1dVVQlSoZWsef31ev0jjzxSUlLi4NqsN+Qoenp6MjMzlyxZsmLFitLSUsdXaKUhB9LQ0LBjx47HHnvs/ffft8lcnI6Unp6+f/9+x6935MfoypUry8rK4uLiCgoK4uLihC7HWoOUbTKZ5syZ093dvWbNGoZh5s+fP5wnyLDm9V+/fv3+/fu5Xt7jSEOOIjk5OT8/f926df7+/r/5zW+G7QwyQw7kwQcfrK6ufvLJJ/Pz80U0f5VGo1GpVElJSXV1dQKsnhnRTpw44efnZzQaGYZhL7a9dOmS0EUNbfCyi4qKgoKC2NtGo1Eul3/55ZfCFDoUa17/vLy82NhYhUJx/PhxAUq0wpCjqK+vd3Z2bmhoYO+mpqb23R5WhhzIuXPnPD09TSYTwzCVlZXOzs46nU6YWjn6wx/+cO+993p5eb377ruOX/sI3xqtrKyMiYmRSqVEFBAQ4OPjY5OZCOxt8LJ9fX1zcnLY2x0dHVqtNiwsTIgyhzbk63/t2rWXXnopJyeH61X5jjTkKMrLy5VK5cWLFzdt2pSXl7dp06aAgACBih3MkAOJiIjw8PBIT08vLy/fvn17TEwMp8vABKRSqUpKSqKiogRZ+wiP0erqan9//767SqWS62QHghi87LCwsNjYWCKqra1dsmTJwoULJ0+eLECVVhh8IBaLZeXKlWlpaUFBQUJUZ60hP0XV1dUdHR2JiYnu7u6vvfZaRETE8NxBMeRAXF1dX3zxxY0bN95///179uyJj493dIniNMJjNCAgoL29ve+uwWDgeim3IKwpe+vWrdHR0Y8//vjn7A/WD0uDDyQ9PX306NExMTHXr183m80tLS1arVaIMocw5Nuh1Wrb29sPHz6ckpJy7ty5zs7OL7/80uFlDm3IgRw6dCglJaWoqKizs/Ozzz579tlnuc6Gd2ca4TEaFRV16dIl9rbZbFar1VOnThW2JGsMXjbDMM8//3xBQUFZWVliYiKn6SEcbPCBXL58+cKFCw888MADDzyg0WjWrFmTmZkpUKWDGfJTFBoaOn78eHbeDYlEolQqh+ebMuRAjh8/vmDBgpkzZ0okkmXLloWGhnKdSfYO5fjdsY7U2dnp5uZ29uxZs9n89ttvT5o0SeiKrDJQ2fv27bty5crx48dDQkKampp0/4s9JjAMDT6Qny85adKkYXuIachRaDQamUx25MgRhmFKS0vHjBnT3d0tZMUDGHIg2dnZPj4+arWaYZjy8nJvb++Ojg4hK+Zo3rx5ghxiGuExyjDMnj17ZDJZaGioXC4vKysTuhxr9Vu2j4/PJ598smXLlpv+L9y5c6ew1Q5ikIH8fLHhHKOMFaP46KOPPD09lUqlm5tbVlaWcJUOYfCBmEymhIQEPz+/iIgIV1fX7du3C1osZ0LF6B0xUZ5Op6urqwsPDx+eX7UGItKybzUyBjLkKHQ6XW1trUKh8PLycnBtnAw5EJPJdOHChbvuusvX19fBtYnUHRGjAAD2I+KtAwCA4QAxCgDAC2IUAIAXxCgAAC+IUQAAXhCjAAC8IEYBAHhBjAIA8IIYBQDgBTEKAMALYhQAgBfEKAAAL4hRAABeEKMAALwgRgEAeEGMAgDw8j91GdiwM3diPwAAAABJRU5ErkJggg==",
      "text/plain": [
       "Winston.FramedPlot(...)"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(Winston.plot(t, abs(g),\"g^\",tnew,abs(itc[tnew]),\"r.\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAAEsCAIAAADfNCTgAAAABmJLR0QA/wD/AP+gvaeTAAAgAElEQVR4nO3deUAUZ54+8G83tDSnXK2gnApERIgGiUZGjVGT6MQZ19EMCvzM/MYjUTTBYzYJEdG4mt1EjYNOXC9MwE0M6MzEBFA88EANMaIgKN2AIBCRo0WO5mq69o/OsAQ1UF1FVwHP5y+6ut73/VYIj3W/EoZhCAAADCUVugAAgL4NMQoAwAliFACAE8QoAAAniFEAAE4QowAAnCBGAQA4QYwCAHCCGAUA4AQxCgDACWIUAIATxCgAACeIUQAATkyFLuAX7t27d+3aNaGrAICBwtfX19fXl2Mn4tob/c1vfsO2ySeffHL27FlWTcrKyo4dO8Z2oF27drFtEhsb297ezqrJvn37kpKSWDWpra09fPgwqyZk0ObExcXV1dWxavL1118fOHCAVROtVrtnzx5WTcigzUlMTPzpp59YNUlLS9u+fTvbgQyozYAmKSkpSqWSVZPr16/HxMSwHciA2nbv3s32D+HChQvXr19n1aSsrGz9+vWsmhDRzJkz2TZ5AkZMPDw82DaZP3/+wYMHWTXJy8vbsmUL24FCQ0PZNlm8eHFbWxurJm+++eaHH37IqklFRcWaNWtYNWEM2px33nmnsrKSVZOYmJiVK1eyatLS0vKnP/2JVRPGoM3ZvHnznTt3WDXZt2/f66+/znYgA2ozoElsbOyVK1dYNfnuu++mT5/OdiADanvjjTdaW1tZNUlISEhOTmbVJDc397nnnmPVhDEocx4nrr1RAIA+x8SAvfrek56eHhISwqpJQUFBQECAh4dHz5tYW1uPHj3aysqK1UDjxo2zs7Nj1WTs2LH29vasmpSVlXl4eLA6WSOXy/39/W1sbFgNZMDm+Pn5KRQKqZTFP71VVVWDBw8ODAzseROpVDp27FhbW1tWtRmwOT4+Ps7OzqamLC4PPHr0iGEYtqeeDKjNgCaenp4uLi6DBg3qeRONRlNTU8P2qNbgPwSJRNLzJi4uLh4eHmZmZj1v0tbWplQq58yZw6o2AzLnceK6xMQ22ojI2tra3NycVRO5XO7k5MR2IE9PT7ZNWIW7npWVlYWFBasmpqamLi4ubAcyYHNcXV3ZNrG0tGxra2PVRCKRGPDfzYDNcXZ2ZtvE3NzcgP9FDajNgCYKhYJtEzMzM2tra7atjPOHwHb/g4hMTU0N+O0Y0ORxOKgHgH6iSdskyLhG3Rtds2aNra1tdHQ0j31OmzZtyJAhPHYorAkTJsjlcqGr4M2YMWNaW1uFroI3I0aMYHvyRMycnJxmz54tdBW8+VH94w+2PwgytJH2RjMyMqKjo2NjYxsbG/ntecyYMf0pRr29vQ04dhYtNzc3Ly8voavgzdChQ/38/ISugje2trbPPfec0FXwQ6vTrj+7/keTH/9x5x/GH91IMXrq1CmVSmXA2SgAgKe533C/WlNNRHuv7c2tylVYKNadWtfS3mLkMowUo5s2bfryyy8nTpxonOEAYCB47/R70eei1U3qjekbl8teSJctff5C4c4rO41chriu1AMA9FBmeeYXN7+QSqTFtcV1LXUfPvRX7NyaIJXYBm5dPHaxs5Xxjn3FFaNKpVJ/O4WFhcW4ceP0C+Pi4mQymaB1AYC4MAzzdurbztbOre2tKQUpq55fpciwISKpRNqma3v/zPtxv497vFVERERtbS0R3b59W61WExHb+wufSFwx6uPjk5mZKXQVACB2R3KOXC27+vnczw9cP3Dx3sURdiNuLgowmR1IEsmcsv/54uYXK4JWBA0L6tJq9+7dXZaEhYVxL0ZcMQoA0K3GtsZ3T787YfgEazPri/cuElHkycgu67yd8nbG/89g9eiUwRCjANDHfHTpo5/qf0p6PclCZpEWnna39u6yE8vmj56/PHB559WatE0WMh6O2btl1Bj9+uuvjTkcAPQ/xbXFn1z+xNvBO/tBdsdCD1uP1ILUXa/uGmY9zPglYW8UAMSoTdcmkz7h2vL6tPXN2mZljXL5t8u7fPXemfc+n/u5Uar7BcQoAIhOeX356pTVx15/wuvV983Zt2/Ovie2kpAxzoQ+DjEKAKLz3un3jt8+fkJ5Yo5P1xff2cnZvabPCPCGJwAQl6tlVxOyE2RS2dqTa1vb+8CrbRCjACAi+vvqh1kPS3o9SaVW/fX7vwpdUfcQowAgIl9kfyH5PvOIdMHv7tvM8Zmz+fzm+w33hS6qGxKGYYSu4f+89tpr+teXKBSK5cu7XoYDgP6tobXhmd3PfBXfNPnGQxo7tvBskt8ev/Bnw/fP2c/vQEePHlWpVER04cKFU6dOcexNXJeYLC0t58+fT0Ss5mABgP5h68Wt9+vvjxnyG6KLRDTSbuTqCau3X9m+PHD5+GHjeRwoODj42WefJaKbN29y701cMSqTyUaNGiV0FQAggLu1d3de3ckQE+h7cfAI0shuKDf9fAPTmpNrzr9xnscnOzumL+Nlj01cMQoAA1bpo9KYF2Oe9m2VpmqIpUjnuUCMAoAoTHGfMsV9itBVGAJX6gEAOEGMAgBwghgFAOAEMQoAvaKsrkzoEowEMQoA/GtobQg+FJxblSt0IcaAGAUA/m29uPXeo3vrT60XuhBjEFeMMgzT0NDQ0NDQ1NQkdC0AYKCih0U7r+50tXGtP5uSGbeFbt0SuqKumpub9VGj0+m49yau+0ZramrWrVtHRMOGDYuOjha6HAAwxLpT60ylphl/zpCO8Bp+aEN72G2T+CNCF/ULhw4dys7OJqL793l474m4YtTR0XHv3r1CVwEAhjt79+zf7/x92/RtrjauTRb2VFORX50/WuiqulixYoX+B14mWBbXQT0A9GntTHvkycgRdiPemfgOEZn/47v1GyeFjsmvaKgQurRehBgFAN7s+3Ff9oPs3z/z+0v3Lp0uOn3aVu3529AbVg0b0zcKXVovEtdBPQD0XQ+bH0afiyainVd37ry6s/NXB68ffGv8W2OdxgpUWu9CjAIAb1LDUp/2lcJSYcxKjAkxCgD8sJPbBToHCl2FAHBuFACAE8QoAAAniFEAAE4QowAAnIjrElN1dfWSJUuIaPjw4Zs2bRK6HAAgIiKdji5fJiLy9KThw4WuhgexsbH6OUHLynh4m5+4YtTBwWH37t1EJJViNxlANJqaaPJkIqL/+i9a3x9e2vTmm2+2t7cT0Z///GfuvYkrRiUSiVwuF7oKAOjnZDKZTCYjIl4mbRZXjAKAGJmbU1oaEZG3t9CliBFiFGCga2fazxeff8nzpaeuIZXSjBlGrKiPwSlIgIHuwPUDi44tqmupE7qQvgoxCjCg1TbXbji74UHjgy0XtghdS1/FW4wqlcqQkBAXF5fZs2crlcrHV8jMzFy8ePGwYcMWLlyYkZHB17gAwMXm85vVTerpntN3fb9LpVYJXU6fxE+MarXa0NBQd3f3jIyM4ODgGTNmaDSazivodLqFCxf6+fldunTJx8dn4cKFbW1tvAwNAAa7U31nd+buZYHL/ucP/2Nuar7u1DqhK+qbGD6kpaU5Ojq2t7frP9ra2iYmJnZeIScnRyaTaTQahmHq6uqkUmleXt7j/YSGhvJSDwD0xKyEWXYf2VU1VjEM83HGxxRDKaoUoYsyKl4yh5+90fz8/KCgoI575n18fPLy8jqv4Ofn5+XltXLlyuTk5OXLl/v7+/v6+vIyNAAY5jvVdykFKRtf3Oho4UhEqyes9nHwWXNyTZsOR4rs8HPDU2FhoZ2dXcdHDw+PurpfXPWTSCTh4eFRUVEXLlwoLCx8++23n9hPY2PjP/7xDyIaNGiQm5ubfqGfnx8vt8gCQIc2Xdvak2tNJCb3Ht179/S7+oXDrYefKz6399reVc+vEra8XpKfn68/nXj//v3GxkYievToEfdu+YlRuVxeU1PT8VGj0Xj/8jbdEydObNq0KT8/39vbW6VSjRkzZsGCBcHBwV36UavV+plBraysxo8fr184atQoU1Pc3wrAh8ZGunOHiP774cn8mnwi2nFlR5dVYtJjFvkvcjB3EKC8XpaWltbQ0EBEt27dqq6uJiK1Ws29W37iydPTMzExseNjUVFRQEBA5xUuX748adIkfbZ6e3tPmTLl4sWLj8eoq6trQkICLyUBwBPk5NALLxBRREpKxEZG6GqMLSIiossSEU2wPHfu3NLS0tTUVIZhDh48WFFRoY/RpKQklUpFRH5+flevXr1z5w4RlZWVXblyZd68ebwMDQAgLH72RhUKRXx8/Lx58+zt7RsaGpKSkiwtLYlo1apVMTEx3t7eoaGh169fnzp1qo2Nzd27d1euXOnj48PL0ADAwqhR9M03RESBA3HSpF4iYRjeduybm5tVKpWvr+/TTmXqdLqbN286Ozs7OTk9cYWwsDAc1AOA0fCSOXxeupHL5f7+/r+yglQqHTduHI8jAgAIDs/UAwBwghgFAOAEMQoAwAliFACAE8QoAAAn4nrIsqqqSv9QgYuLy0cffSR0OQDQP23fvj0rK4uISktLufcmrhhVKBS4bxSgR+rqqKKCiGjECMJLJ1hau3at/gcRPQwKAMZ29Cg98ww98ww9eCB0KQMdYhQAgBMcCwCIlFanrW+tt5PbPfnrl16iI0eIiOztjVkVPA4xCiBSe6/tzavK+9tv//bkr0eOpJEjjVsRPBkO6gHESN2kjkmP2ffjvuwH2ULXAt1AjAKIUUx6zKOWRzZmNu+kviN0LdANxCiA6Nyuvr332t43x7+5dfrWc8Xnjt8+LnRF8GsQowCiE5kaaTXIauPUjUsDlz479Nl1p9Y1a5uFLgqeCjEKIC7/zP/nycKTm6dtdrRwNJGY7Hx1593au49PPAfiIa4r9RqN5ttvvyUia2vrqVOnCl0OgLG1treuP7V+qOXQoOFBP97/kYhszGyeH/78tkvb3hj7xjDrYUIX2E/88MMPDx48IKIuU8EbRlwxyjCMfvpTExMToWsBEMCu73ep1CoimnhgYpev3j/z/uG5hwWoqT9qbm7WR41Op+Pem7hi1NLSMiQkROgqAIzr4UOqqyMTE3Jxmew2OS087YlrSSQShmEkEomRq+uXJk+erP9Bf/jLkbhiFGAg2riRYmPJzo7U6okuXXdCQfxwiQkAgBPsjQIILTSUxo0jMzOh6wADIUYBhDZhAk2YIHQRYDgc1AMAcIIYBQDgBDEKAMAJYhQAgBPEKAAAJ+K6Uo8JlgHACDDBMgAAJ/xOsCyuGAXo2+rrqa2NZDKytha6FDAenBsF4M/vfkcODvTaa0LXAUaFGAUA4AQH9QD8WbWK/u3fyMVF6DrAqHjbG1UqlSEhIS4uLrNnz1YqlY+v8PDhw2XLljk7O7u5ue3YgRkRoD+aN49Wr6Z584SuA4yKnxjVarWhoaHu7u4ZGRnBwcEzZszQaDRd1omMjGxpacnJydm5c2d0dDQv9xkAAAiOn4P69PT04uLibdu2SaXSqKioTz75JDk5ef78+R0rFBYWfv7552VlZY6Ojn/4wx/Mzc1bW1t5GRqgV+GF89AtfvZG8/Pzg4KCpNKfe/Px8cnLy+u8wp07d0aOHHn58uU//vGPH3744ejRo0eOHMnL0AC9avuV7ZfuXRK6ChA1fmK0sLDQzs6u46OHh0eX+fYKCgqqqqr27Nkzc+bMmzdv+vr65uTk8DI0QO950Phgy4Ut606tYxhG6FpAvPg5qJfL5TU1NR0fNRqNt7d35xWam5ubmpoSExMVCsWSJUteeOGFr776yt/fv0s/P/300+uvv05ElpaWgYGB+oVvvfUWJgoFQXxw9oO6lrrvy79PyEkIDwgXuhzgKi4urrGxkYhyc3P1kSWih0E9PT0TExM7PhYVFQUEBHRewc3NzdHR0dHRUf/R1dXV3Nz88X6srKxefvllIpLL5V5eXvqFHecKAIwpqyLrUNahtZPWXim9sv7U+t8/83sbMxuhiwJOnn32Wf1VmSFDhuiPmI8dO8ZDvwwfKisrzczMUlJSdDrdgQMH7O3tGxoaGIZJTExUKpUMw1RXV9vY2Bw5coRhmPz8fFtb27t37z7eT2hoKC/1AHA3JW7KkI+H1DbXXvvpmnSTNOpMlNAVAf94yRx+dvQUCkV8fPy8efNcXV3Xrl179OhRS0tLIlq1atXZs2eJyMHB4fDhw2+99Zabm9vo0aNXrVrl4eHBy9AAvSHxxy+uFVz4aOqHg80GBzoHhgWEfXL5kwJ1gdB1gRhJGP7OnTc3N6tUKl9fX1PTJ58raG5uViqVbm5utra2T1whLCwMb3gCwTVpm5ptrewadczKFZLde4joQeMDn1ifmSNmJr2eJHR1wCdeMofP045yudzf3/9pGapfISAg4GkZCiASH2d8rGN0RHTw+kH7/7S3/097392+mjbNsdvH0ovTha4ORAfP1AP8glanfdj8MC08eFCLtmGk4wK/4Z2/zanMedHjRYFKA5FCjAL8gqnUdOcrO+kVoeuAvgP3EgEAcIIYBQDgBDEKAMAJYhQAgBNxXWLCBMsAYASYYBkAgBN+J1jGQT0AACeIUQAAThCj0K9VVJBEQhIJ7dsndCnQbyFGoV/rePMOXl8PvUZcl5gAeGZtTdHRRETjxwtdCvRbiFHo16ysaNMmoYuAfg4H9QAAnCBGAQA4QYwCAHAirnOjarV6/fr1ROTk5NTxmAEAAL/i4uLy8vKI6P79+9x7E1eMDh48eM2aNUQkk8mErgUA+q25c+e++uqrRBQZGcm9N3HFqImJibOzs9BVAEA/Z2dnp//hV+aO6zmcGwUA4AQxCgDACWIUxKi4trjrorNnycaGbGwoM1OAggCeDjEKotPa3vpKwisXSi78YqlWS/X1VF9PWq1AdQE8mbguMQEQ0a7vdylrlGtOrslcmimV/Otfeg8P0l9UHTZMwNoAHocYBXF50Phgy4UtnraeP97/Me5G3J/H/fnnL3x8aMcOQUsDeDIc1IO4vH/m/WZtc2pY6hT3Ke+fef9RyyOhKwLoBmIURCSrIuvwjcOREyN9HHw+ffXTGk3Nf1z4D6GLAuiGuGJUq9WWl5eXl5dXVlYKXQsYG8Mwb6e8rbBQvDf5PSIa5zTuT+P+pD9PKnRp0N+o1Wp91LS1tXHvTVznRuvq6v76178SkZOTEy8PaUEfcjT36MV7FxeOWXiq8JR+ybNDn21rb1t3at03C78RtjboZ7755pvbt28TUU1NDffeJIyYJlcICwvDBMsDk6ZNM2r3qNK6J08afjLs5MsjXzZySTAQ8JI54tobhQFLq9P+PeTvT/vWwdzBmMUAsIIYBVGwMbMJdA4UugoAQ4jrEhMAQJ+DGAWj+9vfyMGBHByooUHoUgB4wFuMKpXKkJAQFxeX2bNnK5VPvUOlrq7O29v78OHDfI0LfU9TE6nVpFZj7njoH/iJUa1WGxoa6u7unpGRERwcPGPGDI1G88Q1V6xYUVJSIqrbA8DYxo2jiAiKiKBBg4QuBYAH/MRoenp6cXHxtm3b3N3do6Ki6uvrk5OTH18tISFBrVY///zzvAwKfdVLL1FsLMXGkpmZ0KUA8ICfGM3Pzw8KCpJKf+7Nx8dHP11UZ0VFRTExMYcOHZJIJLwMCgAgBvzEaGFhYcfcJkTk4eFRV1fXeQWtVrto0aIdO3Y4OTnxMiIAgEjwc9+oXC7v/EyVRqPx9vbuvMLWrVvt7OxGjBhx69YtjUZTXl5eUlLi7u7epR+lUunp6UlEFhYW48aN0y+Mi4vDRKEAwF1ERERtbS0R3b59W61WE5GFhQX3bvmJUU9Pz8TExI6PRUVFAQEBnVeor6+vqKgIDw8nooKCgv3791dVVe3atatLPz4+PpmYIgIAesfu3bu7LAkLC+PeLT8H9XPnzi0tLU1NTWUY5uDBgxUVFfoYTUpKUqlURPTxxx9n/UtQUFBMTMzjGQoA0BfxszeqUCji4+PnzZtnb2/f0NCQlJRkaWlJRKtWrYqJielygA8A0J/w9kz9ggUL5syZo1KpfH19TU1/7vb+/fuPr3n27Fm+BgUAEByfryaRy+X+/v48dggAIH54ph44+MtfaMQImjJF6DoAhIQX5QEHlZV09y7pdELXASAkxChwEBxMOh0pFELXASAkxChwsHQpLV0qdBEAAsO5UQAATsS1N9rS0pKRkUG/fBgUAIBfHQ+DPu2VnqyILkazs7OJyN7eHjEKAL2kuLi4uLiYiJqamrj3Jq4YtbGxeeutt4SuAgD6uVmzZul/0B/+coRzowAAnCBGB7r61nqhSwDo2xCjA5pWp53xxYySRyVCFwLQhyFGB7S//fC3zPLM98+8L3QhAH0YYnTgqmmqiUmPsTGz+TLny+o502nMGIqIELoogL4HMTpwbTi7ob61/vT/O+022K06+yrl5lIJju4BWEOMDlC5Vbn7r++PeD4iaFjQthnbvnHVFE8PpOBgoesC6HsQowNUZGqkjZnNB1M+IKKFYxZ+t2TKhFdLH0Xipl0A1hCjA9Hf7/w9rSjtzfFvPmp+VPSwqOhhUcTzEZWNlVsvbhW6NIC+R8IwjNA1/J9XXnlFoVAQkYuLy0cffSR0Of1TS3vL6D2jix4WPf6VmYnZrRW3vOy9jF8VgDFt3749KyuLiEpLS8+fP8+xN3E9DKpQKBISEoSuop972PTwoxlP/Seqtb3VmMUACGLt2rX6H3iZYFlcMQpG4GTltGD0AqGrAOg/cG4UAIATxCgAACeIUQAAThCjAACcIEb7nQkT6LnnaM8eoesAGChwpb7fycqitjYqLxe6DoCBAjHa78yZQ21t5OsrdB0AAwVitN85dkzoCgAGFnHFaENDw5EjR4jI1tb2t7/9rdDlAED/dP78+bKyMiKqra3l3pu4YlQqldrb2xORtbW10LUAQL9lZWWljxoTExPuvYkrRi0sLDomPgUA6CWBgYH6H/SHvxzhhicAAE4QowAAnCBGAQA4QYwCAHDCW4wqlcqQkBAXF5fZs2crlcrHV7hx48bUqVNtbW29vb0/++wzvsYFABAWPzGq1WpDQ0Pd3d0zMjKCg4NnzJih0Wg6r9Da2jp9+vRZs2bdunUrKioqIiKioKCAl6H7rdJSmjaNpk2js2eFLgUAfg0/MZqenl5cXLxt2zZ3d/eoqKj6+vrk5OTOK2RlZclksnfffdfFxeWNN9544YUXjh49ysvQ/ZZGQ+nplJ5ODx4IXQoA/Bp+YjQ/Pz8oKEgq/bk3Hx+fvLy8zit4eXmdOXNG/3NjY2Nubu7EiRN5GbrfsrCgKVNoyhQaMkToUgDg1/Bz+31hYaGdnV3HRw8Pj7q6us4rODg4ODg4EFFubm54ePj48eNfeuklXobut1xdifOEhQBgBPzEqFwur6mp6fio0Wi8vb27rMMwzAcffPDZZ59t2LDh7bfflkgkj/ejVCo9PT2JyMLCYty4cfqFcXFxMpmMlzoBYCCLiIjQP0R/+/ZttVpNRBYWFty75SdGPT09ExMTOz4WFRUFBAR0XkGn0y1fvvzevXt5eXlOTk5P68fHxyczM5OXkgAAuti9e3eXJbxMsMzPudG5c+eWlpampqYyDHPw4MGKigp9jCYlJalUKiK6dOnSN998s3//frlcXltbW1tb29TUxMvQAADC4idGFQpFfHz8vHnzXF1d165de/ToUUtLSyJatWrV2bNniejixYuVlZXu7u52//LBBx/wMjQAgLAkDMPw1Vdzc7NKpfL19TU1NfBcQVhYWEJCAl/1iERNU42NmY1MitO7AKLDS+bw+TCoXC739/c3OEP7q/dOv7cnExPMAfRbiLzedaPixqGsQ9Zm1mEBYY4WjkKXAwD8w6tJetc7qe/Yym0bWxs3nNsgdC0A0CsQo73o69yvz5ecj1MsuXli+O8j//v2j6lCVwQA/EOM9pYmbdO/n/73sU5jf2v1nO8Pxa+qmE9SNwpdFADwDzHaWz7O+Li4tvjTVz+VOjhSUFDlaPf0qsxjtzH7MUB/w+cNT9y99tprkyZNIiKFQrF06VKhyzFceX35M7HPeDt4fzzzY/2SdqY9/Hi41SCrvJV5clO5sOUBDHCJiYn6J4PS09NPnTrFsTdxXam3tLScO3cuEcnlfTto3j39bmNb442KGzPjZ3ZeXqWp2nFlx/uT3xeqMAAgookTJ/r5+RFRVlYW997EFaMymWz06NFCV8GVjtG9MfaNxc8ufuK3loMsjVwPAHTh6uqq/8HMzIx7b+KK0f5BKpFO95wudBUAYCS4xAQAwAliFACAE8QoAAAniFEAAE4QowAAnCBG2Th4kCZNouBgoesAABHBDU9s3LtHV64IXQQAiIu4YpRhGI1GQ0RSqVSMDzJ5edHMmfSkOU0BoA9paWlpb28nIp1Ox703cT1T/+qrr7q5uRHRsGHDYmJihC4HAPqn3bt3Z2dnE5FKpTp37hzH3sS1N+ro6Lhv3z6hqwCAfi4iIkL/g4gmWAYAGLAQowAAnCBGAQA4QYwCAHCCGAUA4AQxCgDACWIUAICTgRqj8+eTVEpjxwpdBwD0eQM1RnU6Yhji4zkwABjgxPUUU3V19bJly8gID4O+/jqNGUNOTr04BACIVcfDoOXl5dx7E1eMOjg4fPrpp0QklfbybnJISO/2DwAitnTpUv2rSZYsWcK9N3HFqEQisbCwELoKAOjnOuZV5mWPbaCeGwUA4AliFACAE8QoAAAniFEAAE6MFKNKpTIkJMTFxWX27NlKpdI4gwIAGIExYlSr1YaGhrq7u2dkZAQHB8+YMUM/4RIAQD9gjBhNT08vLi7etm2bu7t7VFRUfX19cnIyb73Hx1NsLJ06xVuHAABsGOO+0fz8/KCgoI77s3x8fPLy8njrffNmKiigRYvo5Zd56xMAoMeMEaOFhYV2dnYdHz08POrq6p64pkajOXHiBBGZmZm5uLjoF/r6+kowpzEAcKZUKrVaLRFVVFQ0NjYS0dOyiBVjxKhcLq+pqen4qNFovL29n7hmTU1NbGwsEVlaWqAThF4AAAfoSURBVAYFBekX+vj4mJo+vc4ffqD2dvrXMwkAAE+Tmpra0NBARLdu3aquriai2tpa7t0aI0Y9PT0TExM7PhYVFQUEBDxxTVdX14SEBHa929pyqQ0ABo7Vq1d3WdJnJlieO3duaWlpamoqwzAHDx6sqKh4WowCAPQ5xtgbVSgU8fHx8+bNs7e3b2hoSEpKsrS0NMK4AABGYKTb7xcsWKBWq1NSUqqrq2fMmMFjz8nJySqViscOhXX+/PmsrCyhq+DN9evXL168KHQVvMnPz09JSRG6Ct6Ul5d3PtvW1z169CguLk6QoY33ojy5XO7v7897t4WFhQ4ODrx3K5TS0lJHR0ehq+BNRUUFL6fwRUKtVhcVFQldBW/q6+vz8/OFroI3LS0tOTk5ggyNZ+oBADgRV4yWlZWxbXLv3j39jQs9V1NTk5GRwXagf/7zn2ybnDhxgmEYVk3KysoqKipYNWlsbDx9+jSrJmTQ5qSlpbF9ivf+/ftsf6c6ne7bb79l1YQM2pxLly6p1WpWTaqqqu7du8d2IANqM6DJtWvX2M6HUVtbW1JSwnYgw/4QdCznPbt161ZhYSGrJk1NTQb8dgzInMeJK0YN+KUWFxc/ePCAVZPKysr09HS2AxlwFunYsWP6iQp6rqSkhO3vtaGhwYATdgZsTnJysv6O5Z4rKytj+3+2Vqs9fvw4qyZk0OacO3euqqqKVZMHDx4UFxezHciA2gxocvXq1dLSUlZN1Gr13bt32Q5kQG3Hjx9n+4dw8+ZNtu8wamxsNGBzDMicx4krRgEA+hwJ26POXhUcHOzn58eqybVr14YOHerq6trzJo2NjRUVFSNHjmQ1UHZ2NtvbXbOzs/39/Vk9yZqdnW1hYeHl5dXzJq2trUVFRaNGjWJbG9vNuXPnzogRIwYNGtTzJkqlsrW1dcyYMT1votPpcnNz2V6NNGBzCgsLnZycWN17V1paWllZGRgY2Nu1GdCkpKRk8ODBtmyeRqmqqlKpVJMmTert2nJycsaMGcPqD+Gnn36SyWQKhaLnTRobGzMzM6dNm8aqtqysrB9++IFVk8eJK0bb29tbWlqErgIABgqZTCaTyTh2Iq4YBQDoc3BuFACAE8QoAAAniFEAAE4QowAAnCBGAQA46cMx2tcnbe62/hs3bkydOtXW1tbb2/uzzz4zfoWs9PDXUVdX5+3tffjwYSOWZohuN+fhw4fLli1zdnZ2c3PbsWOH8StkpdvNyczMXLx48bBhwxYuXGjAo9JisGbNms2bNwsydF+N0b4+aXO39be2tk6fPn3WrFm3bt2KioqKiIgoKCgQqtpu9fzXsWLFipKSEpHfZteTzYmMjNS/Umjnzp3R0dFsH8Q0pm43R6fTLVy40M/P79KlSz4+PgsXLmxraxOqWgNkZGRER0fHxsayfViZN0zflJaW5ujo2N7erv9oa2ubmJgobEmsdFv/1atXhw4d2vExODh4y5YtRi2RjR7+OuLj42fNmhUcHHzo0CHjFshOt5uj/yetrKxM//G7774rKCgwdpU91u3m5OTkyGQyjUbDMExdXZ1UKs3LyxOgUENFR0eHhIS4urr+5S9/EaSAvro32ruTNve+buv38vI6c+aM/ufGxsbc3NyJEycau8oe68mvo6ioKCYm5tChQ+Kf57Xbzblz587IkSMvX778xz/+8cMPPxw9ejTbZ4uNqdvN8fPz8/LyWrlyZXJy8vLly/39/X19fYWo1ECbNm368ssvBfwD6asx2vNJm8Wp2/odHBz0rxfIzc2dPHny+PHjX3rpJWNX2WPdbo5Wq120aNGOHTucnJyMXh1r3W5OQUFBVVXVnj17Zs6cefPmTV9fX6FeGNwT3W6ORCIJDw8/fPjw6tWrv/zyyxdffNHYJfZxfTVGez5pszj1pH6GYaKioiZPnhweHn7y5Ekx78R1uzlbt261s7MbMWLErVu3NBpNeXk5Ly8o6yXdbk5zc3NTU1NiYuKSJUuSkpLGjh371VdfGb3Mnup2c06cOLFp06b8/PyCggKlUvnZZ5/10atMQumrMerp6dn5ra6/MmmzOHVbv06nW7Zs2bVr1/Ly8iIjIzuOyMSp282pr6+vqKgIDw8PDw9XKpX79+8X89XtbjfHzc3N0dGxY7oXV1dXc3Nzo5bIRrebc/ny5UmTJumz1dvbe8qUKf1pBi1jEOSMLHeVlZVmZmYpKSk6ne7AgQP6OUeFLoqFp9WfmJioVCoZhjl//vyQIUNKSkoe/ov+CoA4dbs5nU2bNk3kl5i63Zzq6mobG5sjR44wDJOfn29ra3v37l1ha/4V3W5OfHy8ubn57du3GYYpLS21tLTMz88XuGj2FixYINQlpr4aowzDfP311+bm5sOHDx88eHBaWprQ5bD2xPqdnJz27t3LMMyWLVu6/IO3Zs0aQevtxq9vTmfij1GmB5tz/PhxGxsbV1dXExOTDRs2CFdpj/z65uh0usjIyCFDhnh5eZmYmKxevVrQYg0kYIz27RflNTc3q1QqX19fU1PjTXHKo75efxcDbXOam5uVSqWbmxurlyULpdvN0el0N2/edHZ27hOXAUWlb8coAIDgRH3hAgBA/BCjAACcIEYBADhBjAIAcIIYBQDgBDEKAMAJYhQAgBPEKAAAJ4hRAABOEKMAAJwgRgEAOEGMAgBwghgFAOAEMQoAwAliFACAE8QoAAAn/wvYgTU6fYaI4gAAAABJRU5ErkJggg==",
      "text/plain": [
       "Winston.FramedPlot(...)"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(Winston.plot(t, angle(g),\"g^\",tnew,angle(itc[tnew]),\"r.\"))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Julia 0.4.2",
   "language": "julia",
   "name": "julia-0.4"
  },
  "language_info": {
   "file_extension": ".jl",
   "mimetype": "application/julia",
   "name": "julia",
   "version": "0.4.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}