{"id":6200,"date":"2023-04-20T12:06:16","date_gmt":"2023-04-20T03:06:16","guid":{"rendered":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/?page_id=6200"},"modified":"2024-03-15T14:51:50","modified_gmt":"2024-03-15T05:51:50","slug":"sympy-%e3%81%a7%e5%a4%9a%e9%87%8d%e7%a9%8d%e5%88%86","status":"publish","type":"page","link":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/%e7%90%86%e5%b7%a5%e7%b3%bb%e3%81%ae%e6%95%b0%e5%ad%a6c\/sympy-%e3%81%a7%e7%90%86%e5%b7%a5%e7%b3%bb%e3%81%ae%e6%95%b0%e5%ad%a6c\/sympy-%e3%81%a7%e5%a4%9a%e9%87%8d%e7%a9%8d%e5%88%86\/","title":{"rendered":"SymPy \u3067\u591a\u91cd\u7a4d\u5206"},"content":{"rendered":"<p><!--more--><\/p>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<p>\u5fc5\u8981\u306a\u30d1\u30c3\u30b1\u30fc\u30b8\u3092 import \u3057\u307e\u3059\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[1]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"kn\">from <span class=\"nn\">sympy.abc<\/span> import <span class=\"o\">*<\/span> \r\nfrom<\/span> <span class=\"nn\">sympy<\/span> <span class=\"kn\">import<\/span> <span class=\"o\">*<\/span>\r\n\r\n<span class=\"c1\"># SymPy Plotting Backends (SPB): \u30b0\u30e9\u30d5\u3092\u63cf\u304f\u969b\u306b\u5229\u7528<\/span>\r\n<span class=\"kn\">from<\/span> <span class=\"nn\">spb<\/span> <span class=\"kn\">import<\/span> <span class=\"o\">*<\/span>\r\n\r\n<span class=\"c1\"># \u30b0\u30e9\u30d5\u3092 SVG \u3067 Notebook \u306b\u30a4\u30f3\u30e9\u30a4\u30f3\u8868\u793a<\/span>\r\n<span class=\"o\">%<\/span><span class=\"k\">config<\/span> InlineBackend.figure_formats = ['svg']\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<p>\u591a\u5909\u6570\u95a2\u6570\u306e\u7a4d\u5206\u3067\u3042\u308b\u591a\u91cd\u7a4d\u5206\u306e\u3046\u3061\uff0c\u3082\u3063\u3068\u3082\u7c21\u5358\u306a2\u5909\u6570\u95a2\u6570\u306e\u7a4d\u5206\u3064\u307e\u308a2\u91cd\u7a4d\u5206\u306b\u3064\u3044\u3066\u307e\u3068\u3081\u308b\u3002<\/p>\n<h3 id=\"2\u91cd\u7a4d\u5206\u306f\u7d2f\u6b21\u7a4d\u5206\u3067\">2\u91cd\u7a4d\u5206\u306f\u7d2f\u6b21\u7a4d\u5206\u3067<\/h3>\n<p>2\u91cd\u7a4d\u5206\u306e\u8a08\u7b97\u306f\uff0c\u57fa\u672c\u7684\u306b\u7d2f\u6b21\u7a4d\u5206\u3067\u3084\u308a\u307e\u3059\u3002<\/p>\n<p>$\\displaystyle \\iint_D f(x, y)\\, dx dy$ \u306e\u9818\u57df $D$ \u304c\u4f55\u3068\u4f55\u3067\u56f2\u307e\u308c\u3066\u3044\u308b\u304b\u3092\u660e\u3089\u304b\u306b\u3057\u3066\u7a4d\u5206\u3059\u308b\u3002<\/p>\n<h4 id=\"\u7d2f\u6b21\u7a4d\u5206-1\">\u7d2f\u6b21\u7a4d\u5206 1<\/h4>\n<p>\u9818\u57df $D$ \u304c $x_1 \\leq x \\leq x_2,\\ y_1 \\leq y \\leq y_2$ \u3067\u3042\u308b\u5834\u5408\uff1a<\/p>\n<p>$$\\iint_D f(x, y) \\,dx\\,dy = \\int_{x_1}^{x_2} \\left\\{\\int_{y_1}^{y_2} f(x, y)\\,dy \\right\\} dx$$<\/p>\n<p>\u307e\u305f\u306f<\/p>\n<p>$$\\iint_D f(x, y) \\,dx\\,dy = \\int_{y_1}^{y_2} \\left\\{\\int_{x_1}^{x_2} f(x, y)\\,dx \\right\\} dy$$<\/p>\n<h4 id=\"\u7d2f\u6b21\u7a4d\u5206-2\">\u7d2f\u6b21\u7a4d\u5206 2<\/h4>\n<p>\u9818\u57df $D$ \u304c $x_1 \\leq x \\leq x_1,\\ \\phi_1(x) \\leq y \\leq \\phi_2(x)$ \u3067\u3042\u308b\u5834\u5408\uff1a<\/p>\n<p>$$ \\iint_D f(x, y) \\,dx\\,dy = \\int_{x_1}^{x_2} \\left\\{\\int_{\\phi_1(x)}^{\\phi_2(x)} f(x, y)\\,dy \\right\\} dx$$<\/p>\n<h4 id=\"\u7d2f\u6b21\u7a4d\u5206-3\">\u7d2f\u6b21\u7a4d\u5206 3<\/h4>\n<p>\u9818\u57df $D$ \u304c $\\psi_1(y)\\leq x \\leq \\psi_2(y),\\ y_1 \\leq y \\leq y_2$ \u3067\u3042\u308b\u5834\u5408\uff1a<\/p>\n<p>$$ \\iint_D f(x, y) \\,dx\\,dy = \\int_{y_1}^{y_2} \\left\\{\\int_{\\psi_1(y)}^{\\psi_2(y)} f(x, y)\\,dx \\right\\} dy$$<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h4 id=\"\u4f8b\u984c\">\u4f8b\u984c<\/h4>\n<p>$$I = \\iint_D y \\,dx dy, \\quad D: 0\\leq x\\leq 2, \\ 0 \\leq y \\leq 2$$<\/p>\n<p>\u3053\u308c\u306f\uff0c\u30b1\u30fc\u30b9 1\u3002\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u7d2f\u6b21\u7a4d\u5206\u306b\u3057\u3066&#8230;<br \/>\n\\begin{eqnarray}<br \/>\nf(x, y) &amp;=&amp; y \\\\<br \/>\nx_1 &amp;=&amp; 0\\\\<br \/>\nx_2 &amp;=&amp; 2\\\\<br \/>\ny_1 &amp;=&amp; 0\\\\<br \/>\ny_2 &amp;=&amp; 2<br \/>\n\\end{eqnarray}<\/p>\n<p>SymPy \u3067\u306f2\u91cd\u7a4d\u5206\u306f\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u3002<\/p>\n<p>$\\displaystyle \\int_{y_1}^{y_2} \\left\\{\\int_{x_1}^{x_2} f(x, y) dx \\right\\} dy =<br \/>\n$<\/p>\n<p><code>integrate(f(x, y), (x, x1, x2), (y, y1, y2))<\/code><\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[2]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">f<\/span> <span class=\"o\">=<\/span> <span class=\"n\">y<\/span>\r\n<span class=\"n\">x1<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"n\">x2<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">2<\/span>\r\n<span class=\"n\">y1<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"n\">y2<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">2<\/span>\r\n\r\n<span class=\"c1\"># \u624b\u3063\u53d6\u308a\u65e9\u304f\u7b54\u3048\u3060\u3051\u3092\u8868\u793a\u3055\u305b\u305f\u3044\u306a\u3089<\/span>\r\n<span class=\"n\">integrate<\/span><span class=\"p\">(<\/span><span class=\"n\">f<\/span><span class=\"p\">,<\/span> <span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">x1<\/span><span class=\"p\">,<\/span> <span class=\"n\">x2<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">y<\/span><span class=\"p\">,<\/span> <span class=\"n\">y1<\/span><span class=\"p\">,<\/span> <span class=\"n\">y2<\/span><span class=\"p\">))<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[2]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle 4$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[3]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"c1\"># \u4f55\u3092\u8a08\u7b97\u3059\u308b\u304b\u8868\u793a\u3055\u305b\u3066\u304b\u3089\u7b54\u3048\u3092\u8868\u793a\u3055\u305b\u308b\u306a\u3089<\/span>\r\n\r\n<span class=\"n\">Eq<\/span><span class=\"p\">(<\/span><span class=\"n\">Integral<\/span><span class=\"p\">(<\/span><span class=\"n\">y<\/span><span class=\"p\">,<\/span> <span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">x1<\/span><span class=\"p\">,<\/span> <span class=\"n\">x2<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">y<\/span><span class=\"p\">,<\/span> <span class=\"n\">y1<\/span><span class=\"p\">,<\/span> <span class=\"n\">y2<\/span><span class=\"p\">)),<\/span> \r\n   <span class=\"n\">Integral<\/span><span class=\"p\">(<\/span><span class=\"n\">y<\/span><span class=\"p\">,<\/span> <span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">x1<\/span><span class=\"p\">,<\/span> <span class=\"n\">x2<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">y<\/span><span class=\"p\">,<\/span> <span class=\"n\">y1<\/span><span class=\"p\">,<\/span> <span class=\"n\">y2<\/span><span class=\"p\">))<\/span><span class=\"o\">.<\/span><span class=\"n\">doit<\/span><span class=\"p\">())<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[3]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\int\\limits_{0}^{2}\\int\\limits_{0}^{2} y\\, dx\\, dy = 4$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h4 id=\"\u53c2\u8003\uff1a\u9818\u57df-$D$-\u306e\u5857\u308a\u3064\u3076\u3057\">\u53c2\u8003\uff1a\u9818\u57df $D$ \u306e\u5857\u308a\u3064\u3076\u3057<\/h4>\n<p>\u3064\u3044\u3067\u306b\uff0c\u9818\u57df $D$ \u306e\u90e8\u5206\u3092\u5857\u308a\u3064\u3076\u3057\u3066\u63cf\u3044\u3066\u307f\u307e\u3059\u3002<\/p>\n<p>$x_1 \\leq x \\leq x_2$ \u306e\u7bc4\u56f2\u3067 $y = y1$ \u3068 $y = y2$ \u306b\u56f2\u307e\u308c\u305f\u9818\u57df\u3092\u5857\u308a\u3064\u3076\u3059\u65b9\u6cd5\u306f\uff0cSymPy Plotting Backends \u3067\u306f\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u3057\u307e\u3059\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[4]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"c1\"># NumPy \u306b\u983c\u3089\u306a\u3044\u4f8b<\/span>\r\n<span class=\"c1\"># x1 ~ x2 \u3092 100 \u7b49\u5206\u3057\u305f\u6570\u5024\u30ea\u30b9\u30c8\u3092\u4f5c\u6210<\/span>\r\n<span class=\"n\">x_list<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[(<\/span><span class=\"n\">x1<\/span> <span class=\"o\">+<\/span> <span class=\"p\">(<\/span><span class=\"n\">x2<\/span><span class=\"o\">-<\/span><span class=\"n\">x1<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"n\">i<\/span><span class=\"o\">\/<\/span><span class=\"mi\">100<\/span><span class=\"p\">)<\/span> <span class=\"k\">for<\/span> <span class=\"n\">i<\/span> <span class=\"ow\">in<\/span> <span class=\"nb\">range<\/span><span class=\"p\">(<\/span><span class=\"mi\">101<\/span><span class=\"p\">)]<\/span>\r\n<span class=\"c1\"># x_array \u306b\u5bfe\u5fdc\u3057\u305f y1 \u306e\u6570\u5024\u30ea\u30b9\u30c8\u3092\u4f5c\u6210<\/span>\r\n<span class=\"n\">y1_list<\/span> <span class=\"o\">=<\/span> <span class=\"n\">lambdify<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">y1<\/span><span class=\"p\">)(<\/span><span class=\"n\">x_list<\/span><span class=\"p\">)<\/span>\r\n<span class=\"c1\"># x_array \u306b\u5bfe\u5fdc\u3057\u305f y2 \u306e\u6570\u5024\u30ea\u30b9\u30c8\u3092\u4f5c\u6210<\/span>\r\n<span class=\"n\">y2_list<\/span> <span class=\"o\">=<\/span> <span class=\"n\">lambdify<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">y2<\/span><span class=\"p\">)(<\/span><span class=\"n\">x_list<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plot<\/span><span class=\"p\">(<\/span><span class=\"n\">xlim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mf\">0.5<\/span><span class=\"p\">,<\/span> <span class=\"mf\">2.5<\/span><span class=\"p\">),<\/span> \r\n     <span class=\"n\">ylim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mf\">0.5<\/span><span class=\"p\">,<\/span> <span class=\"mf\">2.5<\/span><span class=\"p\">),<\/span> \r\n     <span class=\"n\">aspect<\/span> <span class=\"o\">=<\/span> <span class=\"s2\">\"equal\"<\/span><span class=\"p\">,<\/span>\r\n     <span class=\"n\">title<\/span> <span class=\"o\">=<\/span> <span class=\"s2\">\"\u9818\u57df $D$\"<\/span><span class=\"p\">,<\/span>\r\n     <span class=\"n\">xlabel<\/span> <span class=\"o\">=<\/span> <span class=\"s2\">\"x\"<\/span><span class=\"p\">,<\/span> <span class=\"n\">ylabel<\/span> <span class=\"o\">=<\/span> <span class=\"s2\">\"y\"<\/span><span class=\"p\">,<\/span>\r\n     <span class=\"n\">fill<\/span><span class=\"o\">=<\/span><span class=\"p\">{<\/span><span class=\"s1\">'x'<\/span><span class=\"p\">:<\/span> <span class=\"n\">x_list<\/span><span class=\"p\">,<\/span>\r\n           <span class=\"s1\">'y1'<\/span><span class=\"p\">:<\/span><span class=\"n\">y1_list<\/span><span class=\"p\">,<\/span>\r\n           <span class=\"s1\">'y2'<\/span><span class=\"p\">:<\/span><span class=\"n\">y2_list<\/span><span class=\"p\">,<\/span>\r\n           <span class=\"s1\">'color'<\/span><span class=\"p\">:<\/span><span class=\"s1\">'yellow'<\/span><span class=\"p\">})<\/span><span class=\"p\">;<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_svg output_subarea \">\n<p><!--?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?--><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-6201\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/pmathc301.svg\" alt=\"\" width=\"640\" height=\"427\" \/><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<p>\u3064\u3044\u3067\u306b\uff0c3\u6b21\u5143\u7684\u306a\u30b0\u30e9\u30d5\u3092\u63cf\u3044\u3066\u307f\u307e\u3059\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[5]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">plot3d<\/span><span class=\"p\">((<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">x1<\/span><span class=\"p\">,<\/span> <span class=\"n\">x2<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">y<\/span><span class=\"p\">,<\/span> <span class=\"n\">y1<\/span><span class=\"p\">,<\/span> <span class=\"n\">y2<\/span><span class=\"p\">)),<\/span> \r\n       <span class=\"p\">(<\/span><span class=\"n\">f<\/span><span class=\"p\">,<\/span> <span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">x1<\/span><span class=\"p\">,<\/span> <span class=\"n\">x2<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">y<\/span><span class=\"p\">,<\/span> <span class=\"n\">y1<\/span><span class=\"p\">,<\/span> <span class=\"n\">y2<\/span><span class=\"p\">)),<\/span>\r\n       <span class=\"n\">xlim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mf\">0.5<\/span><span class=\"p\">,<\/span> <span class=\"mf\">2.5<\/span><span class=\"p\">),<\/span> <span class=\"n\">ylim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mf\">0.5<\/span><span class=\"p\">,<\/span> <span class=\"mf\">2.5<\/span><span class=\"p\">),<\/span> \r\n       <span class=\"n\">zlim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mi\">3<\/span><span class=\"p\">),<\/span> <span class=\"n\">n<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">10<\/span><span class=\"p\">,<\/span>\r\n       <span class=\"n\">rendering_kw<\/span><span class=\"o\">=<\/span><span class=\"p\">[{<\/span><span class=\"s2\">\"color\"<\/span><span class=\"p\">:<\/span><span class=\"s2\">\"yellow\"<\/span><span class=\"p\">},<\/span>\r\n                     <span class=\"p\">{<\/span><span class=\"s2\">\"color\"<\/span><span class=\"p\">:<\/span><span class=\"s2\">\"lightblue\"<\/span><span class=\"p\">}])<\/span><span class=\"p\">;<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_svg output_subarea \">\n<p><!--?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?--><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-6202\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/pmathc302.svg\" alt=\"\" width=\"640\" height=\"427\" \/><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<p>2\u91cd\u7a4d\u5206\u306e\u7d50\u679c\u306f\uff0c\u9ec4\u8272\u3067\u5857\u308a\u3064\u3076\u3055\u308c\u305f\u9818\u57df $D$ \u3068\u8584\u9752\u8272\u3067\u5857\u3089\u308c\u305f\u5e73\u9762\u3067\u631f\u307e\u308c\u305f\u90e8\u5206\u306e\u4f53\u7a4d\u3092\u8868\u3059\u3002<\/p>\n<p>\u7e26\u6a2a\u9ad8\u3055\u3068\u3082\u306b $2$ \u306e\u7acb\u65b9\u4f53\u3092\u659c\u3081\u306b2\u3064\u306b\u5207\u3063\u305f\u7acb\u4f53\u3067\u3042\u308b\u306e\u3067\u4f53\u7a4d\u306f<\/p>\n<p>$$\\frac{2\\cdot 2\\cdot 2}{2} = 4$$<\/p>\n<p>\u306b\u306a\u308b\u306f\u305a\u3067\uff0c\u7b54\u3048\u3082\u78ba\u304b\u306b\u305d\u3046\u306a\u3063\u3066\u3044\u308b\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h4 id=\"\u4f8b\u984c\">\u4f8b\u984c<\/h4>\n<p>$$I = \\iint_D dx dy, \\quad D: x^2 + y^2 \\leq 1$$<\/p>\n<p>\u7d2f\u6b21\u7a4d\u5206\u306e\u5f62\u306b\u3059\u308b\u3002\u307e\u305a\u5148\u306b $-\\sqrt{1-x^2} &lt; y &lt; \\sqrt{1-x^2}$ \u306e\u7bc4\u56f2\u3067 $y$ \u3067\u5b9a\u7a4d\u5206\u3057\uff0c\u305d\u306e\u3042\u3068\u306b $ -1 &lt; x &lt; 1$ \u306e\u7bc4\u56f2\u3067 $x$ \u3067\u5b9a\u7a4d\u5206\u3059\u308b\u3002<br \/>\n\u88ab\u7a4d\u5206\u95a2\u6570\u306f $1$\u3002<\/p>\n<p>\u30b1\u30fc\u30b9 2 \u306e\u7d2f\u6b21\u7a4d\u5206\u306e\u5f62\uff1a<br \/>\n\\begin{eqnarray}<br \/>\nf(x, y) &amp;=&amp; 1 \\\\<br \/>\nx_1 &amp;=&amp; -1\\\\<br \/>\nx_2 &amp;=&amp; 1 \\\\<br \/>\n\\phi_1(x) &amp;=&amp; -\\sqrt{1-x^2}\\\\<br \/>\n\\phi_2(x) &amp;=&amp; \\sqrt{1-x^2}<br \/>\n\\end{eqnarray}<\/p>\n<p>\\begin{eqnarray}<br \/>\n\\iint_D dx dy &amp;=&amp; \\int_{x_1}^{x_2} \\left\\{\\int_{\\phi_1(x)}^{\\phi_2(x) } f(x, y)\\,dy \\right\\} dx \\\\<br \/>\n&amp;=&amp; \\cdots<br \/>\n\\end{eqnarray}<\/p>\n<p>$I$ \u306f\u534a\u5f84 $1$ \u306e\u5186\u306e\u9762\u7a4d\u3067\u3042\u308b\u304b\u3089\uff0c$\\pi r^2 = \\pi \\ (\\because r=1)$ \u3068\u306a\u308b\u306f\u305a\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[6]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">f<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">1<\/span>\r\n<span class=\"n\">x1<\/span> <span class=\"o\">=<\/span> <span class=\"o\">-<\/span><span class=\"mi\">1<\/span>\r\n<span class=\"n\">x2<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">1<\/span>\r\n<span class=\"n\">phi1<\/span> <span class=\"o\">=<\/span> <span class=\"o\">-<\/span><span class=\"n\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"mi\">1<\/span><span class=\"o\">-<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">phi2<\/span> <span class=\"o\">=<\/span> <span class=\"n\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"mi\">1<\/span><span class=\"o\">-<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">Eq<\/span><span class=\"p\">(<\/span><span class=\"n\">Integral<\/span><span class=\"p\">(<\/span><span class=\"n\">f<\/span><span class=\"p\">,<\/span> <span class=\"p\">(<\/span><span class=\"n\">y<\/span><span class=\"p\">,<\/span> <span class=\"n\">phi1<\/span><span class=\"p\">,<\/span> <span class=\"n\">phi2<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">x1<\/span><span class=\"p\">,<\/span> <span class=\"n\">x2<\/span><span class=\"p\">)),<\/span> \r\n   <span class=\"n\">Integral<\/span><span class=\"p\">(<\/span><span class=\"n\">f<\/span><span class=\"p\">,<\/span> <span class=\"p\">(<\/span><span class=\"n\">y<\/span><span class=\"p\">,<\/span> <span class=\"n\">phi1<\/span><span class=\"p\">,<\/span> <span class=\"n\">phi2<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">x1<\/span><span class=\"p\">,<\/span> <span class=\"n\">x2<\/span><span class=\"p\">))<\/span><span class=\"o\">.<\/span><span class=\"n\">doit<\/span><span class=\"p\">())<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[6]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\int\\limits_{-1}^{1}\\int\\limits_{- \\sqrt{1 &#8211; x^{2}}}^{\\sqrt{1 &#8211; x^{2}}} 1\\, dy\\, dx = \\pi$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h4 id=\"\u53c2\u8003\uff1a\u9818\u57df-$D$-\u306e\u5857\u308a\u3064\u3076\u3057\">\u53c2\u8003\uff1a\u9818\u57df $D$ \u306e\u5857\u308a\u3064\u3076\u3057<\/h4>\n<p>\u3064\u3044\u3067\u306b\uff0c\u9818\u57df $D$ \u306e\u90e8\u5206\u3092\u5857\u308a\u3064\u3076\u3057\u3066\u63cf\u3044\u3066\u307f\u307e\u3059\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[7]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"c1\"># NumPy \u3092 import \u3057\u3066\u4f7f\u3046\u5834\u5408<\/span>\r\n<span class=\"kn\">import<\/span> <span class=\"nn\">numpy<\/span> <span class=\"k\">as<\/span> <span class=\"nn\">np<\/span>\r\n<span class=\"n\">x_array<\/span> <span class=\"o\">=<\/span> <span class=\"n\">np<\/span><span class=\"o\">.<\/span><span class=\"n\">linspace<\/span><span class=\"p\">(<\/span><span class=\"n\">x1<\/span><span class=\"p\">,<\/span> <span class=\"n\">x2<\/span><span class=\"p\">,<\/span> <span class=\"mi\">100<\/span><span class=\"p\">)<\/span>\r\n<span class=\"k\">def<\/span> <span class=\"nf\">y1<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"k\">return<\/span> <span class=\"o\">-<\/span><span class=\"n\">np<\/span><span class=\"o\">.<\/span><span class=\"n\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"mi\">1<\/span><span class=\"o\">-<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"k\">def<\/span> <span class=\"nf\">y2<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"k\">return<\/span> <span class=\"n\">np<\/span><span class=\"o\">.<\/span><span class=\"n\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"mi\">1<\/span><span class=\"o\">-<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">y1_array<\/span> <span class=\"o\">=<\/span> <span class=\"n\">y1<\/span><span class=\"p\">(<\/span><span class=\"n\">x_array<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">y2_array<\/span> <span class=\"o\">=<\/span> <span class=\"n\">y2<\/span><span class=\"p\">(<\/span><span class=\"n\">x_array<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plot<\/span><span class=\"p\">(<\/span><span class=\"n\">xlim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mf\">1.5<\/span><span class=\"p\">,<\/span> <span class=\"mf\">1.5<\/span><span class=\"p\">),<\/span> \r\n     <span class=\"n\">ylim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mf\">1.5<\/span><span class=\"p\">,<\/span> <span class=\"mf\">1.5<\/span><span class=\"p\">),<\/span> \r\n     <span class=\"n\">aspect<\/span> <span class=\"o\">=<\/span> <span class=\"s2\">\"equal\"<\/span><span class=\"p\">,<\/span>\r\n     <span class=\"n\">title<\/span> <span class=\"o\">=<\/span> <span class=\"s2\">\"\u9818\u57df $D$\"<\/span><span class=\"p\">,<\/span>\r\n     <span class=\"n\">xlabel<\/span> <span class=\"o\">=<\/span> <span class=\"s2\">\"x\"<\/span><span class=\"p\">,<\/span> <span class=\"n\">ylabel<\/span> <span class=\"o\">=<\/span> <span class=\"s2\">\"y\"<\/span><span class=\"p\">,<\/span>\r\n     <span class=\"n\">fill<\/span><span class=\"o\">=<\/span><span class=\"p\">{<\/span><span class=\"s1\">'x'<\/span><span class=\"p\">:<\/span> <span class=\"n\">x_array<\/span><span class=\"p\">,<\/span>\r\n           <span class=\"s1\">'y1'<\/span><span class=\"p\">:<\/span><span class=\"n\">y1_array<\/span><span class=\"p\">,<\/span>\r\n           <span class=\"s1\">'y2'<\/span><span class=\"p\">:<\/span><span class=\"n\">y2_array<\/span><span class=\"p\">,<\/span>\r\n           <span class=\"s1\">'color'<\/span><span class=\"p\">:<\/span><span class=\"s1\">'yellow'<\/span><span class=\"p\">})<\/span><span class=\"p\">;<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_svg output_subarea \">\n<p><!--?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?--><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-6203\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/pmathc303.svg\" alt=\"\" width=\"640\" height=\"427\" \/><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h3 id=\"\u30ac\u30a6\u30b9\u7a4d\u5206\u306e\u6e96\u5099\u3068\u3057\u3066\u306e2\u91cd\u7a4d\u5206\">\u30ac\u30a6\u30b9\u7a4d\u5206\u306e\u6e96\u5099\u3068\u3057\u3066\u306e2\u91cd\u7a4d\u5206<\/h3>\n<p>$$ S = \\int_{-\\infty}^{\\infty} \\int_{-\\infty}^{\\infty} e^{-x^2 &#8211; y^2} dx\\,dy$$<\/p>\n<h4 id=\"\u6975\u5ea7\u6a19\u306b\u5909\u63db\u3057\u3066\u7a4d\u5206\">\u6975\u5ea7\u6a19\u306b\u5909\u63db\u3057\u3066\u7a4d\u5206<\/h4>\n<p>\u6975\u5ea7\u6a19\u306b\u5909\u63db\u3057\u3066\u8a08\u7b97\u3059\u308b\u3002<\/p>\n<p>\u30c7\u30ab\u30eb\u30c8\u5ea7\u6a19 $x, y$ \u30682\u6b21\u5143\u6975\u5ea7\u6a19 $r, \\theta$ \u3068\u306e\u9593\u306e\u95a2\u4fc2\u306f<br \/>\n\\begin{eqnarray}<br \/>\nx &amp;=&amp; r\\cos\\theta\\\\<br \/>\ny &amp;=&amp; r\\sin\\theta<br \/>\n\\end{eqnarray}<br \/>\n\u3067\u3042\u308a\uff0c\u5fae\u5c0f\u9762\u7a4d\u8981\u7d20 $dx\\,dy$ \u306f\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u5909\u63db\u3055\u308c\u308b\u3002<\/p>\n<p>\\begin{eqnarray}<br \/>\ndx\\,dy &amp;=&amp; \\frac{\\partial(x, y)}{\\partial(r,\\theta)} \\,dr\\, d\\theta<br \/>\n\\end{eqnarray}<\/p>\n<p>\u3053\u3053\u3067\uff0c\u3053\u3053\u3067\uff0c$$ \\frac{\\partial(x,y)}{\\partial(r,\\theta)} \\equiv<br \/>\n\\begin{vmatrix}<br \/>\n\\frac{\\partial x}{\\partial r} &amp; \\frac{\\partial x}{\\partial \\theta}\\\\<br \/>\n\\frac{\\partial y}{\\partial r} &amp; \\frac{\\partial y}{\\partial \\theta}\\\\<br \/>\n\\end{vmatrix}<br \/>\n= \\frac{\\partial x}{\\partial r} \\frac{\\partial y}{\\partial \\theta} &#8211; \\frac{\\partial y}{\\partial r}\\frac{\\partial x}{\\partial \\theta}$$ \u3092\u30e4\u30b3\u30d3\u30a2\u30f3\u3068\u3044\u3046\u3002<\/p>\n<h5 id=\"\u30e4\u30b3\u30d3\u884c\u4f8b\">\u30e4\u30b3\u30d3\u884c\u4f8b<\/h5>\n<p>SymPy \u3067\u30e4\u30b3\u30d3\u30a2\u30f3\u3092\u8a08\u7b97\u3059\u308b\u4f8b\uff1a<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[8]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">X<\/span> <span class=\"o\">=<\/span> <span class=\"n\">Matrix<\/span><span class=\"p\">([<\/span><span class=\"n\">r<\/span><span class=\"o\">*<\/span><span class=\"n\">cos<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">),<\/span> <span class=\"n\">r<\/span><span class=\"o\">*<\/span><span class=\"n\">sin<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">)])<\/span>\r\n<span class=\"n\">Y<\/span> <span class=\"o\">=<\/span> <span class=\"n\">Matrix<\/span><span class=\"p\">([<\/span><span class=\"n\">r<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">])<\/span>\r\n\r\n<span class=\"n\">jaco<\/span> <span class=\"o\">=<\/span> <span class=\"n\">X<\/span><span class=\"o\">.<\/span><span class=\"n\">jacobian<\/span><span class=\"p\">(<\/span><span class=\"n\">Y<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">jaco<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[8]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\left[\\begin{matrix}\\cos{\\left(\\theta \\right)} &amp; &#8211; r \\sin{\\left(\\theta \\right)}\\\\\\sin{\\left(\\theta \\right)} &amp; r \\cos{\\left(\\theta \\right)}\\end{matrix}\\right]$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h5 id=\"\u30e4\u30b3\u30d3\u30a2\u30f3\uff08\u30e4\u30b3\u30d3\u884c\u5217\u5f0f\uff09\">\u30e4\u30b3\u30d3\u30a2\u30f3\uff08\u30e4\u30b3\u30d3\u884c\u5217\u5f0f\uff09<\/h5>\n<p><code>jacobian()<\/code> \u306f\u30e4\u30b3\u30d3\u884c\u5217\u3002\u3053\u306e\u884c\u5217\u306e\u884c\u5217\u5f0f\u304c\uff0c\u6211\u3005\u304c\u6c42\u3081\u305f\u3044\u30e4\u30b3\u30d3\u30a2\u30f3\u306f\u30e4\u30b3\u30d3\u884c\u5217\u306e\u884c\u5217\u5f0f\u3002<code>det()<\/code> \u3067\u884c\u5217\u5f0f\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[9]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">det<\/span><span class=\"p\">(<\/span><span class=\"n\">jaco<\/span><span class=\"p\">)<\/span><span class=\"o\">.<\/span><span class=\"n\">simplify<\/span><span class=\"p\">()<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[9]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle r$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<p>\u30e4\u30b3\u30d3\u30a2\u30f3\u304c $r$ \u3068\u306a\u308b\u3053\u3068\u304c\u308f\u304b\u3063\u305f\u3002<\/p>\n<p>\\begin{eqnarray}<br \/>\n\\therefore\\ \\ S &amp;=&amp; \\int_{-\\infty}^{\\infty} \\int_{-\\infty}^{\\infty} e^{-x^2 &#8211; y^2} dx\\,dy \\\\<br \/>\n&amp;=&amp; \\int_0^{2\\pi} \\int_0^{\\infty} e^{-r^2} r \\,dr \\,d\\theta\\\\<br \/>\n&amp;=&amp; 2 \\pi \\int_0^{\\infty} e^{-r^2} r \\,dr = \\cdots<br \/>\n\\end{eqnarray}<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[10]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">Eq<\/span><span class=\"p\">(<\/span><span class=\"n\">Integral<\/span><span class=\"p\">(<\/span><span class=\"n\">exp<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"n\">r<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"n\">r<\/span><span class=\"p\">,<\/span> <span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">pi<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">r<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"n\">oo<\/span><span class=\"p\">)),<\/span>\r\n   <span class=\"n\">Integral<\/span><span class=\"p\">(<\/span><span class=\"n\">exp<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"n\">r<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"n\">r<\/span><span class=\"p\">,<\/span> <span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">pi<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">r<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"n\">oo<\/span><span class=\"p\">))<\/span><span class=\"o\">.<\/span><span class=\"n\">doit<\/span><span class=\"p\">())<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[10]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\int\\limits_{0}^{\\infty}\\int\\limits_{0}^{2 \\pi} r e^{- r^{2}}\\, d\\theta\\, dr = \\pi$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h3 id=\"\u30ac\u30a6\u30b9\u7a4d\u5206\">\u30ac\u30a6\u30b9\u7a4d\u5206<\/h3>\n<p>$$ I = \\int_{-\\infty}^{\\infty} e^{-x^2 } dx$$<\/p>\n<p>\u3068\u3059\u308b\u3068\uff0c\u4e0a\u3067<\/p>\n<p>\\begin{eqnarray}<br \/>\nS &amp;=&amp; \\int_{-\\infty}^{\\infty} \\int_{-\\infty}^{\\infty} e^{-x^2 &#8211; y^2} dx\\,dy \\\\<br \/>\n&amp;=&amp; I^2 \\\\<br \/>\n&amp;=&amp; \\pi<br \/>\n\\end{eqnarray}<\/p>\n<p>\u3068\u3044\u3046\u3053\u3068\u304c\u308f\u304b\u3063\u305f\u306e\u3067\uff0c<\/p>\n<p>$$ I = \\int_{-\\infty}^{\\infty} e^{-x^2 } dx = \\sqrt{S} = \\sqrt{\\pi}$$<\/p>\n<p>\u3068\u3044\u3046\u3053\u3068\u306b\u306a\u308b\u3002SymPy \u304c\u30ac\u30a6\u30b9\u7a4d\u5206 $I$ \u3092\u77e5\u3063\u3066\u3044\u308b\u3053\u3068\u3092\u78ba\u8a8d\u3059\u308b\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[11]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">Eq<\/span><span class=\"p\">(<\/span><span class=\"n\">Integral<\/span><span class=\"p\">(<\/span><span class=\"n\">exp<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"n\">oo<\/span><span class=\"p\">,<\/span> <span class=\"n\">oo<\/span><span class=\"p\">)),<\/span>\r\n   <span class=\"n\">Integral<\/span><span class=\"p\">(<\/span><span class=\"n\">exp<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"n\">oo<\/span><span class=\"p\">,<\/span> <span class=\"n\">oo<\/span><span class=\"p\">))<\/span><span class=\"o\">.<\/span><span class=\"n\">doit<\/span><span class=\"p\">())<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[11]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\int\\limits_{-\\infty}^{\\infty} e^{- x^{2}}\\, dx = \\sqrt{\\pi}$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<p>\u88ab\u7a4d\u5206\u95a2\u6570\u306f\u5076\u95a2\u6570\u3067\u3042\u308b\u304b\u3089\uff0c\u7a4d\u5206\u7bc4\u56f2\u3092 $0$ \u304b\u3089\u306b\u3059\u308b\u3068\u7b54\u3048\u306f\u534a\u5206\u306b\u306a\u308b\u306f\u305a\u3067\uff0c\u3053\u308c\u3082\u78ba\u8a8d\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[12]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">Eq<\/span><span class=\"p\">(<\/span><span class=\"n\">Integral<\/span><span class=\"p\">(<\/span><span class=\"n\">exp<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"n\">oo<\/span><span class=\"p\">)),<\/span>\r\n   <span class=\"n\">Integral<\/span><span class=\"p\">(<\/span><span class=\"n\">exp<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"n\">oo<\/span><span class=\"p\">))<\/span><span class=\"o\">.<\/span><span class=\"n\">doit<\/span><span class=\"p\">())<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[12]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\int\\limits_{0}^{\\infty} e^{- x^{2}}\\, dx = \\frac{\\sqrt{\\pi}}{2}$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<p>\u30ac\u30a6\u30b9\u7a4d\u5206\u306f\uff0c\u7a4d\u5206\u306e\u4e0a\u9650\u304c $\\infty$ \u306e\u3068\u304d\u306b\u3060\u3051\u89e3\u6790\u7684\u306b\u7a4d\u5206\u3067\u304d\u308b\u4f8b\u3002<\/p>\n<p>\u7a4d\u5206\u7bc4\u56f2\u304c $0$ \u304b\u3089\u4e00\u822c\u306b $x$ \u307e\u3067\u3060\u3068\u7a4d\u5206\u3067\u304d\u306a\u3044\uff08\u7d50\u679c\u304c\u521d\u7b49\u95a2\u6570\u3067\u8868\u305b\u306a\u3044\uff09\u3002\u5b9f\u969b\u306b\u3084\u3063\u3066\u307f\u308b\u3068&#8230;<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[13]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">Eq<\/span><span class=\"p\">(<\/span><span class=\"n\">Integral<\/span><span class=\"p\">(<\/span><span class=\"n\">exp<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"n\">t<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">t<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"n\">x<\/span><span class=\"p\">)),<\/span>\r\n   <span class=\"n\">Integral<\/span><span class=\"p\">(<\/span><span class=\"n\">exp<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"n\">t<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">t<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"n\">x<\/span><span class=\"p\">))<\/span><span class=\"o\">.<\/span><span class=\"n\">doit<\/span><span class=\"p\">())<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[13]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\int\\limits_{0}^{x} e^{- t^{2}}\\, dt = \\frac{\\sqrt{\\pi} \\operatorname{erf}{\\left(x \\right)}}{2}$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h4 id=\"\u53c2\u8003\uff1a\u8aa4\u5dee\u95a2\u6570\">\u53c2\u8003\uff1a\u8aa4\u5dee\u95a2\u6570<\/h4>\n<p>$\\operatorname{erf} (x)$ \u306f\u300c\u8aa4\u5dee\u95a2\u6570\u300d\u3068\u547c\u3070\u308c\uff0c\u3053\u308c\u304c\u5b9a\u7fa9\uff1a<\/p>\n<p>\\begin{eqnarray}<br \/>\n\\operatorname{erf} (x) &amp;\\equiv&amp; \\frac{2}{\\sqrt{\\pi}} \\int_0^x e^{-t^2} \\, dt<br \/>\n\\end{eqnarray}<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h3 id=\"\u30ac\u30a6\u30b9\u7a4d\u5206\u306b\u95a2\u9023\u3057\u305f\u7a4d\u5206\">\u30ac\u30a6\u30b9\u7a4d\u5206\u306b\u95a2\u9023\u3057\u305f\u7a4d\u5206<\/h3>\n<p>$$I(m) \\equiv \\int_{0}^{\\infty} x^m e^{-x^2 } dx$$<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[14]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"k\">def<\/span> <span class=\"nf\">I<\/span><span class=\"p\">(<\/span><span class=\"n\">m<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"k\">return<\/span> <span class=\"n\">integrate<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"n\">m<\/span> <span class=\"o\">*<\/span> <span class=\"n\">exp<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"n\">oo<\/span><span class=\"p\">))<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[15]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[15]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\frac{\\sqrt{\\pi}}{2}$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[16]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[16]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\frac{\\sqrt{\\pi}}{4}$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[17]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">4<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[17]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\frac{3 \\sqrt{\\pi}}{8}$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[18]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">6<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[18]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\frac{15 \\sqrt{\\pi}}{16}$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<p>$m = 2 n$ ($n$ \u306f\u6574\u6570) \u306e\u5834\u5408\u3092\u63a8\u6e2c\u3059\u308b\u305f\u3081\u306b\uff0c$I(0)$ \u306e\u5024\u3067\u5272\u3063\u3066\u307f\u308b\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[19]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[19]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\frac{1}{2}$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[20]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">4<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[20]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\frac{3}{4}$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[21]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">6<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[21]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\frac{15}{8}$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<p>\u3068\u3044\u3046\u3053\u3068\u306f\uff0c<\/p>\n<p>\\begin{eqnarray}<br \/>\nI(2) &amp;=&amp; 1\\times \\frac{1}{2} \\times I(0), \\\\<br \/>\nI(4) &amp;=&amp; \\frac{3\\cdot 1}{2^2} \\times I(0), \\\\<br \/>\nI(6) &amp;=&amp; \\frac{5\\cdot3\\cdot 1}{2^3} \\times I(0), \\dots<br \/>\n\\end{eqnarray}<\/p>\n<p>\u3060\u304b\u3089\uff0c<\/p>\n<p>$$\\displaystyle I(2 n) = \\frac{(2n-1)!!}{2^{n}} \\frac{\\sqrt{\\pi}}{2}$$<\/p>\n<p>\u3068\u63a8\u6e2c\u3055\u308c\u308b\u3002\u5b9f\u969b\u78ba\u304b\u3081\u3066\u307f\u308b\u3068&#8230;<\/p>\n<p>\uff08SymPy \u3067\u306f $n! = $ <code>factorial(n)<\/code>\uff0c$n!! = $ <code>factorial2(n)<\/code>\uff09<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[22]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"c1\"># m &gt;= 2 \u306e\u5076\u6570\u306b\u5bfe\u3057\u3066\uff0c\u4e0a\u8a18\u306e 2n \u3092 m \u3068\u304a\u3044\u3066<\/span>\r\n\r\n<span class=\"k\">def<\/span> <span class=\"nf\">EvenI<\/span><span class=\"p\">(<\/span><span class=\"n\">m<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"k\">return<\/span> <span class=\"n\">factorial2<\/span><span class=\"p\">(<\/span><span class=\"n\">m<\/span><span class=\"o\">-<\/span><span class=\"mi\">1<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"mi\">2<\/span><span class=\"o\">**<\/span><span class=\"p\">(<\/span><span class=\"n\">m<\/span><span class=\"o\">\/<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span> <span class=\"o\">*<\/span> <span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[23]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"c1\"># \u4f8b\u3048\u3070 m = 8 \u306e\u3068\u304d<\/span>\r\n<span class=\"n\">EvenI<\/span><span class=\"p\">(<\/span><span class=\"mi\">8<\/span><span class=\"p\">)<\/span> <span class=\"o\">-<\/span> <span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">8<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[23]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle 0$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<p>\u3053\u308c\u3067 $m \\geq 2$ \u304c\u5076\u6570\u306e\u5834\u5408\u306e $I(m)$ \u306f\u6c42\u307e\u3063\u305f\u3002<\/p>\n<p>\u6b21\u306b\uff0c$m$ \u304c\u5947\u6570\uff08$m = 2 n +1$\uff09\u306e\u5834\u5408\u306e $I(m)$ \u306e\u4e00\u822c\u5f62\u3092\u4e88\u60f3\u3057\u3066\u307f\u308b\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[24]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">1<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[24]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\frac{1}{2}$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[25]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">3<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[25]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\frac{1}{2}$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[26]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">5<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[26]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle 1$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[27]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">7<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[27]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle 3$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[28]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">3<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">1<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[28]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle 1$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[29]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">5<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">1<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[29]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle 2$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[30]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">7<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"n\">I<\/span><span class=\"p\">(<\/span><span class=\"mi\">1<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[30]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle 6$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<p>\u4ee5\u4e0a\u306e\u3053\u3068\u304b\u3089\u6b21\u306e\u3088\u3046\u306a\u63a8\u6e2c\u304c\u3067\u304d\u308b\u3002<\/p>\n<p>$$ I(2n+1) = n! \\times I(1) $$<\/p>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"excerpt":{"rendered":"","protected":false},"author":33,"featured_media":0,"parent":6176,"menu_order":30,"comment_status":"closed","ping_status":"closed","template":"","meta":{"inline_featured_image":false,"footnotes":""},"class_list":["post-6200","page","type-page","status-publish","hentry","nodate","item-wrap"],"aioseo_notices":[],"_links":{"self":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/6200","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/users\/33"}],"replies":[{"embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/comments?post=6200"}],"version-history":[{"count":2,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/6200\/revisions"}],"predecessor-version":[{"id":8093,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/6200\/revisions\/8093"}],"up":[{"embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/6176"}],"wp:attachment":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/media?parent=6200"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}