{"id":3702,"date":"2022-09-05T13:09:12","date_gmt":"2022-09-05T04:09:12","guid":{"rendered":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/?p=3702"},"modified":"2023-03-28T17:58:58","modified_gmt":"2023-03-28T08:58:58","slug":"gnuplot-%e3%81%a7%e7%a3%81%e5%8a%9b%e7%b7%9a%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\/3702\/","title":{"rendered":"gnuplot \u3067\u78c1\u529b\u7dda\u3092\u6570\u5024\u7684\u306b\u89e3\u3044\u3066\u63cf\u304f"},"content":{"rendered":"<p>gnuplot \u3060\u3051\u3067\uff0c\u78c1\u529b\u7dda\u3092\u6570\u5024\u7684\u306b\u89e3\u304d\uff0c\u78c1\u529b\u7dda\u306b\u6cbf\u3063\u305f\u78c1\u5834\u30d9\u30af\u30c8\u30eb\u3092\u63cf\u304f\u3002\u78c1\u5834\u30d9\u30af\u30c8\u30eb\u304c\u78c1\u529b\u7dda\u306e\u63a5\u30d9\u30af\u30c8\u30eb\u3067\u3042\u308b\u3053\u3068\u304b\u3089\uff0c\u5fae\u5206\u65b9\u7a0b\u5f0f\u3092\u30eb\u30f3\u30b2\u30fb\u30af\u30c3\u30bf\u6cd5\u3067\u89e3\u304d\uff0c\u6570\u5024\u89e3\u3092\u30d5\u30a1\u30a4\u30eb\u306b\u4fdd\u5b58\u3057\uff0c\u305d\u308c\u3092\u8aad\u307f\u8fbc\u3093\u3067\u78c1\u5834\u30d9\u30af\u30c8\u30eb\u306e\u59cb\u70b9\u306e\u5ea7\u6a19\u3068\u3057\u3066\u30d9\u30af\u30c8\u30eb\u5834\u3092\u63cf\u304f\u3002\u3053\u308c\u3092\u5168\u3066 gnuplot \u3060\u3051\u3067\u3084\u308b\u3002<!--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>\u300c<a href=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/3659\/\" target=\"_blank\" rel=\"noopener\">gnuplot \u3067\u96fb\u6c17\u529b\u7dda\u3092\u6570\u5024\u7684\u306b\u89e3\u3044\u3066\u63cf\u304f<\/a>\u300d<\/li>\n<\/ul>\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>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<h3 id=\"\u5e73\u884c2\u76f4\u7dda\u96fb\u6d41\u304c\u3064\u304f\u308b\u78c1\u5834\">\u5e73\u884c2\u76f4\u7dda\u96fb\u6d41\u304c\u3064\u304f\u308b\u78c1\u5834<\/h3>\n<p>$$\\boldsymbol{B} = \\frac{1}{2\\pi\\varepsilon_0 c^2} \\frac{\\boldsymbol{I}_1\\times \\left(\\boldsymbol{\\rho} &#8211; \\boldsymbol{r}_1\\right)}{|\\boldsymbol{\\rho} &#8211; \\boldsymbol{r}_1|^2}<br \/>\n+ \\frac{1}{2\\pi\\varepsilon_0 c^2} \\frac{\\boldsymbol{I}_2\\times \\left(\\boldsymbol{\\rho} &#8211; \\boldsymbol{r}_2\\right)}{|\\boldsymbol{\\rho} &#8211; \\boldsymbol{r}_2|^2}$$\\begin{eqnarray}<br \/>\n\\boldsymbol{I}_1 &amp;=&amp; (0, 0, \\pm I)\\\\<br \/>\n\\boldsymbol{I}_2 &amp;=&amp; (0, 0, I) \\\\<br \/>\n\\boldsymbol{\\rho} &amp;=&amp; (x, y, 0) \\\\<br \/>\n\\boldsymbol{r}_1 &amp;=&amp; (-1, 0, 0) \\\\<br \/>\n\\boldsymbol{r}_2 &amp;=&amp; (1, 0, 0) \\\\<br \/>\n\\end{eqnarray}<\/p>\n<p>\u3068\u3057\uff0c\u9069\u5b9c\u5b9a\u6570\u90e8\u5206\u3092\u898f\u683c\u5316\u3057\u3066<br \/>\n\\begin{eqnarray}<br \/>\nB_x(x,y) &amp;=&amp; \\pm \\frac{-y}{\\left(x-x_1\\right)^2 + \\left(y-y_1\\right)^2 }<br \/>\n+\\frac{-y}{\\left(x-x_2\\right)^2 + \\left(y-y_2\\right)^2} \\\\<br \/>\nB_y(x,y) &amp;=&amp; \\pm \\frac{x}{\\left(x-x_1\\right)^2 + \\left(y-y_1\\right)^2 }<br \/>\n+\\frac{x}{\\left(x-x_2\\right)^2 + \\left(y-y_2\\right)^2} \\\\<br \/>\n\\end{eqnarray}<\/p>\n<p>\u898f\u683c\u5316\u3057\u305f $\\hat{\\boldsymbol{B}}$ \u306f<\/p>\n<p>\\begin{eqnarray}<br \/>\nB(x,y) &amp;=&amp; \\sqrt{\\left\\{B_x(x,y) \\right\\}^2 + \\left\\{B_y(x,y) \\right\\}^2}\\\\<br \/>\n\\hat{B}_x(x,y) &amp;=&amp; \\frac{B_x(x,y)}{B(x,y)} \\\\<br \/>\n\\hat{B}_y(x,y) &amp;=&amp; \\frac{B_y(x,y)}{B(x,y)} \\\\<br \/>\n\\end{eqnarray}<\/p>\n<ul>\n<li>\u53c2\u8003\uff1a<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\/#i-2\">\u76f4\u7dda\u96fb\u6d41\u306b\u3088\u308b\u78c1\u5834<\/a><\/li>\n<\/ul>\n<h4 id=\"\u5e73\u884c2\u76f4\u7dda\u96fb\u6d41\u306b\u3088\u308b\u78c1\u529b\u7dda\">\u5e73\u884c2\u76f4\u7dda\u96fb\u6d41\u306b\u3088\u308b\u78c1\u529b\u7dda<\/h4>\n<p>\u78c1\u5834\u30d9\u30af\u30c8\u30eb $\\boldsymbol{B}$ \u304c\u78c1\u529b\u7dda\u306e\u63a5\u30d9\u30af\u30c8\u30eb\u3067\u3042\u308b\u304b\u3089\uff0c\u5fae\u5206\u65b9\u7a0b\u5f0f<\/p>\n<p>$$\\frac{d\\boldsymbol{r}}{ds} = \\hat{\\boldsymbol{B}}$$<\/p>\n<p>\u3092\uff08\u30eb\u30f3\u30b2\u30fb\u30af\u30c3\u30bf\u6cd5\u3067\uff09\u6570\u5024\u7684\u306b\u89e3\u3044\u3066\uff0c\u78c1\u529b\u7dda $\\boldsymbol{r}(s)$ \u3092\u6c42\u3081\u308b\u3002<\/p>\n<p>\\begin{eqnarray}<br \/>\n\\frac{dx}{ds} &amp;=&amp; \\hat{B}_x(x, y) \\equiv F_1(x,y) \\\\<br \/>\n\\frac{dy}{ds} &amp;=&amp; \\hat{B}_y(x, y) \\equiv F_2(x,y)<br \/>\n\\end{eqnarray}<\/p>\n<h4 id=\"\u540c\u3058\u5411\u304d\u306b\u96fb\u6d41-$I$-\u304c\u6d41\u308c\u3066\u3044\u308b\u5834\u5408\uff1aplot-\u7248\">\u540c\u3058\u5411\u304d\u306b\u96fb\u6d41 $I$ \u304c\u6d41\u308c\u3066\u3044\u308b\u5834\u5408\uff1aplot \u7248<\/h4>\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\"># \u78c1\u5834 B<\/span>\r\n<span class=\"nf\">Bx<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"o\">-<\/span> <span class=\"n\">y<\/span><span class=\"o\">\/<\/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\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">By<\/span>(x, y) <span class=\"o\">=<\/span>   <span class=\"n\">x<\/span><span class=\"o\">\/<\/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\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"c\"># \u5e73\u884c2\u76f4\u7dda\u96fb\u6d41\u304c\u3064\u304f\u308b\u78c1\u5834<\/span>\r\n<span class=\"c\"># \u96fb\u7dda\u306e\u4f4d\u7f6e<\/span>\r\n<span class=\"nv\">x1<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">-1<\/span>\r\n<span class=\"nv\">y1<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"nv\">x2<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">1<\/span>\r\n<span class=\"nv\">y2<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n\r\n<span class=\"nf\">B2x<\/span>(x, y) <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\">x1<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y1<\/span><span class=\"p\">)<\/span> <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\">x2<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">B2y<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">By<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">-<\/span><span class=\"n\">x1<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y1<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"nf\">By<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">-<\/span><span class=\"n\">x2<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y2<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"c\"># \u78c1\u5834 B \u306e\u5927\u304d\u3055<\/span>\r\n<span class=\"nf\">B<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"nf\">B2x<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"nf\">B2y<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"c\"># \u898f\u683c\u5316\u3055\u308c\u305f hat B<\/span>\r\n<span class=\"nf\">F1<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">B2x<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/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\">y<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">F2<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">B2y<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/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\">y<\/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.02<\/span>\r\n<span class=\"nv\">N<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">300<\/span>\r\n<span class=\"err\">array<\/span> <span class=\"err\">X[N]<\/span>\r\n<span class=\"err\">array<\/span> <span class=\"err\">Y[N]<\/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\">para<\/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\">para<\/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\">2<\/span><span class=\"o\">:<\/span><span class=\"mi\">11<\/span><span class=\"p\">]{<\/span>\r\n  <span class=\"nv\">x0<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">4.<\/span><span class=\"o\">\/<\/span><span class=\"n\">j<\/span> <span class=\"o\">+<\/span> <span class=\"mi\">1<\/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\">N<\/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\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\">\"para-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\">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=\"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=\"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\"># \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=\"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=\"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\"># \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=\"k\">set<\/span> <span class=\"nb\">print<\/span>  \r\n  \r\n  <span class=\"k\">if<\/span> <span class=\"p\">(<\/span><span class=\"n\">j<\/span><span class=\"o\">%<\/span><span class=\"mi\">2<\/span>  <span class=\"o\">==<\/span><span class=\"mi\">0<\/span><span class=\"p\">)<\/span> <span class=\"err\">{<\/span>\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\">\"para.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\">30<\/span><span class=\"o\">:<\/span><span class=\"mi\">30<\/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=\"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=\"n\">iend<\/span><span class=\"mi\">-20<\/span><span class=\"o\">:<\/span><span class=\"n\">iend<\/span><span class=\"mi\">-20<\/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=\"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=\"mi\">30<\/span><span class=\"o\">:<\/span><span class=\"mi\">30<\/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=\"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=\"mi\">-20<\/span><span class=\"o\">:<\/span><span class=\"n\">iend<\/span><span class=\"mi\">-20<\/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\">set<\/span> <span class=\"nb\">print<\/span>\r\n  <span class=\"err\">}<\/span>\r\n<span class=\"err\">}<\/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\">\"\u540c\u3058\u5411\u304d\u306b\u6d41\u308c\u308b\u5e73\u884c2\u76f4\u7dda\u96fb\u6d41\u306b\u3088\u308b\u78c1\u5834\"<\/span>\r\n<span class=\"k\">unset<\/span> <span class=\"nb\">colorbox<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">palette<\/span> <span class=\"n\">cubehelix<\/span> <span class=\"n\">start<\/span> <span class=\"mi\">2<\/span> <span class=\"n\">cycles<\/span> <span class=\"mi\">0<\/span> <span class=\"n\">saturation<\/span> <span class=\"mi\">3<\/span> <span class=\"n\">negative<\/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\">2<\/span>\r\n\r\n<span class=\"c\"># \u77e2\u5370\u306e\u9577\u3055\u8abf\u6574\u7528<\/span>\r\n<span class=\"nv\">scaling<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.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\">2<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">at<\/span> <span class=\"n\">x1<\/span><span class=\"o\">,<\/span> <span class=\"n\">y1<\/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\">x1<\/span><span class=\"o\">,<\/span><span class=\"n\">y1<\/span> <span class=\"s\">\"\u30fb\"<\/span> <span class=\"n\">tc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">font<\/span> <span class=\"s\">\",24\"<\/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\">2<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">at<\/span> <span class=\"n\">x2<\/span><span class=\"o\">,<\/span> <span class=\"n\">y2<\/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=\"n\">x2<\/span><span class=\"o\">,<\/span><span class=\"n\">y2<\/span> <span class=\"s\">\"\u30fb\"<\/span> <span class=\"n\">tc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">font<\/span> <span class=\"s\">\",24\"<\/span>\r\n\r\n<span class=\"k\">plot<\/span> <span class=\"s\">\"para-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\">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\">\"para.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\">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\">6<\/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-3847\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/vec-para1.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\">\".\/vec-para1.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<h4 id=\"\u540c\u3058\u5411\u304d\u306b\u96fb\u6d41-$I$-\u304c\u6d41\u308c\u3066\u3044\u308b\u5834\u5408\uff1asplot-\u7248\">\u540c\u3058\u5411\u304d\u306b\u96fb\u6d41 $I$ \u304c\u6d41\u308c\u3066\u3044\u308b\u5834\u5408\uff1asplot \u7248<\/h4>\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\">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[5]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-gnuplot\">\n<pre><span class=\"c\"># \u78c1\u5834 B<\/span>\r\n<span class=\"nf\">Bx<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"o\">-<\/span> <span class=\"n\">y<\/span><span class=\"o\">\/<\/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\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">By<\/span>(x, y) <span class=\"o\">=<\/span>   <span class=\"n\">x<\/span><span class=\"o\">\/<\/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\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"c\"># \u5e73\u884c2\u76f4\u7dda\u96fb\u6d41\u304c\u3064\u304f\u308b\u78c1\u5834<\/span>\r\n<span class=\"c\"># \u96fb\u7dda\u306e\u4f4d\u7f6e<\/span>\r\n<span class=\"nv\">x1<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">-1<\/span>\r\n<span class=\"nv\">y1<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"nv\">x2<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">1<\/span>\r\n<span class=\"nv\">y2<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n\r\n<span class=\"nf\">B2x<\/span>(x, y) <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\">x1<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y1<\/span><span class=\"p\">)<\/span> <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\">x2<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">B2y<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">By<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">-<\/span><span class=\"n\">x1<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y1<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"nf\">By<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">-<\/span><span class=\"n\">x2<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y2<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"c\"># \u78c1\u5834 B \u306e\u5927\u304d\u3055<\/span>\r\n<span class=\"nf\">B<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"nf\">B2x<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"nf\">B2y<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"c\"># \u898f\u683c\u5316\u3055\u308c\u305f hat B<\/span>\r\n<span class=\"nf\">F1<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">B2x<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/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\">y<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">F2<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">B2y<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/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\">y<\/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.02<\/span>\r\n<span class=\"nv\">N<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">300<\/span>\r\n<span class=\"err\">array<\/span> <span class=\"err\">X[N]<\/span>\r\n<span class=\"err\">array<\/span> <span class=\"err\">Y[N]<\/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\">para<\/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\">para<\/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=\"err\">do<\/span> <span class=\"err\">for<\/span> <span class=\"err\">[<\/span><span class=\"nv\">k<\/span><span class=\"o\">=<\/span><span class=\"mi\">1<\/span><span class=\"o\">:<\/span><span class=\"mi\">3<\/span><span class=\"o\">:<\/span><span class=\"mi\">2<\/span><span class=\"p\">]<\/span> <span class=\"p\">{<\/span>\r\n  <span class=\"nv\">z<\/span> <span class=\"o\">=<\/span> <span class=\"n\">k<\/span>\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\">2<\/span><span class=\"o\">:<\/span><span class=\"mi\">10<\/span><span class=\"o\">:<\/span><span class=\"mi\">2<\/span><span class=\"p\">]{<\/span>\r\n  <span class=\"nv\">x0<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">4.<\/span><span class=\"o\">\/<\/span><span class=\"n\">j<\/span> <span class=\"o\">+<\/span> <span class=\"mi\">1<\/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\">N<\/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\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\">\"para-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\">if<\/span> <span class=\"p\">(<\/span><span class=\"nf\">abs<\/span><span class=\"p\">(<\/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\">x2<\/span><span class=\"p\">)<\/span> <span class=\"err\">&lt;<\/span> <span class=\"err\">0.07)<\/span> <span class=\"err\">{<\/span><span class=\"k\">print<\/span> <span class=\"s\">\"       \"<\/span><span class=\"p\">}<\/span>\r\n    <span class=\"k\">else<\/span> <span class=\"err\">{<\/span>\r\n      <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/span><span class=\"p\">)<\/span>\r\n      <span class=\"err\">}<\/span>\r\n  <span class=\"err\">}<\/span>\r\n  <span class=\"k\">print<\/span> <span class=\"s\">\"         \"<\/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=\"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\">if<\/span> <span class=\"p\">(<\/span><span class=\"nf\">abs<\/span><span class=\"p\">(<\/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\">x2<\/span><span class=\"p\">)<\/span> <span class=\"err\">&lt;<\/span> <span class=\"err\">0.07)<\/span> <span class=\"err\">{<\/span><span class=\"k\">print<\/span> <span class=\"s\">\"       \"<\/span><span class=\"p\">}<\/span>\r\n    <span class=\"k\">else<\/span> <span class=\"err\">{<\/span>\r\n    <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/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\">print<\/span> <span class=\"s\">\"         \"<\/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=\"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 %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/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\">print<\/span> <span class=\"s\">\"         \"<\/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 %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/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\">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\">if<\/span> <span class=\"p\">((<\/span><span class=\"n\">j<\/span><span class=\"o\">==<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"err\">||(<\/span><span class=\"nv\">j<\/span><span class=\"o\">==<\/span><span class=\"mi\">6<\/span><span class=\"p\">)<\/span><span class=\"o\">||<\/span><span class=\"p\">(<\/span><span class=\"n\">j<\/span><span class=\"o\">==<\/span><span class=\"mi\">10<\/span><span class=\"p\">))<\/span> <span class=\"p\">{<\/span>\r\n  <span class=\"k\">set<\/span> <span class=\"nb\">print<\/span> <span class=\"s\">\"para.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\">12<\/span><span class=\"o\">:<\/span><span class=\"mi\">12<\/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 %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/span><span class=\"p\">)<\/span>\r\n  <span class=\"err\">}<\/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=\"n\">iend<\/span><span class=\"mi\">-20<\/span><span class=\"o\">:<\/span><span class=\"n\">iend<\/span><span class=\"mi\">-20<\/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 %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/span><span class=\"p\">)<\/span>\r\n  <span class=\"err\">}<\/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=\"mi\">30<\/span><span class=\"o\">:<\/span><span class=\"mi\">30<\/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 %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/span><span class=\"p\">)<\/span>\r\n  <span class=\"err\">}<\/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=\"mi\">-20<\/span><span class=\"o\">:<\/span><span class=\"n\">iend<\/span><span class=\"mi\">-20<\/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 %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><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  <span class=\"err\">}<\/span>\r\n<span class=\"err\">}<\/span>\r\n<span class=\"err\">}<\/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\">size<\/span> <span class=\"n\">ratio<\/span> <span class=\"mi\">1<\/span>\r\n\r\n<span class=\"k\">set<\/span> <span class=\"nb\">title<\/span> <span class=\"s\">\"\u540c\u3058\u5411\u304d\u306b\u6d41\u308c\u308b\u5e73\u884c2\u76f4\u7dda\u96fb\u6d41\u306b\u3088\u308b\u78c1\u5834\"<\/span>\r\n<span class=\"k\">unset<\/span> <span class=\"nb\">colorbox<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">palette<\/span> <span class=\"n\">cubehelix<\/span> <span class=\"n\">start<\/span> <span class=\"mi\">2<\/span> <span class=\"n\">cycles<\/span> <span class=\"mi\">0<\/span> <span class=\"n\">saturation<\/span> <span class=\"mi\">3<\/span> <span class=\"n\">negative<\/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\">2<\/span>\r\n<span class=\"c\"># \u77e2\u5370\u306e\u9577\u3055\u8abf\u6574\u7528<\/span>\r\n<span class=\"nv\">scaling<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.1<\/span>\r\n<span class=\"c\"># \u96fb\u7dda<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">arrow<\/span> <span class=\"mi\">1<\/span> <span class=\"n\">from<\/span> <span class=\"n\">x1<\/span><span class=\"o\">,<\/span><span class=\"n\">y1<\/span><span class=\"o\">,<\/span><span class=\"mi\">0<\/span> <span class=\"n\">to<\/span> <span class=\"n\">x1<\/span><span class=\"o\">,<\/span><span class=\"n\">y1<\/span><span class=\"o\">,<\/span><span class=\"mf\">4.3<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">4<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">lt<\/span> <span class=\"mi\">3<\/span> <span class=\"n\">filled<\/span> <span class=\"n\">head<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">arrow<\/span> <span class=\"mi\">2<\/span> <span class=\"n\">from<\/span> <span class=\"n\">x2<\/span><span class=\"o\">,<\/span><span class=\"n\">y2<\/span><span class=\"o\">,<\/span><span class=\"mi\">0<\/span> <span class=\"n\">to<\/span> <span class=\"n\">x2<\/span><span class=\"o\">,<\/span><span class=\"n\">y2<\/span><span class=\"o\">,<\/span><span class=\"mf\">4.3<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">4<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">lt<\/span> <span class=\"mi\">3<\/span> <span class=\"n\">filled<\/span> <span class=\"n\">head<\/span>\r\n\r\n<span class=\"k\">set<\/span> <span class=\"nb\">view<\/span> <span class=\"mi\">45<\/span><span class=\"o\">,<\/span><span class=\"mi\">0<\/span><span class=\"o\">,,<\/span>\r\n\r\n<span class=\"k\">splot<\/span> <span class=\"s\">\"para-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=\"mi\">3<\/span><span class=\"o\">:<\/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\">\"para.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=\"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=\"mi\">0<\/span><span class=\"p\">)<\/span><span class=\"o\">:<\/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\">lw<\/span> <span class=\"mi\">6<\/span> <span class=\"n\">lc<\/span> <span class=\"n\">palette<\/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-3848\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/vec-para2.svg\" alt=\"\" width=\"640\" 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[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\">\".\/vec-para2.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<h4 id=\"\u53cd\u5bfe\u5411\u304d\u306b\u96fb\u6d41-$I$-\u304c\u6d41\u308c\u3066\u3044\u308b\u5834\u5408\uff1aplot-\u7248\">\u53cd\u5bfe\u5411\u304d\u306b\u96fb\u6d41 $I$ \u304c\u6d41\u308c\u3066\u3044\u308b\u5834\u5408\uff1aplot \u7248<\/h4>\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-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[8]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-gnuplot\">\n<pre><span class=\"c\"># \u78c1\u5834 B<\/span>\r\n<span class=\"nf\">Bx<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"o\">-<\/span> <span class=\"n\">y<\/span><span class=\"o\">\/<\/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\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">By<\/span>(x, y) <span class=\"o\">=<\/span>   <span class=\"n\">x<\/span><span class=\"o\">\/<\/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\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"c\"># \u5e73\u884c2\u76f4\u7dda\u96fb\u6d41\u304c\u3064\u304f\u308b\u78c1\u5834<\/span>\r\n<span class=\"c\"># \u96fb\u7dda\u306e\u4f4d\u7f6e<\/span>\r\n<span class=\"nv\">x1<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">-1<\/span>\r\n<span class=\"nv\">y1<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"nv\">x2<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">1<\/span>\r\n<span class=\"nv\">y2<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n\r\n<span class=\"nf\">B2x<\/span>(x, y) <span class=\"o\">=<\/span> <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\">x1<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y1<\/span><span class=\"p\">)<\/span> <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\">x2<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">B2y<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"o\">-<\/span><span class=\"nf\">By<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">-<\/span><span class=\"n\">x1<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y1<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"nf\">By<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">-<\/span><span class=\"n\">x2<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y2<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"c\"># \u78c1\u5834 B \u306e\u5927\u304d\u3055<\/span>\r\n<span class=\"nf\">B<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"nf\">B2x<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"nf\">B2y<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"c\"># \u898f\u683c\u5316\u3055\u308c\u305f hat B<\/span>\r\n<span class=\"nf\">F1<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">B2x<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/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\">y<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">F2<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">B2y<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/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\">y<\/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.02<\/span>\r\n<span class=\"nv\">N<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">300<\/span>\r\n<span class=\"err\">array<\/span> <span class=\"err\">X[N]<\/span>\r\n<span class=\"err\">array<\/span> <span class=\"err\">Y[N]<\/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\">para<\/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\">para<\/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\">2<\/span><span class=\"o\">:<\/span><span class=\"mi\">8<\/span><span class=\"p\">]{<\/span>\r\n  <span class=\"nv\">x0<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">4.<\/span><span class=\"o\">\/<\/span><span class=\"n\">j<\/span> <span class=\"o\">+<\/span> <span class=\"mi\">1<\/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\">N<\/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\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\">\"para-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=\"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=\"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=\"k\">if<\/span> <span class=\"p\">(<\/span><span class=\"n\">j<\/span><span class=\"o\">%<\/span><span class=\"mi\">2<\/span>  <span class=\"o\">==<\/span><span class=\"mi\">1<\/span><span class=\"p\">)<\/span> <span class=\"err\">{<\/span>\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\">\"para.txt\"<\/span> <span class=\"n\">append<\/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=\"mi\">25<\/span><span class=\"o\">:<\/span><span class=\"mi\">25<\/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=\"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=\"mi\">25<\/span><span class=\"o\">:<\/span><span class=\"mi\">25<\/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\">set<\/span> <span class=\"nb\">print<\/span>\r\n  <span class=\"err\">}<\/span>\r\n<span class=\"err\">}<\/span>\r\n\r\n<span class=\"c\"># \u5185\u5074\u306e\u78c1\u529b\u7dda\u3092\u8ffd\u52a0\u3067\u6570\u5024\u8a08\u7b97<\/span>\r\n<span class=\"k\">h<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.2<\/span>\r\n<span class=\"nv\">N<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">30<\/span>\r\n<span class=\"err\">array<\/span> <span class=\"err\">X[N]<\/span>\r\n<span class=\"err\">array<\/span> <span class=\"err\">Y[N]<\/span>\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\">0<\/span><span class=\"o\">:<\/span><span class=\"mi\">2<\/span><span class=\"p\">]{<\/span>\r\n  <span class=\"nv\">x0<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.1<\/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\">N<\/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=\"nf\">abs<\/span><span class=\"p\">(<\/span><span class=\"n\">Y1<\/span><span class=\"p\">)<\/span> <span class=\"err\">&lt;<\/span> <span class=\"err\">ymin)<\/span> <span class=\"err\">{break}<\/span>\r\n    <span class=\"k\">if<\/span> <span class=\"p\">(<\/span><span class=\"nf\">abs<\/span><span class=\"p\">(<\/span><span class=\"n\">Y1<\/span><span class=\"p\">)<\/span> <span class=\"err\">&gt;<\/span> <span class=\"err\">ymax)<\/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\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\">\"para-line.txt\"<\/span> <span class=\"n\">append<\/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=\"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=\"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=\"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\"># \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=\"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=\"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=\"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=\"k\">if<\/span> <span class=\"p\">(<\/span><span class=\"n\">j<\/span><span class=\"o\">%<\/span><span class=\"mi\">2<\/span>  <span class=\"o\">==<\/span><span class=\"mi\">0<\/span><span class=\"p\">)<\/span> <span class=\"err\">{<\/span>\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\">\"para.txt\"<\/span> <span class=\"n\">append<\/span>\r\n  <span class=\"c\"># 1\/4<\/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\">3<\/span><span class=\"o\">:<\/span><span class=\"mi\">3<\/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=\"c\"># 3\/4<\/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\">3<\/span><span class=\"o\">:<\/span><span class=\"mi\">3<\/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\">set<\/span> <span class=\"nb\">print<\/span>\r\n  <span class=\"err\">}<\/span>\r\n<span class=\"err\">}<\/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\">\"\u53cd\u5bfe\u5411\u304d\u306b\u6d41\u308c\u308b\u5e73\u884c2\u76f4\u7dda\u96fb\u6d41\u306b\u3088\u308b\u78c1\u5834\"<\/span>\r\n<span class=\"k\">unset<\/span> <span class=\"nb\">colorbox<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">palette<\/span> <span class=\"n\">cubehelix<\/span> <span class=\"n\">start<\/span> <span class=\"mi\">2<\/span> <span class=\"n\">cycles<\/span> <span class=\"mi\">0<\/span> <span class=\"n\">saturation<\/span> <span class=\"mi\">3<\/span> <span class=\"n\">negative<\/span>\r\n\r\n<span class=\"c\"># \u77e2\u5370\u306e\u9577\u3055\u8abf\u6574\u7528<\/span>\r\n<span class=\"nv\">scaling<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.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\">2<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">at<\/span> <span class=\"n\">x1<\/span><span class=\"o\">,<\/span> <span class=\"n\">y1<\/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\">x1<\/span><span class=\"o\">,<\/span><span class=\"n\">y1<\/span> <span class=\"s\">\"\u00d7\"<\/span> <span class=\"n\">tc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">font<\/span> <span class=\"s\">\",24\"<\/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\">2<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">at<\/span> <span class=\"n\">x2<\/span><span class=\"o\">,<\/span> <span class=\"n\">y2<\/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=\"n\">x2<\/span><span class=\"o\">,<\/span><span class=\"n\">y2<\/span> <span class=\"s\">\"\u30fb\"<\/span> <span class=\"n\">tc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">font<\/span> <span class=\"s\">\",24\"<\/span>\r\n\r\n<span class=\"k\">plot<\/span> <span class=\"s\">\"para-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\">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\">\"para.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\">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\">6<\/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-3849\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/vec-para3.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[9]:<\/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\">\".\/vec-para3.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<h4 id=\"\u53cd\u5bfe\u5411\u304d\u306b\u96fb\u6d41-$I$-\u304c\u6d41\u308c\u3066\u3044\u308b\u5834\u5408\uff1asplot-\u7248\">\u53cd\u5bfe\u5411\u304d\u306b\u96fb\u6d41 $I$ \u304c\u6d41\u308c\u3066\u3044\u308b\u5834\u5408\uff1asplot \u7248<\/h4>\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-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\">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[11]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-gnuplot\">\n<pre><span class=\"c\"># \u78c1\u5834 B<\/span>\r\n<span class=\"nf\">Bx<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"o\">-<\/span> <span class=\"n\">y<\/span><span class=\"o\">\/<\/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\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">By<\/span>(x, y) <span class=\"o\">=<\/span>   <span class=\"n\">x<\/span><span class=\"o\">\/<\/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\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"c\"># \u5e73\u884c2\u76f4\u7dda\u96fb\u6d41\u304c\u3064\u304f\u308b\u78c1\u5834<\/span>\r\n<span class=\"c\"># \u96fb\u7dda\u306e\u4f4d\u7f6e<\/span>\r\n<span class=\"nv\">x1<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">-1<\/span>\r\n<span class=\"nv\">y1<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"nv\">x2<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">1<\/span>\r\n<span class=\"nv\">y2<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n\r\n<span class=\"nf\">B2x<\/span>(x, y) <span class=\"o\">=<\/span> <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\">x1<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y1<\/span><span class=\"p\">)<\/span> <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\">x2<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">B2y<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"o\">-<\/span><span class=\"nf\">By<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">-<\/span><span class=\"n\">x1<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y1<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"nf\">By<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">-<\/span><span class=\"n\">x2<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y2<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"c\"># \u78c1\u5834 B \u306e\u5927\u304d\u3055<\/span>\r\n<span class=\"nf\">B<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"nf\">B2x<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"nf\">B2y<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"c\"># \u898f\u683c\u5316\u3055\u308c\u305f hat B<\/span>\r\n<span class=\"nf\">F1<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">B2x<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/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\">y<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">F2<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">B2y<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">,<\/span> <span class=\"n\">y<\/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\">y<\/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.02<\/span>\r\n<span class=\"nv\">N<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">300<\/span>\r\n<span class=\"err\">array<\/span> <span class=\"err\">X[N]<\/span>\r\n<span class=\"err\">array<\/span> <span class=\"err\">Y[N]<\/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\">para<\/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\">para<\/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=\"err\">do<\/span> <span class=\"err\">for<\/span> <span class=\"err\">[<\/span><span class=\"nv\">k<\/span><span class=\"o\">=<\/span><span class=\"mi\">1<\/span><span class=\"o\">:<\/span><span class=\"mi\">3<\/span><span class=\"o\">:<\/span><span class=\"mi\">2<\/span><span class=\"p\">]<\/span> <span class=\"p\">{<\/span>\r\n  <span class=\"nv\">z<\/span> <span class=\"o\">=<\/span> <span class=\"n\">k<\/span>\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\">2<\/span><span class=\"o\">:<\/span><span class=\"mi\">10<\/span><span class=\"o\">:<\/span><span class=\"mi\">2<\/span><span class=\"p\">]{<\/span>\r\n  <span class=\"nv\">x0<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">4.<\/span><span class=\"o\">\/<\/span><span class=\"n\">j<\/span> <span class=\"o\">+<\/span> <span class=\"mi\">1<\/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\">N<\/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\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\">\"para-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\">if<\/span> <span class=\"p\">(<\/span><span class=\"nf\">abs<\/span><span class=\"p\">(<\/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\">x2<\/span><span class=\"p\">)<\/span> <span class=\"err\">&lt;<\/span> <span class=\"err\">0.07)<\/span> <span class=\"err\">{<\/span><span class=\"k\">print<\/span> <span class=\"s\">\"       \"<\/span><span class=\"p\">}<\/span>\r\n    <span class=\"k\">else<\/span> <span class=\"err\">{<\/span>\r\n      <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/span><span class=\"p\">)}<\/span>\r\n  <span class=\"err\">}<\/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 %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/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\">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\">if<\/span> <span class=\"p\">(<\/span><span class=\"nf\">abs<\/span><span class=\"p\">(<\/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\">x2<\/span><span class=\"p\">)<\/span> <span class=\"err\">&lt;<\/span> <span class=\"err\">0.07)<\/span> <span class=\"err\">{<\/span><span class=\"k\">print<\/span> <span class=\"s\">\"       \"<\/span><span class=\"p\">}<\/span>\r\n    <span class=\"k\">else<\/span> <span class=\"err\">{<\/span>\r\n    <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/span><span class=\"p\">)}<\/span>\r\n  <span class=\"err\">}<\/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 %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/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\">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\">if<\/span> <span class=\"p\">((<\/span><span class=\"n\">j<\/span><span class=\"o\">==<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"err\">||(<\/span><span class=\"nv\">j<\/span><span class=\"o\">==<\/span><span class=\"mi\">6<\/span><span class=\"p\">))<\/span> <span class=\"p\">{<\/span>\r\n  <span class=\"k\">set<\/span> <span class=\"nb\">print<\/span> <span class=\"s\">\"para.txt\"<\/span> <span class=\"n\">append<\/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=\"mi\">-55<\/span><span class=\"o\">:<\/span><span class=\"n\">iend<\/span><span class=\"mi\">-55<\/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 %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/span><span class=\"p\">)<\/span>\r\n  <span class=\"err\">}<\/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=\"mi\">-55<\/span><span class=\"o\">:<\/span><span class=\"n\">iend<\/span><span class=\"mi\">-55<\/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 %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><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  <span class=\"err\">}<\/span>\r\n<span class=\"err\">}<\/span>\r\n<span class=\"err\">}<\/span>\r\n\r\n<span class=\"c\"># \u5185\u5074\u306e\u78c1\u529b\u7dda\u3092\u8ffd\u52a0\u3067\u6570\u5024\u8a08\u7b97<\/span>\r\n<span class=\"k\">h<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.02<\/span>\r\n<span class=\"nv\">N<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">300<\/span>\r\n<span class=\"err\">array<\/span> <span class=\"err\">X[N]<\/span>\r\n<span class=\"err\">array<\/span> <span class=\"err\">Y[N]<\/span>\r\n\r\n<span class=\"err\">do<\/span> <span class=\"err\">for<\/span> <span class=\"err\">[<\/span><span class=\"nv\">k<\/span><span class=\"o\">=<\/span><span class=\"mi\">1<\/span><span class=\"o\">:<\/span><span class=\"mi\">3<\/span><span class=\"o\">:<\/span><span class=\"mi\">2<\/span><span class=\"p\">]<\/span> <span class=\"p\">{<\/span>\r\n  <span class=\"nv\">z<\/span> <span class=\"o\">=<\/span> <span class=\"n\">k<\/span>\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\">0<\/span><span class=\"o\">:<\/span><span class=\"mi\">2<\/span><span class=\"p\">]{<\/span>\r\n  <span class=\"nv\">x0<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.1<\/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\">N<\/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=\"o\">-<\/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=\"o\">-<\/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\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\">\"para-line.txt\"<\/span> <span class=\"n\">append<\/span>\r\n  <span class=\"c\"># 1\/4<\/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 %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/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\">print<\/span> <span class=\"s\">\"        \"<\/span>\r\n  <span class=\"c\"># 2\/4<\/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\">if<\/span> <span class=\"p\">(<\/span><span class=\"nf\">abs<\/span><span class=\"p\">(<\/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\">x2<\/span><span class=\"p\">)<\/span> <span class=\"err\">&lt;<\/span> <span class=\"err\">0.1)<\/span> <span class=\"err\">{<\/span><span class=\"k\">print<\/span> <span class=\"s\">\"        \"<\/span><span class=\"p\">}<\/span>\r\n    <span class=\"k\">else<\/span> <span class=\"err\">{<\/span>\r\n    <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/span><span class=\"p\">)<\/span>\r\n    <span class=\"err\">}<\/span>\r\n  <span class=\"err\">}<\/span>\r\n  <span class=\"k\">print<\/span> <span class=\"s\">\"        \"<\/span>\r\n  <span class=\"k\">print<\/span> <span class=\"s\">\"        \"<\/span>\r\n  <span class=\"c\"># 3\/4<\/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 %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/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\">print<\/span> <span class=\"s\">\"        \"<\/span>\r\n  <span class=\"c\"># 4\/4<\/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\">if<\/span> <span class=\"p\">(<\/span><span class=\"nf\">abs<\/span><span class=\"p\">(<\/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\">x2<\/span><span class=\"p\">)<\/span> <span class=\"err\">&lt;<\/span> <span class=\"err\">0.1)<\/span> <span class=\"err\">{<\/span><span class=\"k\">print<\/span> <span class=\"s\">\"        \"<\/span><span class=\"p\">}<\/span>\r\n    <span class=\"k\">else<\/span> <span class=\"err\">{<\/span>\r\n   <span class=\"k\">print<\/span> <span class=\"nf\">sprintf<\/span><span class=\"p\">(<\/span><span class=\"s\">\"%8.4f %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/span><span class=\"p\">)<\/span>\r\n   <span class=\"err\">}<\/span>\r\n  <span class=\"err\">}<\/span>\r\n  <span class=\"k\">print<\/span> <span class=\"s\">\"        \"<\/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=\"k\">if<\/span> <span class=\"p\">((<\/span><span class=\"n\">j<\/span><span class=\"o\">%<\/span><span class=\"mi\">2<\/span>  <span class=\"o\">==<\/span><span class=\"mi\">0<\/span><span class=\"p\">)<\/span><span class=\"err\">)<\/span> <span class=\"err\">{<\/span>\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\">\"para.txt\"<\/span> <span class=\"n\">append<\/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=\"mi\">16<\/span><span class=\"o\">:<\/span><span class=\"mi\">16<\/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 %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><span class=\"o\">,<\/span> <span class=\"n\">z<\/span><span class=\"p\">)<\/span>\r\n  <span class=\"err\">}<\/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=\"mi\">16<\/span><span class=\"o\">:<\/span><span class=\"mi\">16<\/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 %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><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  <span class=\"err\">}<\/span>\r\n<span class=\"err\">}<\/span>\r\n<span class=\"err\">}<\/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\">size<\/span> <span class=\"n\">ratio<\/span> <span class=\"mi\">1<\/span>\r\n\r\n<span class=\"k\">set<\/span> <span class=\"nb\">title<\/span> <span class=\"s\">\"\u53cd\u5bfe\u5411\u304d\u306b\u6d41\u308c\u308b\u5e73\u884c2\u76f4\u7dda\u96fb\u6d41\u306b\u3088\u308b\u78c1\u5834\"<\/span>\r\n<span class=\"k\">unset<\/span> <span class=\"nb\">colorbox<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">palette<\/span> <span class=\"n\">cubehelix<\/span> <span class=\"n\">start<\/span> <span class=\"mi\">2<\/span> <span class=\"n\">cycles<\/span> <span class=\"mi\">0<\/span> <span class=\"n\">saturation<\/span> <span class=\"mi\">3<\/span> <span class=\"n\">negative<\/span>\r\n\r\n<span class=\"c\"># \u77e2\u5370\u306e\u9577\u3055\u8abf\u6574\u7528<\/span>\r\n<span class=\"nv\">scaling<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.2<\/span>\r\n<span class=\"c\"># \u96fb\u7dda<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">arrow<\/span> <span class=\"mi\">1<\/span> <span class=\"n\">from<\/span> <span class=\"n\">x1<\/span><span class=\"o\">,<\/span><span class=\"n\">y1<\/span><span class=\"o\">,<\/span><span class=\"mi\">4<\/span> <span class=\"n\">to<\/span> <span class=\"n\">x1<\/span><span class=\"o\">,<\/span><span class=\"n\">y1<\/span><span class=\"o\">,<\/span><span class=\"mi\">0<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">4<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">lt<\/span> <span class=\"mi\">3<\/span> <span class=\"n\">filled<\/span> <span class=\"n\">head<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">arrow<\/span> <span class=\"mi\">2<\/span> <span class=\"n\">from<\/span> <span class=\"n\">x2<\/span><span class=\"o\">,<\/span><span class=\"n\">y2<\/span><span class=\"o\">,<\/span><span class=\"mi\">0<\/span> <span class=\"n\">to<\/span> <span class=\"n\">x2<\/span><span class=\"o\">,<\/span><span class=\"n\">y2<\/span><span class=\"o\">,<\/span><span class=\"mi\">4<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">4<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">lt<\/span> <span class=\"mi\">3<\/span> <span class=\"n\">filled<\/span> <span class=\"n\">head<\/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\">2<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">view<\/span> <span class=\"mi\">45<\/span><span class=\"o\">,<\/span><span class=\"mi\">0<\/span><span class=\"o\">,,<\/span>\r\n\r\n<span class=\"k\">splot<\/span> <span class=\"s\">\"para-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=\"mi\">3<\/span><span class=\"o\">:<\/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\">palette<\/span> <span class=\"nb\">notitle<\/span><span class=\"o\">,<\/span>\\\r\n  <span class=\"s\">\"para.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=\"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=\"mi\">0<\/span><span class=\"p\">)<\/span><span class=\"o\">:<\/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\">lw<\/span> <span class=\"mi\">6<\/span> <span class=\"n\">lc<\/span> <span class=\"n\">palette<\/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-3850\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/vec-para4.svg\" alt=\"\" width=\"640\" 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[12]:<\/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\">\".\/vec-para4.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>\n<\/div>\n<\/div>\n","protected":false},"excerpt":{"rendered":"<p>gnuplot \u3060\u3051\u3067\uff0c\u78c1\u529b\u7dda\u3092\u6570\u5024\u7684\u306b\u89e3\u304d\uff0c\u78c1\u529b\u7dda\u306b\u6cbf\u3063\u305f\u78c1\u5834\u30d9\u30af\u30c8\u30eb\u3092\u63cf\u304f\u3002\u78c1\u5834\u30d9\u30af\u30c8\u30eb\u304c\u78c1\u529b\u7dda\u306e\u63a5\u30d9\u30af\u30c8\u30eb\u3067\u3042\u308b\u3053\u3068\u304b\u3089\uff0c\u5fae\u5206\u65b9\u7a0b\u5f0f\u3092\u30eb\u30f3\u30b2\u30fb\u30af\u30c3\u30bf\u6cd5\u3067\u89e3\u304d\uff0c\u6570\u5024\u89e3\u3092\u30d5\u30a1\u30a4\u30eb\u306b\u4fdd\u5b58\u3057\uff0c\u305d\u308c\u3092\u8aad\u307f\u8fbc\u3093\u3067\u78c1\u5834\u30d9\u30af\u30c8\u30eb\u306e\u59cb\u70b9\u306e\u5ea7\u6a19\u3068\u3057\u3066\u30d9\u30af\u30c8\u30eb\u5834\u3092\u63cf\u304f\u3002\u3053\u308c\u3092\u5168\u3066 gnuplot \u3060\u3051\u3067\u3084\u308b\u3002<\/p><p><a class=\"more-link btn\" href=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/3702\/\">\u7d9a\u304d\u3092\u8aad\u3080<\/a><\/p>\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":[15,19],"tags":[],"class_list":["post-3702","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\/3702","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=3702"}],"version-history":[{"count":26,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/posts\/3702\/revisions"}],"predecessor-version":[{"id":3912,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/posts\/3702\/revisions\/3912"}],"wp:attachment":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/media?parent=3702"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/categories?post=3702"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/tags?post=3702"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}