{"id":3873,"date":"2022-09-12T17:24:40","date_gmt":"2022-09-12T08:24:40","guid":{"rendered":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/?p=3873"},"modified":"2023-03-28T17:58:20","modified_gmt":"2023-03-28T08:58:20","slug":"gnuplot-%e3%81%a7%e5%86%86%e9%9b%bb%e6%b5%81%e3%81%ab%e3%82%88%e3%82%8b%e7%a3%81%e5%a0%b4%e3%82%92%e6%95%b0%e5%80%a4%e7%9a%84%e3%81%ab%e8%a7%a3%e3%81%84%e3%81%a6%e6%8f%8f%e3%81%8f","status":"publish","type":"post","link":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/3873\/","title":{"rendered":"gnuplot \u3067\u5186\u96fb\u6d41\u306b\u3088\u308b\u78c1\u5834\u3092\u6570\u5024\u7684\u306b\u89e3\u3044\u3066\u63cf\u304f"},"content":{"rendered":"<p><!--more--><\/p>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<p>\u53c2\u8003\uff1a<\/p>\n<ul>\n<li><a href=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/%E9%9B%BB%E7%A3%81%E6%B0%97%E5%AD%A6-i\/%E9%9D%99%E7%A3%81%E5%A0%B4%EF%BC%9A%E9%9B%BB%E6%B5%81%E5%AF%86%E5%BA%A6%E3%81%8B%E3%82%89%E7%9B%B4%E6%8E%A5%E9%9D%99%E7%A3%81%E5%A0%B4%E3%82%92%E6%B1%82%E3%82%81%E3%82%8B\/#y_0_xz\">\u9759\u78c1\u5834\uff1a\u96fb\u6d41\u5bc6\u5ea6\u304b\u3089\u76f4\u63a5\u9759\u78c1\u5834\u3092\u6c42\u3081\u308b<\/a><\/li>\n<\/ul>\n<p>gnuplot \u3060\u3051\u3067\u30d7\u30ed\u30b0\u30e9\u30df\u30f3\u30b0\u3057\u305f\u308a\uff0c\u6570\u5024\u8a08\u7b97\u3057\u305f\u308a\u3059\u308b\u969b\u306e\u53c2\u8003\u3068\u3057\u3066<\/p>\n<ul>\n<li>\u300c<a href=\"https:\/\/home.hirosaki-u.ac.jp\/jupyter\/gnuplot-prog\/\">\u306f\u3057\u3099\u3081\u3066\u306e gnuplot \u30d5\u309a\u30ed\u30af\u3099\u30e9\u30df\u30f3\u30af\u3099 &#8211; \u5f18\u524d\u5927\u5b66 Home Sweet Home<\/a>\u300d<\/li>\n<li>\u300c<a href=\"https:\/\/home.hirosaki-u.ac.jp\/jupyter\/gnuplot-sk\/\">gnuplot \u3066\u3099\u6570\u5024\u89e3\u6790 \u2013 \u5f18\u524d\u5927\u5b66 Home Sweet Home<\/a>\u300d<\/li>\n<\/ul>\n<hr \/>\n<\/div>\n<h3 id=\"\u5186\u96fb\u6d41\u304c\u3064\u304f\u308b-xz-\u5e73\u9762\u5185\u306e\u78c1\u5834\">\u5186\u96fb\u6d41\u304c\u3064\u304f\u308b xz \u5e73\u9762\u5185\u306e\u78c1\u5834<\/h3>\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<p>\\begin{eqnarray}<br \/>\nB_x(x,0,z) &amp;=&amp; \\frac{I}{4\\pi \\varepsilon_0 c^2}<br \/>\n\\int_0^{2\\pi} \\frac{a z \\cos \\phi&#8217; }{\\left\\{ x^2 + z^2 + a^2 &#8211; 2 a x \\cos\\phi&#8217; \\right\\}^{\\frac{3}{2}}}\\,d\\phi&#8217; \\\\ \\ \\\\<br \/>\nB_y(x,0,z) &amp;=&amp; 0 \\\\ \\ \\\\<br \/>\nB_z(x,0,z) &amp;=&amp; \\frac{I }{4\\pi \\varepsilon_0 c^2}<br \/>\n\\int_0^{2\\pi} \\frac{a^2 &#8211; a x \\cos\\phi&#8217;}{\\left\\{x^2 + z^2 + a^2 &#8211; 2 a x \\cos\\phi&#8217; \\right\\}^{\\frac{3}{2}}}\\, d\\phi&#8217;<br \/>\n\\end{eqnarray}<\/p>\n<p>\u9069\u5b9c\uff0c\u5b9a\u6570\u90e8\u5206\u3092\u898f\u683c\u5316\u3057\u3066<\/p>\n<p>\\begin{eqnarray}<br \/>\nB_x(x,z) &amp;\\Rightarrow&amp;<br \/>\n\\int_0^{2\\pi} \\frac{a z \\cos \\phi&#8217; }{\\left\\{ x^2 + z^2 + a^2 &#8211; 2 a x \\cos\\phi&#8217; \\right\\}^{\\frac{3}{2}}}\\,d\\phi&#8217; \\\\ \\ \\\\<br \/>\nB_y(x,z) &amp;=&amp; 0 \\\\ \\ \\\\<br \/>\nB_z(x,z) &amp;\\Rightarrow&amp;<br \/>\n\\int_0^{2\\pi} \\frac{a^2 &#8211; a x \\cos\\phi&#8217;}{\\left\\{x^2 + z^2 + a^2 &#8211; 2 a x \\cos\\phi&#8217; \\right\\}^{\\frac{3}{2}}}\\, d\\phi&#8217;<br \/>\n\\end{eqnarray}<\/p>\n<p>\u78c1\u5834\u30d9\u30af\u30c8\u30eb\u304c\u78c1\u529b\u7dda\u306e\u63a5\u30d9\u30af\u30c8\u30eb\u3067\u3042\u308b\u3053\u3068\u304b\u3089\uff0c\u5b9a\u7a4d\u5206\u3092\u542b\u3080\u5fae\u5206\u65b9\u7a0b\u5f0f<\/p>\n<p>\\begin{eqnarray}<br \/>\nB_x(x,z) &amp;=&amp; \\int_0^{2\\pi} f_1(x,z,\\phi)\\,d\\phi \\\\<br \/>\nB_z(x,z) &amp;=&amp;\\int_0^{2\\pi} f_2(x,z,\\phi)\\,d\\phi \\\\<br \/>\nf_1(x,z,\\phi) &amp;\\equiv&amp; \\frac{a z \\cos \\phi}{\\left\\{ x^2 + z^2 + a^2 &#8211; 2 a x \\cos\\phi\\right\\}^{\\frac{3}{2}}} \\\\<br \/>\nf_2(x,z,\\phi) &amp;\\equiv&amp; \\frac{a^2 &#8211; a x \\cos\\phi}{\\left\\{x^2 + z^2 + a^2 &#8211; 2 a x \\cos\\phi \\right\\}^{\\frac{3}{2}}}\\ \\\\<br \/>\nB(x, z) &amp;\\equiv&amp; \\sqrt{B_x(x,z)^2 + B_z(x,z)^2} \\\\<br \/>\nF_1(x, z) &amp;\\equiv&amp; \\frac{B_x(x,z)}{B(x,z)} \\\\<br \/>\nF_2(x, z) &amp;\\equiv&amp; \\frac{B_z(x,z)}{B(x,z)} \\\\<br \/>\n\\frac{dx}{ds} &amp;=&amp; F_1(x, z)\\\\<br \/>\n\\frac{dz}{ds} &amp;=&amp; F_2(x, z)<br \/>\n\\end{eqnarray}<\/p>\n<p>\u3092\u6570\u5024\u7684\u306b\u89e3\u3044\u3066\u78c1\u529b\u7dda\u3092\u6c42\u3081\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<h3 id=\"\u5186\u96fb\u6d41\u306b\u3088\u308b\u78c1\u5834\uff1aplot-\u7248\">\u5186\u96fb\u6d41\u306b\u3088\u308b\u78c1\u5834\uff1aplot \u7248<\/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[1]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-gnuplot\">\n<pre><span class=\"err\">%gnu<\/span><span class=\"k\">plot<\/span> <span class=\"n\">inline<\/span> <span class=\"n\">svg<\/span> <span class=\"n\">size<\/span> <span class=\"mi\">480<\/span><span class=\"o\">,<\/span><span class=\"mi\">480<\/span> <span class=\"n\">fixed<\/span> <span class=\"n\">enhanced<\/span> <span class=\"n\">font<\/span> <span class=\"s\">'Noto Sans CJK JP,14'<\/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-gnuplot\">\n<pre><span class=\"c\"># \u88ab\u7a4d\u5206\u95a2\u6570<\/span>\r\n<span class=\"nf\">f1<\/span>(x,z,ph) <span class=\"o\">=<\/span> <span class=\"n\">a<\/span><span class=\"o\">*<\/span><span class=\"n\">z<\/span><span class=\"o\">*<\/span><span class=\"nf\">cos<\/span><span class=\"p\">(<\/span><span class=\"n\">ph<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"o\">+<\/span><span class=\"n\">z<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"o\">+<\/span><span class=\"n\">a<\/span><span class=\"o\">**<\/span><span class=\"mi\">2-2<\/span><span class=\"o\">*<\/span><span class=\"n\">a<\/span><span class=\"o\">*<\/span><span class=\"n\">x<\/span><span class=\"o\">*<\/span><span class=\"nf\">cos<\/span><span class=\"p\">(<\/span><span class=\"n\">ph<\/span><span class=\"p\">))<\/span><span class=\"o\">**<\/span><span class=\"mi\">3<\/span>\r\n<span class=\"nf\">f2<\/span>(x,z,ph) <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"n\">a<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"o\">-<\/span><span class=\"n\">a<\/span><span class=\"o\">*<\/span><span class=\"n\">x<\/span><span class=\"o\">*<\/span><span class=\"nf\">cos<\/span><span class=\"p\">(<\/span><span class=\"n\">ph<\/span><span class=\"p\">))<\/span><span class=\"o\">\/<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"o\">+<\/span><span class=\"n\">z<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"o\">+<\/span><span class=\"n\">a<\/span><span class=\"o\">**<\/span><span class=\"mi\">2-2<\/span><span class=\"o\">*<\/span><span class=\"n\">a<\/span><span class=\"o\">*<\/span><span class=\"n\">x<\/span><span class=\"o\">*<\/span><span class=\"nf\">cos<\/span><span class=\"p\">(<\/span><span class=\"n\">ph<\/span><span class=\"p\">))<\/span><span class=\"o\">**<\/span><span class=\"mi\">3<\/span>\r\n\r\n<span class=\"c\"># 2*pi \u7a4d\u5206\u533a\u9593\u3002<\/span>\r\n<span class=\"c\"># N  \u5206\u5272\u6570 N \u306f\u5076\u6570\u3068\u3059\u308b\u3053\u3068\uff08\u30b7\u30f3\u30d7\u30bd\u30f3\u6cd5\u306e\u6c7a\u3081\u4e8b\uff09<\/span>\r\n<span class=\"c\"># \u4e09\u9805\u6f14\u7b97\u5b50\u306b\u3088\u308b\u518d\u5e30\u7684\u5b9a\u7fa9\u95a2\u6570\u3067\u548c\u3092\u3068\u308b\u3002<\/span>\r\n\r\n<span class=\"nf\">simp_f1<\/span>(x, z, N, i) <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"n\">i<\/span> <span class=\"o\">==<\/span> <span class=\"mi\">2<\/span> <span class=\"o\">?<\/span> \\\r\n  <span class=\"p\">(<\/span><span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">pi<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"p\">(<\/span><span class=\"n\">N<\/span><span class=\"o\">*<\/span><span class=\"mf\">3.0<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span> \\\r\n    <span class=\"p\">(<\/span><span class=\"nf\">f1<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"o\">,<\/span><span class=\"p\">(<\/span><span class=\"n\">i<\/span><span class=\"mi\">-2<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">pi<\/span><span class=\"o\">\/<\/span><span class=\"n\">N<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"mf\">4.0<\/span><span class=\"o\">*<\/span><span class=\"nf\">f1<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"o\">,<\/span><span class=\"p\">(<\/span><span class=\"n\">i<\/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=\"n\">pi<\/span><span class=\"o\">\/<\/span><span class=\"n\">N<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"nf\">f1<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"o\">,<\/span><span class=\"n\">i<\/span><span class=\"o\">*<\/span><span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">pi<\/span><span class=\"o\">\/<\/span><span class=\"n\">N<\/span><span class=\"p\">))<\/span> <span class=\"o\">:<\/span> \\\r\n  <span class=\"nf\">simp_f1<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">z<\/span><span class=\"o\">,<\/span> <span class=\"n\">N<\/span><span class=\"o\">,<\/span> <span class=\"n\">i<\/span><span class=\"mi\">-2<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> \\\r\n  <span class=\"p\">(<\/span><span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">pi<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"p\">(<\/span><span class=\"n\">N<\/span><span class=\"o\">*<\/span><span class=\"mf\">3.0<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span> \\\r\n    <span class=\"p\">(<\/span><span class=\"nf\">f1<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"o\">,<\/span><span class=\"p\">(<\/span><span class=\"n\">i<\/span><span class=\"mi\">-2<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">pi<\/span><span class=\"o\">\/<\/span><span class=\"n\">N<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"mf\">4.0<\/span><span class=\"o\">*<\/span><span class=\"nf\">f1<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"o\">,<\/span><span class=\"p\">(<\/span><span class=\"n\">i<\/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=\"n\">pi<\/span><span class=\"o\">\/<\/span><span class=\"n\">N<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"nf\">f1<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"o\">,<\/span><span class=\"n\">i<\/span><span class=\"o\">*<\/span><span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">pi<\/span><span class=\"o\">\/<\/span><span class=\"n\">N<\/span><span class=\"p\">)))<\/span>\r\n\r\n<span class=\"nf\">simp_f2<\/span>(x, z, N, i) <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"n\">i<\/span> <span class=\"o\">==<\/span> <span class=\"mi\">2<\/span> <span class=\"o\">?<\/span> \\\r\n  <span class=\"p\">(<\/span><span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">pi<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"p\">(<\/span><span class=\"n\">N<\/span><span class=\"o\">*<\/span><span class=\"mf\">3.0<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span> \\\r\n    <span class=\"p\">(<\/span><span class=\"nf\">f2<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"o\">,<\/span><span class=\"p\">(<\/span><span class=\"n\">i<\/span><span class=\"mi\">-2<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">pi<\/span><span class=\"o\">\/<\/span><span class=\"n\">N<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"mf\">4.0<\/span><span class=\"o\">*<\/span><span class=\"nf\">f2<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"o\">,<\/span><span class=\"p\">(<\/span><span class=\"n\">i<\/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=\"n\">pi<\/span><span class=\"o\">\/<\/span><span class=\"n\">N<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"nf\">f2<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"o\">,<\/span><span class=\"n\">i<\/span><span class=\"o\">*<\/span><span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">pi<\/span><span class=\"o\">\/<\/span><span class=\"n\">N<\/span><span class=\"p\">))<\/span> <span class=\"o\">:<\/span> \\\r\n  <span class=\"nf\">simp_f2<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">z<\/span><span class=\"o\">,<\/span> <span class=\"n\">N<\/span><span class=\"o\">,<\/span> <span class=\"n\">i<\/span><span class=\"mi\">-2<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> \\\r\n  <span class=\"p\">(<\/span><span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">pi<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"p\">(<\/span><span class=\"n\">N<\/span><span class=\"o\">*<\/span><span class=\"mf\">3.0<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span> \\\r\n    <span class=\"p\">(<\/span><span class=\"nf\">f2<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"o\">,<\/span><span class=\"p\">(<\/span><span class=\"n\">i<\/span><span class=\"mi\">-2<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">pi<\/span><span class=\"o\">\/<\/span><span class=\"n\">N<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"mf\">4.0<\/span><span class=\"o\">*<\/span><span class=\"nf\">f2<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"o\">,<\/span><span class=\"p\">(<\/span><span class=\"n\">i<\/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=\"n\">pi<\/span><span class=\"o\">\/<\/span><span class=\"n\">N<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"nf\">f2<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"o\">,<\/span><span class=\"n\">i<\/span><span class=\"o\">*<\/span><span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">pi<\/span><span class=\"o\">\/<\/span><span class=\"n\">N<\/span><span class=\"p\">)))<\/span>\r\n\r\n<span class=\"c\"># \u5186\u96fb\u6d41\u306e\u534a\u5f84<\/span>\r\n<span class=\"nv\">a<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">1.<\/span>\r\n\r\n<span class=\"c\"># \u5b9f\u969b\u306e\u8a08\u7b97\u306e\u969b\u306f\uff0csimp_f(a, b, N, N) \u3068\u540c\u3058 N \u30922\u56de\u66f8\u304f<\/span>\r\n<span class=\"c\"># N = 40 \u304f\u3089\u3044\u3067\u3088\u3055\u305d\u3046<\/span>\r\n\r\n<span class=\"nf\">Bx<\/span>(x,z) <span class=\"o\">=<\/span> <span class=\"nf\">simp_f1<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">z<\/span><span class=\"o\">,<\/span> <span class=\"mi\">40<\/span><span class=\"o\">,<\/span> <span class=\"mi\">40<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">Bz<\/span>(x,z) <span class=\"o\">=<\/span> <span class=\"nf\">simp_f2<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">z<\/span><span class=\"o\">,<\/span> <span class=\"mi\">40<\/span><span class=\"o\">,<\/span> <span class=\"mi\">40<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"nf\">B<\/span>(x,z) <span class=\"o\">=<\/span> <span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"nf\">Bx<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"nf\">Bz<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"nf\">F1<\/span>(x,z) <span class=\"o\">=<\/span> <span class=\"nf\">Bx<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"nf\">B<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">F2<\/span>(x,z) <span class=\"o\">=<\/span> <span class=\"nf\">Bz<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"nf\">B<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span><span class=\"n\">z<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"c\"># \u8a08\u7b97\u7bc4\u56f2<\/span>\r\n<span class=\"nv\">xmin<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"nv\">xmax<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">3<\/span>\r\n<span class=\"nv\">ymin<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"nv\">ymax<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">3<\/span>\r\n\r\n<span class=\"c\"># \u30eb\u30f3\u30b2\u30fb\u30af\u30c3\u30bf\u6cd5\u306e\u30d1\u30e9\u30e1\u30fc\u30bf<\/span>\r\n<span class=\"k\">h<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.05<\/span>\r\n<span class=\"nv\">Nend<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">100<\/span>\r\n<span class=\"err\">array<\/span> <span class=\"err\">X[Nend]<\/span>\r\n<span class=\"err\">array<\/span> <span class=\"err\">Y[Nend]<\/span>\r\n\r\n<span class=\"c\"># \u30c7\u30fc\u30bf\u30d5\u30a1\u30a4\u30eb\u306e\u524a\u9664<\/span>\r\n<span class=\"c\"># \u3042\u3068\u3067 append \u3067\u958b\u304f\u306e\u3067<\/span>\r\n<span class=\"err\">!<\/span> <span class=\"err\">rm<\/span> <span class=\"err\">-<\/span><span class=\"k\">f<\/span> <span class=\"n\">en<\/span><span class=\"o\">.<\/span><span class=\"n\">txt<\/span>\r\n<span class=\"err\">!<\/span> <span class=\"err\">rm<\/span> <span class=\"err\">-<\/span><span class=\"k\">f<\/span> <span class=\"n\">en<\/span><span class=\"o\">-<\/span><span class=\"n\">line<\/span><span class=\"o\">.<\/span><span class=\"n\">txt<\/span>\r\n\r\n<span class=\"c\"># \u521d\u671f\u6761\u4ef6<\/span>\r\n<span class=\"c\"># x \u8ef8\u4e0a\u304b\u3089\u306f\u3058\u3081\u308b<\/span>\r\n<span class=\"err\">do<\/span> <span class=\"err\">for<\/span> <span class=\"err\">[<\/span><span class=\"nv\">j<\/span><span class=\"o\">=<\/span><span class=\"mi\">1<\/span><span class=\"o\">:<\/span><span class=\"mi\">4<\/span><span class=\"p\">]{<\/span>\r\n  <span class=\"nv\">x0<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.2<\/span><span class=\"o\">*<\/span><span class=\"n\">j<\/span>\r\n  <span class=\"nv\">y0<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n\r\n  <span class=\"c\"># 4\u6b21\u306e\u30eb\u30f3\u30b2\u30fb\u30af\u30c3\u30bf\u6cd5<\/span>\r\n  <span class=\"err\">do<\/span> <span class=\"err\">for<\/span> <span class=\"err\">[<\/span><span class=\"nv\">i<\/span><span class=\"o\">=<\/span><span class=\"mi\">1<\/span><span class=\"o\">:<\/span> <span class=\"n\">Nend<\/span><span class=\"p\">]{<\/span>\r\n    <span class=\"err\">X[i]<\/span> <span class=\"err\">=<\/span> <span class=\"err\">x0<\/span>\r\n    <span class=\"err\">Y[i]<\/span> <span class=\"err\">=<\/span> <span class=\"err\">y0<\/span>\r\n    <span class=\"nv\">k1<\/span> <span class=\"o\">=<\/span> <span class=\"n\">h<\/span><span class=\"o\">*<\/span><span class=\"nf\">F1<\/span><span class=\"p\">(<\/span><span class=\"n\">x0<\/span><span class=\"o\">,<\/span> <span class=\"n\">y0<\/span><span class=\"p\">)<\/span>\r\n    <span class=\"nv\">m1<\/span> <span class=\"o\">=<\/span> <span class=\"n\">h<\/span><span class=\"o\">*<\/span><span class=\"nf\">F2<\/span><span class=\"p\">(<\/span><span class=\"n\">x0<\/span><span class=\"o\">,<\/span> <span class=\"n\">y0<\/span><span class=\"p\">)<\/span>\r\n    <span class=\"nv\">k2<\/span> <span class=\"o\">=<\/span> <span class=\"n\">h<\/span><span class=\"o\">*<\/span><span class=\"nf\">F1<\/span><span class=\"p\">(<\/span><span class=\"n\">x0<\/span><span class=\"o\">+<\/span><span class=\"mf\">0.5<\/span><span class=\"o\">*<\/span><span class=\"n\">k1<\/span><span class=\"o\">,<\/span> <span class=\"n\">y0<\/span><span class=\"o\">+<\/span><span class=\"mf\">0.5<\/span><span class=\"o\">*<\/span><span class=\"n\">m1<\/span><span class=\"p\">)<\/span>\r\n    <span class=\"nv\">m2<\/span> <span class=\"o\">=<\/span> <span class=\"n\">h<\/span><span class=\"o\">*<\/span><span class=\"nf\">F2<\/span><span class=\"p\">(<\/span><span class=\"n\">x0<\/span><span class=\"o\">+<\/span><span class=\"mf\">0.5<\/span><span class=\"o\">*<\/span><span class=\"n\">k1<\/span><span class=\"o\">,<\/span> <span class=\"n\">y0<\/span><span class=\"o\">+<\/span><span class=\"mf\">0.5<\/span><span class=\"o\">*<\/span><span class=\"n\">m1<\/span><span class=\"p\">)<\/span>\r\n    <span class=\"nv\">k3<\/span> <span class=\"o\">=<\/span> <span class=\"n\">h<\/span><span class=\"o\">*<\/span><span class=\"nf\">F1<\/span><span class=\"p\">(<\/span><span class=\"n\">x0<\/span><span class=\"o\">+<\/span><span class=\"mf\">0.5<\/span><span class=\"o\">*<\/span><span class=\"n\">k2<\/span><span class=\"o\">,<\/span> <span class=\"n\">y0<\/span><span class=\"o\">+<\/span><span class=\"mf\">0.5<\/span><span class=\"o\">*<\/span><span class=\"n\">m2<\/span><span class=\"p\">)<\/span>\r\n    <span class=\"nv\">m3<\/span> <span class=\"o\">=<\/span> <span class=\"n\">h<\/span><span class=\"o\">*<\/span><span class=\"nf\">F2<\/span><span class=\"p\">(<\/span><span class=\"n\">x0<\/span><span class=\"o\">+<\/span><span class=\"mf\">0.5<\/span><span class=\"o\">*<\/span><span class=\"n\">k2<\/span><span class=\"o\">,<\/span> <span class=\"n\">y0<\/span><span class=\"o\">+<\/span><span class=\"mf\">0.5<\/span><span class=\"o\">*<\/span><span class=\"n\">m2<\/span><span class=\"p\">)<\/span>\r\n    <span class=\"nv\">k4<\/span> <span class=\"o\">=<\/span> <span class=\"n\">h<\/span><span class=\"o\">*<\/span><span class=\"nf\">F1<\/span><span class=\"p\">(<\/span><span class=\"n\">x0<\/span><span class=\"o\">+<\/span><span class=\"n\">k3<\/span><span class=\"o\">,<\/span> <span class=\"n\">y0<\/span><span class=\"o\">+<\/span><span class=\"n\">m3<\/span><span class=\"p\">)<\/span>\r\n    <span class=\"nv\">m4<\/span> <span class=\"o\">=<\/span> <span class=\"n\">h<\/span><span class=\"o\">*<\/span><span class=\"nf\">F2<\/span><span class=\"p\">(<\/span><span class=\"n\">x0<\/span><span class=\"o\">+<\/span><span class=\"n\">k3<\/span><span class=\"o\">,<\/span> <span class=\"n\">y0<\/span><span class=\"o\">+<\/span><span class=\"n\">m3<\/span><span class=\"p\">)<\/span>\r\n    <span class=\"nv\">X1<\/span> <span class=\"o\">=<\/span> <span class=\"n\">x0<\/span> <span class=\"o\">+<\/span> <span class=\"p\">(<\/span><span class=\"n\">k1<\/span> <span class=\"o\">+<\/span> <span class=\"mf\">2.<\/span><span class=\"o\">*<\/span><span class=\"n\">k2<\/span> <span class=\"o\">+<\/span> <span class=\"mf\">2.<\/span><span class=\"o\">*<\/span><span class=\"n\">k3<\/span> <span class=\"o\">+<\/span> <span class=\"n\">k4<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"mi\">6<\/span>\r\n    <span class=\"nv\">Y1<\/span> <span class=\"o\">=<\/span> <span class=\"n\">y0<\/span> <span class=\"o\">+<\/span> <span class=\"p\">(<\/span><span class=\"n\">m1<\/span> <span class=\"o\">+<\/span> <span class=\"mf\">2.<\/span><span class=\"o\">*<\/span><span class=\"n\">m2<\/span> <span class=\"o\">+<\/span> <span class=\"mf\">2.<\/span><span class=\"o\">*<\/span><span class=\"n\">m3<\/span> <span class=\"o\">+<\/span> <span class=\"n\">m4<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"mi\">6<\/span>\r\n    <span class=\"nv\">iend<\/span> <span class=\"o\">=<\/span> <span class=\"n\">i<\/span>\r\n    <span class=\"c\"># \u8a08\u7b97\u7bc4\u56f2\u3092\u9038\u8131\u3057\u305f\u3089 do loop \u304b\u3089 break<\/span>\r\n    <span class=\"k\">if<\/span> <span class=\"p\">(<\/span><span class=\"n\">X1<\/span> <span class=\"o\">&lt;<\/span> <span class=\"n\">xmin<\/span><span class=\"p\">)<\/span> <span class=\"err\">{break}<\/span>\r\n    <span class=\"k\">if<\/span> <span class=\"p\">(<\/span><span class=\"n\">X1<\/span> <span class=\"o\">&gt;<\/span> <span class=\"n\">xmax<\/span><span class=\"p\">)<\/span> <span class=\"err\">{break}<\/span>\r\n    <span class=\"k\">if<\/span> <span class=\"p\">(<\/span><span class=\"n\">Y1<\/span> <span class=\"o\">&lt;<\/span> <span class=\"n\">ymin<\/span><span class=\"p\">)<\/span> <span class=\"err\">{break}<\/span>\r\n    <span class=\"k\">if<\/span> <span class=\"p\">(<\/span><span class=\"n\">Y1<\/span> <span class=\"o\">&gt;<\/span> <span class=\"n\">ymax<\/span><span class=\"p\">)<\/span> <span class=\"err\">{break}<\/span>\r\n    <span class=\"nv\">x0<\/span> <span class=\"o\">=<\/span> <span class=\"n\">X1<\/span>\r\n    <span class=\"nv\">y0<\/span> <span class=\"o\">=<\/span> <span class=\"n\">Y1<\/span>\r\n  <span class=\"err\">}<\/span>\r\n  <span class=\"c\"># \u78c1\u529b\u7dda\u306e\u5ea7\u6a19\u30c7\u30fc\u30bf\u30d5\u30a1\u30a4\u30eb<\/span>\r\n  <span class=\"k\">set<\/span> <span class=\"nb\">print<\/span> <span class=\"s\">\"en-line.txt\"<\/span> <span class=\"n\">append<\/span>\r\n  <span class=\"c\"># \u7b2c1\u8c61\u9650<\/span>\r\n  <span class=\"err\">do<\/span> <span class=\"err\">for<\/span> <span class=\"err\">[<\/span><span class=\"nv\">i<\/span><span class=\"o\">=<\/span><span class=\"mi\">1<\/span><span class=\"o\">:<\/span><span class=\"n\">iend<\/span><span class=\"p\">]{<\/span>\r\n    <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %8.4f\"<\/span><span class=\"o\">,<\/span> <span class=\"n\">X<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">]<\/span><span class=\"o\">,<\/span> <span class=\"n\">Y<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">])<\/span>\r\n  <span class=\"err\">}<\/span>\r\n  <span class=\"k\">if<\/span> <span class=\"p\">(<\/span><span class=\"n\">Y<\/span><span class=\"p\">[<\/span><span class=\"n\">iend<\/span><span class=\"p\">]<\/span> <span class=\"o\">&gt;<\/span> <span class=\"mf\">0.5<\/span><span class=\"p\">)<\/span> <span class=\"err\">{<\/span><span class=\"k\">print<\/span> <span class=\"s\">\"         \"<\/span><span class=\"p\">}<\/span>\r\n  <span class=\"c\"># \u7b2c4\u8c61\u9650<\/span>\r\n  <span class=\"err\">do<\/span> <span class=\"err\">for<\/span> <span class=\"err\">[<\/span><span class=\"nv\">i<\/span><span class=\"o\">=<\/span><span class=\"n\">iend<\/span><span class=\"o\">:<\/span><span class=\"mi\">1<\/span><span class=\"o\">:<\/span><span class=\"mi\">-1<\/span><span class=\"p\">]{<\/span>\r\n    <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %8.4f\"<\/span><span class=\"o\">,<\/span> <span class=\"n\">X<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">]<\/span><span class=\"o\">,<\/span> <span class=\"o\">-<\/span><span class=\"n\">Y<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">])<\/span>\r\n  <span class=\"err\">}<\/span>\r\n  <span class=\"k\">print<\/span> <span class=\"s\">\"         \"<\/span>\r\n  <span class=\"c\"># \u7b2c2\u8c61\u9650<\/span>\r\n  <span class=\"err\">do<\/span> <span class=\"err\">for<\/span> <span class=\"err\">[<\/span><span class=\"nv\">i<\/span><span class=\"o\">=<\/span><span class=\"mi\">1<\/span><span class=\"o\">:<\/span><span class=\"n\">iend<\/span><span class=\"p\">]{<\/span>\r\n    <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %8.4f\"<\/span><span class=\"o\">,<\/span> <span class=\"o\">-<\/span><span class=\"n\">X<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">]<\/span><span class=\"o\">,<\/span> <span class=\"n\">Y<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">])<\/span>\r\n  <span class=\"err\">}<\/span>\r\n  <span class=\"k\">if<\/span> <span class=\"p\">(<\/span><span class=\"n\">Y<\/span><span class=\"p\">[<\/span><span class=\"n\">iend<\/span><span class=\"p\">]<\/span> <span class=\"o\">&gt;<\/span> <span class=\"mf\">0.5<\/span><span class=\"p\">)<\/span> <span class=\"err\">{<\/span><span class=\"k\">print<\/span> <span class=\"s\">\"         \"<\/span><span class=\"p\">}<\/span>\r\n  <span class=\"c\"># \u7b2c3\u8c61\u9650<\/span>\r\n  <span class=\"err\">do<\/span> <span class=\"err\">for<\/span> <span class=\"err\">[<\/span><span class=\"nv\">i<\/span><span class=\"o\">=<\/span><span class=\"n\">iend<\/span><span class=\"o\">:<\/span><span class=\"mi\">1<\/span><span class=\"o\">:<\/span><span class=\"mi\">-1<\/span><span class=\"p\">]{<\/span>\r\n    <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %8.4f\"<\/span><span class=\"o\">,<\/span> <span class=\"o\">-<\/span><span class=\"n\">X<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">]<\/span><span class=\"o\">,<\/span> <span class=\"o\">-<\/span><span class=\"n\">Y<\/span><span class=\"p\">[<\/span><span class=\"n\">i<\/span><span class=\"p\">])<\/span>\r\n  <span class=\"err\">}<\/span>\r\n  <span class=\"k\">print<\/span> <span class=\"s\">\"         \"<\/span>\r\n  <span class=\"k\">set<\/span> <span class=\"nb\">print<\/span>\r\n\r\n  <span class=\"c\"># \u30d9\u30af\u30c8\u30eb\u59cb\u70b9\u306e\u5ea7\u6a19\u30c7\u30fc\u30bf\u30d5\u30a1\u30a4\u30eb<\/span>\r\n  <span class=\"c\"># \u9069\u5b9c\u9593\u5f15\u304f<\/span>\r\n  <span class=\"k\">set<\/span> <span class=\"nb\">print<\/span> <span class=\"s\">\"en.txt\"<\/span> <span class=\"n\">append<\/span>\r\n  <span class=\"c\"># \u7b2c1\u8c61\u9650<\/span>\r\n    <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %8.4f\"<\/span><span class=\"o\">,<\/span> <span class=\"n\">X<\/span><span class=\"p\">[<\/span><span class=\"mi\">15<\/span><span class=\"p\">]<\/span><span class=\"o\">,<\/span> <span class=\"n\">Y<\/span><span class=\"p\">[<\/span><span class=\"mi\">15<\/span><span class=\"p\">])<\/span>\r\n  <span class=\"c\"># \u7b2c2\u8c61\u9650<\/span>\r\n    <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %8.4f\"<\/span><span class=\"o\">,<\/span> <span class=\"o\">-<\/span><span class=\"n\">X<\/span><span class=\"p\">[<\/span><span class=\"mi\">15<\/span><span class=\"p\">]<\/span><span class=\"o\">,<\/span> <span class=\"n\">Y<\/span><span class=\"p\">[<\/span><span class=\"mi\">15<\/span><span class=\"p\">])<\/span>\r\n  <span class=\"c\"># \u7b2c3\u8c61\u9650<\/span>\r\n    <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %8.4f\"<\/span><span class=\"o\">,<\/span> <span class=\"o\">-<\/span><span class=\"n\">X<\/span><span class=\"p\">[<\/span><span class=\"mi\">20<\/span><span class=\"p\">]<\/span><span class=\"o\">,<\/span> <span class=\"o\">-<\/span><span class=\"n\">Y<\/span><span class=\"p\">[<\/span><span class=\"mi\">20<\/span><span class=\"p\">])<\/span>\r\n  <span class=\"c\"># \u7b2c4\u8c61\u9650<\/span>\r\n    <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %8.4f\"<\/span><span class=\"o\">,<\/span> <span class=\"n\">X<\/span><span class=\"p\">[<\/span><span class=\"mi\">20<\/span><span class=\"p\">]<\/span><span class=\"o\">,<\/span> <span class=\"o\">-<\/span><span class=\"n\">Y<\/span><span class=\"p\">[<\/span><span class=\"mi\">20<\/span><span class=\"p\">])<\/span>\r\n  <span class=\"k\">set<\/span> <span class=\"nb\">print<\/span>\r\n<span class=\"err\">}<\/span>\r\n\r\n<span class=\"c\"># \u78c1\u529b\u7dda\u306e\u5ea7\u6a19\u30c7\u30fc\u30bf\u30d5\u30a1\u30a4\u30eb<\/span>\r\n<span class=\"c\"># z \u8ef8\u4e0a\u306e\u30c7\u30fc\u30bf\u8ffd\u52a0<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">print<\/span> <span class=\"s\">\"en-line.txt\"<\/span> <span class=\"n\">append<\/span>\r\n  <span class=\"err\">do<\/span> <span class=\"err\">for<\/span> <span class=\"err\">[<\/span><span class=\"nv\">i<\/span><span class=\"o\">=<\/span><span class=\"mi\">0<\/span><span class=\"o\">:<\/span><span class=\"mi\">100<\/span><span class=\"p\">]<\/span> <span class=\"p\">{<\/span>\r\n    <span class=\"nv\">z<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">-3<\/span> <span class=\"o\">+<\/span> <span class=\"mf\">6.<\/span><span class=\"o\">\/<\/span><span class=\"mi\">100<\/span><span class=\"o\">*<\/span><span class=\"n\">i<\/span>\r\n    <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %8.4f\"<\/span><span class=\"o\">,<\/span> <span class=\"mi\">0<\/span><span class=\"o\">,<\/span> <span class=\"n\">z<\/span><span class=\"p\">)<\/span>\r\n  <span class=\"err\">}<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">print<\/span>\r\n\r\n<span class=\"k\">set<\/span> <span class=\"nb\">print<\/span> <span class=\"s\">\"en.txt\"<\/span> <span class=\"n\">append<\/span>\r\n  <span class=\"err\">do<\/span> <span class=\"err\">for<\/span> <span class=\"err\">[<\/span><span class=\"nv\">i<\/span><span class=\"o\">=<\/span><span class=\"mi\">0<\/span><span class=\"o\">:<\/span><span class=\"mi\">10<\/span><span class=\"o\">:<\/span><span class=\"mi\">5<\/span><span class=\"p\">]<\/span> <span class=\"p\">{<\/span>\r\n    <span class=\"nv\">z<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">-2<\/span> <span class=\"o\">+<\/span> <span class=\"mf\">4.<\/span><span class=\"o\">\/<\/span><span class=\"mi\">10<\/span><span class=\"o\">*<\/span><span class=\"n\">i<\/span>\r\n    <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %8.4f\"<\/span><span class=\"o\">,<\/span> <span class=\"mi\">0<\/span><span class=\"o\">,<\/span> <span class=\"n\">z<\/span><span class=\"p\">)<\/span>\r\n  <span class=\"err\">}<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">print<\/span>\r\n\r\n<span class=\"k\">reset<\/span>\r\n<span class=\"k\">unset<\/span> <span class=\"nb\">xtics<\/span>\r\n<span class=\"k\">unset<\/span> <span class=\"nb\">ytics<\/span>\r\n<span class=\"k\">unset<\/span> <span class=\"nb\">border<\/span>\r\n\r\n<span class=\"k\">set<\/span> <span class=\"nb\">xrange<\/span> <span class=\"p\">[<\/span><span class=\"mi\">-3<\/span><span class=\"o\">:<\/span><span class=\"mi\">3<\/span><span class=\"p\">]<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">yrange<\/span> <span class=\"p\">[<\/span><span class=\"mi\">-3<\/span><span class=\"o\">:<\/span><span class=\"mi\">3<\/span><span class=\"p\">]<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">size<\/span> <span class=\"n\">ratio<\/span> <span class=\"mi\">1<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">zeroaxis<\/span>\r\n\r\n<span class=\"k\">set<\/span> <span class=\"nb\">title<\/span> <span class=\"s\">\"\u5186\u96fb\u6d41\u306b\u3088\u308b\u78c1\u5834\uff08xz\u9762\uff09\"<\/span>\r\n<span class=\"k\">unset<\/span> <span class=\"nb\">colorbox<\/span>\r\n<span class=\"c\"># green.pal \u6539<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">palette<\/span> <span class=\"nf\">defined<\/span> <span class=\"p\">(<\/span>\\\r\n  <span class=\"mi\">0<\/span> <span class=\"s\">'#E5F5E0'<\/span><span class=\"o\">,<\/span>\\\r\n  <span class=\"mi\">1<\/span> <span class=\"s\">'#C7E9C0'<\/span><span class=\"o\">,<\/span>\\\r\n  <span class=\"mi\">2<\/span> <span class=\"s\">'#A1D99B'<\/span><span class=\"o\">,<\/span>\\\r\n  <span class=\"mi\">3<\/span> <span class=\"s\">'#74C476'<\/span><span class=\"o\">,<\/span>\\\r\n  <span class=\"mi\">4<\/span> <span class=\"s\">'#41AB5D'<\/span><span class=\"o\">,<\/span>\\\r\n  <span class=\"mi\">5<\/span> <span class=\"s\">'#238B45'<\/span><span class=\"o\">,<\/span>\\\r\n  <span class=\"mi\">6<\/span> <span class=\"s\">'#005A32'<\/span><span class=\"o\">,<\/span>\\\r\n  <span class=\"mi\">7<\/span> <span class=\"s\">'#002000'<\/span> \\\r\n<span class=\"p\">)<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">key<\/span> <span class=\"n\">inside<\/span> <span class=\"n\">sample<\/span> <span class=\"mi\">1<\/span>\r\n<span class=\"c\"># \u96fb\u7dda\u306e\u4f4d\u7f6e<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">label<\/span> <span class=\"mi\">1<\/span> <span class=\"n\">point<\/span> <span class=\"n\">pt<\/span> <span class=\"mi\">6<\/span> <span class=\"n\">ps<\/span> <span class=\"mi\">1<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">at<\/span> <span class=\"n\">a<\/span><span class=\"o\">,<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">label<\/span> <span class=\"mi\">2<\/span> <span class=\"n\">center<\/span> <span class=\"n\">at<\/span> <span class=\"n\">first<\/span> <span class=\"n\">a<\/span><span class=\"o\">,<\/span> <span class=\"mi\">0<\/span> <span class=\"s\">\"\u00d7\"<\/span> <span class=\"n\">tc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">font<\/span> <span class=\"s\">\",14\"<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">label<\/span> <span class=\"mi\">3<\/span> <span class=\"n\">point<\/span> <span class=\"n\">pt<\/span> <span class=\"mi\">6<\/span> <span class=\"n\">ps<\/span> <span class=\"mi\">1<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">at<\/span> <span class=\"o\">-<\/span><span class=\"n\">a<\/span><span class=\"o\">,<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">label<\/span> <span class=\"mi\">4<\/span> <span class=\"n\">center<\/span> <span class=\"n\">at<\/span> <span class=\"n\">first<\/span> <span class=\"o\">-<\/span><span class=\"n\">a<\/span><span class=\"o\">,<\/span> <span class=\"mi\">0<\/span> <span class=\"s\">\"\u30fb\"<\/span> <span class=\"n\">tc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">font<\/span> <span class=\"s\">\",14\"<\/span>\r\n<span class=\"nv\">scaling<\/span><span class=\"o\">=<\/span><span class=\"mf\">0.05<\/span>\r\n<span class=\"k\">plot<\/span> <span class=\"s\">\"en-line.txt\"<\/span> \\\r\n  <span class=\"nb\">u<\/span> <span class=\"mi\">1<\/span><span class=\"o\">:<\/span><span class=\"mi\">2<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"nf\">log<\/span><span class=\"p\">(<\/span><span class=\"nf\">B<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">)))<\/span> <span class=\"nb\">w<\/span> <span class=\"n\">l<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">3<\/span> <span class=\"n\">lc<\/span> <span class=\"n\">palette<\/span> <span class=\"nb\">notitle<\/span><span class=\"o\">,<\/span> \\\r\n  <span class=\"s\">\"en.txt\"<\/span>   \\\r\n  <span class=\"nb\">u<\/span> <span class=\"mi\">1<\/span><span class=\"o\">:<\/span><span class=\"mi\">2<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"n\">scaling<\/span><span class=\"o\">*<\/span><span class=\"nf\">F1<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">))<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"n\">scaling<\/span><span class=\"o\">*<\/span><span class=\"nf\">F2<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">))<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"nf\">log<\/span><span class=\"p\">(<\/span><span class=\"nf\">B<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">)))<\/span> \\\r\n  <span class=\"nb\">w<\/span> <span class=\"n\">vec<\/span> <span class=\"n\">lc<\/span> <span class=\"n\">palette<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">8<\/span> <span class=\"n\">filled<\/span> <span class=\"n\">head<\/span> <span class=\"nb\">title<\/span> <span class=\"s\">\"\u78c1\u5834\"<\/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><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-full wp-image-3880\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/en-jiba2.svg\" alt=\"\" width=\"480\" height=\"480\" \/><\/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[3]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-gnuplot\">\n<pre><span class=\"k\">set<\/span> <span class=\"nb\">output<\/span> <span class=\"s\">\".\/en-jiba2.svg\"<\/span>\r\n<span class=\"k\">replot<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">output<\/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=\"\u5186\u96fb\u6d41\u56de\u8def\u306b\u3088\u308b\u78c1\u5834\uff1asplot-\u7248\">\u5186\u96fb\u6d41\u56de\u8def\u306b\u3088\u308b\u78c1\u5834\uff1asplot \u7248<\/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-gnuplot\">\n<pre><span class=\"err\">%gnu<\/span><span class=\"k\">plot<\/span> <span class=\"n\">inline<\/span> <span class=\"n\">svg<\/span> <span class=\"n\">size<\/span> <span class=\"mi\">640<\/span><span class=\"o\">,<\/span><span class=\"mi\">640<\/span> <span class=\"n\">fixed<\/span> <span class=\"n\">enhanced<\/span> <span class=\"n\">font<\/span> <span class=\"s\">'Noto Sans CJK JP,14'<\/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-gnuplot\">\n<pre><span class=\"c\"># \u534a\u5186\u306e\u5ea7\u6a19\u30c7\u30fc\u30bf\u30d5\u30a1\u30a4\u30eb<\/span>\r\n<span class=\"c\"># theta \u306e\u5206\u5272\u6570<\/span>\r\n<span class=\"nv\">Nth<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">100<\/span>\r\n<span class=\"c\"># \u534a\u5f84<\/span>\r\n<span class=\"nv\">a<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">1<\/span>\r\n\r\n<span class=\"k\">set<\/span> <span class=\"nb\">print<\/span> <span class=\"s\">\"han-en.txt\"<\/span>\r\n  <span class=\"err\">do<\/span> <span class=\"err\">for<\/span> <span class=\"err\">[<\/span><span class=\"nv\">i<\/span><span class=\"o\">=<\/span><span class=\"mi\">0<\/span><span class=\"o\">:<\/span><span class=\"n\">Nth<\/span><span class=\"p\">]{<\/span>\r\n    <span class=\"nv\">th<\/span> <span class=\"o\">=<\/span> <span class=\"n\">pi<\/span><span class=\"o\">\/<\/span><span class=\"n\">Nth<\/span> <span class=\"o\">*<\/span> <span class=\"n\">i<\/span>\r\n    <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %8.4f\"<\/span><span class=\"o\">,<\/span> <span class=\"n\">a<\/span><span class=\"o\">*<\/span><span class=\"nf\">cos<\/span><span class=\"p\">(<\/span><span class=\"n\">th<\/span><span class=\"p\">)<\/span><span class=\"o\">,<\/span> <span class=\"n\">a<\/span><span class=\"o\">*<\/span><span class=\"nf\">sin<\/span><span class=\"p\">(<\/span><span class=\"n\">th<\/span><span class=\"p\">))<\/span>\r\n  <span class=\"err\">}<\/span>\r\n  <span class=\"k\">print<\/span> <span class=\"s\">\"      \"<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">print<\/span>\r\n\r\n<span class=\"c\"># \u534a\u5186\u67f1\u306e\u96fb\u6d41\u30d9\u30af\u30c8\u30eb\u30c7\u30fc\u30bf\u30d5\u30a1\u30a4\u30eb<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">print<\/span> <span class=\"s\">\"han-en-vec2.txt\"<\/span>\r\n<span class=\"c\">#  print 0, 1, 0, -1, 0, 0<\/span>\r\n  <span class=\"k\">print<\/span> <span class=\"mi\">0<\/span><span class=\"o\">,<\/span> <span class=\"mi\">-1<\/span><span class=\"o\">,<\/span> <span class=\"mi\">0<\/span><span class=\"o\">,<\/span> <span class=\"mi\">1<\/span><span class=\"o\">,<\/span>  <span class=\"mi\">0<\/span><span class=\"o\">,<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">print<\/span>\r\n\r\n<span class=\"k\">reset<\/span>\r\n<span class=\"k\">unset<\/span> <span class=\"nb\">xtics<\/span>\r\n<span class=\"k\">unset<\/span> <span class=\"nb\">ytics<\/span>\r\n<span class=\"k\">unset<\/span> <span class=\"nb\">ztics<\/span>\r\n<span class=\"k\">unset<\/span> <span class=\"nb\">border<\/span>\r\n\r\n<span class=\"k\">set<\/span> <span class=\"nb\">xrange<\/span> <span class=\"p\">[<\/span><span class=\"mi\">-3<\/span><span class=\"o\">:<\/span><span class=\"mi\">3<\/span><span class=\"p\">]<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">yrange<\/span> <span class=\"p\">[<\/span><span class=\"mi\">-3<\/span><span class=\"o\">:<\/span><span class=\"mi\">3<\/span><span class=\"p\">]<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">zrange<\/span> <span class=\"p\">[<\/span><span class=\"mi\">-3<\/span><span class=\"o\">:<\/span><span class=\"mi\">3<\/span><span class=\"p\">]<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">view<\/span> <span class=\"mi\">60<\/span><span class=\"o\">,<\/span><span class=\"mi\">15<\/span><span class=\"o\">,,<\/span>\r\n\r\n<span class=\"k\">set<\/span> <span class=\"nb\">key<\/span> <span class=\"n\">inside<\/span> <span class=\"n\">sample<\/span> <span class=\"mi\">1<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">title<\/span> <span class=\"s\">\"\u5186\u96fb\u6d41\u306b\u3088\u308b\u78c1\u5834\"<\/span>\r\n<span class=\"k\">unset<\/span> <span class=\"nb\">colorbox<\/span>\r\n<span class=\"c\"># green.pal \u6539<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">palette<\/span> <span class=\"nf\">defined<\/span> <span class=\"p\">(<\/span>\\\r\n  <span class=\"mi\">0<\/span> <span class=\"s\">'#E5F5E0'<\/span><span class=\"o\">,<\/span>\\\r\n  <span class=\"mi\">1<\/span> <span class=\"s\">'#C7E9C0'<\/span><span class=\"o\">,<\/span>\\\r\n  <span class=\"mi\">2<\/span> <span class=\"s\">'#A1D99B'<\/span><span class=\"o\">,<\/span>\\\r\n  <span class=\"mi\">3<\/span> <span class=\"s\">'#74C476'<\/span><span class=\"o\">,<\/span>\\\r\n  <span class=\"mi\">4<\/span> <span class=\"s\">'#41AB5D'<\/span><span class=\"o\">,<\/span>\\\r\n  <span class=\"mi\">5<\/span> <span class=\"s\">'#238B45'<\/span><span class=\"o\">,<\/span>\\\r\n  <span class=\"mi\">6<\/span> <span class=\"s\">'#005A32'<\/span><span class=\"o\">,<\/span>\\\r\n  <span class=\"mi\">7<\/span> <span class=\"s\">'#002000'<\/span> \\\r\n<span class=\"p\">)<\/span>\r\n<span class=\"nv\">scaling<\/span><span class=\"o\">=<\/span><span class=\"mf\">0.05<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">parametric<\/span>\r\n<span class=\"k\">splot<\/span> <span class=\"p\">[<\/span><span class=\"nb\">u<\/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=\"nb\">u<\/span><span class=\"o\">,<\/span><span class=\"mi\">0<\/span><span class=\"o\">,<\/span><span class=\"mi\">0<\/span> <span class=\"nf\">dt<\/span><span class=\"p\">(<\/span><span class=\"mi\">5<\/span><span class=\"o\">,<\/span><span class=\"mi\">10<\/span><span class=\"p\">)<\/span> <span class=\"n\">lw<\/span> <span class=\"mf\">0.5<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">'black'<\/span> <span class=\"nb\">notitle<\/span><span class=\"o\">,<\/span> \\\r\n               <span class=\"mi\">0<\/span><span class=\"o\">,<\/span><span class=\"nb\">u<\/span><span class=\"o\">,<\/span><span class=\"mi\">0<\/span> <span class=\"nf\">dt<\/span><span class=\"p\">(<\/span><span class=\"mi\">5<\/span><span class=\"o\">,<\/span><span class=\"mi\">10<\/span><span class=\"p\">)<\/span> <span class=\"n\">lw<\/span> <span class=\"mf\">0.5<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">'black'<\/span> <span class=\"nb\">notitle<\/span><span class=\"o\">,<\/span> \\\r\n  <span class=\"s\">\"han-en.txt\"<\/span> \\\r\n    <span class=\"nb\">u<\/span> <span class=\"mi\">1<\/span><span class=\"o\">:<\/span><span class=\"mi\">2<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">)<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span> <span class=\"nb\">w<\/span> <span class=\"n\">l<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">4<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">'red'<\/span> <span class=\"nb\">notitle<\/span><span class=\"o\">,<\/span> \\\r\n  <span class=\"s\">\"en-line.txt\"<\/span> \\\r\n    <span class=\"nb\">u<\/span> <span class=\"mi\">1<\/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\">2<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"nf\">log<\/span><span class=\"p\">(<\/span><span class=\"nf\">B<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">)))<\/span> <span class=\"nb\">w<\/span> <span class=\"n\">l<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">3<\/span> <span class=\"n\">lc<\/span> <span class=\"n\">palette<\/span> <span class=\"nb\">notitle<\/span><span class=\"o\">,<\/span> \\\r\n  <span class=\"s\">\"en.txt\"<\/span> \\\r\n    <span class=\"nb\">u<\/span> <span class=\"mi\">1<\/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\">2<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"n\">scaling<\/span><span class=\"o\">*<\/span><span class=\"nf\">F1<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">))<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">)<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"n\">scaling<\/span><span class=\"o\">*<\/span><span class=\"nf\">F2<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">))<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"nf\">log<\/span><span class=\"p\">(<\/span><span class=\"nf\">B<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">)))<\/span> \\\r\n    <span class=\"nb\">w<\/span> <span class=\"n\">vec<\/span> <span class=\"n\">filled<\/span> <span class=\"n\">head<\/span> <span class=\"n\">lc<\/span> <span class=\"n\">palette<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">6<\/span> <span class=\"nb\">title<\/span> <span class=\"s\">\"\u78c1\u5834\"<\/span><span class=\"o\">,<\/span> \\\r\n  <span class=\"s\">\"han-en-vec2.txt\"<\/span> \\\r\n    <span class=\"nb\">u<\/span> <span class=\"mi\">1<\/span><span class=\"o\">:<\/span><span class=\"mi\">2<\/span><span class=\"o\">:<\/span><span class=\"mi\">3<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"n\">scaling<\/span><span class=\"o\">*<\/span><span class=\"err\">$<\/span><span class=\"mi\">4<\/span><span class=\"p\">)<\/span><span class=\"o\">:<\/span><span class=\"mi\">5<\/span><span class=\"o\">:<\/span><span class=\"mi\">6<\/span> <span class=\"nb\">w<\/span> <span class=\"n\">vec<\/span> <span class=\"n\">filled<\/span> <span class=\"n\">head<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">8<\/span>  <span class=\"n\">lc<\/span> <span class=\"s\">'red'<\/span> <span class=\"nb\">title<\/span> <span class=\"s\">\"\u96fb\u6d41\"<\/span>  <span class=\"o\">,<\/span> \\\r\n  <span class=\"s\">\"han-en.txt\"<\/span> \\\r\n    <span class=\"nb\">u<\/span> <span class=\"mi\">1<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">)<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span> <span class=\"nb\">w<\/span> <span class=\"n\">l<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">4<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">'red'<\/span> <span class=\"nb\">notitle<\/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><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-full wp-image-3881\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/en-jiba3.svg\" alt=\"\" width=\"640\" height=\"640\" \/><\/p>\n<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[5]:<\/div>\n<div class=\"output_text output_subarea output_execute_result\">\n<pre>\tdummy variable is t for curves, u\/v for surfaces\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[6]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-gnuplot\">\n<pre><span class=\"k\">set<\/span> <span class=\"nb\">output<\/span> <span class=\"s\">\".\/en-jiba3.svg\"<\/span>\r\n<span class=\"k\">replot<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">output<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"excerpt":{"rendered":"","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":[15,19],"tags":[],"class_list":["post-3873","post","type-post","status-publish","format-standard","hentry","category-gnuplot","category-19","nodate","item-wrap"],"aioseo_notices":[],"_links":{"self":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/posts\/3873","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=3873"}],"version-history":[{"count":9,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/posts\/3873\/revisions"}],"predecessor-version":[{"id":3896,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/posts\/3873\/revisions\/3896"}],"wp:attachment":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/media?parent=3873"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/categories?post=3873"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/tags?post=3873"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}