{"id":3619,"date":"2022-08-31T10:37:22","date_gmt":"2022-08-31T01:37:22","guid":{"rendered":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/?p=3619"},"modified":"2023-03-28T17:59:35","modified_gmt":"2023-03-28T08:59:35","slug":"palette-%e3%81%a7%e3%83%99%e3%82%af%e3%83%88%e3%83%ab%e3%81%ae%e8%89%b2%e3%82%92%e5%a4%89%e3%81%88%e3%81%a6-gnuplot-%e3%81%a7%e3%83%99%e3%82%af%e3%83%88%e3%83%ab%e5%a0%b4%ef%bc%88%e9%9b%bb%e5%a0%b4","status":"publish","type":"post","link":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/3619\/","title":{"rendered":"palette \u3067\u30d9\u30af\u30c8\u30eb\u306e\u8272\u3092\u5909\u3048\u3066 gnuplot \u3067\u96fb\u5834\u30fb\u78c1\u5834\u3092\u63cf\u304f"},"content":{"rendered":"<p>\u5148\u65e5\uff0c\u300c<a href=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/3595\/\">gnuplot \u3067\u30d9\u30af\u30c8\u30eb\u5834\uff08\u96fb\u5834\u30fb\u78c1\u5834\uff09\u3092\u63cf\u304f<\/a>\u300d\u3067\u63cf\u3044\u3066\u307f\u305f\u304c\uff0c\u3044\u308d\u3044\u308d\u691c\u7d22\u3059\u308b\u3068 palette \u3067\u30d9\u30af\u30c8\u30eb\u306e\u8272\u3092\uff08\u30d9\u30af\u30c8\u30eb\u306e\u5927\u304d\u3055\u306b\u5408\u308f\u305b\u3066\uff09\u53ef\u5909\u306b\u3067\u304d\u308b\u3088\u3046\u306a\u306e\u3067\uff0c\u5099\u5fd8\u9332\u3068\u3057\u3066\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<h3 id=\"\u70b9\u96fb\u8377\u304c\u3064\u304f\u308b\u96fb\u5834-$\\boldsymbol{E}$\">\u70b9\u96fb\u8377\u304c\u3064\u304f\u308b\u96fb\u5834 $\\boldsymbol{E}$<\/h3>\n<p>\u4f4d\u7f6e $\\boldsymbol{r}_1$ \u306b\u3042\u308b\u70b9\u96fb\u8377 $q_1$ \u304c\u3064\u304f\u308b\u96fb\u5834\u30d9\u30af\u30c8\u30eb $\\boldsymbol{E}$ \u306f<\/p>\n<p>\\begin{eqnarray}<br \/>\n\\boldsymbol{E}<br \/>\n&amp;=&amp; \\frac{q_1}{4\\pi \\varepsilon_0} \\frac{\\boldsymbol{r} &#8211; \\boldsymbol{r}_1}{|\\boldsymbol{r} &#8211; \\boldsymbol{r}_1|^3}\\end{eqnarray}<\/p>\n<p>$z=0$ \u306e $xy$ \u5e73\u9762\u30672\u6b21\u5143\u30d9\u30af\u30c8\u30eb\u5834\u3092\u63cf\u304f\u3002<\/p>\n<p>$$\\boldsymbol{r} = (x, y, 0), \\quad \\boldsymbol{r}_1 = (x_1, y_1, 0)$$$$E_x = s \\frac{x-x_1}{\\left((x-x_1)^2 + (y-y_1)^2\\right)^{3\/2}}$$$$E_y = s \\frac{y-y_1}{\\left((x-x_1)^2 + (y-y_1)^2\\right)^{3\/2}}$$<\/p>\n<p>$s$ \u306f $\\mbox{scaling}$ \u4fc2\u6570\u3002\u5b9f\u969b\u306e\u63cf\u753b\u306e\u969b\u306b\u30d9\u30af\u30c8\u30eb\u306e\u9577\u3055\u3092\u5fae\u8abf\u6574\u3059\u308b\u4fc2\u6570\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h4 id=\"\u6b63\u96fb\u8377\u304c\u3064\u304f\u308b\u96fb\u5834\">\u6b63\u96fb\u8377\u304c\u3064\u304f\u308b\u96fb\u5834<\/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 [1]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-gnuplot\">\n<pre><span class=\"c\"># \u96fb\u5834 E<\/span>\r\n<span class=\"nf\">Ex<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"n\">scaling<\/span> <span class=\"o\">*<\/span> <span class=\"n\">x<\/span><span class=\"o\">\/<\/span><span class=\"p\">(<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"n\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">3<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">Ey<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"n\">scaling<\/span> <span class=\"o\">*<\/span> <span class=\"n\">y<\/span><span class=\"o\">\/<\/span><span class=\"p\">(<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"n\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">3<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"c\"># \u96fb\u5834 E \u306e\u5927\u304d\u3055\u306f\u8ddd\u96e2\u306e2\u4e57\u306b\u53cd\u6bd4\u4f8b<\/span>\r\n<span class=\"nf\">E<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"mf\">1.<\/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\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=\"c\"># special filename \"++\" \u306e\u3068\u304d\u306f<\/span>\r\n<span class=\"c\"># urange, vrange \u3067\u63cf\u753b\u7bc4\u56f2\u306e\u8a2d\u5b9a<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">urange<\/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\">vrange<\/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=\"c\"># x \u8ef8\u65b9\u5411\u306e\u30b5\u30f3\u30d7\u30eb\u6570<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">samples<\/span> <span class=\"mi\">13<\/span>\r\n<span class=\"c\"># y \u8ef8\u65b9\u5411\u306e\u30b5\u30f3\u30d7\u30eb\u6570<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">isosamples<\/span> <span class=\"mi\">13<\/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\r\n<span class=\"c\"># \u6b63\u96fb\u8377\u306e\u4f4d\u7f6e\u306b\u9752\u4e38\u3068 +<\/span>\r\n<span class=\"nv\">x1<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"nv\">y1<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">1<\/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=\"mf\">1.5<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">\"blue\"<\/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\">\"+\"<\/span> <span class=\"n\">tc<\/span> <span class=\"s\">\"blue\"<\/span>\r\n\r\n<span class=\"k\">set<\/span> <span class=\"nb\">title<\/span> <span class=\"s\">\"\u6b63\u96fb\u8377\u304c\u3064\u304f\u308b\u96fb\u5834\"<\/span>\r\n\r\n<span class=\"c\"># colorbox \u3092\u975e\u8868\u793a\u306b<\/span>\r\n<span class=\"k\">unset<\/span> <span class=\"nb\">colorbox<\/span>\r\n<span class=\"c\"># palette \u306e cublehelix \u3067\u306e\u8a2d\u5b9a\u4f8b\uff1a <\/span>\r\n<span class=\"c\"># start 0: blue, 1: red, 2: green<\/span>\r\n<span class=\"c\"># cycles 0 \u3067\u540c\u4e00\u8272\u306e\u6fc3\u6de1\u306e\u307f\u306b<\/span>\r\n<span class=\"c\"># saturation 3 \u304c\u597d\u307f<\/span>\r\n<span class=\"c\"># negative \u3067\u5f37\u3044\u3068\u3053\u308d\u304c\u6fc3\u304f\u306a\u308b<\/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\">0<\/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\"># using \u59cb\u70b9\u306ex : \u59cb\u70b9\u306ey : \u30d9\u30af\u30c8\u30eb\u306ex\u6210\u5206 : \u30d9\u30af\u30c8\u30eb\u306ey\u6210\u5206 : variable color \u60c5\u5831<\/span>\r\n<span class=\"k\">plot<\/span> <span class=\"s\">\"++\"<\/span> <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\">Ex<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">-<\/span><span class=\"n\">x1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"o\">-<\/span><span class=\"n\">y1<\/span><span class=\"p\">))<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"nf\">Ey<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">-<\/span><span class=\"n\">x1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"o\">-<\/span><span class=\"n\">y1<\/span><span class=\"p\">))<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"nf\">log<\/span><span class=\"p\">(<\/span><span class=\"nf\">E<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">-<\/span><span class=\"n\">x1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"o\">-<\/span><span class=\"n\">y1<\/span><span class=\"p\">)))<\/span> \\\r\n     <span class=\"nb\">w<\/span> <span class=\"n\">vec<\/span> <span class=\"nb\">notitle<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">4<\/span> <span class=\"n\">lc<\/span> <span class=\"n\">palette<\/span> <span class=\"n\">filled<\/span> <span class=\"n\">head<\/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 wp-image-3630\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/vec-fig01a.svg\" alt=\"\" width=\"480\" height=\"480\" \/><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h3 id=\"$\\boldsymbol{E}$-\u306e\u5411\u304d\u3092\u8868\u3059\u5358\u4f4d\u30d9\u30af\u30c8\u30eb\u5834\">$\\boldsymbol{E}$ \u306e\u5411\u304d\u3092\u8868\u3059\u5358\u4f4d\u30d9\u30af\u30c8\u30eb\u5834<\/h3>\n<p>\u96fb\u5834 $\\boldsymbol{E}$ \u306e\u5927\u304d\u3055\u306f\u96fb\u8377\u304b\u3089\u306e\u8ddd\u96e2\u306e2\u4e57\u306b\u53cd\u6bd4\u4f8b\u3059\u308b\u306e\u3067\uff0c\u4e0a\u56f3\u306e\u3088\u3046\u306b\u30d9\u30af\u30c8\u30eb\u306e\u9577\u3055\u304c\u5927\u304d\u304f\u5909\u5316\u3059\u308b\u3002\u3053\u308c\u306f\u3053\u308c\u3067\u6b63\u3057\u3044\u306e\u3060\u304c\uff0c\u898b\u3084\u3059\u304f\u3059\u308b\u305f\u3081\u306b\uff0c\u4ee5\u4e0b\u306e\u3088\u3046\u306b $\\boldsymbol{E}$ \u306e\u5411\u304d\u3092\u8868\u3059\u5358\u4f4d\u30d9\u30af\u30c8\u30eb\u5834\uff08\u5927\u304d\u3055\u304c $1$\uff09\u3092\u5b9a\u7fa9\u3057\u3066\uff0c\u30d9\u30af\u30c8\u30eb\u306e\u5411\u304d\u306e\u307f\u3092\u6b63\u3057\u304f\u63cf\u753b\u3059\u308b\u3002<\/p>\n<p>$$\\hat{\\boldsymbol{E}} \\equiv \\frac{\\boldsymbol{E}}{\\sqrt{\\boldsymbol{E}\\cdot\\boldsymbol{E}}} =<br \/>\n\\frac{\\boldsymbol{r} &#8211; \\boldsymbol{r}_1}{|\\boldsymbol{r} &#8211; \\boldsymbol{r}_1|}$$$$\\hat{E}_x = s \\frac{x-x_1}{\\sqrt{(x-x_1)^2 + (y-y_1)^2}}$$$$\\hat{E}_y = s \\frac{y-y_1}{\\sqrt{(x-x_1)^2 + (y-y_1)^2}}$$<\/p>\n<p>$s$ \u306f $\\mbox{scaling}$ \u4fc2\u6570\u3002\u5b9f\u969b\u306e\u63cf\u753b\u306e\u969b\u306b\u30d9\u30af\u30c8\u30eb\u306e\u9577\u3055\u3092\u5fae\u8abf\u6574\u3059\u308b\u4fc2\u6570\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h4 id=\"\u6b63\u96fb\u8377\u304c\u3064\u304f\u308b\u96fb\u5834\u306e\u5411\u304d\">\u6b63\u96fb\u8377\u304c\u3064\u304f\u308b\u96fb\u5834\u306e\u5411\u304d<\/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 [2]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-gnuplot\">\n<pre><span class=\"c\"># \u96fb\u5834 E<\/span>\r\n<span class=\"nf\">Ex<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"n\">x<\/span><span class=\"o\">\/<\/span><span class=\"p\">(<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"n\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">3<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">Ey<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"n\">y<\/span><span class=\"o\">\/<\/span><span class=\"p\">(<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"n\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">3<\/span><span class=\"p\">)<\/span>\r\n<span class=\"c\"># \u96fb\u5834 E \u306e\u5927\u304d\u3055<\/span>\r\n<span class=\"nf\">E<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"nf\">Ex<\/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\">Ey<\/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 E<\/span>\r\n<span class=\"nf\">hEx<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"n\">scaling<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">Ex<\/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\">E<\/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\">hEy<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"n\">scaling<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">Ey<\/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\">E<\/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=\"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=\"c\"># special filename \"++\" \u306e\u3068\u304d\u306f<\/span>\r\n<span class=\"c\"># urange, vrange \u3067\u63cf\u753b\u7bc4\u56f2\u306e\u8a2d\u5b9a<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">urange<\/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\">vrange<\/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=\"c\"># x \u8ef8\u65b9\u5411\u306e\u30b5\u30f3\u30d7\u30eb\u6570<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">samples<\/span> <span class=\"mi\">19<\/span>\r\n<span class=\"c\"># y \u8ef8\u65b9\u5411\u306e\u30b5\u30f3\u30d7\u30eb\u6570<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">isosamples<\/span> <span class=\"mi\">19<\/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\r\n<span class=\"c\"># \u6b63\u96fb\u8377\u306e\u4f4d\u7f6e\u306b\u9752\u4e38\u3068 +<\/span>\r\n<span class=\"nv\">x1<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"nv\">y1<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">1<\/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=\"mf\">1.5<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">\"blue\"<\/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\">\"+\"<\/span> <span class=\"n\">tc<\/span> <span class=\"s\">\"blue\"<\/span>\r\n\r\n<span class=\"k\">set<\/span> <span class=\"nb\">title<\/span> <span class=\"s\">\"\u6b63\u96fb\u8377\u304c\u3064\u304f\u308b\u96fb\u5834\u306e\u5411\u304d\"<\/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\">0<\/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=\"k\">plot<\/span> <span class=\"s\">\"++\"<\/span> <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\">hEx<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">-<\/span><span class=\"n\">x1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"o\">-<\/span><span class=\"n\">y1<\/span><span class=\"p\">))<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"nf\">hEy<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">-<\/span><span class=\"n\">x1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"o\">-<\/span><span class=\"n\">y1<\/span><span class=\"p\">))<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"nf\">log<\/span><span class=\"p\">(<\/span><span class=\"nf\">E<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">-<\/span><span class=\"n\">x1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"o\">-<\/span><span class=\"n\">y1<\/span><span class=\"p\">)))<\/span> \\\r\n     <span class=\"nb\">w<\/span> <span class=\"n\">vec<\/span> <span class=\"nb\">notitle<\/span>  <span class=\"n\">lw<\/span> <span class=\"mi\">4<\/span> <span class=\"n\">lc<\/span> <span class=\"n\">palette<\/span> <span class=\"n\">filled<\/span> <span class=\"n\">head<\/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 wp-image-3631\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/vec-fig02a.svg\" alt=\"\" width=\"480\" height=\"480\" \/><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h4 id=\"\u8ca0\u96fb\u8377\u304c\u3064\u304f\u308b\u96fb\u5834\u306e\u5411\u304d\">\u8ca0\u96fb\u8377\u304c\u3064\u304f\u308b\u96fb\u5834\u306e\u5411\u304d<\/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 [3]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-gnuplot\">\n<pre><span class=\"c\"># \u96fb\u5834 E<\/span>\r\n<span class=\"nf\">Ex<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"o\">-<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">-<\/span><span class=\"n\">x2<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"p\">(<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">((<\/span><span class=\"n\">x<\/span><span class=\"o\">-<\/span><span class=\"n\">x2<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"p\">(<\/span><span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y2<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">3<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nf\">Ey<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"o\">-<\/span><span class=\"p\">(<\/span><span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y2<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"p\">(<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">((<\/span><span class=\"n\">x<\/span><span class=\"o\">-<\/span><span class=\"n\">x2<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"p\">(<\/span><span class=\"n\">y<\/span><span class=\"o\">-<\/span><span class=\"n\">y2<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">3<\/span><span class=\"p\">)<\/span>\r\n<span class=\"c\"># \u96fb\u5834 E \u306e\u5927\u304d\u3055<\/span>\r\n<span class=\"nf\">E<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"nf\">Ex<\/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\">Ey<\/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 E<\/span>\r\n<span class=\"nf\">hEx<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"n\">scaling<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">Ex<\/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\">E<\/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\">hEy<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"n\">scaling<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">Ey<\/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\">E<\/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=\"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=\"c\"># special filename \"++\" \u306e\u3068\u304d\u306f<\/span>\r\n<span class=\"c\"># urange, vrange \u3067\u63cf\u753b\u7bc4\u56f2\u306e\u8a2d\u5b9a<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">urange<\/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\">vrange<\/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=\"c\"># x \u8ef8\u65b9\u5411\u306e\u30b5\u30f3\u30d7\u30eb\u6570<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">samples<\/span> <span class=\"mi\">19<\/span>\r\n<span class=\"c\"># y \u8ef8\u65b9\u5411\u306e\u30b5\u30f3\u30d7\u30eb\u6570<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">isosamples<\/span> <span class=\"mi\">19<\/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\r\n<span class=\"c\"># \u8ca0\u96fb\u8377\u306e\u4f4d\u7f6e\u306b\u8d64\u4e38\u3068 -<\/span>\r\n<span class=\"nv\">x2<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"nv\">y2<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">-1<\/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=\"mf\">1.5<\/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\">\"-\"<\/span> <span class=\"n\">tc<\/span> <span class=\"s\">\"red\"<\/span>\r\n\r\n<span class=\"k\">set<\/span> <span class=\"nb\">title<\/span> <span class=\"s\">\"\u8ca0\u96fb\u8377\u304c\u3064\u304f\u308b\u96fb\u5834\u306e\u5411\u304d\"<\/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\">1<\/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=\"k\">plot<\/span> <span class=\"s\">\"++\"<\/span> <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\">hEx<\/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\">hEy<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">))<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"nf\">log<\/span><span class=\"p\">(<\/span><span class=\"nf\">E<\/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=\"nb\">notitle<\/span> <span class=\"n\">lc<\/span> <span class=\"n\">palette<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">4<\/span> <span class=\"n\">filled<\/span> <span class=\"n\">head<\/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 wp-image-3632\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/vec-fig03a.svg\" alt=\"\" width=\"480\" height=\"480\" \/><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h4 id=\"\u6b63\u96fb\u8377\u3068\u8ca0\u96fb\u8377\u304c\u3064\u304f\u308b\u96fb\u5834\u306e\u5411\u304d\">\u6b63\u96fb\u8377\u3068\u8ca0\u96fb\u8377\u304c\u3064\u304f\u308b\u96fb\u5834\u306e\u5411\u304d<\/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 [4]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-gnuplot\">\n<pre><span class=\"c\"># \u96fb\u5834 E<\/span>\r\n<span class=\"nf\">Ex<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"n\">x<\/span><span class=\"o\">\/<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"n\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">3<\/span>\r\n<span class=\"nf\">Ey<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"n\">y<\/span><span class=\"o\">\/<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"n\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">3<\/span>\r\n\r\n<span class=\"c\"># \u6b63\u96fb\u8377\u3068\u8ca0\u96fb\u8377\u304c\u3064\u304f\u308b\u96fb\u5834<\/span>\r\n<span class=\"nf\">E2x<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">Ex<\/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\">Ex<\/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\">E2y<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">Ey<\/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\">Ey<\/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=\"c\"># \u96fb\u5834\u306e\u5927\u304d\u3055<\/span>\r\n<span class=\"nf\">E<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"nf\">E2x<\/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\">E2y<\/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 E<\/span>\r\n<span class=\"nf\">hE2x<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"n\">scaling<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">E2x<\/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\">E<\/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\">hE2y<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"n\">scaling<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">E2y<\/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\">E<\/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=\"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=\"c\"># special filename \"++\" \u306e\u3068\u304d\u306f<\/span>\r\n<span class=\"c\"># urange, vrange \u3067\u63cf\u753b\u7bc4\u56f2\u306e\u8a2d\u5b9a<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">urange<\/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\">vrange<\/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=\"c\"># x \u8ef8\u65b9\u5411\u306e\u30b5\u30f3\u30d7\u30eb\u6570<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">samples<\/span> <span class=\"mi\">19<\/span>\r\n<span class=\"c\"># y \u8ef8\u65b9\u5411\u306e\u30b5\u30f3\u30d7\u30eb\u6570<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">isosamples<\/span> <span class=\"mi\">19<\/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\r\n<span class=\"c\"># \u6b63\u96fb\u8377\u306e\u4f4d\u7f6e\u306b\u9752\u4e38\u3068 +<\/span>\r\n<span class=\"nv\">x1<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"nv\">y1<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">1<\/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=\"mf\">1.5<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">\"blue\"<\/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\">\"+\"<\/span> <span class=\"n\">tc<\/span> <span class=\"s\">\"blue\"<\/span>\r\n\r\n<span class=\"c\"># \u8ca0\u96fb\u8377\u306e\u4f4d\u7f6e\u306b\u8d64\u4e38\u3068 -<\/span>\r\n<span class=\"nv\">x2<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"nv\">y2<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">-1<\/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=\"mf\">1.5<\/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\">\"-\"<\/span> <span class=\"n\">tc<\/span> <span class=\"s\">\"red\"<\/span>\r\n\r\n<span class=\"k\">set<\/span> <span class=\"nb\">title<\/span> <span class=\"s\">\"\u6b63\u96fb\u8377\u3068\u8ca0\u96fb\u8377\u304c\u3064\u304f\u308b\u96fb\u5834\u306e\u5411\u304d\"<\/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\">plot<\/span> <span class=\"s\">\"++\"<\/span> <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\">hE2x<\/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\">hE2y<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">))<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"nf\">log<\/span><span class=\"p\">(<\/span><span class=\"nf\">E<\/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=\"nb\">notitle<\/span> <span class=\"n\">lc<\/span> <span class=\"n\">palette<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">4<\/span> <span class=\"n\">filled<\/span> <span class=\"n\">head<\/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 wp-image-3633\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/vec-fig04a.svg\" alt=\"\" width=\"480\" height=\"480\" \/><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h3 id=\"\u96fb\u6c17\u53cc\u6975\u5b50\u304c\u3064\u304f\u308b\u96fb\u5834\">\u96fb\u6c17\u53cc\u6975\u5b50\u304c\u3064\u304f\u308b\u96fb\u5834<\/h3>\n<p>$$\\boldsymbol{E} = 3 \\frac{\\boldsymbol{r}\\cdot\\boldsymbol{p}}{r^5} \\boldsymbol{r}<br \/>\n&#8211; \\frac{\\boldsymbol{p}}{r^3}$$<\/p>\n<p>$\\boldsymbol{p} = (0, 1, 0), \\quad z = 0$ \u3068\u3057\u3066<\/p>\n<p>\\begin{eqnarray}<br \/>\nE_x &amp;=&amp; \\frac{3}{\\left( x^2 + y^2 \\right)^{5\/2}} x y \\\\<br \/>\nE_y &amp;=&amp; \\frac{3}{\\left( x^2 + y^2 \\right)^{5\/2}} y^2 &#8211; \\frac{1}{\\left( x^2 + y^2 \\right)^{3\/2}} \\\\<br \/>\n\\end{eqnarray}<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h4 id=\"\u96fb\u6c17\u53cc\u6975\u5b50\u304c\u3064\u304f\u308b\u96fb\u5834\u306e\u5411\u304d\">\u96fb\u6c17\u53cc\u6975\u5b50\u304c\u3064\u304f\u308b\u96fb\u5834\u306e\u5411\u304d<\/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 [5]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-gnuplot\">\n<pre><span class=\"c\"># \u96fb\u5834 E<\/span>\r\n<span class=\"nf\">Ex<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"mf\">3.<\/span><span class=\"o\">*<\/span> <span class=\"n\">x<\/span><span class=\"o\">*<\/span><span class=\"n\">y<\/span><span class=\"o\">\/<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"n\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">5<\/span> \r\n<span class=\"nf\">Ey<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"mf\">3.<\/span><span class=\"o\">*<\/span> <span class=\"n\">y<\/span><span class=\"o\">*<\/span><span class=\"n\">y<\/span><span class=\"o\">\/<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"n\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">5<\/span> <span class=\"o\">-<\/span> <span class=\"mf\">1.<\/span><span class=\"o\">\/<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"n\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">3<\/span>\r\n\r\n<span class=\"c\"># \u898f\u683c\u5316\u3055\u308c\u305f hat E<\/span>\r\n<span class=\"nf\">E<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"nf\">Ex<\/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\">Ey<\/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=\"nf\">hEx<\/span>(x,y) <span class=\"o\">=<\/span> <span class=\"n\">scaling<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">Ex<\/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\">E<\/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\">hEy<\/span>(x,y) <span class=\"o\">=<\/span> <span class=\"n\">scaling<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">Ey<\/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\">E<\/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=\"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<span class=\"k\">set<\/span> <span class=\"nb\">zeroaxis<\/span>\r\n\r\n<span class=\"c\"># special filename \"++\" \u306e\u3068\u304d\u306f<\/span>\r\n<span class=\"c\"># urange, vrange \u3067\u8a2d\u5b9a<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">urange<\/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\">vrange<\/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=\"c\"># x \u8ef8\u65b9\u5411\u306e\u30b5\u30f3\u30d7\u30eb\u6570<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">samples<\/span> <span class=\"mi\">19<\/span>\r\n<span class=\"c\"># y \u8ef8\u65b9\u5411\u306e\u30b5\u30f3\u30d7\u30eb\u6570<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">isosamples<\/span> <span class=\"mi\">19<\/span>\r\n\r\n<span class=\"nv\">scaling<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.15<\/span>\r\n\r\n<span class=\"c\"># \u96fb\u6c17\u53cc\u6975\u5b50 p<\/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=\"mi\">0<\/span><span class=\"o\">,<\/span> <span class=\"mi\">-0<\/span><span class=\"mf\">.2<\/span> <span class=\"n\">to<\/span> <span class=\"mi\">0<\/span><span class=\"o\">,<\/span> <span class=\"mf\">0.2<\/span> <span class=\"n\">filled<\/span> <span class=\"n\">head<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">5<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">\"red\"<\/span> <span class=\"n\">lt<\/span> <span class=\"mi\">1<\/span>\r\n\r\n<span class=\"k\">set<\/span> <span class=\"nb\">title<\/span> <span class=\"s\">\"\u96fb\u6c17\u53cc\u6975\u5b50\u304c\u3064\u304f\u308b\u96fb\u5834\u306e\u5411\u304d\"<\/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\">0<\/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=\"k\">plot<\/span> <span class=\"s\">\"++\"<\/span> <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\">hEx<\/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\">hEy<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">))<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"nf\">log<\/span><span class=\"p\">(<\/span><span class=\"nf\">E<\/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=\"nb\">notitle<\/span> <span class=\"n\">lc<\/span> <span class=\"n\">palette<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">4<\/span> <span class=\"n\">filled<\/span> <span class=\"n\">head<\/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 wp-image-3634\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/vec-fig05a.svg\" alt=\"\" width=\"480\" height=\"480\" \/><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h3 id=\"\u76f4\u7dda\u96fb\u6d41\u304c\u3064\u304f\u308b\u78c1\u5834-$\\boldsymbol{B}$\">\u76f4\u7dda\u96fb\u6d41\u304c\u3064\u304f\u308b\u78c1\u5834 $\\boldsymbol{B}$<\/h3>\n<p>$$\\boldsymbol{B} = \\frac{\\boldsymbol{I}\\times\\boldsymbol{\\rho}}{\\rho^2}$$$$ \\boldsymbol{I} = (0, 0, I), \\quad \\boldsymbol{\\rho} = (x, y, 0)$$\\begin{eqnarray}<br \/>\nB_x &amp;=&amp; &#8211; \\frac{ I y}{x^2 + y^2} \\\\<br \/>\nB_y &amp;=&amp; \\frac{ I x}{x^2 + y^2}<br \/>\n\\end{eqnarray}<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h4 id=\"\u76f4\u7dda\u96fb\u6d41\u304c\u3064\u304f\u308b\u78c1\u5834\">\u76f4\u7dda\u96fb\u6d41\u304c\u3064\u304f\u308b\u78c1\u5834<\/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 [6]:<\/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\">scaling<\/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\">scaling<\/span> <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\"># \u78c1\u5834 B \u306e\u5927\u304d\u3055<\/span>\r\n<span class=\"nf\">B<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"mf\">1.<\/span><span class=\"o\">\/<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"n\">x<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"n\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/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<span class=\"k\">set<\/span> <span class=\"nb\">zeroaxis<\/span>\r\n\r\n<span class=\"c\"># special filename \"++\" \u306e\u3068\u304d\u306f<\/span>\r\n<span class=\"c\"># urange, vrange \u3067\u8a2d\u5b9a<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">urange<\/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\">vrange<\/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=\"c\"># x \u8ef8\u65b9\u5411\u306e\u30b5\u30f3\u30d7\u30eb\u6570<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">samples<\/span> <span class=\"mi\">19<\/span>\r\n<span class=\"c\"># y \u8ef8\u65b9\u5411\u306e\u30b5\u30f3\u30d7\u30eb\u6570<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">isosamples<\/span> <span class=\"mi\">19<\/span>\r\n\r\n<span class=\"nv\">scaling<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.25<\/span>\r\n\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=\"mf\">1.5<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">\"blue\"<\/span> <span class=\"n\">at<\/span> <span class=\"mi\">0<\/span><span class=\"o\">,<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">label<\/span> <span class=\"mi\">2<\/span> <span class=\"n\">center<\/span> <span class=\"n\">at<\/span> <span class=\"n\">first<\/span> <span class=\"mi\">0<\/span><span class=\"o\">,<\/span><span class=\"mi\">0<\/span> <span class=\"s\">\"+\"<\/span> <span class=\"n\">tc<\/span> <span class=\"s\">\"blue\"<\/span>\r\n\r\n<span class=\"k\">set<\/span> <span class=\"nb\">title<\/span> <span class=\"s\">\"\u76f4\u7dda\u96fb\u6d41\u304c\u3064\u304f\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\">1<\/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=\"k\">plot<\/span> <span class=\"s\">\"++\"<\/span> <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\">Bx<\/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\">By<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">))<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"nf\">log<\/span><span class=\"p\">(<\/span><span class=\"nf\">B<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">)))<\/span>\\\r\n     <span class=\"nb\">w<\/span> <span class=\"n\">vec<\/span> <span class=\"nb\">notitle<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">4<\/span> <span class=\"n\">lc<\/span> <span class=\"n\">palette<\/span> <span class=\"n\">filled<\/span> <span class=\"n\">head<\/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 wp-image-3635\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/vec-fig06a.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\"><\/div>\n<\/div>\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<h4 id=\"\u76f4\u7dda\u96fb\u6d41\u304c\u3064\u304f\u308b\u78c1\u5834\u306e\u5411\u304d\">\u76f4\u7dda\u96fb\u6d41\u304c\u3064\u304f\u308b\u78c1\u5834\u306e\u5411\u304d<\/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 [7]:<\/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\"># \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\">Bx<\/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\">By<\/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\">hBx<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"n\">scaling<\/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\">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\">hBy<\/span>(x, y) <span class=\"o\">=<\/span> <span class=\"n\">scaling<\/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\">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=\"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<span class=\"k\">set<\/span> <span class=\"nb\">zeroaxis<\/span>\r\n<span class=\"c\"># \u8868\u793a\u7bc4\u56f2\u306e\u8a2d\u5b9a<\/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=\"c\"># special filename \"++\" \u306e\u3068\u304d\u306f<\/span>\r\n<span class=\"c\"># urange, vrange \u3067\u8a2d\u5b9a<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">urange<\/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\">vrange<\/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\r\n<span class=\"k\">set<\/span> <span class=\"nb\">samples<\/span> <span class=\"mi\">19<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">isosamples<\/span> <span class=\"mi\">19<\/span>\r\n\r\n<span class=\"nv\">scaling<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.2<\/span>\r\n\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=\"mf\">1.5<\/span> <span class=\"n\">lc<\/span> <span class=\"s\">\"blue\"<\/span> <span class=\"n\">at<\/span> <span class=\"mi\">0<\/span><span class=\"o\">,<\/span> <span class=\"mi\">0<\/span>\r\n<span class=\"k\">set<\/span> <span class=\"nb\">label<\/span> <span class=\"mi\">2<\/span> <span class=\"n\">center<\/span> <span class=\"n\">at<\/span> <span class=\"n\">first<\/span> <span class=\"mi\">0<\/span><span class=\"o\">,<\/span><span class=\"mi\">0<\/span> <span class=\"s\">\"+\"<\/span> <span class=\"n\">tc<\/span> <span class=\"s\">\"blue\"<\/span>\r\n\r\n<span class=\"k\">set<\/span> <span class=\"nb\">title<\/span> <span class=\"s\">\"\u76f4\u7dda\u96fb\u6d41\u304c\u3064\u304f\u308b\u78c1\u5834\u306e\u5411\u304d\"<\/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\">1<\/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=\"k\">plot<\/span> <span class=\"s\">\"++\"<\/span> <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\">hBx<\/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\">hBy<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">))<\/span><span class=\"o\">:<\/span><span class=\"p\">(<\/span><span class=\"nf\">log<\/span><span class=\"p\">(<\/span><span class=\"nf\">B<\/span><span class=\"p\">(<\/span><span class=\"err\">$<\/span><span class=\"mi\">1<\/span><span class=\"o\">,<\/span><span class=\"err\">$<\/span><span class=\"mi\">2<\/span><span class=\"p\">)))<\/span> \\\r\n     <span class=\"nb\">w<\/span> <span class=\"n\">vec<\/span> <span class=\"nb\">notitle<\/span> <span class=\"n\">lw<\/span> <span class=\"mi\">4<\/span> <span class=\"n\">lc<\/span> <span class=\"n\">palette<\/span> <span class=\"n\">filled<\/span> <span class=\"n\">head<\/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 wp-image-3636\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/vec-fig07a.svg\" alt=\"\" width=\"480\" height=\"480\" \/><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"excerpt":{"rendered":"<p>\u5148\u65e5\uff0c\u300cgnuplot \u3067\u30d9\u30af\u30c8\u30eb\u5834\uff08\u96fb\u5834\u30fb\u78c1\u5834\uff09\u3092\u63cf\u304f\u300d\u3067\u63cf\u3044\u3066\u307f\u305f\u304c\uff0c\u3044\u308d\u3044\u308d\u691c\u7d22\u3059\u308b\u3068 palette \u3067\u30d9\u30af\u30c8\u30eb\u306e\u8272\u3092\uff08\u30d9\u30af\u30c8\u30eb\u306e\u5927\u304d\u3055\u306b\u5408\u308f\u305b\u3066\uff09\u53ef\u5909\u306b\u3067\u304d\u308b\u3088\u3046\u306a\u306e\u3067\uff0c\u5099\u5fd8\u9332\u3068\u3057\u3066\u3002<\/p><p><a class=\"more-link btn\" href=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/3619\/\">\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-3619","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\/3619","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=3619"}],"version-history":[{"count":10,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/posts\/3619\/revisions"}],"predecessor-version":[{"id":3675,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/posts\/3619\/revisions\/3675"}],"wp:attachment":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/media?parent=3619"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/categories?post=3619"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/tags?post=3619"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}