Skip to content

Instantly share code, notes, and snippets.

@mkatsura
Last active December 20, 2015 01:49
Show Gist options
  • Save mkatsura/6052119 to your computer and use it in GitHub Desktop.
Save mkatsura/6052119 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"metadata": {
"name": "SympyIntro"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# \u8907\u96d1\u306a\u6570\u5f0f\u306e\u53d6\u308a\u6271\u3044\uff08sympy\uff09 #\n",
"\u6842\u8aa0\n",
"\n",
"[INDEX](https://sites.google.com/site/makotokatsura4ipythonnotebook/ipythonnb4jpnexp)\n",
"\n",
"## Sympy\u30e2\u30b8\u30e5\u30fc\u30eb\u306b\u3064\u3044\u3066 ##\n",
"[sympy](http://sympy.org/en/index.html)\n",
"\u3068\u3044\u3046\u30e2\u30b8\u30e5\u30fc\u30eb\u3092\u7528\u3044\u308b\u3068\u3001\u6570\u5f0f\u3092TeX\u5f62\u5f0f\u3067\u7dba\u9e97\u306b\u51fa\u529b\u3057\u305f\u308a\u3001\u4ee3\u6570\u8a08\u7b97\u3092\u884c\u3046\u3053\u3068\u3082\u3067\u304d\u307e\u3059\u3002\n",
"\u3053\u306e\u30e2\u30b8\u30e5\u30fc\u30eb\u306f\u306a\u305c\u304b\u3001\u79c1\u306ePC\u306e\u4e00\u3064\u306b\u306f\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb\u3055\u308c\u3066\u3044\u306a\u304b\u3063\u305f\u3088\u3046\u3067\u3001\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb\u3059\u308b\u5fc5\u8981\u304c\u3042\u308a\u307e\u3057\u305f\u3002\n",
"pythonxy\u306b\u306f\u5165\u3063\u3066\u3044\u307e\u3059\u306e\u3067\u3001\u30a4\u30f3\u30b9\u30c8\u30fc\u30eb\u3057\u5fd8\u308c\u3066\u3044\u305f\u306e\u304b\u3082\u3057\u308c\u307e\u305b\u3093\u3002\n",
"\n",
"\u6b63\u5f0f\u3067\u306a\u3044[64\u30d3\u30c3\u30c8\u7248](http://www.lfd.uci.edu/~gohlke/pythonlibs/#sympy)\n",
"\u306f\u898b\u3064\u304b\u308a\u307e\u3059\u304c\u3001\u79c1\u81ea\u8eab\u306f32\u30d3\u30c3\u30c8\u306eWindows PC\u3092\u5229\u7528\u3057\u3066\u3044\u308b\u306e\u3067\u3001\u3053\u3061\u3089\u306e\u52d5\u4f5c\u78ba\u8a8d\u306f\u3057\u3066\u3044\u307e\u305b\u3093\u3002\n",
"[winpython](http://code.google.com/p/winpython/)\u306e64\u30d3\u30c3\u30c8\u7248\u306b\u3082\u3053\u308c\u304c\u5165\u3063\u3066\u3044\u308b\u306e\u3067\u306f\u306a\u3044\u304b\u3068\u63a8\u5bdf\u3057\u307e\u3059\u3002\n",
"\n",
"## \u6570\u5f0f\u306e\u8868\u793a\u65b9\u6cd5 ##"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"%load_ext sympy.interactive.ipythonprinting\n"
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 1
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"\u3053\u306e\u30b3\u30de\u30f3\u30c9\u3067sympy\u3068\u3044\u3046\u30e2\u30b8\u30e5\u30fc\u30eb\u3092\u7528\u3044\u3066\u66f8\u304b\u308c\u305f\u6570\u5f0f\u306e\u753b\u9762\u51fa\u529b\u304c\u5f97\u3089\u308c\u307e\u3059\u3002"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import sympy\n",
"from sympy import var,Eq\n",
"var('a:z')#\u666e\u901a\u306e\u5c0f\u6587\u5b57\u306ea\u304b\u3089z\u306f\u5168\u3066sympy\u306e\u30b7\u30f3\u30dc\u30eb\uff08\u5909\u6570\uff09\u3068\u3057\u3066\u6271\u3046\u3053\u3068\u306b\u3057\u307e\u3057\u305f\u3002\n",
"eq1=Eq(a*x**2+b*x+c,0)\n",
"eq1\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"latex": [
"$$a x^{2} + b x + c = 0$$"
],
"output_type": "pyout",
"png": "iVBORw0KGgoAAAANSUhEUgAAAIYAAAAZCAYAAAD0SVYHAAAABHNCSVQICAgIfAhkiAAABENJREFU\naIHt2WuoVFUUwPGfmlaaaab0sFKvdk0kTYoUxQfWh5IokCyCsMAoIqLIIgqhF0WYFWFpENiUX4Io\nyMeHsoge9LIoohdKaWWCIRWGVpbWh3WmOTMdc05zz8yVzh8ud/ba+5y91p6111p7DyUlGfQp+P1T\nMR1DMAP34PWC5yzp5RyF+1PtS7AHIzujTklvYRL2Y2zSPhp/Cgcp+R/TR6SSarqaKBxjSsc0KumV\nrMaDbZxvNt4Rzvh0G+c9FBiP5/EwHhLrc1wnFFmEpYovdhsZhD9wdZvn7c0MwTZcnpLdjk8woJ2K\nXCAcA47A6BbeNQmH5Rh/rogYE1qYMy95dWw39+J79ToOw++4tiroW7ASs0WIWo/jcR5OaOF9N+Gk\nHONnYie+aGHOvOTVsd0swLsiklb5QazRgqqgSM/uwjpxbE0zpMA5G5mFN0XUKGEwTsWGjL7tmFZt\nZDnGXFyDrzECr+BGnJUacyYWYp9IDVclzwwV9xR34MtEkU7RX5yKVqkVvZOELZ+mxjVrS5GcjsXY\nhV/wG+7Drz08z6jk/66Mvt3iSuHwZP46FmGHWigcnSj3YmpMFx5TS0MVbBI3nDPE3cXiFpT/Nyqa\nr1Gmi0ixVjgJEeY/VyuCi7Alj47Evc52TE7as/AzLsoYuwof5fybk3q+uiZ3Zbx7ddL3j9PJZFGA\nXNYg34klqfYK9ZHgWZGz4GQsw7EZE/cEFc0v+q0idw5MyeYL4ycm7SJsyaPjFLE7F6ZkF+IDxdQp\nU4X9d2b0PZP0ndjYsU4sZP+UbEIyeHZKNqbhuW2i0m0HFc0v+nq80CC7QdhzdtIuwpaK5nXcgG/Q\nr8U5m6XLgR1jbdI3mFqNMVScGJ4TUaPKHOxV20WwJfV5vMjDr7aucx1PqYXWNKeIL3VvRt8isdOq\nTBN5WoNsHzYn7VZsaVXHETgHKxOd2sEO8eUfk9E3CD+JNPa3Y4wTXvt2w+A52CjqjC581dA/VyzA\nWynZWK0Xa1ccQF4R3r71IM8PF2fzjSlZP3GvsQE/ZjyT15ZWdewStc77BxmX5gn5f1JYjNeSz7vx\noUiTjYwTNUkd3cKT5qdkRwoPW5q0H0lkS0UVTVyrphe/Lx7PqXgeKpoL0wNF4didkl0qdmY1jRRl\nS7M6jhFrfnFG3yic/x/mboa78Z36W+ixiS7XVQXVanwzPlYzqD8eFTeVW8UO3Il5uEUUb6eJM3H6\naLNEhNhOswcvq914jsRy3Iz3ElmnbdmCl8QpJM1McUR+o6B5V4q0kb4Svx6fiYiEeq/pFj+qbBJh\nd7k4sl0pHOe2ZNwDwkmqx54VItXsxRrZlyc9RUVzYZrI9cvE7hgnjF6T6h+uGFvy6DhUrPl+UfgP\nEBv0yURWFGeIInuTKDaHicL82wLnLJSK1n5raQcVvV/Hg1L0byU9zS49fxvY0xwKOpaUlJSUlJSU\ndIS/ABU//PrKyBZhAAAAAElFTkSuQmCC\n",
"prompt_number": 2,
"text": [
" 2 \n",
"a\u22c5x + b\u22c5x + c = 0"
]
}
],
"prompt_number": 2
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"\u4e0b\u8a18\u306e\u3088\u3046\u306b\u3001\u65b9\u7a0b\u5f0f\u3092\u89e3\u304f\u4e8b\u3082\u3067\u304d\u307e\u3059\u3002"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"eq2=eq1.subs({a:1,b:0,c:1})# \u5404\u4fc2\u6570\u306b\u5024\u3092\u4ee3\u5165\u3002\n",
"from sympy import solve\n",
"ans=solve(eq2)\n",
"ans\n",
"\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"latex": [
"$$\\begin{bmatrix}- \\mathbf{\\imath}, & \\mathbf{\\imath}\\end{bmatrix}$$"
],
"output_type": "pyout",
"prompt_number": 3,
"text": [
"[-\u2148, \u2148]"
]
}
],
"prompt_number": 3
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"\u3053\u306e\u3088\u3046\u306b\u3001i\u306e\u4e0a\u306e\u70b9\u304c\u3042\u308a\u307e\u305b\u3093\u3002\n",
"\u3053\u308c\u306f\u865a\u6570\u5358\u4f4d\u306e$i$\u306b\\imath\u3068\u3044\u3046TeX\u30b3\u30de\u30f3\u30c9\u304c\u4f7f\u7528\u3055\u308c\u3066\u3044\u308b\u306e\u304c\u539f\u56e0\u306e\u3088\u3046\u3067\u3059\u3002\n",
"\u3053\u308c\u306f\u79c1\u306b\u306f\u6c17\u6301\u3061\u60aa\u3044\u306e\u3067\u3001\u4ee5\u4e0b\u306e\u30b3\u30de\u30f3\u30c9\u3092\u5b9a\u7fa9\u3057\u3066\u307f\u307e\u3057\u305f\u3002"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"def print_eq(eq):\n",
" from IPython.display import display, Math\n",
" from sympy import latex\n",
" display(Math(latex(eq).replace('\\imath','i')))\n",
" \n",
" "
],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 4
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"print_eq(ans)\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"latex": [
"$$\\begin{bmatrix}- \\mathbf{i}, & \\mathbf{i}\\end{bmatrix}$$"
],
"output_type": "display_data",
"text": [
"<IPython.core.display.Math object at 0x0488BAF0>"
]
}
],
"prompt_number": 5
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"\u3053\u3063\u3061\u306a\u3089\u79c1\u306f\u8a31\u5bb9\u3067\u304d\u307e\u3059\u3002"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"print ans\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"[-I, I]\n"
]
}
],
"prompt_number": 6
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"\u5358\u306bprint\u3068\u3059\u308c\u3070\u3001Sympy\u3067\u306e\u30ad\u30e3\u30e9\u30af\u30bf\u304c\u51fa\u529b\u3055\u308c\u307e\u3059\u3002\u4ee5\u4e0b\u3082\u540c\u69d8\u3067\u3059\u3002"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"print eq1\n",
"print eq2\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"a*x**2 + b*x + c == 0\n",
"x**2 + 1 == 0\n"
]
}
],
"prompt_number": 7
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## \u8907\u96d1\u306a\u5f0f\u306e\u66f8\u304d\u65b9 ##\n",
"\u5c11\u3057\u3060\u3051\u8907\u96d1\u306a\u4f8b\u3092\u793a\u3057\u307e\u3059\u3002"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from sympy import atan, pi\n",
"y=2*x+1\n",
"exp0=atan(y**2/(2*y+10)**2)\n",
"exp1=sympy.integrals.Integral(exp0,(x,-10.0,pi))\n",
"print_eq(exp1)\n",
"print(exp1)\n",
"print(exp0)\n",
"\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"latex": [
"$$\\int_{-10.0}^{\\pi} \\operatorname{atan}{\\left (\\frac{\\left(2 x + 1\\right)^{2}}{\\left(4 x + 12\\right)^{2}} \\right )}\\, dx$$"
],
"output_type": "display_data",
"text": [
"<IPython.core.display.Math object at 0x059B5290>"
]
},
{
"output_type": "stream",
"stream": "stdout",
"text": [
"Integral(atan((2*x + 1)**2/(4*x + 12)**2), (x, -10.0, pi))\n",
"atan((2*x + 1)**2/(4*x + 12)**2)\n"
]
}
],
"prompt_number": 8
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"\u81ea\u52d5\u7684\u306by\u3092\u4ee3\u5165\u3057\u3066\u304f\u308c\u308b\u306e\u3067\u3001\u8907\u96d1\u306a\u5f0f\u3092\u66f8\u304f\u306b\u306f\u4fbf\u5229\u3067\u3059\u3002\n",
"\n",
"## \u4ed6\u306e\u30e2\u30b8\u30e5\u30fc\u30eb\u306e\u5229\u7528 ##\n",
"\n",
"\u4eca\u5ea6\u306f\u3001\u3053\u306e\u5f0f\u3092Numpy\u3067\u6271\u3063\u3066\u307f\u307e\u3057\u3087\u3046\u3002\n",
"\u624b\u3067\u5c11\u3057\u66f8\u304d\u63db\u3048\u307e\u3059\u3002\n",
"\u4e0a\u306eprint\u6587\u306e\u51fa\u529b\u3092\u30b3\u30d4\u30da\u3057\u3066\u7c21\u5358\u306a\u95a2\u6570\u3092\u5b9a\u7fa9\u3057\u307e\u3059\u3002"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"def exp0_func(x):\n",
" from numpy import arctan # numpy\u3067\u306farctan, sympy\u3067\u306fatan\n",
" return arctan((2*x + 1)**2/(4*x + 12)**2)\n",
"_x=np.arange(-10.0,np.pi,0.1)\n",
"plt.plot(_x,exp0_func(_x))\n",
"\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "pyout",
"prompt_number": 9,
"text": [
"[Line2D(_line0)]"
]
},
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAD9CAYAAAChtfywAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XtclFXiBvBnFNxMLc0LGjOKMCgoAq6UmlmTmbCWtJnb\nomUliuTmlt3WrrtQidp22ZRPG9pamWZmq+H2Q8zbaGqAt8S84oUY0EwSi7IcGM7vjxMoKsMMzMyZ\nd+b5fj7zkYl3Zp5t5elw3vO+RyeEECAiIs1qoToAERE1D4uciEjjWORERBrHIici0jgWORGRxrHI\niYg0zm6RJycnIygoCP369WvwGLPZjP79+yMqKgomk8nV+YiIqBE6e+vIv/jiC7Rt2xb3338/9uzZ\nc8n3z5w5gyFDhmD16tXQ6/UoLy9Hp06d3BqYiIjqszsiHzp0KDp06NDg9z/88EPcfffd0Ov1AMAS\nJyJSIKA5Ly4qKkJVVRVuueUWVFZW4tFHH8X48eMvOU6n0zXnY4iI/JYjF98362RnVVUVdu7ciZyc\nHKxevRovvfQSioqKGgyj1cc//vEP5Rn8Nb+WszO/+ofW8zuqWSNyg8GATp06oXXr1mjdujVuuukm\n7N69G+Hh4c15WyIickKzRuR33nknNm/eDJvNhrNnzyI/Px99+vRxVTYiInKA3RH52LFjsXHjRpSX\nl8NgMCA9PR1VVVUAgNTUVERERCAhIQHR0dFo0aIFUlJSfLLItb6sUsv53ZVdCCA/H/joI2D3biAp\nCRg3DmjXzrWfo+V/9wDza4Xd5Ycu+xCdzqn5HiJ3+uknYPRo4JtvZHlHRwMffACYzcCyZcCtt6pO\nSCQ52p0scvIrp08DI0cCUVFAVhbQsuX5761dCzzwALBzJxAUpC4jUS1Hu5OX6JPf+PlnYNgw4MYb\ngfnz65c4AAwfDkyYANx/P1BToyYjUVNwRE5+IzUVOHsWWLgQaOjShupqwGSSUy+PP+7ReESX4NQK\n0QVWrACefBLYtQu46ir7x379NZCQABQXAwHNWqBL1DycWiH6zfHjwEMPAYsXN17igJw/1+uB1avd\nn43IFVjk5POeegqYNAkYNMjx10yaBLzzjvsyEbkSp1bIp+XlAWPGAAcOAG3bOv66ykqge3dg/36g\na1f35SOyh1Mr5PeEAKZNAzIynCtxQF4YNHq0PDFK5O1Y5OSzliwBbDbgvvua9vra6RX+MknejkVO\nPqmqCnjuOeD114EWTfxbPmgQcO4ccPCga7MRuRqLnHzSBx8AYWHA0KFNfw+dDrj5ZmDjRtflInIH\nFjn5nOpqYMYM4B//aP573XwzsGlT89+HyJ1Y5ORzFi+WK06aMxqvddNNckTOeXLyZixy8imuHI0D\ngNEo77ty9Khr3o/IHVjk5FNWrAA6d5ZTIq7AeXLSAhY5+QwhgFdflfdUceV+3yxy8nYscvIZW7YA\n338PJCa69n15wpO8HYucfMZrrwGPPXbpfcabKyJC3su8pMS170vkKixy8glFRcDmzcCDD7r+vXU6\nuXqFo3LyVnaLPDk5GUFBQejXr5/dN9m2bRsCAgKwfPlyl4YjctSbbwKTJwNt2rjn/QcMAL76yj3v\nTdRcdot8woQJyM3NtfsGNpsN06dPR0JCAu9wSEr8+CPw4YfAlCnu+4yoKGDvXve9P1Fz2C3yoUOH\nokOHDnbfYO7cuRgzZgw6d+7s0mBEjnr/fbnfpl7vvs/o21fuHETkjZq1kVVZWRmys7Oxfv16bNu2\nDTo7a77S0tLqvjaZTDCZTM35aCIA8mKdzEz3bwIREgJUVAA//ABcfbV7P4v8l9lshtlsdvp1zSry\nadOmYdasWXU3P7c3tXJhkRO5ypo1QOvWwI03uvdzWrQAIiPl9MoNN7j3s8h/XTzITU9Pd+h1zSry\nHTt2ICkpCQBQXl6OVatWITAwEImuXshL1IDMTOCvf3XtBUANqZ0nZ5GTt2lWkR+94AYUEyZMwKhR\no1ji5DHFxcCXXwJLl3rm8/r25QlP8k52i3zs2LHYuHEjysvLYTAYkJ6ejqqqKgBAamqqRwISNWT+\nfGD8eODKKz3zeVFRQCOLuIiU4ObLpElWq7xVrdksr7z0BIsFuP564MQJz3weETdfJp+WnS1PPnqq\nxAG5vPHsWXk/FyJvwiInTXr7beChhzz7mTod58nJO7HISXMOHpQX59x1l+c/OyqKFwaR92GRk+bM\nmwckJwOtWnn+szkiJ2/UrOWHRJ72yy/AwoVAfr6az4+KAnhvOPI2HJGTpixbBsTFAaGhaj4/PBw4\nfFjNZxM1hEVOmqLiJOeFgoOB8nLg3Dl1GYguxiInzdi9W+7Sc/vt6jK0bCmXIX7zjboMRBdjkZNm\nZGUBKSlAgOIzOz17ytsDEHkLnuwkTaisBD76CNizR3USWeTHjqlOQXQeR+SkCUuWACaTnKNWLSSE\nRU7ehUVOXk8I9Sc5L8SpFfI2LHLyetu2yZ15hg9XnUTiiJy8DYucvN7bbwOpqXKXHm/AOXLyNryN\nLXm1igpZnIcOAV26qE4jCSHvgV5eDrRpozoN+TLexpZ8wnvvyXXj3lLigLwLYo8enCcn78EiJ69V\nUwP8+9/AX/6iOsmleMKTvAmLnLzWunVA69beudkx58nJm7DIyWu99Rbw8MNyKsPbcOUKeZNGizw5\nORlBQUHo16/fZb+/ePFixMTEIDo6GkOGDEFhYaHLQ5L/sViATZuAceNUJ7k8Tq2QN2m0yCdMmIBc\nO1uHh4aGYtOmTSgsLMQLL7yAyZMnuzQg+aesLODee4G2bVUnuTyOyMmbNHqvlaFDh6LYztBj8ODB\ndV8PHDgQpaWllz0uLS2t7muTyQSTyeRwSPIvVivwzjvAhg2qkzSMc+TkDmazGWaz2enXObSOvLi4\nGKNGjcKeRu5Y9Oqrr+LQoUOYN29e/Q/hOnJywkcfye3c1q9XnaRhQgBXXSWngNq3V52GfJWj3emy\nux9u2LABCxYswJYtW1z1luSn3noLeOQR1Sns0+nk9EpxMRAbqzoN+TuXrFopLCxESkoKVq5ciQ4d\nOrjiLclP7dkDHDkC3Hmn6iSN0+uBsjLVKYhcUOQlJSUYPXo0Fi1aBKPR6IpM5MfeeguYPBkIDFSd\npHHBwcDx46pTEDkwtTJ27Fhs3LgR5eXlMBgMSE9PR1VVFQAgNTUVL774IioqKjBlyhQAQGBgIAoK\nCtybmnxSRYWcH9+7V3USx1x7LUfk5B140yzyGq+8AhQWAosWqU7imKwsYPt2YP581UnIV/GmWaQp\nVVXA3LnAY4+pTuI4Tq2Qt2CRk1f473/l2uwBA1QncRynVshbsMhJOSGAN97Q1mgckCNyFjl5AxY5\nKffll3KThsRE1Umc07mz3ILu3DnVScjfschJuTfekBcAtWypOolzWrQAunYFTpxQnYT8HYuclCou\nlpfiJyerTtI0nF4hb8AiJ6XmzgUefBBo1051kqbhyhXyBi671wqRsyor5Z6cO3aoTtJ0HJGTN+CI\nnJRZsAAYNkzefEqruASRvAFH5KSE1Qq8/jrw8ceqkzRPcLC8GpVIJY7ISYlFi4DwcGDgQNVJmodT\nK+QNOCInj7PZgFmz5L1KtI5TK+QNOCInj1u2DOjUCfCF3f5qV63wnnCkEoucPKqmBsjIAJ5/Xu6y\no3Vt2wIBAfIKTyJVWOTkUZ99JovvD39QncR1OL1CqrHIyWOEAF5+GXj2Wd8YjdfiRUGkGoucPGbt\nWuCnn4DRo1UncS2uXCHVWOTkMTNmAM88I2825Us4tUKq2f2RSk5ORlBQEPr169fgMY888gjCw8MR\nExODXbt2uTwg+YYtW4CSEmDsWNVJXI8jclLNbpFPmDABubm5DX4/JycHhw8fRlFREebNm1e3ATPR\nxV54Qc6NB/jglQvdugHffqs6Bfkzu0U+dOhQdOjQocHvr1y5Eg888AAAYODAgThz5gxOnjzp2oSk\neevWARYL8NtfFZ/TtSuLnNRq1viorKwMBoOh7rler0dpaSmCgoIuOTYtLa3ua5PJBJMvXA1CjRIC\neO45ID0dCAxUncY9unXj5hLkGmazGWaz2enXNfsXXXHRJW26BtaVXVjk5D8++wz4+WcgKUl1Evep\nHZEL4VvLKsnzLh7kpqenO/S6Zq0fCA4OhsViqXteWlqK4ODg5rwl+RCbTV7B+dJLvrdS5UJXXgm0\nasWrO0mdZv14JSYmYuHChQCAvLw8tG/f/rLTKuSfPvhAXsJ+552qk7gf58lJJbtTK2PHjsXGjRtR\nXl4Og8GA9PR0VFVVAQBSU1MxcuRI5OTkwGg0ok2bNnj33Xc9Epq83y+/yJUqS5f6x3RD7Tx5RITq\nJOSPdOLiSW53fIhOd8lcOvm2mTPlFm6ffKI6iWckJcnfPHxxnTyp42h3+uCqXlLt1CngtdeAL79U\nncRzuHKFVPLhU1CkygsvAPfeK3cA8hddu7LISR2OyMmlvvoKWLECOHBAdRLP6tYN+Ppr1SnIX3FE\nTi4jBPDoo8CLLwJ2Lgj2SRyRk0oscnKZZcvkWupJk1Qn8Tzeb4VU4tQKuURlJfDEE8DixUDLlqrT\neB5H5KQSlx+SSzz+OFBRAfjrpQQ1NcAVV8j/oP3ud6rTkK/g8kPymF275Ejcn0/2tWgBdOkCnDwJ\ndO+uOg35G86RU7PYbMCUKXL3n86dVadRi/PkpAqLnJpl7lx5w6jkZNVJ1OM8OanCqRVqsqNHgZdf\nlldw+vLdDR3FETmpwh8/ahIhgJQU4Omn/esKTns4IidVWOTUJPPmyRUa06apTuI9OCInVVjk5LTD\nh+WGEe+/75ubKTcVR+SkCoucnFJdDYwfL2+MFRmpOo134YicVGGRk1Nmz5a7/kydqjqJ9+GInFTh\nlZ3ksK1bgdGjge3bAb1edRrv8+uvwNVXyz/9YVckcj9Hu5MjcnLI6dPAuHHA/Pks8YZccYXciPn0\nadVJyN+wyKlRQgATJwJ//CMwapTqNN6NOwWRCo0WeW5uLiIiIhAeHo7Zs2df8v3y8nIkJCQgNjYW\nUVFReO+999yRkxT6178Ai0XOj5N9XbvyhCd5nt0it9lsmDp1KnJzc7Fv3z4sWbIE+/fvr3dMZmYm\n+vfvj6+++gpmsxlPPPEEqqur3RqaPGfTJlngn3zCu/o5giNyUsFukRcUFMBoNCIkJASBgYFISkpC\ndnZ2vWO6deuGH3/8EQDw448/omPHjgjg4mKfcPy43B3+/feBkBDVabSBI3JSwW7jlpWVwWAw1D3X\n6/XIz8+vd0xKSgqGDRuGa6+9FpWVlfj4448v+15paWl1X5tMJphMpqanJrc7e1bOiT/8MBAfrzqN\ndnTrJv8DSNQUZrMZZrPZ6dfZLXKdA2uoMjIyEBsbC7PZjCNHjuC2227D7t270a5du3rHXVjk5N1q\naoD77wd69waefVZ1Gm3p2hXYuVN1CtKqiwe56enpDr3O7tRKcHAwLBZL3XOLxQL9RWvPtm7dij/9\n6U8AgLCwMPTs2RMHDx50NDd5oWeeAb77DnjnHa6HdhbnyEkFu0UeFxeHoqIiFBcXw2q1YunSpUhM\nTKx3TEREBNauXQsAOHnyJA4ePIjQ0FD3JSa3mj8fWLFCPnhy03mcIycV7E6tBAQEIDMzE/Hx8bDZ\nbJg4cSIiIyORlZUFAEhNTcWzzz6LCRMmICYmBjU1NXjllVdwzTXXeCQ8udaaNfIeKl98AXTsqDqN\nNnFETirwEn0CAOzZA9x6K/Df/wJDh6pOo11CyCs8z5wBWrdWnYa0jpfok8MOHpQrU+bOZYk3l07H\n6RXyPBa5nztyBBg+HJg5E/jzn1Wn8Q0scvI0FrkfKymR0ynPPQc88IDqNL6D8+TkaSxyP3X8ODBs\nmNyq7aGHVKfxLRyRk6exyP3QiRNyOmXiRO656Q4ckZOnscj9zJEjwI03AvfeKy/8Idfjlm/kaSxy\nP7J7N3DTTcBTT8l5cXIPbvlGnsbbFPqJzZvlNm2ZmcA996hO49s4IidPY5H7gZwcuSpl8WJgxAjV\naXwfR+Tkabyy04cJIS/yycgAPv0UGDRIdSL/YLUCbdrITZhbtlSdhrTM0e7kiNxHWa3yXuJ5ecDW\nrQDvY+Y5rVoBV18NfP890KWL6jTkD3iy0wedOiWXF373HUtcFS5BJE9ikfuYwkLg+uvlPVNWrAAu\n2t+DPIQXBZEnsch9hBDyXuK33irnxGfMAFrw/11luOUbeRLnyH3ADz8AqanAvn1y1/vISNWJSK8H\nyspUpyB/wTGbxm3bBvz+98A11wD5+Sxxb6HXA6WlqlOQv2CRa1R1NTB7NnD77fLPt97iRgbehEVO\nnsSpFQ3auxeYMEGeyCwoAEJCVCeii7HIyZM4IteQ6mq5AYTJBEyaBKxdyxL3Vixy8qRGizw3NxcR\nEREIDw/H7NmzL3uM2WxG//79ERUVBZPJ5OqMBGDXLnll5oYNwPbtwOTJclsx8k6dOgGVlcAvv6hO\nQv7A7iX6NpsNvXv3xtq1axEcHIzrrrsOS5YsQeQFZ9TOnDmDIUOGYPXq1dDr9SgvL0enTp3qfwgv\n0W+yigrg+eflpsgZGXJKhQWuDaGhwOefA0aj6iSkVS7ZfLmgoABGoxEhISEIDAxEUlISsrOz6x3z\n4Ycf4u6774ZerweAS0qcmqamBvjPf+QqFCHk0sLkZJa4lnB6hTzF7snOsrIyGAyGuud6vR75+fn1\njikqKkJVVRVuueUWVFZW4tFHH8X48eMvea+0tLS6r00mE6dg7NiyBXjiCVnaOTlyeSFpD4ucnGU2\nm2E2m51+nd0i1zkw/KuqqsLOnTuxbt06nD17FoMHD8agQYMQHh5e77gLi5wub/9+4Omn5Xz4Sy8B\n48fz6kwtMxhY5OSciwe56enpDr3Obk0EBwfDYrHUPbdYLHVTKLUMBgNGjBiB1q1bo2PHjrjpppuw\ne/duJ6LT8eNASorcvWfoUODQIXn/cJa4tnFETp5ityri4uJQVFSE4uJiWK1WLF26FImJifWOufPO\nO7F582bYbDacPXsW+fn56NOnj1tD+4rSUuCRR4CoKHll5qFDwJNPAldcoToZuQKLnDzF7tRKQEAA\nMjMzER8fD5vNhokTJyIyMhJZWVkAgNTUVERERCAhIQHR0dFo0aIFUlJSWOSNKCkBZs0CPvpI7mS/\nb5+8Wx75Fr0euOAXWiK34Q5BHrR3L/D663K3nsmTgccfBzp3Vp2K3OXECSA2Fjh5UnUS0iruEOQl\nhADWrQNee02exHz4YTmF0rGj6mTkbl26yOsAzp0Dfvc71WnIl7HI3eTsWTl1MmcOUFUlR98rVnD+\n25+0bHn+vuQ9e6pOQ76MRe5iBw4Ab78NLFoEDB4s58Lj43khj7+qPeHJIid3YpG7wM8/y9H2ggXy\nxOXEicCOHUCPHqqTkWpcuUKewCJvopoaYONGYOFCefJy8GC5S89dd8ld1IkAFjl5BovcSYcOyfL+\n4AOgfXt54c7MmVw+SJdnMADHjqlOQb6ORe6AQ4fk3Qc/+UTuwzhuHLByJRATozoZeTu9HvjiC9Up\nyNexyC9DCLnmu7a8y8uB0aOBV1+Vl9AH8N8aOchgAL75RnUK8nW8IOg31dXA1q3yboOffiqXD959\nt3zccAPve0JNc+oUEBEBfP+96iSkRY52p18X+XffAatWyfJes0Zum3b77cCoUcB113HJIDWfEMBV\nV8lL9du3V52GtIZXdl7GuXNAXp680jI3V859Dx8OjBwJvPEGcO21qhOSr9Hp5E5BR4/yvvLkPj5d\n5NXVwM6dwPr18vHll3LHnWHDgNmzgSFDuFSQ3C8sjEVO7uVTRf7rr3Jj4q1bgc2b5WoBvV4W98MP\nAx9/zF9vyfNCQ4EjR1SnIF+m6SL/7jtZ2lu2yMfu3XLEPWQIcN99wPz5QFCQ6pTk70JD5d9NInfR\nTJGfOwcUFspL3/PyZHGfOiWvqBwyBHj5ZeD664G2bVUnJaovLEzewoHIXbyyyM+dA/bskaW9fbv8\n88ABoFcvYMAAYNAguTlx375cFkjej1Mr5G7Klx+ePi1Le8+e8yPu/fuB8HBZ2gMGAHFxQHQ00Lq1\nu5MSuZ7VCrRrB/z0ExAYqDoNaYnXrSP/9VeBAwfql/aePcCPP8o9K6OjgX79ZHHHxLC0ybeEhMiV\nU6GhqpOQlnjdOvL27eU9mfv1k6U9ZYr8s0cPXnhDvq92eoVFTu7Q6Axzbm4uIiIiEB4ejtmzZzd4\n3LZt2xAQEIDly5df9vsVFfJe3UuXAs89ByQmylEKS5z8Qe1aciJ3sFvkNpsNU6dORW5uLvbt24cl\nS5Zg//79lz1u+vTpSEhIaPDXAG5xRv6s9upOInewW+QFBQUwGo0ICQlBYGAgkpKSkJ2dfclxc+fO\nxZgxY9CZW8ITXRZXrpA72Z0jLysrg8FgqHuu1+uRn59/yTHZ2dlYv349tm3bBl0DcyVpaWl1X5tM\nJphMpqanJtIYTq2QI8xmM8xms9Ovs1vkDZXyhaZNm4ZZs2bVnV1taGrlwiIn8je1I3IheF6IGnbx\nIDc9Pd2h19kt8uDgYFgslrrnFosFer2+3jE7duxAUlISAKC8vByrVq1CYGAgEhMTHc1O5PM6dJAF\nfvo00LGj6jTka+wWeVxcHIqKilBcXIxrr70WS5cuxZIlS+odc/SC3xcnTJiAUaNGscSJLqLTySuT\nDx6UG5UQuZLdk50BAQHIzMxEfHw8+vTpgz//+c+IjIxEVlYWsrKyPJWRyCdERQFff606Bfki5Zfo\nE/mL118HiouBOXNUJyGtcLQ7ecspIg/hiJzchUVO5CH9+rHIyT1Y5EQe0rUrYLMBJ0+qTkK+hkVO\n5CE6HadXyD1Y5EQexCInd2CRE3kQi5zcgUVO5EH9+skNVYhcievIiTyoogLo3h344QfuN0uN4zpy\nIi/UoQNw9dVASYnqJORLWOREHsZ5cnI1FjmRh0VFcZ6cXItFTuRhAwYA27erTkG+hEVO5GE33ABs\n3So3mSByBRY5kYd17y5XrBQXq05CvoJFTuRhOt35UTmRK7DIiRRgkZMrsciJFGCRkyvxyk4iBaxW\n4JprgBMngHbtVKchb8UrO4m8WKtWQP/+QEGB6iTkCxot8tzcXERERCA8PByzZ8++5PuLFy9GTEwM\noqOjMWTIEBQWFrolKJGvueEGYMsW1SnIF9gtcpvNhqlTpyI3Nxf79u3DkiVLsH///nrHhIaGYtOm\nTSgsLMQLL7yAyZMnuzUwka/gPDm5it0iLygogNFoREhICAIDA5GUlITs7Ox6xwwePBhXX301AGDg\nwIEoLS11X1oiHzJ4MJCXB1RVqU5CWhdg75tlZWUwGAx1z/V6PfLz8xs8/j//+Q9Gjhx52e+lpaXV\nfW0ymWAymZxLSuRjunQBwsOBzZuBW25RnYa8gdlshtlsdvp1dotcp9M5/EYbNmzAggULsKWBSb8L\ni5yIpFGjgJUrWeQkXTzITU9Pd+h1dqdWgoODYbFY6p5bLBbo9fpLjissLERKSgpWrlyJDh06OBiZ\niBITgf/9j/ddoeaxW+RxcXEoKipCcXExrFYrli5disTExHrHlJSUYPTo0Vi0aBGMRqNbwxL5mpgY\nuab8wAHVSUjL7E6tBAQEIDMzE/Hx8bDZbJg4cSIiIyORlZUFAEhNTcWLL76IiooKTJkyBQAQGBiI\nAi6OJXKITienV/73PyAyUnUa0ipe2UmkWG4u8PLL8qQn0YUc7U4WOZFiv/4KBAUBR44AnTqpTkPe\nhJfoE2nEFVcAI0YAy5erTkJaxSIn8gITJwLz5qlOQVrFIifyAiNGAN9/z708qWlY5EReoEULICUF\n+G1BGJFTeLKTyEt8+61cglhcDPx2+yLyczzZSaQxXbsCw4cDixerTkJawyIn8iJ/+Qvwr3/xjojk\nHBY5kRcxmYDu3YEFC1QnIS3hHDmRl9m+Xd5Mq6gIaNNGdRpSiVd2EmlYUhIQFQU8/7zqJHQhqxWw\nWICSEuCbb+Sj9uuwMNevOmKRE2nY4cPAoEHA7t1AcLDqNP7DZgNKS4GjR+s/iotlYZeXA926AT16\nyEf37ue/NhqB0FDX5mGRE2nciy8CGzcCn38OtGypOo3vqKw8X9BHjtQv7JISeb+b0FD5CAsDevYE\nQkJkWXfrBgTYvWesa7HIiTTOZgOGDQP+8Afg6adVp9GWs2flOYaDB4FDh+SfRUWyuM+ePV/UtWVd\n+3VIiLz3jbdgkRP5AIsFiIsDPv1UbtZM59lscgRdW9QXlvapU7Kge/cGevU6/2dYmNwr1YldLJVi\nkRP5iJwcIDkZWL1a7ijkb86dk+W8d6987N8vC/vIEaBjx/plXfvo3t03pqNY5EQ+5JNPgL/+FViz\nRq5m8UVWqyzo2sKufXzzjZyn7ttXPiIjZVmHh/v+8kxHu9OD0/ZE1FRjxgDV1fIS/vffB+LjVSdq\nuqoqOV99cWEfOyZPKNYW9j33yD979QJatVKd2rtxRO4As9kMk8mkOkaTaTm/lrMDrs9vNgP33gs8\n+CCQlgYEBrrsrRv4vKbnr66W0x8XF/bhw4Bef76wax+9e7v+RKPW//647KZZubm5iIiIQHh4OGbP\nnn3ZYx555BGEh4cjJiYGu3btcj6tlzObzaojNIuW82s5O+D6/CYTsGuXfERFyV2F3DlGciR/dbWc\nw16+XO49OnYsEB0NtGsH3H478N57cqXIHXcACxcCp0/LEfmnnwIzZgDjxsm5f3esFtH63x9H2Z1a\nsdlsmDp1KtauXYvg4GBcd911SExMROQF233n5OTg8OHDKCoqQn5+PqZMmYK8vDy3ByfyV126AP/3\nf3J9+fTpsjxTUmSBtm/vvs/96Sc5wj58WJ5w3LsX2LdPlnK3bnJU3acPMHIk8NRTQEQEcOWV7stD\n59kt8oKCAhiNRoSEhAAAkpKSkJ2dXa/IV65ciQceeAAAMHDgQJw5cwYnT55EUFCQ+1IT+TmdTs6T\n33abPAG6YIFcaz5gAHDLLcDAgbJUg4MdW2onhLxQ5sQJ+Th+XM5ZZ2cD69bJ8v7hh/PrriMj5fr2\nJ5+UX7Ns4sfUAAAF0ElEQVSwFRN2LFu2TEyaNKnu+QcffCCmTp1a75g77rhDbNmype75rbfeKrZv\n317vGAB88MEHH3w04eEIuyNynYOr5i+ejL/4dVo+0UlE5O3snuwMDg6GxWKpe26xWKDX6+0eU1pa\nimDe5YeIyGPsFnlcXByKiopQXFwMq9WKpUuXIjExsd4xiYmJWLhwIQAgLy8P7du35/w4EZEH2Z1a\nCQgIQGZmJuLj42Gz2TBx4kRERkYi67eb7qampmLkyJHIycmB0WhEmzZt8O6773okOBER/cahmfQm\n+vjjj0WfPn1EixYtxI4dO+p9LyMjQxiNRtG7d2+xevVqd8Zwifz8fHHdddeJ2NhYERcXJwoKClRH\nctqcOXNERESE6Nu3r/jb3/6mOo7TXn31VaHT6cT333+vOopTnnzySRERESGio6PFXXfdJc6cOaM6\nkkNWrVolevfuLYxGo5g1a5bqOE4pKSkRJpNJ9OnTR/Tt21e8+eabqiM1SXV1tYiNjRV33HGH3ePc\nWuT79+8XBw8eFCaTqV6R7927V8TExAir1SqOHTsmwsLChM1mc2eUZrv55ptFbm6uEEKInJwcYTKZ\nFCdyzvr168Xw4cOF1WoVQgjx3XffKU7knJKSEhEfHy9CQkI0V+Sff/553d/v6dOni+nTpytO1Ljq\n6moRFhYmjh07JqxWq4iJiRH79u1THcthJ06cELt27RJCCFFZWSl69eqlqfy1XnvtNTFu3DgxatQo\nu8e5dfPliIgI9OrV65J/np2djbFjxyIwMBAhISEwGo0oKChwZ5Rm69atG3744QcAwJkzZzR3Qvff\n//43nnnmGQT+dk13586dFSdyzuOPP45XXnlFdYwmue2229CihfxRGzhwIEpLSxUnatyF15AEBgbW\nXUOiFV27dkVsbCwAoG3btoiMjMTx48cVp3JOaWkpcnJyMGnSpEZX/rm1yBty/Pjxeqtf9Ho9ysrK\nVERx2KxZs/DEE0+ge/fueOqppzBz5kzVkZxSVFSETZs2YdCgQTCZTNi+fbvqSA7Lzs6GXq9HdHS0\n6ijNtmDBAowcOVJ1jEaVlZXBYDDUPdfCz2hDiouLsWvXLgwcOFB1FKc89thj+Oc//1k3CLCn2Xc/\nvO222/Dtt99e8s8zMjIwatQoh9/H0TXr7tTQ/5YZM2Zgzpw5mDNnDu666y4sW7YMycnJWLNmjYKU\nDbOXv7q6GhUVFcjLy8O2bdtwzz334OjRowpSXp697DNnzsTnn39e988aG52o4MjPwYwZM9CqVSuM\nGzfO0/Gc5g0/j67w008/YcyYMXjzzTfRtm1b1XEc9tlnn6FLly7o37+/Y/eL8cQ8z8Vz5DNnzhQz\nZ86sex4fHy/y8vI8EaXJ2rVrV/d1TU2NuOqqqxSmcV5CQoIwm811z8PCwkR5ebnCRI7Zs2eP6NKl\niwgJCREhISEiICBA9OjRQ5w8eVJ1NKe8++674oYbbhC//PKL6igO+fLLL0V8fHzd84yMDM2d8LRa\nrWLEiBHijTfeUB3Fac8884zQ6/UiJCREdO3aVVx55ZVi/PjxDR7vsSK/8LL92pOd586dE0ePHhWh\noaGipqbGE1GarH///nVFuHbtWhEXF6c4kXPefvtt8fe//10IIcTBgweFwWBQnKhptHiyc9WqVaJP\nnz7i1KlTqqM4rKqqSoSGhopjx46Jc+fOae5kZ01NjRg/fryYNm2a6ijNZjab1a5aWb58udDr9eKK\nK64QQUFBIiEhoe57M2bMEGFhYaJ37951q0G82bZt28T1118vYmJixKBBg8TOnTtVR3KK1WoV9913\nn4iKihK///3vxYYNG1RHapKePXtqrsiNRqPo3r27iI2NFbGxsWLKlCmqIzkkJydH9OrVS4SFhYmM\njAzVcZzyxRdfCJ1OJ2JiYur+va9atUp1rCYxm82NrlrxyMYSRETkPkpWrRARkeuwyImINI5FTkSk\ncSxyIiKNY5ETEWkci5yISOP+H7uX50jMp2DhAAAAAElFTkSuQmCC\n"
}
],
"prompt_number": 9
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"\u3053\u306e\u3088\u3046\u306b\u95a2\u6570\u3092\u308f\u3056\u308f\u3056\u5b9a\u7fa9\u3057\u3066\u3001x\u3001y\u3068\u3044\u3046\u5909\u6570\u3092\u4f7f\u308f\u306a\u304b\u3063\u305f\u306e\u306f\u3001\u666e\u901a\u306e\u30a2\u30eb\u30d5\u30a1\u30d9\u30c3\u30c8\u306ea\u304b\u3089z\u306fsympy\u306e\u30b7\u30f3\u30dc\u30eb\u3068\u3057\u3066\u4e88\u7d04\u3057\u3066\u3057\u307e\u3063\u305f\u304b\u3089\u3067\u3059\u3002\n",
"\n",
"\u6b21\u306b\u3001scipy\u3067\u7a4d\u5206\u3092\u5b9f\u884c\u3057\u3066\u307f\u307e\u3057\u3087\u3046\u3002"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from scipy import integrate\n",
"integrate.quad(exp0_func,-10.0,np.pi)\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"latex": [
"$$\\begin{pmatrix}6.94987255622, & 3.9625733733e-11\\end{pmatrix}$$"
],
"output_type": "pyout",
"prompt_number": 10,
"text": [
"(6.94987255622, 3.9625733733e-11)"
]
}
],
"prompt_number": 10
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"\u6700\u521d\u306e\u5024\u304c\u7a4d\u5206\u7d50\u679c\u30012\u756a\u76ee\u306e\u5024\u306f\u8aa4\u5dee\u306e\u898b\u7a4d\u3082\u308a\u3067\u3059\u3002\n",
"\u3053\u306e\u3088\u3046\u306b\u3001Numpy,Scipy\u3068Sympy\u306e\u9593\u306e\u5f0f\u306e\u5909\u63db\u306f\u624b\u52d5\u3067\u305d\u308c\u307b\u3069\u9762\u5012\u3067\u306f\u3042\u308a\u307e\u305b\u3093\u3002\n",
"\n",
"\u3068\u3053\u308d\u304c\u3001\u5b9f\u306f\u3001"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"exp1.evalf()\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"latex": [
"$$6.94987255622221$$"
],
"output_type": "pyout",
"png": "iVBORw0KGgoAAAANSUhEUgAAAJ0AAAASCAYAAABBwNzbAAAABHNCSVQICAgIfAhkiAAABVhJREFU\naIHtmWlsVUUUx3+vtD67QNgKxWBABRW3WIMkgjZuiUG/YBMTiASXL0LcUHFNwMZCRDDVKIpo1Icx\nKtaFD5igxiD4oY1bJBpJlYho3agRd1tc6of/XN+86dx35xb1fnn/5OXdOfOfc889d+bMOXOhggr+\nZ+Q8sqOBNuA3YMD8twE/Jeg6Ergd6Af+BOqBm4GvE8YtA2qAO4epbzJwKzAI5IE6YDXwvsWZCzwE\n7LSe6y+r/y3gAXM9C7gFqDW63wRWAF849+0B2oFtwC/ATDPuamDXMHgQ5vtQ+7LiRagDuoGTYvr/\nwXHAZ8Bs024C9qAHL4cjgG+BVku2EL3k6jLjpqAX4eoP1TcO2AxMcsb2ANMt2Q1oUsb95hreKcAr\nwGjTbgB2APuAqY6Nro4DwJWeZwzlhfg+1L6seBFORQt50NNXgmr0spZassOBPuDahLGbge+AKktW\ni6LUpWXGPWwMaxumvhtRxHBxG3CX1V6PnqXG0TkHWGe1XwKmObqajY3POPK9xv4XgTXADI8dobxQ\n34falxVvhuEWUJRLnHSXo1U4Oono4BDgdxRuXfQAW2LGtQLzGTrp0uhbj16mi+uBDqu9zsNpALai\nbSDCzyjaTHC4+1HktfG6R6cPIbxQ34falxXPRoGYSWev+gXAbuD7GCVxGItWar+n7weUC7hoAM5n\naPRIq+89YB7wFDDGyPJoKy5YvKs8utainPFXS7YHmIjyRxsDKNL+Vwj1fah9WfFSIYeS1u1AC7AS\nuBdtc80JY0egvMwXmXrRbHfzulXAUebajXRp9OUp5g5focn2JMUcLQ5z0Jbnoh7lUjYOM/q3OfIu\nYDny01rgBVQIuEjipfF9qH1Z8WwUSNhexxvCLmCxJT8TVU7HlxsMPIbCrF0NN6EKcZDSsHwyqjYj\n+HK6NPpGom0yStS3oFVZDu+ggiMEq1H1PNuR70Z5V4SFwDcMfTlJvIP1fZx9WfMKJEy6iYbQz9Bw\n2Ut8XhahEfgIuMK0q9GKfdfoHWfkVcBGlLdF8E26UH2gMv4R4ALgE4pR78QYW88BPkx4ngjTUD6z\n0tNX5bRHoK36vpS8g/F9Ofuy5hVImHQ1hvCBp68b7d35BEPGoBypwxgzBXgbbR1RxFoCnOWM8026\nUH3XoWopQj1wD4qIO2PsfB5F0iTk0RbfkUS0sNf80vCG6/tQ+7LiFQioXvcBb3jk283gSZ6+JHwO\nvGaum4D7PZy4SZekL4eqpxM8vMVGr3s8UYNW6KqE++RQcXJHTP8O/L7qpbQACuWl9X2SfVnzILB6\n7Ub5hYs8Wm19ATey0YhOrjtN+1zgWJQgR78oSs037Vbi4eprREcMH3u4G4AfgUMd+SwUDePK/Ajt\nKMdaYckWWdfNqAJ3MZ7SSBfKS+v7JPuy5gVjAco17Lwih8r4Zx3udIe3FK3WyZZsGVrR7ou3MRV/\npAvRl0OfxE736B1puO62dIm535IyNl2GfyXb1W4npXklFA9Llw+Dl8b3IfZlyYtQICbS2UcZm4Br\n0Ceau43sInRQe5PFa0EHnq8C5xlZA3LaH6bdjL4WXIj/vC1CjfMfIUTfIDqt34CKiE+NfBTwKJq4\nA47eqOo9EGPP2ehYYys6eolQTakD15j7LjJ25lB+2WX60vJCfR9qX1Y8G1FwqKP0LHQIxqIZ2on2\n7qcpnqdFOAaV/A9aslpj1OPAc2hCnlbmPqPQ+c6XFCu3LnTQm1ZfC9qaN6HD5k503ODDPBQ9Zsb0\n7yf++2y7wz3D3OsJc/92/FE9lBfi+1D7suJNAF5GRVHU34fe9cWeZ66gggoqqKCCCv4d/A38FTCb\nmtYBiAAAAABJRU5ErkJggg==\n",
"prompt_number": 11,
"text": [
"6.94987255622221"
]
}
],
"prompt_number": 11
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"\u3053\u308c\u3060\u3051\u3067\u7a4d\u5206\u304c\u5b9f\u884c\u3067\u304d\u3066\u3057\u307e\u3046\u306e\u3067\u3001\u4eca\u306e\u5834\u5408\u306a\u3089Scipy\u3092\u4f7f\u3046\u5fc5\u8981\u306f\u3042\u308a\u307e\u305b\u3093\u3067\u3057\u305f\u3002\n",
"\u30d7\u30ed\u30c3\u30c8\u306b\u95a2\u3057\u3066\u3082\u3001"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"from sympy.plotting import plot\n",
"plot(exp0,(x,-10.0,pi))\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAWkAAAD9CAYAAAB++4avAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlYlWXeB/DvUdDMvVTIczCEw6pshuJyaccxgzHFchnB\nskbNeJ3XymkZs3lrsMmt3bKF5lUbk3ytScMaIDPFNRZX3HJIJQG3aDR1NA8c7vePO04gi8BZ7uec\n8/1c17nkeJ6e53cpfrv5Pfdz3zohhAAREWlSK9UFEBFRwxjSREQaxpAmItIwhjQRkYYxpImINIwh\nTUSkYQxpIjc3bdo0+Pj4ICIiosFjcnJyEBMTg759+8JkMjmvOLohHedJE7m3bdu2oUOHDnjwwQdx\n4MCBOp9fuHABQ4YMwZdffgmDwYDy8nJ069ZNQaVUH46kidzc0KFD0bVr1wY//+ijjzB+/HgYDAYA\nYEBrjJfqAohIraKiIlRUVGD48OG4dOkSHn/8cUyZMqXOcTqdDn/5y1+s700mE1sjTsCQJvJwFRUV\n2LNnD77++mtcuXIFgwYNwsCBAxEUFFTn2NTUVOcX6OEY0kQezs/PD926dUO7du3Qrl07DBs2DPv3\n7683pMn52JMm8nBjx47F9u3bYbFYcOXKFeTl5SE8PFx1WfQLjqSJmqCyEli9GsjPB86cASZNAsaP\nB3Q61ZXdWHJyMrZs2YLy8nL4+flh3rx5qKioAACkpKQgNDQUCQkJiIyMRKtWrTBjxgyGtIZwCh7R\nDVy8CCQlAd27A6NHAxUVwKJF8usXXwRaecjPozqdDowL52NIEzWiqgq4914gOBhYvBho3Vr+/pUr\nwIMPAiNGADNnqq3RWRjSajCkiRrx4otAVhaweTPQpk3tz4qK5Aj788+Bnj3V1OdMDGk1POQHNaLm\n27IFePdd4JNP6gY0AAQFASYT8MYbTi+NPAhH0kT1MJuB8HDg7beB+PiGjystBfr1A44dAzp2dF59\nKnAkrQZH0kT1eOcdICSk8YAGAIMBGDYM+Owz59RFnocjaaLrnD8vA3rzZqBPnxsfn5EBvPaabI+4\nM46k1eBImug6CxYAY8c2LaABICEBuHYNOH3asXWRZ2JIE9Vw7hywYwfQnCUq2rYFbr9dzgIhsjeG\nNFENS5cCERGAXt+8/27sWKCepZqJbMaeNNEvrlwB/P2B7dvlwyvNceqUDPcffnDfJxDZk1bDTb+d\niJpvxQpgyJDmBzQgH2bp1g0oLLR/XeTZGNJEkAsovfYa8PTTLT9HfDywbZv9aiICGNJEAORNv+ho\nYPDglp8jJgb45hv71UQEMKSJAADLl8updLYYNAjYudM+9RBV441D8nhnzwKhocDJk7Y92i0EEBUF\nfPUV4ONjv/q0gjcO1eBImjzeypXAfffZvvaGTifXnN6zxz51EQEMafJwQgDLlgHTp9vnfIMHA0eP\n2udcRABDmjxcdQ/ZlhuGNYWFsS9N9sWQJo9WPYq2116FkZGcK032xZAmj3XxIrBundwGy15CQmRv\n+z//sd85bTVt2jT4+PggIiKi0eMKCgrg5eWFtWvXOqkyagqGNHmsNWvkzir2nInh7S1XxNNSX3rq\n1KnIzs5u9BiLxYI5c+YgISGBMzg0hiFNHsueNwxrCgsDjhyx/3lbaujQoejatWujx7z11luYMGEC\nunfv7qSqqKm8VBdApMKhQ0BJie0PsNQnPBw4fNj+53WUsrIyZGRkYNOmTSgoKICukQZ9ao01XE0m\nE0wmk+ML9HAMafJIy5YBDz0EeDngX0BUlNzVxVXMnj0bixYtsj6s0li7I7U5C22TXTCkyeOYzcCq\nVY6bKhcQAPz5z445tyPs3r0bSUlJAIDy8nJkZWXB29sbiYmJiisjgCFNHig7Gxg9GjAaHXN+oxE4\nfhywWIDWrR1zDXs6fvy49eupU6dizJgxDGgNYUiTx1m50jG96Go33ywfDy8uBgIDHXedpkpOTsaW\nLVtQXl4OPz8/zJs3DxUVFQCAlJQUxdXRjXCBJfIoFy7I/Qi//x7o0sVx1xkzBnj0UeDuux13DWfj\nAktqcAoeeZRPPwVGjHBsQAOAr69seRDZiiFNHiU9Hbj/fsdfJyCAIU32wZAmj1FWBuzbB9xzj+Ov\nFRgIHDvm+OuQ+2NIk8dYvVquG33TTY6/FkfSZC8MafIYzmp1ADKkW7WS61UT2YIhTR7h8GHg3Dng\nzjudc71bbgG++w44f9451yP3xZAmj5CeDiQnO/fhkuqpfkS2YEiT2xMC+Ogj57U6qjGkyR4Y0uT2\ndu4E2rUDoqOde12GNNkDQ5rcXvUNQ3ttkdVUDGmyB67dQW6togL45BMgP9/51w4KknOziWzBkTS5\ntQ0b5MMrvXs7/9o+PsA33zj/uuReGNLk1j78EBg4UM21DQagtFTNtcl9cBU8clsXLwJ+fvLJv1tv\ndf71KyrksqVXrzpmBxhn4yp4anAkTW5r7Vq5G7iKgAbkzuHduwNnzqi5PrkHhjS5rfR04IEH1NbA\nlgfZiiFNbunUKWD3brlNlkp6PUOabMOQJrdUveJdu3Zq6zAYOA2PbMOQJre0apX6VgfAdgfZjiFN\nbufgQaC83Hkr3jVGC+2OadOmwcfHBxEREfV+np6ejqioKERGRmLIkCEoLCx0coXUGIY0uZ3qFe9a\naeC7WwvtjqlTpyI7O7vBzwMCArB161YUFhbiueeewyOPPOLE6uhG3GD2JtGvqqpkSH/xhepKJL0e\nqKxUW8PQoUNRXFzc4OeDBg2yfh0XF4fSRob+qamp1q9NJhNMJpMdKqTGMKTJrezYAfTrB0RGqq5E\n6tkT2LtXLpfq7AWeWmLZsmUYNWpUg5/XDGlyDoY0uZXly4GhQ1VX8av27WXb5fJloGNH1dU0bvPm\nzVi+fDl27NihuhSqQQNdOyL7uHwZWLfO+Yv734iPD3D2rOoqGldYWIgZM2Zg/fr16Nq1q+pyqAaG\nNLmNTz+Vo2hfX9WV1Kb1kD558iTGjRuHVatWwWg0qi6HrsN2B7mNDz4AHn1UdRV1qQ7p5ORkbNmy\nBeXl5fDz88O8efNQUVEBAEhJScELL7yA8+fPY+bMmQAAb29v5KtYgJvqxVXwyC0cPw7Excnpbm3a\nqK6mtpQUuXXXLxnosrgKnhpsd5BbWLlSzo3WWkAD6kfS5NoY0uTyqqqAv/8dmDpVdSX1Y0iTLRjS\n5PK2bgU6dXL+buBNxZAmWzCkyeV98AHw+99r92ERhjTZgiFNLu2nn+Ta0ZMnq66kYQxpsgVDmlxa\nejrQpYsMQq3q0QM4d051FeSqGNLksoQA3ntPTnHTss6dAbNZbkhL1FwMaXJZubky+IYPV11J43Q6\nOZpmy4NagiFNListTY6itbBu9I2wL00t5QLf3kR1nT8PZGTIWR2ugCFNLcWQJpe0ciXw298C3bqp\nrqRpGNLUUgxpcjlCyFbHf/2X6kqajiFNLcWQJpezbZv8VUuL+98IQ5paiiFNLqf6hqFWnzCsD0Oa\nWoohTS7lzBn5euAB1ZU0D0OaWoohTS7lnXeA4GDg1ltVV9I8fOqQWoo7s5DLuHpVtjq2blVdSfNx\nJE0txZE0uYxVq4D+/YGQENWVNN8ttwCXLsnHw4magyFNLkEI4PXXgSeeUF1Jy7RqBXTvrqblMW3a\nNPj4+CAiIqLBYx577DEEBQUhKioKe/fudWJ1dCMMaXIJX34JeHtrf52OxqhqeUydOhXZ2dkNfp6Z\nmYnvvvsORUVFeP/9960b0pI2MKTJJbz2mhxFu9K0u+upCumhQ4eia9euDX6+fv16PPTQQwCAuLg4\nXLhwAWfZQNcM3jgkzTt4UL6SklRXYhut3jwsKyuDn5+f9b3BYEBpaSl86lmkOzU11fq1yWSCyWRy\nQoWejSFNmvf668Af/gC0bau6EttoNaQBQAhR672ugR9ZaoY0OQfbHaRpp08DxcWutU5HQ7Qa0nq9\nHiUlJdb3paWl0Ov1CiuimhjSpGkvvwxERLjOaneN0WpIJyYmYuXKlQCA3NxcdOnSpd5WB6nBdgdp\n1pkzcifwgwdVV2Ifqp46TE5OxpYtW1BeXg4/Pz/MmzcPFRUVAICUlBSMGjUKmZmZMBqNaN++PVas\nWOH8IqlBOnF9M4pII558EqisBJYsUV2JfezfL9ccOXBAdSUto9Pp6vSuyfE4kiZNOnsWWLHCfUbR\ngHbbHaRt7EmTJr38shx19uypuhL76dZN7hz+S6eBqEk4kibNOXcOWL7cddsCDfHyAi5cAH78EfD1\nVV0NuQqOpElzXn4ZmDwZcMdZYD16AD/8oLoKciUcSZOmnDsHLFsGFBaqrsQxundnSFPzcCRNmvLK\nK0ByMmAwqK7EMRjS1FwcSZNmnDwJfP458NVXqitxHFXLlZLr4kiaNGPuXGDiRPcdRQPsSVPzMaRJ\nE3JzgS1bgD/9SXUljsV2BzUXQ5qUEwKYPRuYPx/o0EF1NY7FkKbmYkiTcqtXy8e/p0xRXYnjsSdN\nzcUbh6TUlSvAM88A6elyH0B3x540NZcH/LMgLXv1VSAuDhg6VHUlzsF2BzUXV8EjZU6dkmtFFxQA\nAQGqq3EOi0XuMHPtGtC6tepqmoer4KnBkTQp88ILwNNPe05AAzKYu3SR63cQNQV70qTEP/8pH1px\nt0WUmqJ68f8ePVRXQq6AIU1Od/EiMHOmXC/65ptVV+N87EtTc7DdQU43dy4wciQwYoTqStRgSFNz\ncCRNTrVjB7BuHXDokOpK1GFIU3NwJE1O8/PPwMMPA2++CXTtqroadVRsSJudnY3Q0FAEBQVh8eLF\ndT4vLy9HQkICoqOj0bdvX3zwwQfOLZAaxJAmp1myBIiJAcaPV12JWs4eSVssFsyaNQvZ2dk4fPgw\nVq9ejSNHjtQ6ZunSpYiJicG+ffuQk5ODJ598EpWVlc4rkhrEkCanyMmRIf3664BOp7oatZwd0vn5\n+TAajfD394e3tzeSkpKQkZFR65jbbrsNFy9eBABcvHgRt956K7y82A3VAv4tkMP9+KNcl2P5crlj\ntqdzdkiXlZXBz8/P+t5gMCAvL6/WMTNmzMBvfvMb9OzZE5cuXcLHH39c77lSU1OtX5tMJphMJkeU\nTDUwpMmhhJB96IkTgYQE1dVog7PX79A14UeXBQsWIDo6Gjk5OTh27BhGjhyJ/fv3o2PHjrWOqxnS\n5Bxsd5BDvfce8P33wMKFqivRDmevhKfX61FSUmJ9X1JSAsN1Oyvs3LkTEydOBAAEBgaid+/eOHr0\nqPOKpAYxpMlhDh4EnntOLkXatq3qarTj1luB8+flOh7OEBsbi6KiIhQXF8NsNmPNmjVITEysdUxo\naCg2btwIADh79iyOHj2KAE96Xl/D2O4gh7hwAZgwAXj7bSAkRHU12uLlBXTuDPz733JU7fjreWHp\n0qWIj4+HxWLB9OnTERYWhrS0NABASkoKnn32WUydOhVRUVGoqqrCSy+9hFtuucXxxdENcRU8sjuL\nBRg9GggKknOiqa6wMODTT4HwcNWVNB1XwVOD7Q6yu2efBcxmuVY01Y87tFBTsd1BdpWeDnzyiVwj\n2ttbdTXaxUfDqakY0mQ3u3bJDWU3bZI3x6hhDGlqKrY7yC6+/17O5EhLk7utUOO41yE1FUOabFZe\nDsTHA7/9LTBunOpqXANH0tRUDGmyyeXLwKhRctGkxx5TXY3r4I1DaiqGNLWY2SxHzlFRwIsvqq7G\ntXTrBvz0k+oqyBUwpKlFKirkJrKdOwPvvsuV7ZrLxwcoLVVdBbkCzu6gZquoAJKTgatX5QMZXNGy\n+W67DThzRnUV5Ar4xCE1i9kMTJoknyr85BOuydFSVVXATTcBly65zp8hnzhUg+0OarJr1+R6HADw\nj3+4TrhoUatWarbRItfDkKYmuXgR+P3v5d6EH38MtGmjuiLX5+vLlgfdGEOabujUKeDOO+VNwmXL\n+Li3vTCkqSkY0tSoQ4eAwYOB3/1OzuLgTUL7YUhTU/CfHDUoJ0feJHz1VeCBB1RX434Y0tQUDGmq\n17JlwIcfAh99BIwYoboa9+TrCxw5oroK0jqGNNViNgOPPw5s3gysWycXpyfH8PWVf85EjWFIk1VZ\nmdzVu0cPID8f6NRJdUXuje0OagreOCQAQGam3PJqzBhg7VoGtDMwpKkpGNIe7uefZXtj5kxgyRJg\n7lz5oAU5XnVIO+MhvuzsbISGhiIoKAiLFy+u95icnBzExMSgb9++MJlMji+KmoSPhXuwQ4eAyZPl\nhrF/+5t8UIWcq317GdQdOzruGhaLBSEhIdi4cSP0ej369++P1atXI6zGDYcLFy5gyJAh+PLLL2Ew\nGFBeXo5u3brVOg8fC1eDYyYPVFkJvPSSHD0/+qhcg4MBrYYzWh75+fkwGo3w9/eHt7c3kpKSkJGR\nUeuYjz76COPHj4fBYACAOgFN6vDGoYc5cACYOhXo0gX4+9+B3r1VV+TZqkM6KMhx1ygrK4Ofn5/1\nvcFgQF5eXq1jioqKUFFRgeHDh+PSpUt4/PHHMWXKlDrnSk1NtX5tMpnYFnEChrSHuHoVeOMN+Zo/\nH5g+nWtAa4GvL3D2rGOvoWvCX3RFRQX27NmDr7/+GleuXMGgQYMwcOBABF33f4+aIU3OwZD2AJ9/\nLm8OxsYCe/YAer3qiqiaM9oder0eJSUl1vclJSXWtkY1Pz8/dOvWDe3atUO7du0wbNgw7N+/v05I\nk/OxJ+3Gjh4FHn5Y7qDy/vty9ToGtLY4I6RjY2NRVFSE4uJimM1mrFmzBomJibWOGTt2LLZv3w6L\nxYIrV64gLy8P4eHhji2MmoQh7YZ++AGYNQsYMgSIiQEKC4G77lJdFdXHGSHt5eWFpUuXIj4+HuHh\n4Zg0aRLCwsKQlpaGtLQ0AEBoaCgSEhIQGRmJuLg4zJgxgyGtEZyC50YuX5Yr1b38MpCUBDz/vNzw\nlLTr88+BtDTgiy9UV3JjnIKnBnvSbuDKFeCdd2Q4jxgB7NwJGI2qq6Km4FOHdCMMaRd2+TKQng7M\nmyfXfP76a6BvX9VVUXMwpOlGGNIu6McfgaVLgbffljum/POfsvdMrqd6n8OqKj6OT/Xjt4ULKSoC\nnn1WPvhQWgps3y6fFmRAu662beUj4f/+t+pKSKsY0honBLBpk1ydbvBg+f7gQbnWRnCw6urIHtjy\noMaw3aFRP/wAfPCBnN8cHQ0kJgJr1gA336y6MrK36pDm/QSqD0NaQ6qqgG3bgPfeA7KygHvvBVau\nBAYO5CPc7iwkBDh9WnUVpFUMaQ04ckSGcXq6HFVNmSKn1HFlOs/QqRNQ46ltoloY0oqUlso9BFeu\nlNtWPfCAnKUREaG6MnI2g4Eb0lLDGNJOVFYGfPqpXEPjyBFg7Fhg4UJg+HCgdWvV1ZEqBgPw1Veq\nqyCtYkg72LffAp99Jl/t2gG33y6n0d11F9CmjerqSAv0evmTFVF9uHaHnV27Jucvf/EFkJ0tb/gN\nHy5vAt55J4OZ6jp9Ws7gcfS60rbi2h1qMKRtJARw/DiwYYMM5aNHgc6d5c7b99wj//HxSTJqjMUi\nf8q6dEk+3KJVDGk12O5ogbNngc2bgY0b5XoZt90mnwKcNAkYORLo3l11heRKWreWs3rKyoCAANXV\nkNYwpJvg1CkgN1fe3NmyBTCbgT595IpzTzwBhIVxHjPZJjRUTsNjSNP1GNLXsViAw4eBvDxg61bZ\nX754EUhOlqPlGTOAqCjOxiD76t4dOHlSdRWkRR4f0qdPA7t2yZFybi5QUCDbF/fcAwwbJmdihIRw\npEyO1asXQ5rq51EhfeaM3Ih11y5g3z45Wv75Z6B/f/no9VNPAQMGALfeqrpS8jS33y6/N4mu5zEh\n/frrwIsvAv36ydf99wOvvgr4+3OUTOr16iXn0hNdz2Mmh/33fwPl5fLm3+LFwPjxQO/eDGjShl69\ngO+/d9z5s7OzERoaiqCgICxevLjB4woKCuDl5YW1a9c6rhhqFo8J6TZtGMikXdU9aUdMQ7ZYLJg1\naxays7Nx+PBhrF69GkfqWSzEYrFgzpw5SEhI4HxoDfGYkCbSsk6d5P2Q8nL7nzs/Px9GoxH+/v7w\n9vZGUlISMjIy6hz31ltvYcKECejOif6a4jE9aSKtu3gROHHC/g9DlZWVwc/Pz/reYDAgLy+vzjEZ\nGRnYtGkTCgoKoGvgx87U1FTr1yaTCSaTyb7FUh0MaSKNCAiQSwwMGGDf8zYUuDXNnj0bixYtsj76\n3VC7o2ZIk3MwpIk0ojqk7U2v16Okxq4CJSUlMBgMtY7ZvXs3kpKSAADl5eXIysqCt7c3EhMT7V8Q\nNQtDmkgjAgLkw1T2Fhsbi6KiIhQXF6Nnz55Ys2YNVq9eXeuY4zX+7zB16lSMGTOGAa0RvHFIpBGB\ngcCxY/Y/r5eXF5YuXYr4+HiEh4dj0qRJCAsLQ1paGtLS0ux/QbIrLlVKpBHHj8u1xx05X9oWXKpU\nDYY0kUZUVAAdOsh1pbW4OQRDWg22O4g0wttbjqQdcfOQXBdDmkhDWrcG/vUv1VWQljCkiTQkOFhu\nwUZUjSFNpCEhIRxJU20MaSIN4UiarseQJtKQkBCGNNXGkCbSkJ495SbH58+rroS0giFNpCE6HXD5\nstwMmQhgSBNpTp8+wKFDqqsgrWBIE2kMQ5pqYkgTaUzfvsDBg6qrIK1gSBNpTESE3EaLy2QQwJAm\n0pyePYFTp4DTp1VXQlrAkCbSGJ0OiIkB9u1TXQlpAUOaSIOioxnSJDGkiTQoJgYoLlZdBWkBQ5pI\ng/r1AzZsUF0FaQFDmkiDgoKAn34Czp5VXQmpxpAm0qBWrYABA4C8PPucLzs7G6GhoQgKCsLixYvr\nfJ6eno6oqChERkZiyJAhKCwstM+FyWYMaSKNiouzT0hbLBbMmjUL2dnZOHz4MFavXo0jR47UOiYg\nIABbt25FYWEhnnvuOTzyyCO2X5jsgiFNpFFDhwJlZbafJz8/H0ajEf7+/vD29kZSUhIyMjJqHTNo\n0CB07twZABAXF4fS0lLbL0x24aW6ACKq34ABwLhxwPvv27Z7eFlZGfz8/KzvDQYD8hoZoi9btgyj\nRo2q97PU1FTr1yaTCSaTqeWFUZMwpIk0qnNneQNx925g0KCWn0en0zX52M2bN2P58uXYsWNHvZ/X\nDGlyDrY7iDRs2DAgJ8e2c+j1epSUlFjfl5SUwGAw1DmusLAQM2bMwPr169G1a1fbLkp2w5Am0rCR\nI4Ht2207R2xsLIqKilBcXAyz2Yw1a9YgMTGx1jEnT57EuHHjsGrVKhiNRtsuSHalE4JrbRFp1cWL\ngF4v50vffHPLz5OVlYXZs2fDYrFg+vTpmDt3LtLS0gAAKSkpePjhh7Fu3Tr06tULAODt7Y38/Pxa\n59DpdGBcOB9Dmkjjhg4F/ud/gPh4tXUwpNVgu4NI4yZNAhq4j0cegCNpIo07eBC45x654FIzJmrY\nHUfSanAkTaRxffrIedJcutQzMaSJNE6nAx56CNi4UXUlpALbHUQuYM8eYMIE4NgxdS0PtjvU4Eia\nyAXExADt2vEGoidiSBO5AJ0OmDUL+Owz1ZWQs7HdQeQizp0DQkJky+OWW5x/fbY71OBImshF9OgB\njBkD/O//qq6EnIkjaSIXsn8/cP/9QEGB7FE7E0fSajCkiVzMuHHA4MHAU08597qeFNIWC1BSAnz3\nHXD8eO3XffcBf/6z82phSBO5mEOHgOHDZYB06uS867pbSFdVySAuKpKvf/0L+PFHID9fPt3ZvTtg\nNALh4XKRq4AA+QoKApy5kitDmsgFTZkin0R85hnnXdNVQ/ryZeDoUeDIEfnrt9/KX2+6SW5PFhws\ngzc4WL6MRqB3b+e3kxrCkCZyQWVlQHQ08PXXQGSkc66p9ZD+6Sfg8GH5k8bRo8CBA/J9eTlQvXx2\naKicIRMaKoO5Qwe1NTcFQ5rIRX3wAfDGG/LHc1v2QGwqrYT0zz/L0fCBA/J18KBsXWzbBoSFyZ8w\n+vUDAgPle39/oHVr1VW3HEOayEUJAfzhD/Im17vvOj6IVIT0mTNAYSGwd6+c2VJYKOeJ33233ASh\nb18gIkL+6u8PtHLDScUMaSIX9p//AGPHyptcK1cC3t6Ou5YjQ7qqSt4I3bdPrlOyb598VVbKVkXX\nrrKtExUlR8dt2zqkDE1iSBO5uKtX5eJLbdsC//d/jmt92CukKyvlTbzDh4GdO2Uo798PdOsGJCQA\nPXvKfnt0tJxVoXINbS2w+YeDHFu3MnYSV6jTFWoEWKe92Vpnu3bA2rWy/TFunOzZaoXZLFsV6emy\nNRMXB3TuDEycKDfY7dULmDcP+P57OQf5nXfkVmGjRwMGQ8sC2t3+3hnSGuIKNQKs097sUWfbtsDH\nHwM+PvLR8fXrZWhXy87ORmhoKIKCgrB48eJ6z/HYY48hKCgIUVFR2Lt3b7NruHwZ+OYbGbRPPgnc\ncQfQpQvw4IPArl1yVsUrr8g+87ffAm+9JY/7zW/sO+/Y3f7evRxbBhE5i7e3XNcjIwN4/nnghRfk\nKDU+3oJZs2Zh48aN0Ov16N+/PxITExEWFmb9bzMzM/Hdd9+hqKgIeXl5mDlzJnJzc+u9TlWVfNjj\n4EH5OnFCzqwoKZEPfkRHyyciJ06UfWRbdjknhjSRW9HpgHvvlTfb1q2TD7v85S/5uHrViKwsfwwY\nAEycmISMjIxaIb1+/Xo89NBDAIC4uDicP38B+/adxZUrPjh2TN7UA+S6IRkZchW+vn3la+RI4PHH\n5UjZkTcuPZXNNw51nt7VJyJqoabEr809aSEEX3zxpeHXP/7xDzz88MPW9x9++CFmzZpV65jRo0dj\n+/bt1vcjRozA7t27ldfu7q+mcMOp30RUk16vR0lJifV9SUkJDAZDo8eUlpZCr9c7rUZqGEOayM3F\nxsaiqKgIxcXFMJvNWLNmDRKrF7P4RWJiIlauXAkAyM3NRZcuXeDj46OiXLoObxwSuTkvLy8sXboU\n8fHxsFgsmD59OsLCwpCWlgYASElJwahRo5CZmQmj0Yj27dtjxYoViqsmK9FCH3/8sQgPDxetWrUS\nu3btsv6jP5Q+AAAFhklEQVT+hg0bxB133CEiIiLEHXfcITZt2tTSS9isZo27d++u9dmCBQuE0WgU\nISEh4ssvv1RUYV15eXmif//+Ijo6WsTGxor8/HzVJTXozTffFKGhoaJPnz7iT3/6k+pyGvXKK68I\nnU4nfvzxR9Wl1Oupp54SoaGhIjIyUtx3333iwoULqkuqJSsrS4SEhAij0SgWLVqkupx6nTx5UphM\nJhEeHi769OkjlixZorqkBlVWVoro6GgxevToGx7b4pA+cuSIOHr0qDCZTLUCcO/eveL06dNCCCEO\nHjwo9Hp9Sy9hs4ZqPHTokIiKihJms1mcOHFCBAYGCovFoqzOmu68806RnZ0thBAiMzNTmEwmxRXV\nb9OmTeKuu+4SZrNZCCHEuXPnFFfUsJMnT4r4+Hjh7++v2ZDesGGD9Xtwzpw5Ys6cOYor+lVlZaUI\nDAwUJ06cEGazWURFRYnDhw+rLquO06dPi7179wohhLh06ZIIDg7WZJ1CCPHqq6+KyZMnizFjxtzw\n2Bb3pENDQxEcHFzn96Ojo+Hr6wsACA8Px9WrV1FRUdHyob4NGqoxIyMDycnJ8Pb2hr+/P4xGI/Lz\n8xVUWNdtt92Gn376CQBw4cIFzd68effddzF37lx4/zIxtnv37ooratgTTzyBl156SXUZjRo5ciRa\n/bKEW1xcHEpLSxVX9Kv8/HwYjUb4+/vD29sbSUlynrXW+Pr6Ijo6GgDQoUMHhIWF4dSpU4qrqqu0\ntBSZmZnWGTc34tAbh59++inuuOMO6z9krTh16lStu9sGgwFlZWUKK/rVokWL8OSTT6JXr154+umn\nsXDhQtUl1auoqAhbt27FwIEDYTKZsGvXLtUl1SsjIwMGgwGRzloZ3w6WL1+OUaNGqS7DqqysDH5+\nftb3Wvr30pDi4mLs3bsXcXFxqkup449//CNefvll6/+Ub6TRG4cjR47EmTNn6vz+ggULMGbMmEZP\nfOjQITzzzDP46quvmlRIS9lSY03OfCinoZrnz5+PN998E2+++Sbuu+8+fPLJJ5g2bZrD/wwb0lid\nlZWVOH/+PHJzc1FQUIDf/e53OH78uIIqG69z4cKF2LBhg/X3mjJycZSmfK/Onz8fbdq0weTJk51d\nXoNc7YG1y5cvY8KECViyZAk6aGzrlS+++AI9evRATExM09cYsbW3cn2/VwghSkpKRHBwsNi5c6et\np7eL62tcuHChWLhwofV9fHy8yM3NVVFaHR07drR+XVVVJTp16qSwmoYlJCSInJwc6/vAwEBRXl6u\nsKK6Dhw4IHr06CH8/f2Fv7+/8PLyErfffrs4e/as6tLqtWLFCjF48GBx9epV1aXU8s0334j4+Hjr\n+wULFmj25qHZbBZ33323eP3111WXUq+5c+cKg8Eg/P39ha+vr7j55pvFlClTGv1v7BLSNWd3nD9/\nXkRGRop169bZemq7ub7G6huH165dE8ePHxcBAQGiqqpKYYW/iomJsYbfxo0bRWxsrOKK6vfee++J\n559/XgghxNGjR4Wfn5/iim5MyzcOs7KyRHh4uPjhhx9Ul1JHRUWFCAgIECdOnBDXrl3T7I3Dqqoq\nMWXKFDF79mzVpTRJTk6OY2d3rF27VhgMBnHTTTcJHx8fkZCQIIQQ4q9//ato3769iI6Otr5UfeM1\nVKMQQsyfP18EBgaKkJAQ62wKLSgoKBADBgwQUVFRYuDAgWLPnj2qS6qX2WwWDzzwgOjbt6/o16+f\n2Lx5s+qSbqh3796aDWmj0Sh69epl/Tczc+ZM1SXVkpmZKYKDg0VgYKBYsGCB6nLqtW3bNqHT6URU\nVJT1zzErK0t1WQ3Kyclp0uwO7sxCRKRhfCyciEjDGNJERBrGkCYi0jCGNBGRhjGkiYg0jCFNRKRh\n/w/p2ShZnc0kYQAAAABJRU5ErkJggg==\n"
},
{
"output_type": "pyout",
"prompt_number": 12,
"text": [
" Plot object containing: \n",
"[0]: cartesian line: atan((2*x + 1)**2/(4*x + 12)**2) for x over (-10.0, 3.141\n",
"\n",
" \n",
"592653589793)"
]
}
],
"prompt_number": 12
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"\u3053\u306e\u3088\u3046\u306bMatplotlib\u3092\u4f7f\u308f\u305a\u306bSympy\u5185\u90e8\u306e\u95a2\u6570\u3067\u3082\u7c21\u5358\u306a\u30b0\u30e9\u30d5\u4f5c\u6210\u306f\u53ef\u80fd\u3067\u3059\u3002\n",
"\n",
"## sympy\u3060\u3051\u3067\u5341\u5206\u304b\uff1f\u305d\u308c\u3068\u3082\uff1f ##\n",
"\u3053\u308c\u307e\u3067\u898b\u3066\u304d\u305f\u3088\u3046\u306b\u3001Sympy\u306f\u305d\u308c\u5358\u72ec\u3067\u76f8\u5f53\u306a\u4e8b\u304c\u3067\u304d\u308b\u30e2\u30b8\u30e5\u30fc\u30eb\u3067\u3042\u308b\u3053\u3068\u304c\u60f3\u50cf\u3067\u304d\u307e\u3059\u3002\n",
"\u79c1\u81ea\u8eab\u306f\u3001\u73fe\u6642\u70b9\u3067\u306f\u3001Sympy\u306e\u5168\u8c8c\u304c\u628a\u63e1\u3067\u304d\u3066\u3044\u307e\u305b\u3093\u3002\n",
"\uff08\u4ed6\u306e\u30e2\u30b8\u30e5\u30fc\u30eb\u3060\u3063\u3066\u305d\u3046\u3067\u3059\u304c\u3002\u3002\uff09\n",
"\n",
"\u3057\u304b\u3057\u3001\u8907\u96d1\u306a\u6570\u5f0f\u3092\u6271\u3046\u5fc5\u8981\u304c\u3042\u308b\u5834\u5408\u306b\u306fSympy\u3092\u5229\u7528\u3057\u3001\u305d\u308c\u4ee5\u5916\u306e\u6642\u306b\u306fNumpy,Scipy\u3092\u7528\u3044\u308b\u3068\u3059\u308b\u3068\u3001\n",
"\u79c1\u306e\u5834\u5408\u306b\u306fSympy\u306e\u4f7f\u7528\u983b\u5ea6\u306f\u4eca\u5f8c\u3082\u3042\u307e\u308a\u9ad8\u304f\u306a\u3089\u306a\u3044\u3068\u4e88\u60f3\u3057\u3066\u3044\u307e\u3059\u3002"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Sympy\u306b\u3088\u308b\u4ee3\u6570\u8a08\u7b97 ##\n",
"Sympy\u306e\u4ee3\u6570\u8a08\u7b97\u306e\u6a5f\u80fd\u306b\u3064\u3044\u3066\u306f\u6b86\u3069\u89e6\u308c\u307e\u305b\u3093\u3067\u3057\u305f\u3002\n",
"\n",
"[\u30c1\u30e5\u30fc\u30c8\u30ea\u30a2\u30eb](http://docs.sympy.org/0.7.2/tutorial.html#tutorial)\n",
"\u306b\u3042\u308b\u3082\u306e\u3092\u629c\u7c8b\u3057\u305f\u3060\u3051\u3067\u3082\u3001\u4ee5\u4e0b\u306e\u901a\u308a\u3067\u3059\u3002\n",
"\n",
"* \u30ab\u30c3\u30b3\u306e\u5c55\u958b\n",
"* \u901a\u5206\n",
"* \u90e8\u5206\u5206\u6570\u5c55\u958b\n",
"* \u5fae\u5206\uff08\u504f\u5fae\u5206\u3092\u542b\u3080\uff09\n",
"* \u7a4d\u5206\uff08\u6570\u5024\u7a4d\u5206\u3082\u3001\u4e0d\u5b9a\u7a4d\u5206\u3082\uff09\n",
"* \u30c6\u30fc\u30e9\u30fc\u5c55\u958b\n",
"* \u7d1a\u6570\u306e\u548c\n",
"* \u6975\u9650\n",
"* \u884c\u5217\u6f14\u7b97\n",
"* \u5404\u7a2e\u306e\u7279\u6b8a\u95a2\u6570\n",
"\n",
"\u3053\u308c\u4ee5\u5916\u306e\u91cd\u8981\u306a\u6a5f\u80fd\u3068\u3057\u3066\u306f\u3001\u8a08\u7b97\u5f0f\u306e\u5358\u7d14\u5316(\n",
"[simplify](http://docs.sympy.org/0.7.2/modules/simplify/simplify.html#simplify)\n",
")\u304c\u3042\u308a\u307e\u3059\u3002\n",
"\u3069\u306e\u3088\u3046\u306b\u5f0f\u3092\u5358\u7d14\u5316\u3059\u3079\u304d\u304b\u3001\u3068\u3044\u3046\u306e\u306f\u81ea\u660e\u3067\u306f\u306a\u3044\uff08\u4eba\u9593\u306e\u90fd\u5408\u306b\u3088\u308b\uff09\u306e\u3067\u3001\u8272\u3093\u306a\u8a2d\u5b9a\u3092\u5fc5\u8981\u3068\u3057\u307e\u3059\u3002\n",
"\u624b\u8a08\u7b97\u306e\u88dc\u52a9\u3068\u3044\u3046\u4f4d\u7f6e\u3065\u3051\u3067\u8003\u3048\u3066\u304a\u304f\u3068\u826f\u3044\u3068\u601d\u3044\u307e\u3059\u3002\n",
"\n",
"\u307e\u305f\u3001[physics](http://docs.sympy.org/0.7.2/modules/physics/index.html#module-sympy.physics)\n",
"\u3068\u3044\u3046\u30e2\u30b8\u30e5\u30fc\u30eb\u3082\u3042\u308a\u307e\u3059\u3002\n",
"\u3053\u3053\u306b\u306f\u91cf\u5b50\u529b\u5b66\u3084\u529b\u5b66\u306e\u8a08\u7b97\u306b\u4f7f\u3048\u308b\u30e2\u30b8\u30e5\u30fc\u30eb\u306e\u4ed6\u306b\u3001\u5358\u4f4d\u63db\u7b97\u306b\u5229\u7528\u3067\u304d\u308b\u3088\u3046\u306a\u30e2\u30b8\u30e5\u30fc\u30eb\u3082\u3042\u308a\u3001\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u4fbf\u5229\u306b\u4f7f\u3048\u305d\u3046\u3067\u3059\u3002"
]
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import sympy.physics.units as unit\n",
"unit.find_unit(\"avogad\")[1]\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "pyout",
"prompt_number": 13,
"text": [
"avogadro_number"
]
}
],
"prompt_number": 13
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"unit.avogadro_number\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"latex": [
"$$602214179000000000000000$$"
],
"output_type": "pyout",
"png": "iVBORw0KGgoAAAANSUhEUgAAAPIAAAASCAYAAABhLIMmAAAABHNCSVQICAgIfAhkiAAABORJREFU\neJzt2nmsXVMUx/HPq9bYKkVbUxQ1J2KIIYiIeUhEGhFiJiGGGILU0FChERGCEIJERVEpYmoi/hDx\nh/EPBBFTipjqiZpKn2rrj7Xve/sd596ec/uXZP+Sm3P32t+9770rd5299jqboqKi/70Gamw7Yjb+\nwlC6zsbvGbMTbsFXWIVNcRUWV+baF1djPWyFt3E9vu2T62h9vInde/66EV2Jcek7t53vGNyP9434\nZGXW/w7uTe+3ww1YhhXYADPxQ2XOpv4rXOHacMPaFV/jgNSeikUikDuaiG9wWma7Fh9i7cy2F17G\nRqk9Hq/hR0zrg+toHxE8q7r9iIq2wdLKb2gz3xWpr9vrmMRti58wIxt7mrgBjM1sTf1XuMK14YY1\nFp/gssy2NQZxaWabI4Is/3NOwnJckNkWYnrlM/YUf/75fXC7JHauWD2bBvIDiZ1dsTed7z7hh3EY\nk9kPxD1Z+1n8XGHWE6vzWZmtqf8KV7g23LDOwd9GVsZu+hQv1Ng/wCtZ+w+xuk+ucEvEytWWyzVX\ns0CegZPVB3LT+e6psY3HSyIlJ+6My8WWoKpP8GLWbuq/whWuMZevHqfgc/xSM6ijCdhBBF5V32Hv\nrL0IU8Q+MdeQWKnacm01Hscavar3o4trbLeJvfCfqT1J3DGX1bC/ihoAzf1XuMK14YaX6wEcJFaU\ng3GkCIRpuBHvJm6bdP2tZuKl2BDriCDcP32RvNCzhQjaVzNbU66trtG9uLUmOhBr4a3MNiiCet0a\nfgtsJnzd1H+FK1wbbqgTyJuIP+FkUfCaleyHiMLT/vgoDSRS8LqJidR8cWovrTCXiIrvdZVxTbg2\n2kOk7F/0Ob6X7saJFdsKPInjxU2xk6ZPFYFMrNpN/Ve4wrXhFndS67XSdVs8ksGvitTw1tReka51\n+8lxlbmqmi7S1FvwehemDddNY3A5bu9j7Op0mEj3F9X0zRTFrvNSe6z4He+l9grN/Ve4wrXhhvfI\nP6fr5+JZaa5vcIRYwgdrJuyos8f9vaZvHTwuKsizavrbcr10vihe1d3F1lQXigp3nQaxn1iF7xDF\ntQdFZrFM+Lip/wpXuDbc8B55eRq0pGbAkKjKThIp8yps3GXiX/w3kAfwsKjyXt/jizXlemmq2Brc\n1+f4XhqHo3BXD2aJqCnkmiIyi1Wa+29l4QrXghtVtX5TnBipqlO8GhR5+bviuWpV042kkbluwsdG\nB+cZa8D10uHYWTzT7bwWpr6TU3tG/dDVal/hvG6PxOq0mTiptiC1m/qvcIVrw40K5CdElSx/5DMg\nAuN5/JNsC0UKmR/v3D592FOVDztb3FVuqtgP6pNbneaJbcAJ2eui1Dc/tZ9pOWdHnUMrdY+YiIM0\nP4rA7ehMccx0bmZr6r/CFa4NN6wxeEOcS+7oJLEST8tsm4sl/fTMdqeoaudHxg4Vq9e8ymu+uGm0\n5XLNFynH+l36c+2Q2Dk9mCbzXZWYc7v0z8KXIr0nTqctNnLctaOm/itc4Rpz+bGvlThOFGoWiH3z\ngHj09GXGfS8eS80R56QniP3z0UYXmJ4WpfFT/Vc398FNxqPYErsl21fizOlDeKwydkM8Jw6cE2em\nDxUV+Gf7mO8zUcF/v+Z7ElXyiaLaPiG9P0HcHHM19V/hCteGKyoqKioqKioqKipac/0L7hNToHRj\ndKgAAAAASUVORK5CYII=\n",
"prompt_number": 14,
"text": [
"602214179000000000000000"
]
}
],
"prompt_number": 14
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"unit.find_unit(\"planc\")\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "pyout",
"prompt_number": 15,
"text": [
"[planck]"
]
}
],
"prompt_number": 15
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"unit.planck\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"latex": [
"$$\\frac{41412931}{62500000000000000000000000000000000000000} \\frac{kg m^{2}}{s}$$"
],
"output_type": "pyout",
"png": "iVBORw0KGgoAAAANSUhEUgAAAUYAAAAnCAYAAABzCOG7AAAABHNCSVQICAgIfAhkiAAABjZJREFU\neJzt3XusXFUVx/FPH7y9Iqil1AoFTAPII4GiXpWXUN8gYq2KCiqVWCAKigaBP4z+Y3yAMTEqf/iH\nQWLwEU2M+KyvKPWBEBExVkESrEqrprZgFeT6xzrjnJb77Nxp79z7+yaT2Wd+Z5+zJ9mzZq+1916H\nEEIIOzBvTzcghBB2I/NxGfZrjj+0B9sSQggzgnPw9Kb8JZw82knzd1tzQghhz3MUXteU/6BrJEMI\nYc6yD4aa8jexZLSTMmIMIcwUFuGzuLOP9/g3tuJUfB8bRztpYR8bEEIIU+FB/BAP9Pk+B+JMvL/P\n9wkhhGnhFjy/z/e4DHs1r7P7fK8wg/mox3e2IXwRh01TvWfjSrwP38ZpLe0k3IjL8XGsaGmn4Y1Y\ng5s8vqNOpp1hdjAfG7AAL8NHlNu7oNGfh+vV5MkH8CkcrPrdVzCMi3EFPoaz8Ho1MryqucYF+Cc2\n4x84rs/fKcxQTsdfcUbrszXKgI1g2TTUewI+2DpejYfxNOW2bMLxjbYI9+rGtzfjoqb86qZeJzg+\nmXaG2cMp+C7epPrOQhzRaMvxG92+8UllENeokd8GvKHR9lexxOHm+Nimbggoo/R2FWQ+YxR9LIMz\n1Xon4DG1FAKe2JyzGuc35XY8+8/q35/6xz6gKa/CI7qdf6J2htnFtWoy5Ad4yU7al3F16/jHyhMZ\nUkb0Ty1tGOtbxxfiC1NpSGalZzeX4NO7od5dqjPe2xx31oZtUK4Q7N06f3+c2JR/jYea8vlqhLhV\nmIusxCtVH7i5+eyZzfsL1WiSMobL1ez1VuUyr2td52x8p3V8QXO9gybbkBjG2cu5uFW5FP2uN4Kf\nNu/UP/v1uAPfa96f1WjDqt8d2Kp/kooBPYQbptjeMDs4QHkPt+Me5VVQs8fwFxWSoTyL9cpLoQxj\n2xC2DeNBKqb9Nbxlso2JYZydLMGT1Ghsd9Rrc7Hq1J1g96Oqc5+I12K7MoDt9WO/VIH2n+NHKmYZ\n5hYnqAXXj6rY9u1Yq1xoeAcuVRN1Fyl3u8Mz8K2mPA9PxU+a4+24TcUtv9631oeB4M14rxq5Xa2M\n0I1qpq/NzrG7Xa3X4eXKMMK+rXPayUqeojrrEjxH/Qg6Afajm2uvmuT9wtxkvZqoCaEn/mhqky+7\nUu90ZRQXN6/zdGcFH8Bzm/K1uu7yChUb2qc5fin+o0YAU2lnmN28G69pysfiV/rs7S6Y+JQwwCzF\ndSoGc4gKVP9Ore26VBmmQ5Uh+1kP9Y5ULvAq5UJfpdzmd6pY5Xy1TGeliiVdo+JDG/FfvECNHs/D\nu5RLbRLtDHODV6l+c5jySi7Hlj3aohBCCCGEEEIIIYQQQhiHzjKKkXHPCiGE/jFPbFAIIYQQQhgo\ndl7HeKF6itZytbD3NpVnbzVerBbn3t+84DPqSVvX4BVqa1dnj+MKtXtikVp3dA/+Pgmtl7rRokUb\nPK1XzlH7o49u7nPHNF0XtZewk1NvmdqBsNTYefaoLBiL1eLbNnsr49n5/BTdRbvjab3UjRYt2uBp\nvTKkdsJQe6Rvnabrohq+CYe3PjvS+Hn2KMM4Gitxd+t4HrapPbHjab3UjRYt2uBpvbKf2rp6l0pE\ncsg0XPP/yUOH1eb+w1Uq+5PxDZWifqw8e51GrVXbc85SqabuViPOv7XuM6KbRnzxONp9PdSNFi3a\n4Gn36Y3tOAYvUsmVr7RjQttdomMYO67xCD6Hr6qN/8epPHsd2nn2KCt9i3K7H1Qpgo5RRvbhUb7A\n0ASaHupGixZt8LReWKbSkx2qnvmyl2kaMXYyVHQ2ZP+ied+Gf6kN2x12zrMHn1dGEX6vJm2Ob67X\nTjVF5djbPIFmAj1atGizS+uFzfiwCu1douZEPtHjNdE1jHeq0WJ7lrr9nI6OgXyPShG1TGVD2aLc\nabrW/xH8Vg2fOyxUmXTvn0DTQ91o0aINntYL29Tk8E0qb+gN+rBQfJ3y06nZnS3KxR4rz95SO06+\nvE09oIb64ht1H3l5plrKM5HWS91o0aINnjYjaY8Q16ncd0eppJDXKYs8Vp69TWqUeC5OVTPYb1Vp\n6x9To9C1eLJa43iFCsCOp+mhbrRo0QZPCyGEEEIIIYQQQgghhBBCCCGEEOYY/wNFl2aBASp79gAA\nAABJRU5ErkJggg==\n",
"prompt_number": 16,
"text": [
" 2 \n",
" 41412931\u22c5kg\u22c5m \n",
"\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n",
"62500000000000000000000000000000000000000\u22c5s"
]
}
],
"prompt_number": 16
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"unit.planck.evalf()\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"latex": [
"$$6.62606896 \\times 10^{-34} \\frac{kg m^{2}}{s}$$"
],
"output_type": "pyout",
"png": "iVBORw0KGgoAAAANSUhEUgAAAMUAAAAjCAYAAAAkJc5vAAAABHNCSVQICAgIfAhkiAAAB8pJREFU\neJztnHmMFEUUh79ZFtdlXTlcWJYjct8awRNURIJERY0HSiCgYFRUPBEURRABAxg5EoMIElwFETEg\nCRrFEzDIpXigEkA8URGIqIiACusfvxqnpumG7mb2IFtfMtnpOl4dXa/q1auaBYfDkUaivCvgcJQx\nWcAgINc8P16OdXE4KgSXAw3N9wXA6d4EWWVaHYej/GkK9Dbft5BSEIej0pID5JvvS4B63gRupXBU\nFOoAzwOflHI5+4HdwPnAUuAnb4LsUq6AwxGW7cByYGsZlFUduBAYXQZlOSopnYF+wE3AHKBbQLqJ\nwHmHkTP/CPGZYhBQ1XwOqatzyToywU7gPuA54FrztxCZKUkuQIO+FzJbvGQBG4FWwMVoJh8I1AAO\nmDTnAtcAa4E2powNRvYEE5YPNAIWA3WBlsAfwBNGRh/gaeBvoAoyoz6P23CHI4h2QJ753hP4h9Rm\nFmSu3IWUoUuAjDOBd4D+QH1k2je24lsAX1pypwFD0Wy/GehrwquhfUNH89zG5HM4yo25wHBP2FDk\n9VlKsFIMR5veZcAlPvELgWHW8wq0QtQHfrTCOwKrrOfrgZdD1dzgt9FuAYwC9iKN22ued/ukjZP3\nLNS4XKABsAYYSXrDoshrCYwDvgNKgAL0En6JKa8B8KCRlYNmnvHAeo+sJsAjwD60vOcBDwDbYpab\nSaqhgXFqQHyUPgtLB6ArsAeYbIVfAbyO2n04LgKuQnVfCNQE2gJfmPjups6g1aIF8D5aId615HQD\n3rae+wDTjbxdURqUpA3wPdDJPNcFvkEvMBN5OwBvIjsR4ATkcdiO7MCo8qojb0VfK+whZCMeF0Pe\nScAioMgKa4xs3eaesJ3A1VZYX+BTDp1ojqZP43AmsrlLAuLD9tntyA4P+vQKkH8L8BF6t/XQTJ1k\nKf4rRR7qz2zUP0lz5w4rzVekxsgAtGcA7V9usNIts8qoaeRWRXueyGSjl3+PFdYQ2AHcnaG8rwHN\nPHnboxc4L4a8x5BC2QOxFrJpb4shbyhwp7dxaNBMsJ4XAb+Sfs6Ti1aN/jHKzQStUf8Wo1UiSCnC\n9llYzkErTNL+b2XK7okG74PIMhiGzKMZQA+PjI7AC+Z7Apht6mIfrPVA95T6IeUaYsJXkJrEEkih\nksqdi5TnZtQ/kbkR7chrHCnhUeT9E82adTzhu5BGR5W3idSMYbOe9CU1rLxpwCs+4YOBSeb7cWgA\nrfFJtxF4NUa5maaYYKUI22dhOcPkyzHPl6I2eyc/gG8J3lNEYRVaEUudt4i/Sw+bdz2yLRt7wrch\nWzSKvHz04qf6xC0Bfo9Rv4FG5ly09IJe9jpS9nldk2a5T/41aBaOWq4fCeSaPBwF+FxoI1gpovRZ\nFPoi8+RedMmuuye+AZrl9yGFvCyi/KGkzLU2wGeUwW2MBNr8LUMHMWOBKchMaJ/BvHloUNnUQy/q\nvYjy2pl84ziUhaQ2ylHql0PKHv8Zvew5pHtDqiAF9lsptpq82RHLDeIl0j0uNoXAauRn91KMv1KE\n7bOKxkRgBHofk5GSlToFqEM2ALda4V2Qh6RtKeUFeXYOkNqIhpXXyaR71EfmbBNXGKN++cAbJk8J\nMocKPWlmIXPPPvysCxw0eerEKNePKqYtIzzhRUgp/RQCgpUibJ85UEeUoOUt1xO3lXQ7OZN5m6F9\nxtgY8s426Ub5yJ1n4urFqN8w4Bm0sfua1KpxipWmNrLNB5rnbNOGdSb9STHKDSILKeEY81wfrWYd\nA3MEK0XYPnMgl1UJ/sfdq9A+IGhZjZs3B812kzzhYeU1IfgFLzZx+RHrdy/y4CTJQ8v1QeRutamJ\nzikmIYU4GfgQmUyJiOUeiQS6mvAUUogjbTKL8VeKsH1WqUm65f5BbkK/w439yONSC82YXuLkTQDP\nIjNlZEx5v6CXWNMnXR7wG6nDsTDytpm62CbJHqQoG5FnqjUyhzDyvGZIIfCBqdfR9KmXEnQ2sBod\nhK0NkcePKH1WabF38KuQHewlB73EHYeREzXvGDS4bIWwD3nCyNsDfIz/L6eakX4vP4y82sh1utkn\n3XR0qex4n7gktdEG0L5ScDR9atMKXabrgkyvqcS7zBmlz8qSkgr2+Z/ewF+k278JNHvM9zSiuSdd\nlLwD8L/HPiOGvNHoeog9QJqihg2KKC+BVgu/q8v5aDAmzZ17kOvV9oIMMWlsxYnSL0G0QwPZPnwa\nA8wk2C1ZTPA5Rdg+c6AOXknqpBDgOjSbNbLCOiMbe0mMvF2R12aO5zMPeDGGvCI0wPpZYVPQfRn7\nykJYeb1MXjvsRDT797TCHkYHUUn3cntkmnQinbDlBtEeKURzn7gR6IpDFZ+45Ka5mk9c2D6rtHiX\n4Fpo45iHbOIEGgBbrDQt0cHVAnQ/JkreXQSf7o4l3fUYRh7AaejqwiY0o9dCVyh+iNE2kNIPRuZN\nCRp0U0n/DUAumnELTJnVkSm40qddYcv1koUuvPVDXjA/7kenx1OQG3g28k4l3b070UZ/JqlrFBC+\nz44Fkv+dYzdybswq3+o4HOVLPjrdBu3pXi/HujgcFYJcZMauR945d/DoqPQkkGJciS4ljs+EUPff\nPBzHKo3Q7zaK0H2yqmRopfDzXDgcxwL/mk9D9JuOAuDJcq2Rw+FwOBwOh8PhcDgcjkP5D/1fg3do\nM1lWAAAAAElFTkSuQmCC\n",
"prompt_number": 17,
"text": [
" 2\n",
"6.62606896e-34\u22c5kg\u22c5m \n",
"\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n",
" s "
]
}
],
"prompt_number": 17
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"unit.find_unit(unit.pressure)\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "pyout",
"prompt_number": 18,
"text": [
"[Pa, pa, atm, bar, kPa, psi, bars, mmHg, pascal, pascals, pressure, atmosphere\n",
", atmospheres]"
]
}
],
"prompt_number": 18
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"3.5*unit.psi\n"
],
"language": "python",
"metadata": {},
"outputs": [
{
"latex": [
"$$24131.6505260893 \\frac{kg}{m s^{2}}$$"
],
"output_type": "pyout",
"png": "iVBORw0KGgoAAAANSUhEUgAAALgAAAAlCAYAAADxyyNaAAAABHNCSVQICAgIfAhkiAAACHBJREFU\neJztnHmQVMUdxz+zLq67y66KrCyKxRLBY40pQEmEKInxiDGmKuWttRJMaTyjmIBnUDwwxErhVSkq\nalkoOagYFcv7jJhS8T4wWoBHUGISwAtExYP1j293Tb83/WZ6Zh4LRfWnamrndf/m9359vO5f/7rf\nQiQSiUT6hG2Bm4EX81LYmJeiSCQHlgOPAcvyUtiQl6JIJCcOBO7JS1khL0WRSA40AIuAXYCDgH2B\nk4CtgK+MzHeBw4BngG6gAzi5zy2NRGpgDPAwMBHYHrnQw5z8nYBXgTZzPQs4qw/ti0Tq4gLgXWA+\n8CNP/m3Auc7148DocgrTLsq3jYJmYAjwNHAh8J8Khk0G+gG/zchvARYA38rI3wtNRZ+ae7cAlwMv\n1agvi52AaeY+a83facBqR2YRcCnwD2ANsCeqk18CrzlyO6PyLgV6gYHAFOD/qXuG6qum7kPKEWrf\nEOA8I9OE6nYGsDAl9w3gIuAz5C60AucA/0vJVduWLo8anS2oM28N7Ab8y+R/jNyWZ9Ao/iYwCFgX\noJvRwAPI3wHoj1a0y4GuMr8bihpuWkb+GGNQb0b+KOAuYAsnbRawChhZg74suoG3gXHmuhN4i1K7\ne1Ofz4HTUjJbopV+j5N2PvAKsHkN+qqp+5ByhNq3DTAPGOykDUMP5YhU2krgUCetB3VaNxJXbVu6\ntJp7NJoyvWrST3dkXqdYH8cDd1bQmeBuYHgqbRRqlLllfnedkZmWSt/V6JyNRtusDnmlyTvKSTvE\npF1Tgz4fjajRJjlpOwArgDNTsktRmW4HrjD3TTMddT63cQcAXwCn1KAvtO5DyxFq3xQ0k6Q5H/id\ncz0PeJ9k1K0ZjeYTnbTQtvQxFviz+V4A5hhbt3Nkfozq8Dg02k+uoDPBx2hk2DaV/gF6snwcChyN\nv4O7zCa7Q04APgIOcNKszhk16PPxczRyblVJEFVcJRbjHz0WAo/UoC+07kPLEWrfLPTgpfkVMNN8\n3xw9GE975BahEdtSS1vWygI0m5fFfSLfQv5Ma0pmLXpa0/QHDqb86B7CzWhKfdBJG438vL/Wqdty\nDJrePsxBVxuavt/25L0L7FGDztC6DylHNfa9CPwU+Avyd0F+eA8aREAjfyMardN8hNYOlvXZllMo\nzgzdyE9/rtKP3ClsL1Q57qJhO1Txj3p+ex7Zi8p6GAb8DDiVsIVJJQrA3mgEGo82EvojX+5i4IWU\nfBMwFfmnXwA7osXfYpM/1Pxd5bnXGqDd6FgbqA/C6j60HNXYNxs4AT04+6JOdBCKZrxsfrMC+ISk\nX+3a2IH60Zee/DzbstPY0IMe0oMJXFyWYwZ6+sal0keiDm6px0WxHIKiDS8jH7DcLmuIPstAI/sa\nyQ2B76Oow24p+deRX2vpQZGHTnM9zui72HOvOSZvUBX6skjXfWg5qrWvDbiP4iL4rlQ+wI3IVXKj\nbp2og/VS6lpV05YbjOHIN7wsld4A3ERyNZ5HB7c0oiluAWrUevUNMrKfUepqLSPpQ0JpY2yGRjC7\nSPoO2eWda/LchVElfT58dR9ajmrtOxe4Hi3g3jT5/wV2d2Q60IxzkrluNLY9b+S3yShHSFtuEJrQ\nVDjTk3cKms5c8uzgoCm4F/h7Dvr6GdlXPHkL0FTdVEHHUvMBxYOzynunyWvz5GXpS5NV96HlqMa+\ns1AEx9KKIiHrKHUptkZx8Jmocw8FnkXx7nJHPiq1ZZ9TQIuOSzx5ncC1nvR6OvgulG7YtBv5dcjP\nrEafj+XAPz3p840eGwd+LENuGcVFVqux62qP3CMo8mEJ0edSru4hrByh9hXM92965E42+nwhTZd3\n0Na6pZa2XK/4fKNLkZ93oZM2wfzdHxVinvOxI8DR5trdDKhEO1ocPY8WXxZ7sKaApvR6yZoi7WJr\nhbkehb8RBlIccdcgm3fwyA0neZY5RJ9LubqHsHKE2teBwo1LPHJ/RItU38LS0oF2QW8x16Ftmd74\nWt+fBMfjHz2uK1PQLmofwbdAkYUlKBxl2dPIP1WlPssISsNrn6TSCijc9jcn7RZK/Um74TLVSbsE\nbaG7U/OORs7dpQzVB2F1H1qOEPsKKGqzt+eebWiWsa7bJDR7DHFkJhsZ+xDU2pb10IA2qs42n7L8\nAK2U/5T6zKV8DHMEKsD0MjJ2cdPiybvcGOk2xhwUGcg6SFNO33g0Hd7vpDUAT5Lc+ToSjXhdTtoY\n5CtavQUU232CpJ8+GHWq45y0q9CZCXfxHaovtO5DyxFq31Emzf1tO3owD3fSfgP8m2LkZxSKBKWj\na7W0ZT38hOJMdSuePQjXkA/I3iG7jNIRpx24Ax3qGYymyBfQFu88FDqag4492hDWSrRIuoHitixo\nu/eHKJbaibaFp5KMFYfq2xn5vrei+KtlAFogtaKRpoAa7o1UufYBzkCLp3a0+zedUp95pElfjEa8\nAWi7/J0a9FVT96HlCLVvPNq5XIsGjc2AP5Dc+2hGs8JAo2tL5EY96bF3IpXbMi8moQf2CvN5AvW9\nSGSToIliROh+kuFPIJ8FXCSyofgKnc3ZB7ml925YcyKR/LHuUiSySXIa2gTrh8LYkcgmw7EoXr+S\njE2r+FZ9pK8ZgE4Yfg9F3LrRQrELHSXoRJGwVcDvU7+1YcHVaMS+sU8sjkSq4ATUOZdQfKWuBYUp\nx5rrboqvrFnaKB7h7SAuKCMbKW1oL8N9mXosOoZgmUDxCIClGW02LUQje/pIr5eN8pxuZJNmNbAf\nyVfn9gcecq6PJfmWEWhjbFe0aTSa+P9QIhsxNyE/3DIfvbgB6tQrkRvza5PWBbxH8ZjBESTfto9E\nNioep3hEuYD8bdt5m9Fi80SKx3X7oxczeoBfoNE7BkgikUgkEolEIpFIJNK3fA1Ho5UZrsvQzwAA\nAABJRU5ErkJggg==\n",
"prompt_number": 19,
"text": [
"24131.6505260893\u22c5kg\n",
"\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\u2500\n",
" 2 \n",
" m\u22c5s "
]
}
],
"prompt_number": 19
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"\u3057\u304b\u3057\u3001Scipy\u306b\u3082\u540c\u69d8\u306b\u4fbf\u5229\u306a[\u30e2\u30b8\u30e5\u30fc\u30eb](http://docs.scipy.org/doc/scipy-0.7.x/reference/constants.html#module-scipy.constants)\u304c\u3042\u308a\u307e\u3059\u3002"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## \u6539\u3081\u3066\u3001Sympy\u306e\u4f7f\u3044\u65b9\u306b\u3064\u3044\u3066\u306e\u79c1\u898b ##\n",
"Sympy\u306f\u3053\u308c\u3060\u3051\u9ad8\u6a5f\u80fd\u3067\u3059\u3002\n",
"\u307e\u305f\u3001Sympy\u3067\u5f97\u305f\u6570\u5f0f\u30aa\u30d6\u30b8\u30a7\u30af\u30c8\u3092Numpy\u7b49\u3067\u76f4\u63a5\u5229\u7528\u3059\u308b\u306f\u3053\u3068\u306f\u3067\u304d\u307e\u305b\u3093\u306e\u3067\u3001\u30b3\u30d4\u30fc\uff06\u30da\u30fc\u30b9\u30c8\u3068\u591a\u5c11\u306e\u66f8\u304d\u63db\u3048\u3092\u5fc5\u8981\u3068\u3057\u307e\u3059\u3002\n",
"\n",
"\u3067\u3059\u304b\u3089\u3001\u8907\u96d1\u306a\u6570\u5f0f\u3092\u53d6\u308a\u6271\u3046\u5fc5\u8981\u304c\u751f\u3058\u305f\u5834\u5408\u306b\u306f\u3001\u4e00\u3064\u306eNotebook\u5185\u3067\u306f\u307b\u307cSympy\u306b\u5fb9\u3057\u3001\n",
"\u5b9f\u9a13\u30c7\u30fc\u30bf\u3092\u53d6\u308a\u6271\u3046\u5834\u5408\u306b\u306f\u9006\u306bSympy\u306f\u4f7f\u308f\u306a\u3044\u3068\u3044\u3063\u305f\u4f7f\u3044\u5206\u3051\u3092\u3057\u305f\u65b9\u304c\u826f\u3044\u3088\u3046\u306b\u601d\u3048\u307e\u3059\u3002\n",
"\n",
"\u5168\u3066\u3092Sympy\u3067\u884c\u3046\u53ef\u80fd\u6027\u3001\u307e\u305f\u3001\n",
"[Sage](http://ja.wikipedia.org/wiki/Sage_(%E6%95%B0%E5%BC%8F%E5%87%A6%E7%90%86%E3%82%B7%E3%82%B9%E3%83%86%E3%83%A0))\n",
"\u306a\u3069\u306e\u5de8\u5927\u306a\u30bd\u30d5\u30c8\u30a6\u30a7\u30a2\u3092\u4f7f\u3046\u3053\u3068\u306a\u3069\u3082\u5c06\u6765\u7684\u306b\u306f\u8003\u3048\u3066\u3082\u3088\u3044\u306e\u304b\u3082\u77e5\u308c\u307e\u305b\u3093\u304c\u3001\n",
"Windows\u74b0\u5883\u306b\u304a\u3044\u3066\u3001Numpy(+Scipy)+Matplotlib\u3092\u30e1\u30a4\u30f3\u3068\u3059\u308b\u5b9f\u9a13\u30c7\u30fc\u30bf\u306e\u89e3\u6790\u3068Sympy\u3092\u30e1\u30a4\u30f3\u3068\u3059\u308b\u6570\u5f0f\u51e6\u7406\u3092\u3042\u308b\u7a0b\u5ea6\u5206\u3051\u3066\u884c\u3046\u4e8b\u306f\u3001\n",
"\u591a\u304f\u306e\u5b9f\u9a13\u5c4b\u306b\u3068\u3063\u3066\u306f\u8efd\u5feb\u3067\u5341\u5206\u306a\u3084\u308a\u65b9\u3060\u308d\u3046\u3068\u601d\u3063\u3066\u3044\u307e\u3059\u3002\n",
"\u305d\u3057\u3066\u3001\u305d\u306e\u4e21\u65b9\u304cipython notebook\u3068\u3044\u3046\u4e00\u3064\u306e\u74b0\u5883\u3067\u5b9f\u884c\u3067\u304d\u308b\u306a\u3093\u3066\u3001\u3068\u3066\u3082\u4fbf\u5229\u306a\u3053\u3068\u3060\u3068\u601d\u3044\u307e\u3059\u3002"
]
}
],
"metadata": {}
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment