{"id":3643,"date":"2022-09-26T16:41:22","date_gmt":"2022-09-26T07:41:22","guid":{"rendered":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/?p=3643"},"modified":"2023-03-28T17:57:37","modified_gmt":"2023-03-28T08:57:37","slug":"gnuplot-%e3%81%a7%e5%a7%8b%e7%82%b9%e3%81%ae%e5%ba%a7%e6%a8%99%e3%82%92%e3%83%95%e3%82%a1%e3%82%a4%e3%83%ab%e3%81%8b%e3%82%89%e8%aa%ad%e3%81%bf%e8%be%bc%e3%82%93%e3%81%a7%e3%83%99%e3%82%af%e3%83%88","status":"publish","type":"post","link":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/3643\/","title":{"rendered":"gnuplot \u3067\u5ea7\u6a19\u3092\u30d5\u30a1\u30a4\u30eb\u304b\u3089\u8aad\u307f\u8fbc\u3093\u3067\u96fb\u6c17\u529b\u7dda\u30fb\u78c1\u529b\u7dda\u3092\u63cf\u304f"},"content":{"rendered":"

\u300cgnuplot \u3067\u30d9\u30af\u30c8\u30eb\u5834\uff08\u96fb\u5834\u30fb\u78c1\u5834\uff09\u3092\u63cf\u304f<\/a>\u300d\u3084\u300cpalette \u3067\u30d9\u30af\u30c8\u30eb\u306e\u8272\u3092\u5909\u3048\u3066 gnuplot \u3067\u96fb\u5834\u30fb\u78c1\u5834\u3092\u63cf\u304f<\/a>\u300d\u3067\u63cf\u3044\u3066\u307f\u305f\u304c\uff0c\u3082\u3046\u5c11\u3057\u3002<\/p>\n

\n
\n
\n

gnuplot \u3060\u3051\u3067\u30d9\u30af\u30c8\u30eb\u306e\u59cb\u70b9\u306e\u5ea7\u6a19\u30d5\u30a1\u30a4\u30eb\u3092\u4f5c\u6210\u3057\uff0c\u305d\u308c\u3092\u8aad\u307f\u8fbc\u3093\u3067\u96fb\u6c17\u529b\u7dda\u30fb\u78c1\u529b\u7dda\u3092\u63cf\u304f\u4f8b\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n

<\/p>\n

\n
<\/div>\n
\n
\n

\u70b9\u96fb\u8377\u306b\u3088\u308b\u96fb\u5834 $\\boldsymbol{E}$<\/h3>\n

\u6b63\u306e\u70b9\u96fb\u8377\u306b\u3088\u308b\u96fb\u5834<\/h4>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[20]:<\/div>\n
\n
\n
\n
%gnu<\/span>plot<\/span> inline<\/span> svg<\/span> size<\/span> 480<\/span>,<\/span>480<\/span> fixed<\/span> enhanced<\/span> font<\/span> 'Noto Sans CJK JP,16'<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[27]:<\/div>\n
\n
\n
\n
reset<\/span>\r\nunset<\/span> xtics<\/span>\r\nunset<\/span> ytics<\/span>\r\nunset<\/span> border<\/span>\r\n\r\n# \u6b63\u96fb\u8377\u306e\u4f4d\u7f6e<\/span>\r\nx1<\/span> =<\/span> 0<\/span>\r\ny1<\/span> =<\/span> 1<\/span>\r\n\r\n# \u96fb\u6c17\u529b\u7dda\u7528\u5ea7\u6a19\u30d5\u30a1\u30a4\u30eb<\/span>\r\nset<\/span> print<\/span> \"xy-line.txt\"<\/span>\r\ndo<\/span> for<\/span> [<\/span>j<\/span>=<\/span>0<\/span>:<\/span> 15<\/span>]<\/span> {<\/span>\r\n  do<\/span> for<\/span> [<\/span>i<\/span>=<\/span>1<\/span>:<\/span>60<\/span>]<\/span> {<\/span>\r\n    r<\/span> =<\/span> 0.2<\/span> *<\/span> i<\/span>\r\n    theta<\/span> =<\/span> pi<\/span>\/<\/span>8<\/span> *<\/span> j<\/span>\r\n    x<\/span> =<\/span> x1<\/span> +<\/span> r<\/span> *<\/span> cos<\/span>(<\/span>theta<\/span>)<\/span>\r\n    y<\/span> =<\/span> y1<\/span> +<\/span> r<\/span> *<\/span> sin<\/span>(<\/span>theta<\/span>)<\/span>\r\n    print<\/span> sprintf<\/span>(<\/span>\"%8.4f  %8.4f\"<\/span>,<\/span> x<\/span>,<\/span> y<\/span>)<\/span>\r\n  }<\/span>\r\n  # with lines \u3059\u308b\u306e\u3067\u9069\u5b9c\u7a7a\u884c\u3092\u3044\u308c\u308b<\/span>\r\n  print<\/span> \"\"<\/span>\r\n}<\/span>\r\nset<\/span> print<\/span>\r\n\r\n# \u30d9\u30af\u30c8\u30eb\u306e\u59cb\u70b9\u7528\u5ea7\u6a19\u30d5\u30a1\u30a4\u30eb<\/span>\r\n# \u77e2\u5370\u306e\u6df7\u96d1\u9632\u6b62\u306e\u305f\u3081\u9069\u5b9c\u9593\u5f15\u304f<\/span>\r\nset<\/span> print<\/span> \"xy.txt\"<\/span>\r\ndo<\/span> for<\/span> [<\/span>j<\/span>=<\/span>0<\/span>:<\/span> 15<\/span>]<\/span> {<\/span>\r\n  do<\/span> for<\/span> [<\/span>i<\/span>=<\/span>2<\/span>:<\/span>8<\/span>:<\/span>3<\/span>]<\/span> {<\/span>\r\n    r<\/span> =<\/span> 0.5<\/span> *<\/span> i<\/span>\r\n    theta<\/span> =<\/span> pi<\/span>\/<\/span>8<\/span> *<\/span> j<\/span>\r\n    x<\/span> =<\/span> x1<\/span> +<\/span> r<\/span> *<\/span> cos<\/span>(<\/span>theta<\/span>)<\/span>\r\n    y<\/span> =<\/span> y1<\/span> +<\/span> r<\/span> *<\/span> sin<\/span>(<\/span>theta<\/span>)<\/span>\r\n    print<\/span> sprintf<\/span>(<\/span>\"%8.4f  %8.4f\"<\/span>,<\/span> x<\/span>,<\/span> y<\/span>)<\/span>\r\n  }<\/span>\r\n}<\/span>\r\nset<\/span> print<\/span>\r\n\r\n# \u96fb\u5834 E<\/span>\r\nEx<\/span>(x, y) =<\/span> (<\/span>x<\/span>-<\/span>x1<\/span>)<\/span>\/<\/span>(<\/span>sqrt<\/span>((<\/span>x<\/span>-<\/span>x1<\/span>)<\/span>**<\/span>2<\/span> +<\/span> (<\/span>y<\/span>-<\/span>y1<\/span>)<\/span>**<\/span>2<\/span>)<\/span>**<\/span>3<\/span>)<\/span>\r\nEy<\/span>(x, y) =<\/span> (<\/span>y<\/span>-<\/span>y1<\/span>)<\/span>\/<\/span>(<\/span>sqrt<\/span>((<\/span>x<\/span>-<\/span>x1<\/span>)<\/span>**<\/span>2<\/span> +<\/span> (<\/span>y<\/span>-<\/span>y1<\/span>)<\/span>**<\/span>2<\/span>)<\/span>**<\/span>3<\/span>)<\/span>\r\n# \u96fb\u5834 E \u306e\u5927\u304d\u3055<\/span>\r\nE<\/span>(x, y) =<\/span> sqrt<\/span>(<\/span>Ex<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>**<\/span>2<\/span> +<\/span> Ey<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>**<\/span>2<\/span>)<\/span>\r\n# \u898f\u683c\u5316\u3055\u308c\u305f hat E<\/span>\r\nhEx<\/span>(x, y) =<\/span> Ex<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\/<\/span>E<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\r\nhEy<\/span>(x, y) =<\/span> Ey<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\/<\/span>E<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\r\n\r\nset<\/span> xrange<\/span> [<\/span>-3<\/span>.5<\/span>:<\/span>3.5<\/span>]<\/span>\r\nset<\/span> yrange<\/span> [<\/span>-3<\/span>.5<\/span>:<\/span>3.5<\/span>]<\/span>\r\nset<\/span> size<\/span> ratio<\/span> 1<\/span>\r\n\r\n# \u77e2\u5370\u306e\u9577\u3055\u8abf\u6574\u7528<\/span>\r\nscaling<\/span> =<\/span> 0.1<\/span>\r\n\r\n# \u6b63\u96fb\u8377\u306e\u4f4d\u7f6e<\/span>\r\nset<\/span> label<\/span> 1<\/span> point<\/span> pt<\/span> 6<\/span> ps<\/span> 1.6<\/span> lc<\/span> \"black\"<\/span> at<\/span> x1<\/span>,<\/span> y1<\/span>\r\nset<\/span> label<\/span> 2<\/span> center<\/span> at<\/span> first<\/span> x1<\/span>,<\/span> y1<\/span> \"\uff0b\"<\/span> tc<\/span> \"black\"<\/span> font<\/span> \",22\"<\/span>\r\n\r\nset<\/span> title<\/span> \"\u6b63\u96fb\u8377\u306b\u3088\u308b\u96fb\u5834\"<\/span>\r\nunset<\/span> colorbox<\/span>\r\nset<\/span> palette<\/span> cubehelix<\/span> start<\/span> 0<\/span> cycles<\/span> 0<\/span> saturation<\/span> 3<\/span> negative<\/span>\r\n\r\nplot<\/span> \\\r\n  \"xy-line.txt\"<\/span> \\\r\n    u<\/span> 1<\/span>:<\/span>2<\/span>:<\/span>(<\/span>log<\/span>(<\/span>E<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>)))<\/span> w<\/span> l<\/span> lw<\/span> 3<\/span> lc<\/span> palette<\/span> notitle<\/span>,<\/span>\\\r\n  \"xy.txt\"<\/span> \\\r\n    u<\/span> 1<\/span>:<\/span>2<\/span>:<\/span>(<\/span>scaling<\/span>*<\/span>hEx<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>))<\/span>:<\/span>(<\/span>scaling<\/span>*<\/span>hEy<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>))<\/span>:<\/span>(<\/span>log<\/span>(<\/span>E<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>)))<\/span> \\\r\n    w<\/span> vec<\/span> notitle<\/span> lw<\/span> 8<\/span> lc<\/span> palette<\/span> filled<\/span> head<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
\n
<\/div>\n
\n
\"\"<\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[3]:<\/div>\n
\n
\n
\n
set<\/span> output<\/span> \".\/vec-fig01d.svg\"<\/span>\r\nreplot<\/span>\r\nset<\/span> output<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
<\/div>\n
\n
\n
\u8ca0\u306e\u70b9\u96fb\u8377\u306b\u3088\u308b\u96fb\u5834<\/h4>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[28]:<\/div>\n
\n
\n
\n
reset<\/span>\r\nunset<\/span> xtics<\/span>\r\nunset<\/span> ytics<\/span>\r\nunset<\/span> border<\/span>\r\nset<\/span> xrange<\/span> [<\/span>-3<\/span>.5<\/span>:<\/span>3.5<\/span>]<\/span>\r\nset<\/span> yrange<\/span> [<\/span>-3<\/span>.5<\/span>:<\/span>3.5<\/span>]<\/span>\r\nset<\/span> size<\/span> ratio<\/span> 1<\/span>\r\n\r\n# \u6975\u5ea7\u6a19\u3067 x, y \u306e\u30c7\u30fc\u30bf\u30d5\u30a1\u30a4\u30eb\u4f5c\u6210<\/span>\r\n# \u8ca0\u96fb\u8377\u306e\u4f4d\u7f6e<\/span>\r\nx2<\/span> =<\/span> 0<\/span>\r\ny2<\/span> =<\/span> -1<\/span>\r\n\r\n# \u96fb\u6c17\u529b\u7dda\u7528\u5ea7\u6a19\u30d5\u30a1\u30a4\u30eb<\/span>\r\nset<\/span> print<\/span> \"xy-line.txt\"<\/span>\r\ndo<\/span> for<\/span> [<\/span>j<\/span>=<\/span>0<\/span>:<\/span>15<\/span>]<\/span> {<\/span>\r\n  do<\/span> for<\/span> [<\/span>i<\/span>=<\/span>1<\/span>:<\/span>60<\/span>]<\/span> {<\/span>\r\n    r<\/span> =<\/span> 0.2<\/span> *<\/span> i<\/span>\r\n    theta<\/span> =<\/span> pi<\/span>\/<\/span>8<\/span> *<\/span> j<\/span>\r\n    x<\/span> =<\/span> x2<\/span> +<\/span> r<\/span> *<\/span> cos<\/span>(<\/span>theta<\/span>)<\/span>\r\n    y<\/span> =<\/span> y2<\/span> +<\/span> r<\/span> *<\/span> sin<\/span>(<\/span>theta<\/span>)<\/span>\r\n    print<\/span> sprintf<\/span>(<\/span>\"%8.4f  %8.4f\"<\/span>,<\/span> x<\/span>,<\/span> y<\/span>)<\/span>\r\n  }<\/span>\r\n  # with lines \u3059\u308b\u306e\u3067\u9069\u5b9c\u7a7a\u884c\u3092\u3044\u308c\u308b<\/span>\r\n  print<\/span> \"\"<\/span>\r\n}<\/span>\r\nset<\/span> print<\/span>\r\n\r\n# \u30d9\u30af\u30c8\u30eb\u306e\u59cb\u70b9\u7528\u5ea7\u6a19\u30d5\u30a1\u30a4\u30eb<\/span>\r\n# \u77e2\u5370\u306e\u6df7\u96d1\u9632\u6b62\u306e\u305f\u3081\u9069\u5b9c\u9593\u5f15\u304f<\/span>\r\nset<\/span> print<\/span> \"xy.txt\"<\/span>\r\ndo<\/span> for<\/span> [<\/span>j<\/span>=<\/span>0<\/span>:<\/span>15<\/span>]<\/span> {<\/span>\r\n  do<\/span> for<\/span> [<\/span>i<\/span>=<\/span>3<\/span>:<\/span>9<\/span>:<\/span>3<\/span>]<\/span> {<\/span>\r\n    r<\/span> =<\/span> 0.5<\/span> *<\/span> i<\/span>\r\n    theta<\/span> =<\/span> pi<\/span>\/<\/span>8<\/span> *<\/span> j<\/span>\r\n    x<\/span> =<\/span> x2<\/span> +<\/span> r<\/span> *<\/span> cos<\/span>(<\/span>theta<\/span>)<\/span>\r\n    y<\/span> =<\/span> y2<\/span> +<\/span> r<\/span> *<\/span> sin<\/span>(<\/span>theta<\/span>)<\/span>\r\n    print<\/span> sprintf<\/span>(<\/span>\"%8.4f  %8.4f\"<\/span>,<\/span> x<\/span>,<\/span> y<\/span>)<\/span>\r\n  }<\/span>\r\n}<\/span>\r\nset<\/span> print<\/span>\r\n\r\n# \u96fb\u5834 E<\/span>\r\nEx<\/span>(x, y) =<\/span> -<\/span>(<\/span>x<\/span>-<\/span>x2<\/span>)<\/span>\/<\/span>(<\/span>sqrt<\/span>((<\/span>x<\/span>-<\/span>x2<\/span>)<\/span>**<\/span>2<\/span> +<\/span> (<\/span>y<\/span>-<\/span>y2<\/span>)<\/span>**<\/span>2<\/span>)<\/span>**<\/span>3<\/span>)<\/span>\r\nEy<\/span>(x, y) =<\/span> -<\/span>(<\/span>y<\/span>-<\/span>y2<\/span>)<\/span>\/<\/span>(<\/span>sqrt<\/span>((<\/span>x<\/span>-<\/span>x2<\/span>)<\/span>**<\/span>2<\/span> +<\/span> (<\/span>y<\/span>-<\/span>y2<\/span>)<\/span>**<\/span>2<\/span>)<\/span>**<\/span>3<\/span>)<\/span>\r\n# \u96fb\u5834 E \u306e\u5927\u304d\u3055<\/span>\r\nE<\/span>(x, y) =<\/span> sqrt<\/span>(<\/span>Ex<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>**<\/span>2<\/span> +<\/span> Ey<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>**<\/span>2<\/span>)<\/span>\r\n# \u898f\u683c\u5316\u3055\u308c\u305f hat E<\/span>\r\nhEx<\/span>(x, y) =<\/span> Ex<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\/<\/span>E<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\r\nhEy<\/span>(x, y) =<\/span> Ey<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\/<\/span>E<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\r\n\r\nset<\/span> xrange<\/span> [<\/span>-3<\/span>.5<\/span>:<\/span>3.5<\/span>]<\/span>\r\nset<\/span> yrange<\/span> [<\/span>-3<\/span>.5<\/span>:<\/span>3.5<\/span>]<\/span>\r\nset<\/span> size<\/span> ratio<\/span> 1<\/span>\r\n\r\n# \u77e2\u5370\u306e\u9577\u3055\u8abf\u6574\u7528<\/span>\r\nscaling<\/span> =<\/span> 0.1<\/span>\r\n\r\n# \u8ca0\u96fb\u8377\u306e\u4f4d\u7f6e<\/span>\r\nset<\/span> label<\/span> 3<\/span> point<\/span> pt<\/span> 6<\/span> ps<\/span> 1.6<\/span> lc<\/span> \"black\"<\/span> at<\/span> x2<\/span>,<\/span> y2<\/span>\r\nset<\/span> label<\/span> 4<\/span> center<\/span> at<\/span> first<\/span> x2<\/span>,<\/span> y2<\/span> \"\u30fc\"<\/span> tc<\/span> \"black\"<\/span> font<\/span> \",20\"<\/span>\r\n\r\nset<\/span> title<\/span> \"\u8ca0\u96fb\u8377\u306b\u3088\u308b\u96fb\u5834\"<\/span>\r\nunset<\/span> colorbox<\/span>\r\nset<\/span> palette<\/span> cubehelix<\/span> start<\/span> 0<\/span> cycles<\/span> 0<\/span> saturation<\/span> 3<\/span> negative<\/span>\r\n\r\nplot<\/span> \"xy-line.txt\"<\/span> u<\/span> 1<\/span>:<\/span>2<\/span>:<\/span>(<\/span>log<\/span>(<\/span>E<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>)))<\/span> w<\/span> l<\/span> lw<\/span> 3<\/span> lc<\/span> palette<\/span> notitle<\/span>,<\/span>\\\r\n  \"xy.txt\"<\/span> \\\r\n  u<\/span> 1<\/span>:<\/span>2<\/span>:<\/span>(<\/span>scaling<\/span>*<\/span>hEx<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>))<\/span>:<\/span>(<\/span>scaling<\/span>*<\/span>hEy<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>))<\/span>:<\/span>(<\/span>log<\/span>(<\/span>E<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>)))<\/span> w<\/span> vec<\/span> \\\r\n  notitle<\/span> lw<\/span> 8<\/span> lc<\/span> palette<\/span> filled<\/span> head<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
\n
<\/div>\n
\n
\"\"<\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[5]:<\/div>\n
\n
\n
\n
set<\/span> output<\/span> \".\/vec-fig02d.svg\"<\/span>\r\nreplot<\/span>\r\nset<\/span> output<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
<\/div>\n
\n
\n
\u4e00\u69d8\u306a\u7dda\u96fb\u8377\u306b\u3088\u308b\u96fb\u5834<\/h3>\n
$z$ \u8ef8\u4e0a\u306e\u4e00\u69d8\u306a\u7dda\u96fb\u8377\u5bc6\u5ea6 $\\lambda (\\mbox{C}\/\\mbox{m})$ \u306b\u3088\u308b\u96fb\u5834\u306f\uff0c$\\boldsymbol{r} \\equiv (x, y, 0), \\ r^2 \\equiv \\boldsymbol{r} \\cdot\\boldsymbol{r} $ \u3068\u3057\u3066\u307e\u3068\u3081\u308b\u3068<\/p>\n
$$\\boldsymbol{E} = \\frac{\\lambda}{2\\pi\\varepsilon_0} \\frac{\\boldsymbol{r}}{r^2}$$<\/p>\n
\u5b9a\u6570\u90e8\u5206\u3092\u9069\u5b9c\u898f\u683c\u5316\u3057\u3066<\/p>\n
\\begin{eqnarray}
\nE_x &\\Rightarrow& \\frac{x}{x^2 + y^2} \\\\
\nE_x &\\Rightarrow& \\frac{y}{x^2 + y^2} \\\\
\n\\end{eqnarray}<\/p>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[40]:<\/div>\n
\n
\n
\n
# \u96fb\u6c17\u529b\u7dda<\/span>\r\nset<\/span> print<\/span> \"xyz-line.txt\"<\/span>\r\ndo<\/span> for<\/span> [<\/span>z<\/span>=<\/span>0<\/span>:<\/span>2<\/span>]{<\/span>\r\n  do<\/span> for<\/span> [<\/span>i<\/span>=<\/span>0<\/span>:<\/span>11<\/span>]{<\/span>\r\n    theta<\/span> =<\/span> 2<\/span>*<\/span>pi<\/span>\/<\/span>12<\/span> *<\/span> i<\/span>\r\n    do<\/span> for<\/span> [<\/span>j<\/span>=<\/span>1<\/span>:<\/span>20<\/span>]<\/span> {<\/span>\r\n      r<\/span> =<\/span> 0.2<\/span>*<\/span>j<\/span>\r\n      x<\/span> =<\/span> r<\/span> *<\/span> cos<\/span>(<\/span>theta<\/span>)<\/span>\r\n      y<\/span> =<\/span> r<\/span> *<\/span> sin<\/span>(<\/span>theta<\/span>)<\/span>\r\n      print<\/span> sprintf<\/span>(<\/span>\" %8.4f %8.4f %8.4f\"<\/span>,<\/span> x<\/span>,<\/span> y<\/span>,<\/span> z<\/span>)<\/span>\r\n    }<\/span>\r\n    print<\/span> \"\"<\/span>\r\n    print<\/span> \"\"<\/span>\r\n  }<\/span>\r\n    print<\/span> \"\"<\/span>\r\n    print<\/span> \"\"<\/span>\r\n}<\/span>\r\nset<\/span> print<\/span>\r\n\r\n# \u96fb\u5834\u30d9\u30af\u30c8\u30eb\u306e\u59cb\u70b9<\/span>\r\nset<\/span> print<\/span> \"xyz.txt\"<\/span>\r\ndo<\/span> for<\/span> [<\/span>z<\/span>=<\/span>0<\/span>:<\/span>2<\/span>]{<\/span>\r\n  do<\/span> for<\/span> [<\/span>i<\/span>=<\/span>0<\/span>:<\/span>11<\/span>:<\/span>1<\/span>]{<\/span>\r\n    theta<\/span> =<\/span> 2<\/span>*<\/span>pi<\/span>\/<\/span>12<\/span> *<\/span> i<\/span>\r\n    do<\/span> for<\/span> [<\/span>j<\/span>=<\/span>3<\/span>:<\/span>6<\/span>:<\/span>3<\/span>]<\/span> {<\/span>\r\n      r<\/span> =<\/span> 0.5<\/span>*<\/span>j<\/span>\r\n      x<\/span> =<\/span> r<\/span> *<\/span> cos<\/span>(<\/span>theta<\/span>)<\/span>\r\n      y<\/span> =<\/span> r<\/span> *<\/span> sin<\/span>(<\/span>theta<\/span>)<\/span>\r\n      print<\/span> sprintf<\/span>(<\/span>\" %8.4f %8.4f %8.4f\"<\/span>,<\/span> x<\/span>,<\/span> y<\/span>,<\/span> z<\/span>)<\/span>\r\n    }<\/span>\r\n    print<\/span> \"\"<\/span>\r\n  }<\/span>\r\n    print<\/span> \"\"<\/span>\r\n}<\/span>\r\nset<\/span> print<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[41]:<\/div>\n
\n
\n
\n
%gnu<\/span>plot<\/span> inline<\/span> svg<\/span> size<\/span> 640<\/span>,<\/span>640<\/span> fixed<\/span> enhanced<\/span> font<\/span> 'Noto Sans CJK JP,16'<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[45]:<\/div>\n
\n
\n
\n
reset<\/span>\r\n\r\nEx<\/span>(x,y) =<\/span> x<\/span>\/<\/span>(<\/span>x<\/span>**<\/span>2<\/span> +<\/span> y<\/span>**<\/span>2<\/span>)<\/span>\r\nEy<\/span>(x,y) =<\/span> y<\/span>\/<\/span>(<\/span>x<\/span>**<\/span>2<\/span> +<\/span> y<\/span>**<\/span>2<\/span>)<\/span>\r\nE<\/span>(x, y) =<\/span> sqrt<\/span>(<\/span>Ex<\/span>(<\/span>x<\/span>,<\/span>y<\/span>)<\/span>**<\/span>2<\/span> +<\/span> Ey<\/span>(<\/span>x<\/span>,<\/span>y<\/span>)<\/span>**<\/span>2<\/span>)<\/span>\r\nex<\/span>(x,y)<\/span> =<\/span> scaling*Ex(x,y)\/E(x,<\/span> y)<\/span>\r\ney<\/span>(x,y) =<\/span> scaling<\/span>*<\/span>Ey<\/span>(<\/span>x<\/span>,<\/span>y<\/span>)<\/span>\/<\/span>E<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\r\n\r\nunset<\/span> xtics<\/span>\r\nunset<\/span> ytics<\/span>\r\nunset<\/span> ztics<\/span>\r\n# unset border<\/span>\r\nset<\/span> xyplane<\/span> 0<\/span>\r\nset<\/span> zeroaxis<\/span>\r\nset<\/span> xrange<\/span> [<\/span>-4<\/span>:<\/span>4<\/span>]<\/span>\r\nset<\/span> yrange<\/span> [<\/span>-4<\/span>:<\/span>4<\/span>]<\/span>\r\nunset<\/span> colorbox<\/span>\r\nset<\/span> palette<\/span> cubehelix<\/span> start<\/span> 0<\/span> cycles<\/span> 0<\/span> saturation<\/span> 3<\/span> negative<\/span>\r\n\r\n# \u7dda\u96fb\u8377<\/span>\r\nset<\/span> arrow<\/span> 1<\/span> from<\/span> 0<\/span>,<\/span>0<\/span>,<\/span>-0<\/span>.2<\/span> to<\/span> 0<\/span>,<\/span>0<\/span>,<\/span>2.2<\/span> lw<\/span> 4<\/span> lc<\/span> \"black\"<\/span> lt<\/span> 3<\/span>  nohead<\/span>\r\nset<\/span> key<\/span> inside<\/span>\r\nset<\/span> title<\/span> \"\u4e00\u69d8\u306a\u7dda\u96fb\u8377\u306b\u3088\u308b\u96fb\u5834\"<\/span>\r\nscaling<\/span> =<\/span> 0.1<\/span>\r\nset<\/span> view<\/span> 65<\/span>,<\/span>15<\/span>,,<\/span>\r\nsplot<\/span> \"xyz-line.txt\"<\/span> u<\/span> 1<\/span>:<\/span>2<\/span>:<\/span>3<\/span>:<\/span>(<\/span>log<\/span>(<\/span>E<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>)))<\/span> w<\/span> l<\/span> lw<\/span> 3<\/span> lc<\/span> palette<\/span> notitle<\/span>,<\/span> \\\r\n      \"xyz.txt\"<\/span> u<\/span> 1<\/span>:<\/span>2<\/span>:<\/span>3<\/span>:<\/span>(<\/span>ex<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>))<\/span>:<\/span>(<\/span>ey<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>))<\/span>:<\/span>(<\/span>0<\/span>)<\/span>:<\/span>(<\/span>log<\/span>(<\/span>E<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>)))<\/span> \\\r\n        w<\/span> vec<\/span> filled<\/span> head<\/span> lw<\/span> 6<\/span> lc<\/span> palette<\/span> title<\/span> \"\u96fb\u5834 E\"<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
\n
<\/div>\n
\n
\"\"<\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[9]:<\/div>\n
\n
\n
\n
set<\/span> output<\/span> \".\/vec-fig-sen.svg\"<\/span>\r\nreplot<\/span>\r\nset<\/span> output<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
<\/div>\n
\n
\n
\u4e00\u69d8\u306a\u9762\u96fb\u8377\u306b\u3088\u308b\u96fb\u5834<\/h3>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[10]:<\/div>\n
\n
\n
\n
%gnu<\/span>plot<\/span> inline<\/span> svg<\/span> size<\/span> 640<\/span>,<\/span>480<\/span> fixed<\/span> enhanced<\/span> font<\/span> 'Noto Sans CJK JP,16'<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[11]:<\/div>\n
\n
\n
\n
reset<\/span>\r\nset<\/span> parametric<\/span>\r\nset<\/span> isosamples<\/span> 30<\/span>\r\nset<\/span> hidden3d<\/span>\r\nset<\/span> xrange<\/span> [<\/span>-3<\/span>:<\/span>3<\/span>]<\/span>\r\nset<\/span> xrange<\/span> [<\/span>-3<\/span>:<\/span>3<\/span>]<\/span>\r\nset<\/span> zrange<\/span> [<\/span>-3<\/span>:<\/span>3<\/span>]<\/span>\r\nunset<\/span> xtics<\/span>\r\nunset<\/span> ytics<\/span>\r\nunset<\/span> ztics<\/span>\r\nunset<\/span> border<\/span>\r\nset<\/span> zeroaxis<\/span>\r\nset<\/span> xyplane<\/span> 0<\/span>\r\nset<\/span> isosamples<\/span> 60<\/span>\r\nset<\/span> view<\/span> ,<\/span>20<\/span>,,<\/span>\r\n\r\n# \u4e00\u69d8\u306a\u96fb\u5834\u30d9\u30af\u30c8\u30eb\u3092 arrow \u3067\u63cf\u304f<\/span>\r\ndo<\/span> for<\/span> [<\/span>i<\/span>=<\/span>1<\/span>:<\/span>5<\/span>]<\/span> {<\/span>\r\n  zi<\/span> =<\/span> -3<\/span> +<\/span> i<\/span>\r\n  # x > 0 \u5074<\/span>\r\n  set<\/span> arrow<\/span> i<\/span> from<\/span> 0.5<\/span>,<\/span>0<\/span>,<\/span>zi<\/span> to<\/span> 3<\/span>,<\/span>0<\/span>,<\/span>zi<\/span> filled<\/span> head<\/span> lc<\/span> \"blue\"<\/span> lw<\/span> 4<\/span>\r\n  # x < 0 \u5074<\/span>\r\n  set<\/span> arrow<\/span> i<\/span>+<\/span>10<\/span> from<\/span> -1<\/span>.2<\/span>,<\/span>0<\/span>,<\/span>zi<\/span> to<\/span> -3<\/span>,<\/span>0<\/span>,<\/span>zi<\/span> filled<\/span> head<\/span> lc<\/span> \"blue\"<\/span> lw<\/span> 4<\/span>\r\n}<\/span>\r\n\r\nset<\/span> title<\/span> \"\u4e00\u69d8\u306a\u9762\u96fb\u8377\u306b\u3088\u308b\u96fb\u5834\"<\/span>\r\nsplot<\/span> 0<\/span>,<\/span> u<\/span>,<\/span> v<\/span> lc<\/span> \"black\"<\/span> lw<\/span> 1<\/span> notitle<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
\n
<\/div>\n
\n
\"\"<\/p>\n<\/div>\n<\/div>\n
\n
Out[11]:<\/div>\n
\n
\tdummy variable is t for curves, u\/v for surfaces\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[12]:<\/div>\n
\n
\n
\n
set<\/span> output<\/span> \".\/vec-fig-men.svg\"<\/span>\r\nreplot<\/span>\r\nset<\/span> output<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
<\/div>\n
\n
\n
\u76f4\u7dda\u96fb\u6d41\u306b\u3088\u308b\u78c1\u5834 $\\boldsymbol{B}$<\/h3>\n
plot \u7248<\/h4>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[13]:<\/div>\n
\n
\n
\n
%gnu<\/span>plot<\/span> inline<\/span> svg<\/span> size<\/span> 480<\/span>,<\/span>480<\/span> fixed<\/span> enhanced<\/span> font<\/span> 'Noto Sans CJK JP,16'<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[14]:<\/div>\n
\n
\n
\n
reset<\/span>\r\nunset<\/span> xtics<\/span>\r\nunset<\/span> ytics<\/span>\r\nunset<\/span> border<\/span>\r\nset<\/span> zeroaxis<\/span>\r\nset<\/span> xrange<\/span> [<\/span>-3<\/span>.5<\/span>:<\/span>3.5<\/span>]<\/span>\r\nset<\/span> yrange<\/span> [<\/span>-3<\/span>.5<\/span>:<\/span>3.5<\/span>]<\/span>\r\nset<\/span> size<\/span> ratio<\/span> 1<\/span>\r\n\r\n# \u78c1\u529b\u7dda\u7528\u5ea7\u6a19\u30d5\u30a1\u30a4\u30eb<\/span>\r\nset<\/span> print<\/span> \"xy-line.txt\"<\/span>\r\ndo<\/span> for<\/span> [<\/span>i<\/span>=<\/span>0<\/span>:<\/span>6<\/span>]<\/span> {<\/span>\r\n  # \u96fb\u7dda\u306b\u8fd1\u3044\u65b9\u304c\u5bc6\u3068\u306a\u308b\u3088\u3046\u306b\u3057\u3066\u307f\u308b<\/span>\r\n  r<\/span> =<\/span> 3<\/span>*<\/span>(<\/span>2.<\/span>\/<\/span>3<\/span>)<\/span>**<\/span>i<\/span>\r\n  do<\/span> for<\/span> [<\/span>j<\/span>=<\/span>0<\/span>:<\/span> 64<\/span>]<\/span> {<\/span>\r\n    theta<\/span> =<\/span> pi<\/span>\/<\/span>32<\/span> *<\/span> j<\/span>\r\n    x<\/span> =<\/span> r<\/span> *<\/span> cos<\/span>(<\/span>theta<\/span>)<\/span>\r\n    y<\/span> =<\/span> r<\/span> *<\/span> sin<\/span>(<\/span>theta<\/span>)<\/span>\r\n    print<\/span> sprintf<\/span>(<\/span>\"%8.4f  %8.4f\"<\/span>,<\/span> x<\/span>,<\/span> y<\/span>)<\/span>\r\n  }<\/span>\r\n  # with lines \u3059\u308b\u306e\u3067\u9069\u5b9c\u7a7a\u884c\u3092\u3044\u308c\u308b<\/span>\r\n  print<\/span> \" \"<\/span>\r\n}<\/span>\r\nset<\/span> print<\/span>\r\n\r\n# \u30d9\u30af\u30c8\u30eb\u306e\u59cb\u70b9\u7528\u5ea7\u6a19\u30d5\u30a1\u30a4\u30eb<\/span>\r\n# \u77e2\u5370\u306e\u6df7\u96d1\u9632\u6b62\u306e\u305f\u3081\u9069\u5b9c\u9593\u5f15\u304f<\/span>\r\nset<\/span> print<\/span> \"xy.txt\"<\/span>\r\ndo<\/span> for<\/span> [<\/span>i<\/span>=<\/span>0<\/span>:<\/span>5<\/span>]<\/span> {<\/span>\r\n  r<\/span> =<\/span> 3<\/span>*<\/span>(<\/span>2.<\/span>\/<\/span>3<\/span>)<\/span>**<\/span>i<\/span>\r\n  do<\/span> for<\/span> [<\/span>j<\/span>=<\/span>0<\/span>:<\/span> 3<\/span>]<\/span> {<\/span>\r\n    theta<\/span> =<\/span> pi<\/span> *<\/span> j<\/span> +<\/span> i<\/span>*<\/span>pi<\/span>\/<\/span>2<\/span>\r\n    x<\/span> =<\/span> r<\/span> *<\/span> cos<\/span>(<\/span>theta<\/span>)<\/span>\r\n    y<\/span> =<\/span> r<\/span> *<\/span> sin<\/span>(<\/span>theta<\/span>)<\/span>\r\n    print<\/span> sprintf<\/span>(<\/span>\"%8.4f  %8.4f\"<\/span>,<\/span> x<\/span>,<\/span> y<\/span>)<\/span>\r\n  }<\/span>\r\n}<\/span>\r\nset<\/span> print<\/span>\r\n\r\n# \u78c1\u5834 B<\/span>\r\nBx<\/span>(x, y) =<\/span> -<\/span> y<\/span>\/<\/span>(<\/span>x<\/span>**<\/span>2<\/span> +<\/span> y<\/span>**<\/span>2<\/span>)<\/span>\r\nBy<\/span>(x, y) =<\/span>   x<\/span>\/<\/span>(<\/span>x<\/span>**<\/span>2<\/span> +<\/span> y<\/span>**<\/span>2<\/span>)<\/span>\r\n# \u78c1\u5834 B \u306e\u5927\u304d\u3055<\/span>\r\nB<\/span>(x, y) =<\/span> sqrt<\/span>(<\/span>Bx<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>**<\/span>2<\/span> +<\/span> By<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>**<\/span>2<\/span>)<\/span>\r\n# \u898f\u683c\u5316\u3055\u308c\u305f hat B<\/span>\r\nhBx<\/span>(x, y) =<\/span> Bx<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\/<\/span>B<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\r\nhBy<\/span>(x, y) =<\/span> By<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\/<\/span>B<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\r\n\r\n# \u77e2\u5370\u306e\u9577\u3055\u8abf\u6574\u7528<\/span>\r\nscaling<\/span> =<\/span> 0.1<\/span>\r\n\r\n# \u96fb\u7dda\u306e\u4f4d\u7f6e<\/span>\r\nset<\/span> label<\/span> 1<\/span> point<\/span> pt<\/span> 6<\/span> ps<\/span> 1.5<\/span> lc<\/span> \"red\"<\/span> at<\/span> 0<\/span>,<\/span> 0<\/span>\r\nset<\/span> label<\/span> 2<\/span> center<\/span> at<\/span> first<\/span> 0<\/span>,<\/span>0<\/span> \"\u30fb\"<\/span> tc<\/span> \"red\"<\/span>\r\n\r\nset<\/span> title<\/span> \"\u76f4\u7dda\u96fb\u6d41\u306b\u3088\u308b\u78c1\u529b\u7dda\"<\/span>\r\nunset<\/span> colorbox<\/span>\r\n# green.pal \u6539<\/span>\r\n# https:\/\/github.com\/Gnuplotting\/gnuplot-palettes<\/span>\r\nset<\/span> palette<\/span> defined<\/span> (<\/span>\\\r\n  0<\/span> '#E5F5E0'<\/span>,<\/span>\\\r\n  1<\/span> '#C7E9C0'<\/span>,<\/span>\\\r\n  2<\/span> '#A1D99B'<\/span>,<\/span>\\\r\n  3<\/span> '#74C476'<\/span>,<\/span>\\\r\n  4<\/span> '#41AB5D'<\/span>,<\/span>\\\r\n  5<\/span> '#238B45'<\/span>,<\/span>\\\r\n  6<\/span> '#005A32'<\/span>,<\/span>\\\r\n  7<\/span> '#004000'<\/span> \\\r\n)<\/span>\r\nplot<\/span> \"xy-line.txt\"<\/span> u<\/span> 1<\/span>:<\/span>2<\/span>:<\/span>(<\/span>log<\/span>(<\/span>B<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>)))<\/span> w<\/span> l<\/span> lw<\/span> 3<\/span> lc<\/span> palette<\/span> notitle<\/span>,<\/span>\\\r\n  \"xy.txt\"<\/span> u<\/span> 1<\/span>:<\/span>2<\/span>:<\/span>(<\/span>scaling<\/span>*<\/span>hBx<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>))<\/span>:<\/span>(<\/span>scaling<\/span>*<\/span>hBy<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>))<\/span>:<\/span>(<\/span>log<\/span>(<\/span>B<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>)))<\/span> \\\r\n  w<\/span> vec<\/span> notitle<\/span> lw<\/span> 7<\/span> lc<\/span> palette<\/span> filled<\/span> head<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
\n
<\/div>\n
\n
\"\"<\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[15]:<\/div>\n
\n
\n
\n
set<\/span> output<\/span> \".\/vec-fig07d.svg\"<\/span>\r\nreplot<\/span>\r\nset<\/span> output<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
<\/div>\n
\n
\n
splot \u7248<\/h4>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[16]:<\/div>\n
\n
\n
\n
%gnu<\/span>plot<\/span> inline<\/span> svg<\/span> size<\/span> 640<\/span>,<\/span>640<\/span> fixed<\/span> enhanced<\/span> font<\/span> 'Noto Sans CJK JP,16'<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[38]:<\/div>\n
\n
\n
\n
reset<\/span>\r\nunset<\/span> xtics<\/span>\r\nunset<\/span> ytics<\/span>\r\nunset<\/span> ztics<\/span>\r\n#unset border<\/span>\r\nset<\/span> zeroaxis<\/span>\r\nset<\/span> xrange<\/span> [<\/span>-3<\/span>.5<\/span>:<\/span>3.5<\/span>]<\/span>\r\nset<\/span> yrange<\/span> [<\/span>-3<\/span>.5<\/span>:<\/span>3.5<\/span>]<\/span>\r\n#set size ratio 1<\/span>\r\n\r\n# \u78c1\u529b\u7dda\u7528\u5ea7\u6a19\u30d5\u30a1\u30a4\u30eb<\/span>\r\n# \u96fb\u6d41\u306e\u5f8c\u308d\u5074<\/span>\r\nset<\/span> print<\/span> \"xyz-linep.txt\"<\/span>\r\ndo<\/span> for<\/span> [<\/span>z<\/span>=<\/span>0<\/span>:<\/span>2<\/span>]<\/span> {<\/span>\r\n  do<\/span> for<\/span> [<\/span>i<\/span>=<\/span>0<\/span>:<\/span>6<\/span>]<\/span> {<\/span>\r\n    r<\/span> =<\/span> 3<\/span>*<\/span>(<\/span>2.<\/span>\/<\/span>3<\/span>)<\/span>**<\/span>i<\/span>\r\n    do<\/span> for<\/span> [<\/span>j<\/span>=<\/span>0<\/span>:<\/span> 32<\/span>]<\/span> {<\/span>\r\n    theta<\/span> =<\/span> pi<\/span>\/<\/span>32<\/span> *<\/span> j<\/span>\r\n      x<\/span> =<\/span> r<\/span> *<\/span> cos<\/span>(<\/span>theta<\/span>)<\/span>\r\n      y<\/span> =<\/span> r<\/span> *<\/span> sin<\/span>(<\/span>theta<\/span>)<\/span>\r\n      print<\/span> sprintf<\/span>(<\/span>\"%8.4f  %8.4f  %8.4f\"<\/span>,<\/span> x<\/span>,<\/span> y<\/span>,<\/span> z<\/span>)<\/span>\r\n      # with lines \u3059\u308b\u306e\u3067\u9069\u5b9c\u7a7a\u884c\u3092\u3044\u308c\u308b<\/span>\r\n      print<\/span> \"\"<\/span>\r\n    }<\/span>\r\n    # with lines \u3059\u308b\u306e\u3067\u9069\u5b9c\u7a7a\u884c\u3092\u3044\u308c\u308b<\/span>\r\n    print<\/span> \"\"<\/span>\r\n  }<\/span>\r\n  # with lines \u3059\u308b\u306e\u3067\u9069\u5b9c\u7a7a\u884c\u3092\u3044\u308c\u308b<\/span>\r\n  print<\/span> \"\"<\/span>\r\n}<\/span>\r\nset<\/span> print<\/span>\r\n\r\n# \u78c1\u529b\u7dda\u7528\u5ea7\u6a19\u30d5\u30a1\u30a4\u30eb<\/span>\r\n# \u96fb\u6d41\u306e\u524d\u5074<\/span>\r\nset<\/span> print<\/span> \"xyz-linem.txt\"<\/span>\r\ndo<\/span> for<\/span> [<\/span>k<\/span>=<\/span>0<\/span>:<\/span>2<\/span>]<\/span> {<\/span>\r\n  do<\/span> for<\/span> [<\/span>i<\/span>=<\/span>0<\/span>:<\/span> 6<\/span>]<\/span> {<\/span>\r\n    r<\/span> =<\/span> 3<\/span>*<\/span>(<\/span>2.<\/span>\/<\/span>3<\/span>)<\/span>**<\/span>i<\/span>\r\n    do<\/span> for<\/span> [<\/span>j<\/span>=<\/span>32<\/span>:<\/span> 64<\/span>]<\/span> {<\/span>\r\n    theta<\/span> =<\/span> pi<\/span>\/<\/span>32<\/span> *<\/span> j<\/span>\r\n      x<\/span> =<\/span> r<\/span> *<\/span> cos<\/span>(<\/span>theta<\/span>)<\/span>\r\n      y<\/span> =<\/span> r<\/span> *<\/span> sin<\/span>(<\/span>theta<\/span>)<\/span>\r\n      z<\/span> =<\/span> 1<\/span>*<\/span>k<\/span>\r\n      print<\/span> sprintf<\/span>(<\/span>\"%8.4f  %8.4f  %8.4f\"<\/span>,<\/span> x<\/span>,<\/span> y<\/span>,<\/span> z<\/span>)<\/span>\r\n      # with lines \u3059\u308b\u306e\u3067\u9069\u5b9c\u7a7a\u884c\u3092\u3044\u308c\u308b<\/span>\r\n      print<\/span> \"\"<\/span>\r\n    }<\/span>\r\n    # with lines \u3059\u308b\u306e\u3067\u9069\u5b9c\u7a7a\u884c\u3092\u3044\u308c\u308b<\/span>\r\n    print<\/span> \"\"<\/span>\r\n  }<\/span>\r\n  # with lines \u3059\u308b\u306e\u3067\u9069\u5b9c\u7a7a\u884c\u3092\u3044\u308c\u308b<\/span>\r\n  print<\/span> \"\"<\/span>\r\n}<\/span>\r\nset<\/span> print<\/span>\r\n\r\n# \u30d9\u30af\u30c8\u30eb\u306e\u59cb\u70b9\u7528\u5ea7\u6a19\u30d5\u30a1\u30a4\u30eb<\/span>\r\n# \u77e2\u5370\u306e\u6df7\u96d1\u9632\u6b62\u306e\u305f\u3081\u9069\u5b9c\u9593\u5f15\u304f<\/span>\r\nset<\/span> print<\/span> \"xyz.txt\"<\/span>\r\ndo<\/span> for<\/span> [<\/span>k<\/span>=<\/span>0<\/span>:<\/span>2<\/span>]<\/span> {<\/span>\r\n  do<\/span> for<\/span> [<\/span>i<\/span>=<\/span>0<\/span>:<\/span>5<\/span>:<\/span>2<\/span>]<\/span> {<\/span>\r\n    r<\/span> =<\/span> 3<\/span>*<\/span>(<\/span>2.<\/span>\/<\/span>3<\/span>)<\/span>**<\/span>i<\/span>\r\n    do<\/span> for<\/span> [<\/span>j<\/span>=<\/span>0<\/span>:<\/span>3<\/span>]<\/span> {<\/span>\r\n      theta<\/span> =<\/span> pi<\/span> *<\/span> j<\/span> +<\/span> i<\/span>*<\/span>pi<\/span>\/<\/span>2<\/span> \r\n      x<\/span> =<\/span> r<\/span> *<\/span> cos<\/span>(<\/span>theta<\/span>)<\/span>\r\n      y<\/span> =<\/span> r<\/span> *<\/span> sin<\/span>(<\/span>theta<\/span>)<\/span>\r\n      z<\/span> =<\/span> 1<\/span>*<\/span>k<\/span>\r\n      print<\/span> sprintf<\/span>(<\/span>\"%8.4f  %8.4f  %8.4f\"<\/span>,<\/span> x<\/span>,<\/span> y<\/span>,<\/span> z<\/span>)<\/span>\r\n    }<\/span>\r\n  }<\/span>\r\n}<\/span>\r\nset<\/span> print<\/span>\r\n\r\n# \u78c1\u5834 B<\/span>\r\nBx<\/span>(x, y) =<\/span> -<\/span> y<\/span>\/<\/span>(<\/span>x<\/span>**<\/span>2<\/span> +<\/span> y<\/span>**<\/span>2<\/span>)<\/span>\r\nBy<\/span>(x, y) =<\/span>   x<\/span>\/<\/span>(<\/span>x<\/span>**<\/span>2<\/span> +<\/span> y<\/span>**<\/span>2<\/span>)<\/span>\r\n# \u78c1\u5834 B \u306e\u5927\u304d\u3055<\/span>\r\nB<\/span>(x, y) =<\/span> sqrt<\/span>(<\/span>Bx<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>**<\/span>2<\/span> +<\/span> By<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>**<\/span>2<\/span>)<\/span>\r\n# \u898f\u683c\u5316\u3055\u308c\u305f hat B<\/span>\r\nhBx<\/span>(x, y) =<\/span> Bx<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\/<\/span>B<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\r\nhBy<\/span>(x, y) =<\/span> By<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\/<\/span>B<\/span>(<\/span>x<\/span>,<\/span> y<\/span>)<\/span>\r\n\r\n# \u77e2\u5370\u306e\u9577\u3055\u8abf\u6574\u7528<\/span>\r\nscaling<\/span> =<\/span> 0.1<\/span>\r\n\r\n# \u96fb\u7dda<\/span>\r\n# set arrow 1 from 0,0,0 to 0,0,3.8 lw 3.5 lc \"red\" lt 4 filled head<\/span>\r\nset<\/span> print<\/span> \"densen.txt\"<\/span>\r\n  print<\/span> 0<\/span>,<\/span>0<\/span>,<\/span>0<\/span>,<\/span>0<\/span>,<\/span>0<\/span>,<\/span>2.8<\/span> \r\nset<\/span> print<\/span>\r\n\r\nset<\/span> title<\/span> \"\u76f4\u7dda\u96fb\u6d41\u306b\u3088\u308b\u78c1\u5834\"<\/span>\r\nunset<\/span> colorbox<\/span>\r\n# green.pal \u6539<\/span>\r\n# https:\/\/github.com\/Gnuplotting\/gnuplot-palettes<\/span>\r\nset<\/span> palette<\/span> defined<\/span> (<\/span>\\\r\n  0<\/span> '#E5F5E0'<\/span>,<\/span>\\\r\n  1<\/span> '#C7E9C0'<\/span>,<\/span>\\\r\n  2<\/span> '#A1D99B'<\/span>,<\/span>\\\r\n  3<\/span> '#74C476'<\/span>,<\/span>\\\r\n  4<\/span> '#41AB5D'<\/span>,<\/span>\\\r\n  5<\/span> '#238B45'<\/span>,<\/span>\\\r\n  6<\/span> '#005A32'<\/span>,<\/span>\\\r\n  7<\/span> '#004000'<\/span> \\\r\n)<\/span>\r\nset<\/span> key<\/span> inside<\/span> sample<\/span> 1<\/span>\r\nset<\/span> view<\/span> 70<\/span>,,,<\/span>\r\nset<\/span> xyplane<\/span> 0<\/span>\r\nset<\/span> zrange<\/span> [<\/span>0<\/span>:<\/span>3<\/span>]<\/span>\r\nsplot<\/span> \\\r\n  \"xyz-linep.txt\"<\/span> u<\/span> 1<\/span>:<\/span>2<\/span>:<\/span>3<\/span>:<\/span>(<\/span>log<\/span>(<\/span>B<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>)))<\/span> w<\/span> l<\/span> lw<\/span> 3<\/span> lc<\/span> palette<\/span> notitle<\/span>,<\/span>\\\r\n  \"densen.txt\"<\/span> w<\/span> vec<\/span> lc<\/span> \"red\"<\/span> lw<\/span> 4<\/span> filled<\/span> head<\/span> title<\/span> \"\u96fb\u6d41\"<\/span>,<\/span> \\\r\n  \"xyz.txt\"<\/span> \\\r\n    u<\/span> 1<\/span>:<\/span>2<\/span>:<\/span>3<\/span>:<\/span>(<\/span>scaling<\/span>*<\/span>Bx<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>))<\/span>:<\/span>(<\/span>scaling<\/span>*<\/span>By<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>))<\/span>:<\/span>(<\/span>0<\/span>)<\/span>:<\/span>(<\/span>log<\/span>(<\/span>B<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>)))<\/span> \\\r\n    w<\/span> vec<\/span> title<\/span> \"\u78c1\u5834\"<\/span> lw<\/span> 6<\/span> lc<\/span> palette<\/span> filled<\/span> head<\/span>,<\/span> \\\r\n  \"xyz-linem.txt\"<\/span> u<\/span> 1<\/span>:<\/span>2<\/span>:<\/span>3<\/span>:<\/span>(<\/span>log<\/span>(<\/span>B<\/span>(<\/span>$<\/span>1<\/span>,<\/span>$<\/span>2<\/span>)))<\/span> w<\/span> l<\/span> lw<\/span> 3<\/span> lc<\/span> palette<\/span> notitle<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
\n
<\/div>\n
\n
\"\"<\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
\n
In\u00a0[18]:<\/div>\n
\n
\n
\n
set<\/span> output<\/span> \".\/vec-fig08d.svg\"<\/span>\r\nreplot<\/span>\r\nset<\/span> output<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"excerpt":{"rendered":"
\u300cgnuplot \u3067\u30d9\u30af\u30c8\u30eb\u5834\uff08\u96fb\u5834\u30fb\u78c1\u5834\uff09\u3092\u63cf\u304f\u300d\u3084\u300cpalette \u3067\u30d9\u30af\u30c8\u30eb\u306e\u8272\u3092\u5909\u3048\u3066 gnuplot \u3067\u96fb\u5834\u30fb\u78c1\u5834\u3092\u63cf\u304f\u300d\u3067\u63cf\u3044\u3066\u307f\u305f\u304c\uff0c\u3082\u3046\u5c11\u3057\u3002<\/p>\n
gnuplot \u3060\u3051\u3067\u30d9\u30af\u30c8\u30eb\u306e\u59cb\u70b9\u306e\u5ea7\u6a19\u30d5\u30a1\u30a4\u30eb\u3092\u4f5c\u6210\u3057\uff0c\u305d\u308c\u3092\u8aad\u307f\u8fbc\u3093\u3067\u96fb\u6c17\u529b\u7dda\u30fb\u78c1\u529b\u7dda\u3092\u63cf\u304f\u4f8b\u3002<\/p>
\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-3643","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\/3643","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=3643"}],"version-history":[{"count":17,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/posts\/3643\/revisions"}],"predecessor-version":[{"id":3968,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/posts\/3643\/revisions\/3968"}],"wp:attachment":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/media?parent=3643"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/categories?post=3643"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/tags?post=3643"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}