{"id":5530,"date":"2023-03-06T17:26:53","date_gmt":"2023-03-06T08:26:53","guid":{"rendered":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/?p=5530"},"modified":"2023-03-14T16:13:09","modified_gmt":"2023-03-14T07:13:09","slug":"5530","status":"publish","type":"post","link":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/5530\/","title":{"rendered":"Python \u306e Matplotlib \u3067\u6b63\u898f\u5206\u5e03\u3092\u03c3\u3054\u3068\u306b\u5857\u308a\u308f\u3051\u308b"},"content":{"rendered":"<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>\u300c<a href=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/5467\/\">gnuplot \u3067\u6b63\u898f\u5206\u5e03\u3092\u03c3\u3054\u3068\u306b\u5857\u308a\u308f\u3051\u308b<\/a>\u300d\u306e Python \u7248\u3002<\/p>\n<p>\u4ee5\u4e0b\u3092\u53c2\u8003\u306b\uff0c<a href=\"https:\/\/ja.wikipedia.org\/wiki\/%E3%83%95%E3%82%A1%E3%82%A4%E3%83%AB:Standard_deviation_diagram.svg\">\u6a19\u6e96\u6b63\u898f\u5206\u5e03\u304c\u3082\u3064\u78ba\u7387\u5bc6\u5ea6\u95a2\u6570\u306e\u30b0\u30e9\u30d5<\/a>\u3092\u63cf\u3044\u3066\u307f\u308b\u3002<\/p>\n<ul>\n<li><a href=\"https:\/\/ja.wikipedia.org\/wiki\/%E6%AD%A3%E8%A6%8F%E5%88%86%E5%B8%83\">\u6b63\u898f\u5206\u5e03 &#8211; Wikipedia<\/a><\/li>\n<li><a href=\"https:\/\/ja.wikipedia.org\/wiki\/%E8%AA%A4%E5%B7%AE%E9%96%A2%E6%95%B0\">\u8aa4\u5dee\u95a2\u6570 &#8211; Wikipedia<\/a><\/li>\n<\/ul>\n<\/div>\n<\/div>\n<\/div>\n<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<h3 id=\"\u6b63\u898f\u5206\u5e03\u95a2\u6570\">\u6b63\u898f\u5206\u5e03\u95a2\u6570<\/h3>\n<p>$$N(\\mu, \\sigma^2) \\equiv \\frac{1}{\\sqrt{2\\pi} \\sigma} \\exp\\left\\{ \u2013 \\frac{(x-\\mu)^2}{2\\sigma^2}\\right\\}$$<\/p>\n<p>\u7c21\u5358\u306e\u305f\u3081\u306b\uff0c$\\mu = 0, \\sigma = 1$ \u3068\u3057\u305f\u6a19\u6e96\u6b63\u898f\u5206\u5e03 $f(x)$ \u3092\u4f7f\u3046\u3053\u3068\u306b\u3059\u308b\u3002<\/p>\n<p>$$f(x) \\equiv N(0, 1) = \\frac{1}{\\sqrt{2\\pi}} \\exp\\left( \u2013 \\frac{x^2}{2}\\right)$$<\/p>\n<h3 id=\"\u8aa4\u5dee\u95a2\u6570\">\u8aa4\u5dee\u95a2\u6570<\/h3>\n<p>$$\\mbox{erf}(x) \\equiv \\frac{2}{\\sqrt{\\pi}} \\int_0^x e^{-t^2}\\, dt $$<\/p>\n<h3 id=\"\u6b63\u898f\u5206\u5e03\u306e\u7a4d\u5206\u3092\u8aa4\u5dee\u95a2\u6570\u3067\u8868\u3059\">\u6b63\u898f\u5206\u5e03\u306e\u7a4d\u5206\u3092\u8aa4\u5dee\u95a2\u6570\u3067\u8868\u3059<\/h3>\n<p>\\begin{eqnarray}<br \/>\nF(x_1, x_2) &amp;\\equiv&amp; \\int_{x_1}^{x_2} f(t)\\, dt \\\\<br \/>\n&amp;=&amp; \\frac{1}{\\sqrt{2\\pi}} \\int_{x_1}^{x_2}\\exp\\left( \u2013 \\frac{t^2}{2}\\right)\\, dt \\\\<br \/>\n&amp;&amp; \\quad\\left(\\frac{t}{\\sqrt{2}} \\equiv x, \\ dt = \\sqrt{2} dx \\right) \\\\<br \/>\n&amp;=&amp; \\frac{1}{\\sqrt{\\pi}} \\int_{\\frac{x_1}{\\sqrt{2}}}^{\\frac{x_2}{\\sqrt{2}}} e^{-x^2} \\,dx \\\\<br \/>\n&amp;=&amp; \\frac{1}{2} \\left(\\mbox{erf}\\left(\\frac{x_2}{\\sqrt{2}}\\right)-\\mbox{erf}\\left(\\frac{x_1}{\\sqrt{2}}\\right) \\right)<br \/>\n\\end{eqnarray}<\/p>\n<p>\u8aa4\u5dee\u95a2\u6570 <code>scipy.special.erf(x)<\/code> \u3092\u4f7f\u3044\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\">import<\/span> <span class=\"nn\">matplotlib.pyplot<\/span> <span class=\"k\">as<\/span> <span class=\"nn\">plt<\/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=\"kn\">from<\/span> <span class=\"nn\">scipy<\/span> <span class=\"kn\">import<\/span> <span class=\"n\">special<\/span>\r\n\r\n<span class=\"c1\"># \u4ee5\u4e0b\u306f\u30b0\u30e9\u30d5\u3092 SVG \u3067 Notebook \u306b\u30a4\u30f3\u30e9\u30a4\u30f3\u8868\u793a\u3055\u305b\u308b\u8a2d\u5b9a<\/span>\r\n<span class=\"o\">%<\/span><span class=\"k\">config<\/span> InlineBackend.figure_formats = ['svg']\r\n<span class=\"c1\"># \u30b0\u30e9\u30d5\u306e\u30b5\u30a4\u30ba<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">rcParams<\/span><span class=\"p\">[<\/span><span class=\"s2\">\"figure.figsize\"<\/span><span class=\"p\">]<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"mf\">6.4<\/span><span class=\"p\">,<\/span> <span class=\"mf\">2.4<\/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[2]:<\/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\">f<\/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\">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=\"o\">\/<\/span> <span class=\"mi\">2<\/span><span class=\"p\">)<\/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\">2<\/span> <span class=\"o\">*<\/span> <span class=\"n\">np<\/span><span class=\"o\">.<\/span><span class=\"n\">pi<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"k\">def<\/span> <span class=\"nf\">F<\/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=\"k\">return<\/span> <span class=\"mf\">0.5<\/span><span class=\"o\">*<\/span><span class=\"p\">(<\/span><span class=\"n\">special<\/span><span class=\"o\">.<\/span><span class=\"n\">erf<\/span><span class=\"p\">(<\/span><span class=\"n\">x2<\/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\">2<\/span><span class=\"p\">))<\/span><span class=\"o\">-<\/span><span class=\"n\">special<\/span><span class=\"o\">.<\/span><span class=\"n\">erf<\/span><span class=\"p\">(<\/span><span class=\"n\">x1<\/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\">2<\/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 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=\"\u6a19\u6e96\u6b63\u898f\u5206\u5e03\u3092-plt.plot()\">\u6a19\u6e96\u6b63\u898f\u5206\u5e03\u3092 plt.plot()<\/h3>\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=\"n\">x<\/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=\"o\">-<\/span><span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"mi\">401<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"c1\"># \u6a2a\u8ef8\u7e26\u8ef8\u306e\u8868\u793a\u7bc4\u56f2<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xlim<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">ylim<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">0.41<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">plot<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s1\">'#0000ff'<\/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><!-- Created with matplotlib (https:\/\/matplotlib.org\/) --><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-5566\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/pbunpu01b.svg\" alt=\"\" width=\"640\" height=\"240\" \/><\/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=\"\u03c3\u3054\u3068\u306b\u5857\u308a\u308f\u3051\u308b\">\u03c3\u3054\u3068\u306b\u5857\u308a\u308f\u3051\u308b<\/h3>\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\"># \u5857\u308a\u308f\u3051\u308b\u8272<\/span>\r\n<span class=\"n\">c<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"s1\">'#0000ff'<\/span><span class=\"p\">,<\/span><span class=\"s1\">'#1e90ff'<\/span><span class=\"p\">,<\/span><span class=\"s1\">'#87ceeb'<\/span><span class=\"p\">,<\/span><span class=\"s1\">'#add8e6'<\/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[5]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"c1\"># \u03c3\u3054\u3068\u306b\u5857\u308a\u308f\u3051\u308b<\/span>\r\n<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\">4<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">fill_between<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">where<\/span><span class=\"o\">=<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"n\">i<\/span><span class=\"o\">-<\/span><span class=\"mi\">1<\/span><span class=\"o\">&lt;=<\/span><span class=\"n\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">&amp;<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">&lt;=<\/span>  <span class=\"o\">-<\/span><span class=\"n\">i<\/span><span class=\"p\">),<\/span> <span class=\"n\">fc<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">])<\/span>\r\n    <span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">fill_between<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">where<\/span><span class=\"o\">=<\/span><span class=\"p\">(<\/span>   <span class=\"n\">i<\/span><span class=\"o\">&lt;=<\/span><span class=\"n\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">&amp;<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">&lt;=<\/span> <span class=\"n\">i<\/span><span class=\"o\">+<\/span><span class=\"mi\">1<\/span><span class=\"p\">),<\/span> <span class=\"n\">fc<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">])<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xlim<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">ylim<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">0.41<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">plot<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/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\"><\/div>\n<div class=\"output_svg output_subarea \">\n<p><!--?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?--><\/p>\n<p><!-- Created with matplotlib (https:\/\/matplotlib.org\/) --><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-5567\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/pbunpu02b.svg\" alt=\"\" width=\"640\" height=\"240\" \/><\/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=\"\u5857\u308a\u5206\u3051\u305f\u9818\u57df\u3092\u767d\u7dda\u3067\u533a\u5207\u308b\">\u5857\u308a\u5206\u3051\u305f\u9818\u57df\u3092\u767d\u7dda\u3067\u533a\u5207\u308b<\/h3>\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=\"c1\"># vlines \u3067\u767d\u3044\u7e26\u7dda\u3092\u5f15\u304f<\/span>\r\n<span class=\"n\">xsigma<\/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=\"o\">-<\/span><span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"mi\">7<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">vlines<\/span><span class=\"p\">(<\/span><span class=\"n\">xsigma<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">xsigma<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s1\">'white'<\/span><span class=\"p\">)<\/span>\r\n\r\n<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\">4<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">fill_between<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">where<\/span><span class=\"o\">=<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"n\">i<\/span><span class=\"o\">-<\/span><span class=\"mi\">1<\/span><span class=\"o\">&lt;=<\/span><span class=\"n\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">&amp;<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">&lt;=<\/span>  <span class=\"o\">-<\/span><span class=\"n\">i<\/span><span class=\"p\">),<\/span> <span class=\"n\">fc<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">])<\/span>\r\n    <span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">fill_between<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">where<\/span><span class=\"o\">=<\/span><span class=\"p\">(<\/span>   <span class=\"n\">i<\/span><span class=\"o\">&lt;=<\/span><span class=\"n\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">&amp;<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">&lt;=<\/span> <span class=\"n\">i<\/span><span class=\"o\">+<\/span><span class=\"mi\">1<\/span><span class=\"p\">),<\/span> <span class=\"n\">fc<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">])<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xlim<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">ylim<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">0.41<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">plot<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/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\"><\/div>\n<div class=\"output_svg output_subarea \">\n<p><!--?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?--><\/p>\n<p><!-- Created with matplotlib (https:\/\/matplotlib.org\/) --><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-5568\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/pbunpu03b.svg\" alt=\"\" width=\"640\" height=\"240\" \/><\/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=\"\u6a2a\u8ef8\u306e\u76ee\u76db\u30e9\u30d9\u30eb\u3092\u8a2d\u5b9a\u3059\u308b\">\u6a2a\u8ef8\u306e\u76ee\u76db\u30e9\u30d9\u30eb\u3092\u8a2d\u5b9a\u3059\u308b<\/h3>\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\"># \u6a2a\u8ef8\u306e\u76ee\u76db\u306e\u4f4d\u7f6e\u3068\u30e9\u30d9\u30eb\u306e\u8a2d\u5b9a<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xticks<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"n\">ticks<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"o\">-<\/span><span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"mi\">2<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"mi\">1<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"mi\">2<\/span><span class=\"p\">,<\/span> <span class=\"mi\">1<\/span><span class=\"p\">],<\/span> \r\n  <span class=\"n\">labels<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"s2\">\"\u03bc-4\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc-3\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc-2\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc-\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc\"<\/span><span class=\"p\">,<\/span>\r\n            <span class=\"s2\">\"\u03bc+4\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc+3\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc+2\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc+\u03c3\"<\/span><span class=\"p\">])<\/span>\r\n\r\n<span class=\"n\">xsigma<\/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=\"o\">-<\/span><span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"mi\">7<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">vlines<\/span><span class=\"p\">(<\/span><span class=\"n\">xsigma<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">xsigma<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s1\">'white'<\/span><span class=\"p\">)<\/span>\r\n\r\n<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\">4<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">fill_between<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">where<\/span><span class=\"o\">=<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"n\">i<\/span><span class=\"o\">-<\/span><span class=\"mi\">1<\/span><span class=\"o\">&lt;=<\/span><span class=\"n\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">&amp;<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">&lt;=<\/span>  <span class=\"o\">-<\/span><span class=\"n\">i<\/span><span class=\"p\">),<\/span> <span class=\"n\">fc<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">])<\/span>\r\n    <span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">fill_between<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">where<\/span><span class=\"o\">=<\/span><span class=\"p\">(<\/span>   <span class=\"n\">i<\/span><span class=\"o\">&lt;=<\/span><span class=\"n\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">&amp;<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">&lt;=<\/span> <span class=\"n\">i<\/span><span class=\"o\">+<\/span><span class=\"mi\">1<\/span><span class=\"p\">),<\/span> <span class=\"n\">fc<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">])<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xlim<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">ylim<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">0.41<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">plot<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/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\"><\/div>\n<div class=\"output_svg output_subarea \">\n<p><!--?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?--><\/p>\n<p><!-- Created with matplotlib (https:\/\/matplotlib.org\/) --><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-5569\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/pbunpu04b.svg\" alt=\"\" width=\"640\" height=\"240\" \/><\/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=\"\u5857\u308a\u5206\u3051\u305f\u9818\u57df\u306b\u30e9\u30d9\u30eb\u3092\u3064\u3051\u308b\">\u5857\u308a\u5206\u3051\u305f\u9818\u57df\u306b\u30e9\u30d9\u30eb\u3092\u3064\u3051\u308b<\/h3>\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=\"c1\"># \u5404\u9818\u57df\u306e\u30d1\u30fc\u30bb\u30f3\u30c8\u3092\u8868\u3059\u6587\u5b57\u5217\u3092\u914d\u5217\u3068\u3057\u3066\u4f5c\u6210<\/span>\r\n<span class=\"c1\"># \u30ea\u30b9\u30c8 sig \u306e\u8981\u7d20\u3092\u5168\u3066 0 \u306b<\/span>\r\n<span class=\"n\">sig<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"mi\">0<\/span><span class=\"p\">]<\/span><span class=\"o\">*<\/span><span class=\"mi\">4<\/span>\r\n<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\">0<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">]<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'<\/span><span class=\"si\">%5.2f<\/span><span class=\"s1\">'<\/span> <span class=\"o\">%<\/span> <span class=\"p\">(<\/span><span class=\"n\">F<\/span><span class=\"p\">(<\/span><span class=\"n\">i<\/span><span class=\"p\">,<\/span><span class=\"n\">i<\/span><span class=\"o\">+<\/span><span class=\"mi\">1<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"mi\">100<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"s1\">'%'<\/span>\r\n    <span class=\"nb\">print<\/span><span class=\"p\">(<\/span><span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/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_subarea output_stream output_stdout output_text\">\n<pre>34.13%\r\n13.59%\r\n 2.14%\r\n 0.13%\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[9]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"c1\"># \u9818\u57df\u306b\u30e9\u30d9\u30eb\u3092\u3064\u3051\u308b<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mf\">0.5<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.20<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"mi\">0<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"white\"<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span> <span class=\"mf\">0.5<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.20<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"mi\">0<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"white\"<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mf\">1.5<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.04<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"mi\">1<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"white\"<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span> <span class=\"mf\">1.5<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.04<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"mi\">1<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"white\"<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mf\">2.5<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.05<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"mi\">2<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"black\"<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span> <span class=\"mf\">2.5<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.05<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"mi\">2<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"black\"<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mf\">3.5<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.02<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"mi\">3<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"black\"<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span> <span class=\"mf\">3.5<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.02<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"mi\">3<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"black\"<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xticks<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"n\">ticks<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"o\">-<\/span><span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"mi\">2<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"mi\">1<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"mi\">2<\/span><span class=\"p\">,<\/span> <span class=\"mi\">1<\/span><span class=\"p\">],<\/span> \r\n  <span class=\"n\">labels<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"s2\">\"\u03bc-4\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc-3\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc-2\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc-\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc\"<\/span><span class=\"p\">,<\/span>\r\n            <span class=\"s2\">\"\u03bc+4\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc+3\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc+2\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc+\u03c3\"<\/span><span class=\"p\">])<\/span>\r\n\r\n<span class=\"n\">xsigma<\/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=\"o\">-<\/span><span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"mi\">7<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">vlines<\/span><span class=\"p\">(<\/span><span class=\"n\">xsigma<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">xsigma<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s1\">'white'<\/span><span class=\"p\">)<\/span>\r\n\r\n<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\">4<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">fill_between<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">where<\/span><span class=\"o\">=<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"n\">i<\/span><span class=\"o\">-<\/span><span class=\"mi\">1<\/span><span class=\"o\">&lt;=<\/span><span class=\"n\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">&amp;<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">&lt;=<\/span>  <span class=\"o\">-<\/span><span class=\"n\">i<\/span><span class=\"p\">),<\/span> <span class=\"n\">fc<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">])<\/span>\r\n    <span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">fill_between<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">where<\/span><span class=\"o\">=<\/span><span class=\"p\">(<\/span>   <span class=\"n\">i<\/span><span class=\"o\">&lt;=<\/span><span class=\"n\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">&amp;<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">&lt;=<\/span> <span class=\"n\">i<\/span><span class=\"o\">+<\/span><span class=\"mi\">1<\/span><span class=\"p\">),<\/span> <span class=\"n\">fc<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">])<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xlim<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">ylim<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">0.41<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">plot<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/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\"><\/div>\n<div class=\"output_svg output_subarea \">\n<p><!--?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?--><\/p>\n<p><!-- Created with matplotlib (https:\/\/matplotlib.org\/) --><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-5570\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/pbunpu05b.svg\" alt=\"\" width=\"640\" height=\"240\" \/><\/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=\"\u67a0\u3068\u7e26\u8ef8\u76ee\u76db\u3092\u975e\u8868\u793a\u306b\">\u67a0\u3068\u7e26\u8ef8\u76ee\u76db\u3092\u975e\u8868\u793a\u306b<\/h3>\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=\"c1\"># \u67a0\u3068\u7e26\u8ef8\u76ee\u76db\u3092\u975e\u8868\u793a\u306b<\/span>\r\n<span class=\"n\">fig<\/span><span class=\"p\">,<\/span> <span class=\"n\">ax<\/span> <span class=\"o\">=<\/span> <span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">subplots<\/span><span class=\"p\">()<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">spines<\/span><span class=\"p\">[<\/span><span class=\"s1\">'right'<\/span><span class=\"p\">]<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">spines<\/span><span class=\"p\">[<\/span><span class=\"s1\">'left'<\/span><span class=\"p\">]<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">spines<\/span><span class=\"p\">[<\/span><span class=\"s1\">'top'<\/span><span class=\"p\">]<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">axes<\/span><span class=\"o\">.<\/span><span class=\"n\">yaxis<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mf\">0.5<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.20<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"mi\">0<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"white\"<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span> <span class=\"mf\">0.5<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.20<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"mi\">0<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"white\"<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mf\">1.5<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.04<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"mi\">1<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"white\"<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span> <span class=\"mf\">1.5<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.04<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"mi\">1<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"white\"<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mf\">2.5<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.05<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"mi\">2<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"black\"<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span> <span class=\"mf\">2.5<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.05<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"mi\">2<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"black\"<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mf\">3.5<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.02<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"mi\">3<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"black\"<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span> <span class=\"mf\">3.5<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.02<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig<\/span><span class=\"p\">[<\/span><span class=\"mi\">3<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"black\"<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xticks<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"n\">ticks<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"o\">-<\/span><span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"mi\">2<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"mi\">1<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"mi\">2<\/span><span class=\"p\">,<\/span> <span class=\"mi\">1<\/span><span class=\"p\">],<\/span> \r\n  <span class=\"n\">labels<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"s2\">\"\u03bc-4\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc-3\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc-2\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc-\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc\"<\/span><span class=\"p\">,<\/span>\r\n            <span class=\"s2\">\"\u03bc+4\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc+3\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc+2\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc+\u03c3\"<\/span><span class=\"p\">])<\/span>\r\n\r\n<span class=\"n\">xsigma<\/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=\"o\">-<\/span><span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"mi\">7<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">vlines<\/span><span class=\"p\">(<\/span><span class=\"n\">xsigma<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">xsigma<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s1\">'white'<\/span><span class=\"p\">)<\/span>\r\n\r\n<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\">4<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">fill_between<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">where<\/span><span class=\"o\">=<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"n\">i<\/span><span class=\"o\">-<\/span><span class=\"mi\">1<\/span><span class=\"o\">&lt;=<\/span><span class=\"n\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">&amp;<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">&lt;=<\/span>  <span class=\"o\">-<\/span><span class=\"n\">i<\/span><span class=\"p\">),<\/span> <span class=\"n\">fc<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">])<\/span>\r\n    <span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">fill_between<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">where<\/span><span class=\"o\">=<\/span><span class=\"p\">(<\/span>   <span class=\"n\">i<\/span><span class=\"o\">&lt;=<\/span><span class=\"n\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">&amp;<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">&lt;=<\/span> <span class=\"n\">i<\/span><span class=\"o\">+<\/span><span class=\"mi\">1<\/span><span class=\"p\">),<\/span> <span class=\"n\">fc<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">])<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xlim<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">ylim<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">0.41<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">plot<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/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\"><\/div>\n<div class=\"output_svg output_subarea \">\n<p><!--?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?--><\/p>\n<p><!-- Created with matplotlib (https:\/\/matplotlib.org\/) --><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-5571\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/pbunpu06b.svg\" alt=\"\" width=\"640\" height=\"240\" \/><\/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=\"-n\u03c3-\u304b\u3089-+n\u03c3-\u307e\u3067\u3092\u5857\u308a\u3064\u3076\u3059\">-n\u03c3 \u304b\u3089 +n\u03c3 \u307e\u3067\u3092\u5857\u308a\u3064\u3076\u3059<\/h3>\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=\"c1\"># \u5404\u9818\u57df\u306e\u30d1\u30fc\u30bb\u30f3\u30c8\u3092\u8868\u3059\u6587\u5b57\u5217\u3092\u914d\u5217\u3068\u3057\u3066\u4f5c\u6210<\/span>\r\n<span class=\"c1\"># \u30ea\u30b9\u30c8 sig \u306e\u8981\u7d20\u3092\u5168\u3066 0 \u306b<\/span>\r\n<span class=\"n\">sig2<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"mi\">0<\/span><span class=\"p\">]<\/span><span class=\"o\">*<\/span><span class=\"mi\">4<\/span>\r\n<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\">0<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"n\">sig2<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">]<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'<\/span><span class=\"si\">%5.2f<\/span><span class=\"s1\">'<\/span> <span class=\"o\">%<\/span> <span class=\"p\">(<\/span><span class=\"n\">F<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"n\">i<\/span><span class=\"o\">-<\/span><span class=\"mi\">1<\/span><span class=\"p\">,<\/span><span class=\"n\">i<\/span><span class=\"o\">+<\/span><span class=\"mi\">1<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"mi\">100<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"s1\">'%'<\/span>\r\n    <span class=\"nb\">print<\/span><span class=\"p\">(<\/span><span class=\"n\">sig2<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/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_subarea output_stream output_stdout output_text\">\n<pre>68.27%\r\n95.45%\r\n99.73%\r\n99.99%\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[12]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"c1\"># \u67a0\u3068\u7e26\u8ef8\u76ee\u76db\u3092\u975e\u8868\u793a\u306b<\/span>\r\n<span class=\"n\">fig<\/span><span class=\"p\">,<\/span> <span class=\"n\">ax<\/span> <span class=\"o\">=<\/span> <span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">subplots<\/span><span class=\"p\">()<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">spines<\/span><span class=\"p\">[<\/span><span class=\"s1\">'right'<\/span><span class=\"p\">]<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">spines<\/span><span class=\"p\">[<\/span><span class=\"s1\">'left'<\/span><span class=\"p\">]<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">spines<\/span><span class=\"p\">[<\/span><span class=\"s1\">'top'<\/span><span class=\"p\">]<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">axes<\/span><span class=\"o\">.<\/span><span class=\"n\">yaxis<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xlim<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">ylim<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">0.41<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xticks<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"n\">ticks<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"o\">-<\/span><span class=\"mi\">1<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mi\">1<\/span><span class=\"p\">],<\/span> \r\n  <span class=\"n\">labels<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"s2\">\"\u03bc-\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc\"<\/span><span class=\"p\">,<\/span> \r\n            <span class=\"s2\">\"\u03bc+\u03c3\"<\/span><span class=\"p\">])<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">plot<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s1\">'gray'<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"c1\"># \u9818\u57df\u306b\u30e9\u30d9\u30eb\u3092\u3064\u3051\u308b<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.15<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig2<\/span><span class=\"p\">[<\/span><span class=\"mi\">0<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"white\"<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"c1\"># -1\u03c3 \u304b\u3089 +1\u03c3 \u307e\u3067\u3092\u5857\u308a\u3064\u3076\u3059<\/span>\r\n<span class=\"n\">x1<\/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=\"o\">-<\/span><span class=\"mi\">1<\/span><span class=\"p\">,<\/span> <span class=\"mi\">1<\/span><span class=\"p\">,<\/span> <span class=\"mi\">101<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">fill_between<\/span><span class=\"p\">(<\/span><span class=\"n\">x1<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x1<\/span><span class=\"p\">),<\/span> <span class=\"n\">fc<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/span><span class=\"p\">[<\/span><span class=\"mi\">0<\/span><span class=\"p\">])<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">plot<\/span><span class=\"p\">(<\/span><span class=\"n\">x1<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x1<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/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\"><\/div>\n<div class=\"output_svg output_subarea \">\n<p><!--?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?--><\/p>\n<p><!-- Created with matplotlib (https:\/\/matplotlib.org\/) --><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-5572\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/pbunpu07b.svg\" alt=\"\" width=\"640\" height=\"240\" \/><\/p>\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[13]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"c1\"># \u67a0\u3068\u7e26\u8ef8\u76ee\u76db\u3092\u975e\u8868\u793a\u306b<\/span>\r\n<span class=\"n\">fig<\/span><span class=\"p\">,<\/span> <span class=\"n\">ax<\/span> <span class=\"o\">=<\/span> <span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">subplots<\/span><span class=\"p\">()<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">spines<\/span><span class=\"p\">[<\/span><span class=\"s1\">'right'<\/span><span class=\"p\">]<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">spines<\/span><span class=\"p\">[<\/span><span class=\"s1\">'left'<\/span><span class=\"p\">]<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">spines<\/span><span class=\"p\">[<\/span><span class=\"s1\">'top'<\/span><span class=\"p\">]<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">axes<\/span><span class=\"o\">.<\/span><span class=\"n\">yaxis<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xlim<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">ylim<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">0.41<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xticks<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"n\">ticks<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"o\">-<\/span><span class=\"mi\">2<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"mi\">1<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mi\">2<\/span><span class=\"p\">,<\/span> <span class=\"mi\">1<\/span><span class=\"p\">],<\/span> \r\n  <span class=\"n\">labels<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"s2\">\"\u03bc-2\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc-\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc\"<\/span><span class=\"p\">,<\/span>\r\n            <span class=\"s2\">\"\u03bc+2\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc+\u03c3\"<\/span><span class=\"p\">])<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">plot<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s1\">'gray'<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"c1\"># \u9818\u57df\u306b\u30e9\u30d9\u30eb\u3092\u3064\u3051\u308b<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.15<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig2<\/span><span class=\"p\">[<\/span><span class=\"mi\">1<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"white\"<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"c1\"># -2\u03c3 \u304b\u3089 +2\u03c3 \u307e\u3067\u3092\u5857\u308a\u3064\u3076\u3059<\/span>\r\n<span class=\"n\">x2<\/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=\"o\">-<\/span><span class=\"mi\">2<\/span><span class=\"p\">,<\/span> <span class=\"mi\">2<\/span><span class=\"p\">,<\/span> <span class=\"mi\">101<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">fill_between<\/span><span class=\"p\">(<\/span><span class=\"n\">x2<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x2<\/span><span class=\"p\">),<\/span> <span class=\"n\">fc<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/span><span class=\"p\">[<\/span><span class=\"mi\">1<\/span><span class=\"p\">])<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">plot<\/span><span class=\"p\">(<\/span><span class=\"n\">x2<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x2<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/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\"><\/div>\n<div class=\"output_svg output_subarea \">\n<p><!--?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?--><\/p>\n<p><!-- Created with matplotlib (https:\/\/matplotlib.org\/) --><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-5573\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/pbunpu08b.svg\" alt=\"\" width=\"640\" height=\"240\" \/><\/p>\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[14]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"c1\"># \u67a0\u3068\u7e26\u8ef8\u76ee\u76db\u3092\u975e\u8868\u793a\u306b<\/span>\r\n<span class=\"n\">fig<\/span><span class=\"p\">,<\/span> <span class=\"n\">ax<\/span> <span class=\"o\">=<\/span> <span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">subplots<\/span><span class=\"p\">()<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">spines<\/span><span class=\"p\">[<\/span><span class=\"s1\">'right'<\/span><span class=\"p\">]<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">spines<\/span><span class=\"p\">[<\/span><span class=\"s1\">'left'<\/span><span class=\"p\">]<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">spines<\/span><span class=\"p\">[<\/span><span class=\"s1\">'top'<\/span><span class=\"p\">]<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">axes<\/span><span class=\"o\">.<\/span><span class=\"n\">yaxis<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xlim<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">ylim<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">0.41<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xticks<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"n\">ticks<\/span><span class=\"o\">=<\/span><span class=\"p\">[<\/span><span class=\"o\">-<\/span><span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"mi\">2<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"mi\">1<\/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=\"mi\">2<\/span><span class=\"p\">,<\/span> <span class=\"mi\">1<\/span><span class=\"p\">],<\/span> \r\n  <span class=\"n\">labels<\/span><span class=\"o\">=<\/span><span class=\"p\">[<\/span><span class=\"s2\">\"\u03bc-3\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc-2\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc-\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc\"<\/span><span class=\"p\">,<\/span>\r\n          <span class=\"s2\">\"\u03bc+3\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc+2\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc+\u03c3\"<\/span><span class=\"p\">])<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">plot<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s1\">'gray'<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"c1\"># \u9818\u57df\u306b\u30e9\u30d9\u30eb\u3092\u3064\u3051\u308b<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.15<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig2<\/span><span class=\"p\">[<\/span><span class=\"mi\">2<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"black\"<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"c1\"># -3\u03c3 \u304b\u3089 +3\u03c3 \u307e\u3067\u3092\u5857\u308a\u3064\u3076\u3059<\/span>\r\n<span class=\"n\">x3<\/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=\"o\">-<\/span><span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"mi\">101<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">fill_between<\/span><span class=\"p\">(<\/span><span class=\"n\">x3<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x3<\/span><span class=\"p\">),<\/span> <span class=\"n\">fc<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/span><span class=\"p\">[<\/span><span class=\"mi\">2<\/span><span class=\"p\">])<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">plot<\/span><span class=\"p\">(<\/span><span class=\"n\">x3<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x3<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/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\"><\/div>\n<div class=\"output_svg output_subarea \">\n<p><!--?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?--><\/p>\n<p><!-- Created with matplotlib (https:\/\/matplotlib.org\/) --><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-5574\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/pbunpu09b.svg\" alt=\"\" width=\"640\" height=\"240\" \/><\/p>\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=\"c1\"># \u67a0\u3068\u7e26\u8ef8\u76ee\u76db\u3092\u975e\u8868\u793a\u306b<\/span>\r\n<span class=\"n\">fig<\/span><span class=\"p\">,<\/span> <span class=\"n\">ax<\/span> <span class=\"o\">=<\/span> <span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">subplots<\/span><span class=\"p\">()<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">spines<\/span><span class=\"p\">[<\/span><span class=\"s1\">'right'<\/span><span class=\"p\">]<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">spines<\/span><span class=\"p\">[<\/span><span class=\"s1\">'left'<\/span><span class=\"p\">]<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">spines<\/span><span class=\"p\">[<\/span><span class=\"s1\">'top'<\/span><span class=\"p\">]<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">ax<\/span><span class=\"o\">.<\/span><span class=\"n\">axes<\/span><span class=\"o\">.<\/span><span class=\"n\">yaxis<\/span><span class=\"o\">.<\/span><span class=\"n\">set_visible<\/span><span class=\"p\">(<\/span><span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xlim<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">ylim<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">0.41<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">xticks<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"n\">ticks<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"o\">-<\/span><span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"mi\">2<\/span><span class=\"p\">,<\/span> <span class=\"o\">-<\/span><span class=\"mi\">1<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mi\">4<\/span><span class=\"p\">,<\/span> <span class=\"mi\">3<\/span><span class=\"p\">,<\/span> <span class=\"mi\">2<\/span><span class=\"p\">,<\/span> <span class=\"mi\">1<\/span><span class=\"p\">],<\/span> \r\n  <span class=\"n\">labels<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"s2\">\"\u03bc-4\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc-3\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc-2\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc-\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc\"<\/span><span class=\"p\">,<\/span>\r\n            <span class=\"s2\">\"\u03bc+4\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc+3\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc+2\u03c3\"<\/span><span class=\"p\">,<\/span> <span class=\"s2\">\"\u03bc+\u03c3\"<\/span><span class=\"p\">])<\/span>\r\n\r\n<span class=\"c1\"># \u9818\u57df\u306b\u30e9\u30d9\u30eb\u3092\u3064\u3051\u308b<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">text<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span><span class=\"mf\">0.15<\/span><span class=\"p\">,<\/span> <span class=\"n\">sig2<\/span><span class=\"p\">[<\/span><span class=\"mi\">3<\/span><span class=\"p\">],<\/span> <span class=\"n\">ha<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'center'<\/span><span class=\"p\">,<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"s2\">\"black\"<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"c1\"># -4\u03c3 \u304b\u3089 +4\u03c3 \u307e\u3067\u3092\u5857\u308a\u3064\u3076\u3059<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">fill_between<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">fc<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/span><span class=\"p\">[<\/span><span class=\"mi\">3<\/span><span class=\"p\">])<\/span>\r\n<span class=\"n\">plt<\/span><span class=\"o\">.<\/span><span class=\"n\">plot<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">,<\/span> <span class=\"n\">f<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"p\">),<\/span> <span class=\"n\">color<\/span><span class=\"o\">=<\/span><span class=\"n\">c<\/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\"><\/div>\n<div class=\"output_svg output_subarea \">\n<p><!--?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?--><\/p>\n<p><!-- Created with matplotlib (https:\/\/matplotlib.org\/) --><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-5575\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/pbunpu10b.svg\" alt=\"\" width=\"640\" height=\"240\" \/><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"excerpt":{"rendered":"<p>\u300cgnuplot \u3067\u6b63\u898f\u5206\u5e03\u3092\u03c3\u3054\u3068\u306b\u5857\u308a\u308f\u3051\u308b\u300d\u306e Python \u7248\u3002<\/p>\n<p>\u4ee5\u4e0b\u3092\u53c2\u8003\u306b\uff0c\u6a19\u6e96\u6b63\u898f\u5206\u5e03\u304c\u3082\u3064\u78ba\u7387\u5bc6\u5ea6\u95a2\u6570\u306e\u30b0\u30e9\u30d5\u3092\u63cf\u3044\u3066\u307f\u308b\u3002<\/p><p><a class=\"more-link btn\" href=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/5530\/\">\u7d9a\u304d\u3092\u8aad\u3080<\/a><\/p>\n<ul>\n<li>\u6b63\u898f\u5206\u5e03 &#8211; Wikipedia<\/li>\n<li>\u8aa4\u5dee\u95a2\u6570 &#8211; Wikipedia<\/li>\n<\/ul>\n","protected":false},"author":33,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"inline_featured_image":false,"footnotes":""},"categories":[13],"tags":[],"class_list":["post-5530","post","type-post","status-publish","format-standard","hentry","category-matplotlib","nodate","item-wrap"],"aioseo_notices":[],"_links":{"self":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/posts\/5530","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/types\/post"}],"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=5530"}],"version-history":[{"count":7,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/posts\/5530\/revisions"}],"predecessor-version":[{"id":5605,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/posts\/5530\/revisions\/5605"}],"wp:attachment":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/media?parent=5530"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/categories?post=5530"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/tags?post=5530"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}