{"id":7371,"date":"2025-01-31T12:00:08","date_gmt":"2025-01-31T03:00:08","guid":{"rendered":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/?page_id=7371"},"modified":"2025-01-31T12:55:40","modified_gmt":"2025-01-31T03:55:40","slug":"sympy-%e3%81%a7%e9%ab%98%e3%81%95-h-%e3%81%8b%e3%82%89%e3%81%ae%e6%96%9c%e6%96%b9%e6%8a%95%e5%b0%84%e3%81%ae%e6%9c%80%e5%a4%a7%e5%88%b0%e9%81%94%e8%b7%9d%e9%9b%a2%e3%82%92%e6%b1%82%e3%82%81%e3%82%8b","status":"publish","type":"page","link":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/%e3%82%b3%e3%83%b3%e3%83%94%e3%83%a5%e3%83%bc%e3%82%bf%e6%bc%94%e7%bf%92\/%e9%ab%98%e3%81%95-h-%e3%81%8b%e3%82%89%e3%81%ae%e6%96%9c%e6%96%b9%e6%8a%95%e5%b0%84%e3%81%ae%e6%9c%80%e5%a4%a7%e5%88%b0%e9%81%94%e8%b7%9d%e9%9b%a2%e3%82%92%e6%b1%82%e3%82%81%e3%82%8b%e6%ba%96\/sympy-%e3%81%a7%e9%ab%98%e3%81%95-h-%e3%81%8b%e3%82%89%e3%81%ae%e6%96%9c%e6%96%b9%e6%8a%95%e5%b0%84%e3%81%ae%e6%9c%80%e5%a4%a7%e5%88%b0%e9%81%94%e8%b7%9d%e9%9b%a2%e3%82%92%e6%b1%82%e3%82%81%e3%82%8b\/","title":{"rendered":"SymPy \u3066\u3099\u9ad8\u3055 h \u304b\u3089\u306e\u659c\u65b9\u6295\u5c04\u306e\u6700\u5927\u5230\u9054\u8ddd\u96e2\u3092\u6c42\u3081\u308b"},"content":{"rendered":"<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<p>\u5730\u9762\u304b\u3089\u9ad8\u3055 $h$ \u306e\u5730\u70b9\u304b\u3089\u7a7a\u6c17\u62b5\u6297\u306a\u3057\u306e\u659c\u65b9\u6295\u5c04\u3092\u884c\u3046\u3068\uff0c\u521d\u901f\u5ea6\u306e\u5927\u304d\u3055\u3092\u4e00\u5b9a\u3068\u3057\u305f\u5834\u5408\uff0c\u6c34\u5e73\u65b9\u5411\u306e\u5230\u9054\u8ddd\u96e2\u304c\u6700\u5927\u3068\u306a\u308b\u306e\u306f\u6253\u3061\u4e0a\u3052\u89d2\u5ea6\uff08\u4ef0\u89d2\uff09\u304c\u4f55\u5ea6\u306e\u3068\u304d\u304b\u3092 Python \u306e SymPy \u3092\u4f7f\u3063\u3066\uff08\u304b\u3064 SciPy \u3084 NumPy \u3092\u4f7f\u308f\u305a\u306b\uff09\u6c42\u3081\u308b\u3002<\/p>\n<p style=\"text-align: right;\">update: 2025.01.31<\/p>\n<\/div>\n<\/div>\n<\/div>\n<p><!--more--><\/p>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h3 id=\"\u30e9\u30a4\u30d6\u30e9\u30ea\u306e-import\">\u30e9\u30a4\u30d6\u30e9\u30ea\u306e import<\/h3>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[1]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"kn\">from<\/span> <span class=\"nn\">sympy.abc<\/span> <span class=\"kn\">import<\/span> <span class=\"o\">*<\/span>\r\n<span class=\"kn\">from<\/span> <span class=\"nn\">sympy<\/span> <span class=\"kn\">import<\/span> <span class=\"o\">*<\/span> \r\n\r\n<span class=\"c1\"># SymPy Plotting Backends (SPB) \u3067\u30b0\u30e9\u30d5\u3092\u63cf\u304f<\/span>\r\n<span class=\"kn\">from<\/span> <span class=\"nn\">spb<\/span> <span class=\"kn\">import<\/span> <span class=\"o\">*<\/span>\r\n\r\n<span class=\"c1\"># \u30b0\u30e9\u30d5\u3092 SVG \u3067 Notebook \u306b\u30a4\u30f3\u30e9\u30a4\u30f3\u8868\u793a<\/span>\r\n<span class=\"o\">%<\/span><span class=\"k\">config<\/span> InlineBackend.figure_formats = ['svg']\r\n\r\n<span class=\"n\">init_printing<\/span><span class=\"p\">()<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h3 id=\"\u904b\u52d5\u65b9\u7a0b\u5f0f\u306e\u89e3\">\u904b\u52d5\u65b9\u7a0b\u5f0f\u306e\u89e3<\/h3>\n<p>\u521d\u671f\u6761\u4ef6\u3092 $t = 0$ \u3067<\/p>\n<p>$$x(0) = 0, \\quad y(0) = h, \\quad v_x(0)\u00a0 = v_0 \\cos\\theta, \\quad v_y(0) = v_0 \\sin \\theta$$<\/p>\n<p>\u3068\u3057\u305f\u3068\u304d\u306e\u89e3\u306f<\/p>\n<p>\\begin{eqnarray}<br \/>\nx(t) &amp;=&amp; v_0 \\cos\\theta\\cdot t \\\\<br \/>\ny(t) &amp;=&amp; h + v_0 \\sin\\theta\\cdot t -\\frac{1}{2} g t^2<br \/>\n\\end{eqnarray}<\/p>\n<h3 id=\"\u7121\u6b21\u5143\u5316\u3055\u308c\u305f\u659c\u65b9\u6295\u5c04\u306e\u89e3\">\u7121\u6b21\u5143\u5316\u3055\u308c\u305f\u659c\u65b9\u6295\u5c04\u306e\u89e3<\/h3>\n<p>\u3053\u306e\u7cfb\u306b\u7279\u5fb4\u7684\u306a\u6642\u9593 $\\displaystyle \\tau \\equiv \\frac{v_0}{g}$ \u304a\u3088\u3073\u9577\u3055 $\\displaystyle \\ell \\equiv v_0 \\tau = \\frac{v_0^2}{g}$ \u3067\u89e3\u3092\u7121\u6b21\u5143\u5316\u3059\u308b\u3068\uff0c<\/p>\n<p>\\begin{eqnarray}<br \/>\nT &amp;\\equiv&amp; \\frac{t}{\\tau} \\\\<br \/>\nH &amp;\\equiv&amp; \\frac{h}{\\ell} \\\\<br \/>\nX &amp;\\equiv&amp; \\frac{x}{\\ell} = \\cos\\theta\\cdot T \\\\<br \/>\nY&amp;\\equiv&amp; \\frac{y}{\\ell} = H+ \\sin\\theta\\cdot T -\\frac{1}{2} T^2 \\\\<br \/>\n\\end{eqnarray}<\/p>\n<p>\u8a73\u7d30\u306f\u4ee5\u4e0b\u306e\u30da\u30fc\u30b8<\/p>\n<ul>\n<li><a href=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/%e3%82%b3%e3%83%b3%e3%83%94%e3%83%a5%e3%83%bc%e3%82%bf%e6%bc%94%e7%bf%92\/%e9%ab%98%e3%81%95-h-%e3%81%8b%e3%82%89%e3%81%ae%e6%96%9c%e6%96%b9%e6%8a%95%e5%b0%84%e3%81%ae%e6%9c%80%e5%a4%a7%e5%88%b0%e9%81%94%e8%b7%9d%e9%9b%a2%e3%82%92%e6%b1%82%e3%82%81%e3%82%8b%e6%ba%96\/\">\u9ad8\u3055 h \u304b\u3089\u306e\u659c\u65b9\u6295\u5c04\u306e\u6700\u5927\u5230\u9054\u8ddd\u96e2\u3092\u6c42\u3081\u308b\u6e96\u5099<\/a><\/li>\n<\/ul>\n<p>\u7121\u6b21\u5143\u5316\u3055\u308c\u305f\u89e3 $X(T, \\theta), Y(T, \\theta)$ \u3092\u95a2\u6570\u3068\u3057\u3066\u5b9a\u7fa9\u3057\u307e\u3059\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[2]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"k\">def<\/span> <span class=\"nf\">X<\/span><span class=\"p\">(<\/span><span class=\"n\">T<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"k\">return<\/span> <span class=\"n\">cos<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">)<\/span> <span class=\"o\">*<\/span> <span class=\"n\">T<\/span>\r\n\r\n<span class=\"n\">var<\/span><span class=\"p\">(<\/span><span class=\"s1\">'H'<\/span><span class=\"p\">,<\/span> <span class=\"n\">positive<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">True<\/span><span class=\"p\">)<\/span>\r\n<span class=\"k\">def<\/span> <span class=\"nf\">Y<\/span><span class=\"p\">(<\/span><span class=\"n\">T<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"k\">return<\/span> <span class=\"n\">H<\/span> <span class=\"o\">+<\/span> <span class=\"n\">sin<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"n\">T<\/span> <span class=\"o\">-<\/span> <span class=\"n\">T<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"o\">\/<\/span><span class=\"mi\">2<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h3 id=\"\u7121\u6b21\u5143\u5316\u3055\u308c\u305f\u6ede\u7a7a\u6642\u9593\u3068\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\">\u7121\u6b21\u5143\u5316\u3055\u308c\u305f\u6ede\u7a7a\u6642\u9593\u3068\u6c34\u5e73\u5230\u9054\u8ddd\u96e2<\/h3>\n<p>\u6ede\u7a7a\u6642\u9593 $T_1 \\ (&gt;0)$ \u306f\u4ee5\u4e0b\u306e\u5f0f\u3092\u89e3\u3044\u3066\u6c42\u3081\u3089\u308c\u308b\u3002<\/p>\n<p>$$Y(T_1, \\theta, H) = 0, \\quad\\Rightarrow\\quad T_1 = T_1(\\theta, H)$$<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[3]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">var<\/span><span class=\"p\">(<\/span><span class=\"s1\">'T1'<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">solve<\/span><span class=\"p\">(<\/span><span class=\"n\">Y<\/span><span class=\"p\">(<\/span><span class=\"n\">T1<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">),<\/span> <span class=\"n\">T1<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[3]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\left[ -\\sqrt{2 H + \\sin^{2}{\\left(\\theta \\right)}} + \\sin{\\left(\\theta \\right)}, \\ \\sqrt{2 H + \\sin^{2}{\\left(\\theta \\right)}} + \\sin{\\left(\\theta \\right)}\\right]$<\/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<p>$T_1 &gt; 0$ \u3067\u3042\u308b\u304b\u3089\uff0c\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u95a2\u6570\u3068\u3057\u3066\u5b9a\u7fa9\u3057\u307e\u3059\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[4]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"k\">def<\/span> <span class=\"nf\">T1<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"k\">return<\/span> <span class=\"n\">sin<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"n\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">H<\/span> <span class=\"o\">+<\/span> <span class=\"n\">sin<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">T1<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[4]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\sqrt{2 H + \\sin^{2}{\\left(\\theta \\right)}} + \\sin{\\left(\\theta \\right)}$<\/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<p>\u6ede\u7a7a\u6642\u9593 $T_1(\\theta)$ \u306e\u9593\u306b\u6c34\u5e73\u65b9\u5411\u306b\u9032\u3080\u8ddd\u96e2\u3092\u6c34\u5e73\u5230\u9054\u8ddd\u96e2 $L(\\theta)$ \u306f<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[5]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"k\">def<\/span> <span class=\"nf\">L<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"k\">return<\/span> <span class=\"n\">X<\/span><span class=\"p\">(<\/span><span class=\"n\">T1<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">),<\/span> <span class=\"n\">theta<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">L<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[5]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\left(\\sqrt{2 H + \\sin^{2}{\\left(\\theta \\right)}} + \\sin{\\left(\\theta \\right)}\\right) \\cos{\\left(\\theta \\right)}$<\/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=\"$H-=-0$-\u306e\u5834\u5408\u306e\u8ecc\u9053\">$H = 0$ \u306e\u5834\u5408\u306e\u8ecc\u9053<\/h3>\n<p>\u307e\u305a\u306f\uff0c\u3068\u308a\u3042\u3048\u305a\u306e\u30b0\u30e9\u30d5\u3002$\\theta = 45^{\\circ}$ \u3067\u306e\u659c\u65b9\u6295\u5c04\u306e\u8ecc\u9053\u3092\u30b0\u30e9\u30d5\u306b\u3057\u3066\u307f\u308b\u3002SymPy \u3067\u5a92\u4ecb\u5909\u6570\u8868\u793a\u306e\u66f2\u7dda\u306e\u30b0\u30e9\u30d5\u3092\u63cf\u304f\u65b9\u6cd5\u306f\u6388\u696d\u3067\u3084\u308a\u307e\u3057\u305f\u3002\u4ee5\u4e0b\u3092\u53c2\u8003\u306b\u3002<\/p>\n<ul>\n<li><a href=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/%e3%82%b3%e3%83%b3%e3%83%94%e3%83%a5%e3%83%bc%e3%82%bf%e6%bc%94%e7%bf%92\/python-%e3%81%a7%e3%82%b0%e3%83%a9%e3%83%95%e4%bd%9c%e6%88%90\/sympy-plotting-backends-%e3%81%a7%e3%82%b0%e3%83%a9%e3%83%95%e4%bd%9c%e6%88%90%ef%bc%9athe-graphic-module-%e7%b7%a8\/\">SymPy Plotting Backends \u3067\u30b0\u30e9\u30d5\u4f5c\u6210\uff1aThe Graphic Module \u7de8<\/a>\n<ul>\n<li><a href=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/%e3%82%b3%e3%83%b3%e3%83%94%e3%83%a5%e3%83%bc%e3%82%bf%e6%bc%94%e7%bf%92\/python-%e3%81%a7%e3%82%b0%e3%83%a9%e3%83%95%e4%bd%9c%e6%88%90\/sympy-plotting-backends-%e3%81%a7%e3%82%b0%e3%83%a9%e3%83%95%e4%bd%9c%e6%88%90%ef%bc%9athe-graphic-module-%e7%b7%a8\/#_line_parametric_2d\">\u5a92\u4ecb\u5909\u6570\u8868\u793a\u306e\u30b0\u30e9\u30d5 line_parametric_2d()<\/a><\/li>\n<\/ul>\n<\/li>\n<\/ul>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[6]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"c1\"># \u03b8 = 45\u00b0 \u306e\u30b0\u30e9\u30d5<\/span>\r\n<span class=\"n\">th<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">45<\/span>\r\n<span class=\"n\">theta<\/span> <span class=\"o\">=<\/span> <span class=\"n\">rad<\/span><span class=\"p\">(<\/span><span class=\"n\">th<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">H<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n\r\n<span class=\"n\">graphics<\/span><span class=\"p\">(<\/span>\r\n    <span class=\"n\">line_parametric_2d<\/span><span class=\"p\">(<\/span>\r\n        <span class=\"n\">X<\/span><span class=\"p\">(<\/span><span class=\"n\">T<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">),<\/span> <span class=\"n\">Y<\/span><span class=\"p\">(<\/span><span class=\"n\">T<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">T<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"n\">T1<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">)))<\/span>\r\n<span class=\"p\">)<\/span><span class=\"p\">;<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_svg output_subarea \"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-10023\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/symtosha01.svg\" alt=\"\" width=\"640\" height=\"427\" \/><\/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<p>\u4ee5\u4e0b\u306e\u3088\u3046\u306a\u30aa\u30d7\u30b7\u30e7\u30f3\u3092\u8a2d\u5b9a\u3057\u3066\uff0c\u3082\u3046\u5c11\u3057\u4f53\u88c1\u3092\u6574\u3048\u307e\u3059\u3002<\/p>\n<ul>\n<li>$\\theta$ \u306e\u5024\u3092\u51e1\u4f8b\u306b<\/li>\n<li>\u30ab\u30e9\u30fc\u30de\u30c3\u30d7\u306f\u4e0d\u8981 <code>use_cm = False<\/code> \uff08\u30c7\u30d5\u30a9\u30eb\u30c8\u3067 <code>False<\/code> \u306b\u3057\u3066\u307b\u3057\u3044\u3068\u3053\u308d\uff09<\/li>\n<li>\u5ea7\u6a19\u8ef8\u306e\u30e9\u30d9\u30eb\u3068\u30b0\u30e9\u30d5\u306e\u30bf\u30a4\u30c8\u30eb<\/li>\n<li>\u6a2a\u8ef8\u7e26\u8ef8\u306e\u8868\u793a\u7bc4\u56f2<\/li>\n<li>\u7e26\u8ef8\u6a2a\u8ef8\u306e\u30a2\u30b9\u30da\u30af\u30c8\u6bd4<\/li>\n<\/ul>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[7]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"c1\"># \u03b8 = 45\u00b0 \u306e\u30b0\u30e9\u30d5<\/span>\r\n<span class=\"n\">th<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">45<\/span>\r\n<span class=\"n\">theta<\/span> <span class=\"o\">=<\/span> <span class=\"n\">rad<\/span><span class=\"p\">(<\/span><span class=\"n\">th<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">H<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n\r\n<span class=\"n\">graphics<\/span><span class=\"p\">(<\/span>\r\n    <span class=\"n\">line_parametric_2d<\/span><span class=\"p\">(<\/span>\r\n        <span class=\"n\">X<\/span><span class=\"p\">(<\/span><span class=\"n\">T<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">),<\/span> <span class=\"n\">Y<\/span><span class=\"p\">(<\/span><span class=\"n\">T<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">T<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"n\">T1<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">)),<\/span> \r\n        <span class=\"c1\"># \u51e1\u4f8b\u306b \u03b8 \u306e\u5024\u3092\u5165\u308c\u308b<\/span>\r\n        <span class=\"n\">label<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'\u03b8 = <\/span><span class=\"si\">%2d<\/span><span class=\"s1\">\u00b0'<\/span> <span class=\"o\">%<\/span> <span class=\"n\">th<\/span><span class=\"p\">,<\/span>\r\n        <span class=\"n\">use_cm<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">False<\/span><span class=\"p\">),<\/span> \r\n    <span class=\"c1\"># \u30bf\u30a4\u30c8\u30eb\u306b H \u306e\u5024\u3092\u5165\u308c\u308b<\/span>\r\n    <span class=\"n\">title<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'H = <\/span><span class=\"si\">%.1f<\/span><span class=\"s1\"> \u306e\u5834\u5408\u306e\u659c\u65b9\u6295\u5c04'<\/span> <span class=\"o\">%<\/span> <span class=\"n\">H<\/span><span class=\"p\">,<\/span>\r\n    <span class=\"n\">xlabel<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'X'<\/span><span class=\"p\">,<\/span> <span class=\"n\">ylabel<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'Y'<\/span><span class=\"p\">,<\/span> \r\n    <span class=\"n\">xlim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">1.1<\/span><span class=\"p\">),<\/span> <span class=\"n\">ylim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">0.5<\/span><span class=\"p\">),<\/span> \r\n    <span class=\"n\">aspect<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'equal'<\/span><span class=\"p\">,<\/span> \r\n    <span class=\"n\">legend<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">True<\/span>\r\n<span class=\"p\">)<\/span><span class=\"p\">;<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_svg output_subarea \"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-10024\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/symtosha02.svg\" alt=\"\" width=\"640\" height=\"427\" \/><\/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<p>$\\theta$ \u306e\u5024\u3092 $43^{\\circ}$ \u304b\u3089 $47^{\\circ}$ \u307e\u3067 $1^{\\circ}$ \u523b\u307f\u3067\u5909\u3048\u3066\uff0c\u8907\u6570\u306e\u66f2\u7dda\u3092\u4e00\u5ea6\u306b\u63cf\u304f\u4f8b\u3002<\/p>\n<p>\u306a\u304a\uff0c\u7d30\u304b\u3044\u3053\u3068\u3067\u3059\u304c\uff0c\u30b0\u30e9\u30d5\u306e\u6700\u9ad8\u70b9\u306f $\\theta$ \u306e\u5024\u304c\u5927\u304d\u3044\u3068\u304d\u3067\u3059\u306e\u3067\uff0c\u5148\u306b $47^{\\circ}$ \u304b\u3089\u63cf\u304d\u306f\u3058\u3081\u3066\uff0c$1^{\\circ}$ \u305a\u3064\u6e1b\u3089\u3057\u3066\u307f\u307e\u3059\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[8]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"c1\"># th \u3092\u5909\u3048\u3066 line_parametric_2d() \u3092<\/span>\r\n<span class=\"c1\"># \u4f55\u884c\u3082\u30b3\u30d4\u30da\u3059\u308b\u306e\u306f\u5927\u5909\u306a\u306e\u3067\uff0c<\/span>\r\n<span class=\"c1\"># \u95a2\u6570\u3092\u5b9a\u7fa9\u3057\u3066\u307f\u308b<\/span>\r\n\r\n<span class=\"k\">def<\/span> <span class=\"nf\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"n\">th<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"n\">theta<\/span> <span class=\"o\">=<\/span> <span class=\"n\">rad<\/span><span class=\"p\">(<\/span><span class=\"n\">th<\/span><span class=\"p\">)<\/span>\r\n    <span class=\"k\">return<\/span> <span class=\"n\">line_parametric_2d<\/span><span class=\"p\">(<\/span>\r\n        <span class=\"n\">X<\/span><span class=\"p\">(<\/span><span class=\"n\">T<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">),<\/span> <span class=\"n\">Y<\/span><span class=\"p\">(<\/span><span class=\"n\">T<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">T<\/span><span class=\"p\">,<\/span> <span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"n\">T1<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">)),<\/span> \r\n        <span class=\"c1\"># \u51e1\u4f8b\u306b \u03b8 \u306e\u5024\u3092\u5165\u308c\u308b<\/span>\r\n        <span class=\"n\">label<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'\u03b8 = <\/span><span class=\"si\">%2d<\/span><span class=\"s1\">\u00b0'<\/span> <span class=\"o\">%<\/span> <span class=\"n\">th<\/span><span class=\"p\">,<\/span>\r\n        <span class=\"n\">use_cm<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">False<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[9]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">H<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n\r\n<span class=\"n\">graphics<\/span><span class=\"p\">(<\/span>\r\n    <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">47<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">46<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">45<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">44<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">43<\/span><span class=\"p\">),<\/span>\r\n    <span class=\"c1\"># \u30bf\u30a4\u30c8\u30eb\u306b H \u306e\u5024\u3092\u5165\u308c\u308b<\/span>\r\n    <span class=\"n\">title<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'H = <\/span><span class=\"si\">%.1f<\/span><span class=\"s1\"> \u306e\u5834\u5408\u306e\u659c\u65b9\u6295\u5c04'<\/span> <span class=\"o\">%<\/span> <span class=\"n\">H<\/span><span class=\"p\">,<\/span>\r\n    <span class=\"n\">xlabel<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'X'<\/span><span class=\"p\">,<\/span> <span class=\"n\">ylabel<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'Y'<\/span><span class=\"p\">,<\/span> \r\n    <span class=\"n\">xlim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">1.1<\/span><span class=\"p\">),<\/span> <span class=\"n\">ylim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">0.5<\/span><span class=\"p\">),<\/span> \r\n    <span class=\"n\">aspect<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'equal'<\/span><span class=\"p\">,<\/span> \r\n    <span class=\"n\">legend<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">True<\/span>\r\n<span class=\"p\">)<\/span><span class=\"p\">;<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_svg output_subarea \"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-10025\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/symtosha03.svg\" alt=\"\" width=\"640\" height=\"427\" \/><\/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<p>\u6700\u5927\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b\u89d2\u5ea6 $\\theta$ \u3092\u8abf\u3079\u308b\u305f\u3081\u306b\uff0c\u7740\u5730\u70b9\u4ed8\u8fd1\u3092\u62e1\u5927\u8868\u793a\u3057\u307e\u3059\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[10]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">H<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>\r\n\r\n<span class=\"n\">graphics<\/span><span class=\"p\">(<\/span>\r\n    <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">47<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">46<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">45<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">44<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">43<\/span><span class=\"p\">),<\/span>\r\n    <span class=\"c1\"># \u30bf\u30a4\u30c8\u30eb\u306b H \u306e\u5024\u3092\u5165\u308c\u308b<\/span>\r\n    <span class=\"n\">title<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'H = <\/span><span class=\"si\">%.1f<\/span><span class=\"s1\"> \u306e\u5834\u5408\u306e\u659c\u65b9\u6295\u5c04'<\/span> <span class=\"o\">%<\/span> <span class=\"n\">H<\/span><span class=\"p\">,<\/span>\r\n    <span class=\"n\">xlabel<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'X'<\/span><span class=\"p\">,<\/span> <span class=\"n\">ylabel<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'Y'<\/span><span class=\"p\">,<\/span> \r\n    <span class=\"c1\"># \u7740\u5730\u70b9\u4ed8\u8fd1\u3092\u62e1\u5927\u8868\u793a<\/span>\r\n    <span class=\"n\">xlim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"mf\">0.99<\/span><span class=\"p\">,<\/span> <span class=\"mf\">1.001<\/span><span class=\"p\">),<\/span> <span class=\"n\">ylim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">0.01<\/span><span class=\"p\">),<\/span> \r\n    <span class=\"n\">legend<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">True<\/span>\r\n<span class=\"p\">)<\/span><span class=\"p\">;<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_svg output_subarea \"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-10026\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/symtosha04.svg\" alt=\"\" width=\"640\" height=\"427\" \/><\/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<p>\u4e0a\u306e\u30b0\u30e9\u30d5\u304b\u3089\uff0c$H=0$ \u306e\u5834\u5408\u306b\u306f\u78ba\u304b\u306b $\\theta=45^{\\circ}$ \u306e\u6642\u306b\u6700\u5927\u5230\u9054\u8ddd\u96e2\u306b\u306a\u3063\u3066\u3044\u308b\u3053\u3068\u304c\u308f\u304b\u308a\u307e\u3059\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h3 id=\"$H-=-0.2$-\u306e\u5834\u5408\u306e\u8ecc\u9053\">$H = 0.2$ \u306e\u5834\u5408\u306e\u8ecc\u9053<\/h3>\n<p>\u4eca\u5ea6\u306f\uff0c$H=0.2$ \u306e\u5834\u5408\u306e\u8ecc\u9053\u3092\u30b0\u30e9\u30d5\u306b\u3057\u3066\u307f\u307e\u3059\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[11]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">H<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.2<\/span>\r\n\r\n<span class=\"n\">graphics<\/span><span class=\"p\">(<\/span>\r\n    <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">47<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">46<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">45<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">44<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">43<\/span><span class=\"p\">),<\/span>\r\n    <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">42<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">41<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">40<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">39<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">38<\/span><span class=\"p\">),<\/span>\r\n    <span class=\"c1\"># \u30bf\u30a4\u30c8\u30eb\u306b H \u306e\u5024\u3092\u5165\u308c\u308b<\/span>\r\n    <span class=\"n\">title<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'H = <\/span><span class=\"si\">%.1f<\/span><span class=\"s1\"> \u306e\u5834\u5408\u306e\u659c\u65b9\u6295\u5c04'<\/span> <span class=\"o\">%<\/span> <span class=\"n\">H<\/span><span class=\"p\">,<\/span>\r\n    <span class=\"n\">xlabel<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'X'<\/span><span class=\"p\">,<\/span> <span class=\"n\">ylabel<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'Y'<\/span><span class=\"p\">,<\/span> \r\n    <span class=\"n\">xlim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">1.4<\/span><span class=\"p\">),<\/span> <span class=\"n\">ylim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">0.7<\/span><span class=\"p\">),<\/span>\r\n    <span class=\"n\">aspect<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'equal'<\/span><span class=\"p\">,<\/span>\r\n    <span class=\"n\">legend<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">True<\/span>\r\n<span class=\"p\">)<\/span><span class=\"p\">;<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_svg output_subarea \"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-10027\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/symtosha05.svg\" alt=\"\" width=\"640\" height=\"427\" \/><\/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<p>\u6700\u5927\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b\u89d2\u5ea6 $\\theta$ \u3092\u8abf\u3079\u308b\u305f\u3081\u306b\uff0c\u7740\u5730\u70b9\u4ed8\u8fd1\u3092\u62e1\u5927\u8868\u793a\u3057\u307e\u3059\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[12]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">H<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.2<\/span>\r\n\r\n<span class=\"n\">graphics<\/span><span class=\"p\">(<\/span>\r\n    <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">47<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">46<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">45<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">44<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">43<\/span><span class=\"p\">),<\/span>\r\n    <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">42<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">41<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">40<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">39<\/span><span class=\"p\">),<\/span> <span class=\"n\">lp2d<\/span><span class=\"p\">(<\/span><span class=\"mi\">38<\/span><span class=\"p\">),<\/span>\r\n    <span class=\"c1\"># \u30bf\u30a4\u30c8\u30eb\u306b H \u306e\u5024\u3092\u5165\u308c\u308b<\/span>\r\n    <span class=\"n\">title<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'H = <\/span><span class=\"si\">%.1f<\/span><span class=\"s1\"> \u306e\u5834\u5408\u306e\u659c\u65b9\u6295\u5c04'<\/span> <span class=\"o\">%<\/span> <span class=\"n\">H<\/span><span class=\"p\">,<\/span>\r\n    <span class=\"n\">xlabel<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'X'<\/span><span class=\"p\">,<\/span> <span class=\"n\">ylabel<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'Y'<\/span><span class=\"p\">,<\/span> \r\n    <span class=\"c1\"># \u7740\u5730\u70b9\u4ed8\u8fd1\u3092\u62e1\u5927\u8868\u793a<\/span>\r\n    <span class=\"n\">xlim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"mf\">1.15<\/span><span class=\"p\">,<\/span> <span class=\"mf\">1.19<\/span><span class=\"p\">),<\/span> <span class=\"n\">ylim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"mi\">0<\/span><span class=\"p\">,<\/span> <span class=\"mf\">0.02<\/span><span class=\"p\">),<\/span>\r\n    <span class=\"n\">legend<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">True<\/span>\r\n<span class=\"p\">)<\/span><span class=\"p\">;<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_svg output_subarea \"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-10028\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/symtosha06.svg\" alt=\"\" width=\"640\" height=\"427\" \/><\/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<p>\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u304c\u6700\u5927\u3068\u306a\u308b\u306e\u306f $\\theta=45^{\\circ}$ \u306e\u3068\u304d<span style=\"font-family: helvetica, arial, sans-serif;\"><strong>\u3067\u306f\u306a\u3044\uff01<\/strong><\/span>\u3053\u3068\u304c\u308f\u304b\u308b\u3060\u308d\u3046\u3002<\/p>\n<h3 id=\"$H-=-0.2$-\u306e\u5834\u5408\u306e\u6c34\u5e73\u5230\u9054\u8ddd\u96e2-$L$\">$H = 0.2$ \u306e\u5834\u5408\u306e\u6c34\u5e73\u5230\u9054\u8ddd\u96e2 $L$<\/h3>\n<p>\u5b9a\u91cf\u7684\u306b\u8abf\u3079\u308b\u305f\u3081\u306b\uff0c\u6c34\u5e73\u5230\u9054\u8ddd\u96e2 <code>L(th, H)<\/code> \u3092\u30b0\u30e9\u30d5\u306b\u3057\u3066\u307f\u308b\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[13]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">var<\/span><span class=\"p\">(<\/span><span class=\"s1\">'th'<\/span><span class=\"p\">)<\/span> <span class=\"c1\"># \u5909\u6570 th \u306e\u5ba3\u8a00\u30fb\u521d\u671f\u5316\uff08\u524d\u306b\u4f7f\u3063\u3066\u3044\u308b\u306e\u3067\uff09<\/span>\r\n<span class=\"n\">H<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.2<\/span>\r\n\r\n<span class=\"n\">graphics<\/span><span class=\"p\">(<\/span>\r\n    <span class=\"n\">line<\/span><span class=\"p\">(<\/span><span class=\"n\">L<\/span><span class=\"p\">(<\/span><span class=\"n\">rad<\/span><span class=\"p\">(<\/span><span class=\"n\">th<\/span><span class=\"p\">),<\/span> <span class=\"n\">H<\/span><span class=\"p\">),<\/span> <span class=\"p\">(<\/span><span class=\"n\">th<\/span><span class=\"p\">,<\/span> <span class=\"mi\">38<\/span><span class=\"p\">,<\/span> <span class=\"mi\">47<\/span><span class=\"p\">)),<\/span> \r\n    <span class=\"n\">title<\/span> <span class=\"o\">=<\/span> <span class=\"s2\">\"H = <\/span><span class=\"si\">%.1f<\/span><span class=\"s2\"> \u306e\u5834\u5408\u306e\u6253\u3061\u51fa\u3057\u89d2\u5ea6\u3068\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\"<\/span> <span class=\"o\">%<\/span> <span class=\"n\">H<\/span><span class=\"p\">,<\/span>\r\n    <span class=\"n\">xlabel<\/span> <span class=\"o\">=<\/span> <span class=\"s2\">\"\u6253\u3061\u51fa\u3057\u89d2\u5ea6 (\u00b0)\"<\/span><span class=\"p\">,<\/span> \r\n    <span class=\"n\">ylabel<\/span> <span class=\"o\">=<\/span> <span class=\"s2\">\"\u898f\u683c\u5316\u3055\u308c\u305f\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\"<\/span><span class=\"p\">,<\/span> \r\n<span class=\"p\">)<\/span><span class=\"p\">;<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_svg output_subarea \"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-10029\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/symtosha07.svg\" alt=\"\" width=\"640\" height=\"427\" \/><\/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<p>\u4e0a\u306e\u30b0\u30e9\u30d5\u3092\u307f\u308b\u3068\uff0c$H=0.2$ \u306e\u5834\u5408\u306b\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u304c\u6700\u5927\u3068\u306a\u308b\u306e\u306f $\\theta=45^{\\circ}$ \u306e\u3068\u304d\u3067\u306f\u306a\u304f&#8230; \u3069\u306e\u304f\u3089\u3044\uff1f $40^{\\circ}$ \u306e\u3042\u305f\u308a\uff1f<\/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=\"L-\u304c\u6700\u5927\u3068\u306a\u308b\u89d2\u5ea6-$\\theta_{\\rm-m}$-\u3092\u6570\u5024\u7684\u306b\u6c42\u3081\u308b-nsolve()\">L \u304c\u6700\u5927\u3068\u306a\u308b\u89d2\u5ea6 $\\theta_{\\rm m}$ \u3092\u6570\u5024\u7684\u306b\u6c42\u3081\u308b <code>nsolve()<\/code><\/h4>\n<p>$\\displaystyle \\frac{dL(\\theta)}{d\\theta} = 0$ \u3068\u306a\u308b\u89d2\u5ea6 $\\theta_{\\rm m}$ \u3092\u6570\u5024\u7684\u306b\u6c42\u3081\u3066\u307f\u308b\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[14]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">var<\/span><span class=\"p\">(<\/span><span class=\"s1\">'theta H'<\/span><span class=\"p\">,<\/span> <span class=\"n\">positive<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">True<\/span><span class=\"p\">)<\/span> <span class=\"c1\"># \u5909\u6570\u306e\u521d\u671f\u5316\uff08\u6570\u5024\u3092\u5165\u308c\u3066\u305f\u306e\u3067\uff09<\/span>\r\n\r\n<span class=\"k\">def<\/span>  <span class=\"nf\">dL<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"k\">return<\/span> <span class=\"n\">diff<\/span><span class=\"p\">(<\/span><span class=\"n\">L<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">),<\/span> <span class=\"n\">theta<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">dL<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[14]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle -\\left(\\sqrt{2 H + \\sin^{2}{\\left(\\theta \\right)}} + \\sin{\\left(\\theta \\right)}\\right) \\sin{\\left(\\theta \\right)} + \\left(\\cos{\\left(\\theta \\right)} + \\frac{\\sin{\\left(\\theta \\right)} \\cos{\\left(\\theta \\right)}}{\\sqrt{2 H + \\sin^{2}{\\left(\\theta \\right)}}}\\right) \\cos{\\left(\\theta \\right)}$<\/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<p>\u65b9\u7a0b\u5f0f\u306e\u6570\u5024\u89e3\u3092\u6c42\u3081\u308b\u3002\u6388\u696d\u3067\u3084\u308a\u307e\u3057\u305f\u3002\u4ee5\u4e0b\u3092\u53c2\u7167\uff1a<\/p>\n<ul>\n<li><a href=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/%e3%82%b3%e3%83%b3%e3%83%94%e3%83%a5%e3%83%bc%e3%82%bf%e6%bc%94%e7%bf%92\/python-%e3%81%a7%e6%95%b0%e5%80%a4%e8%a7%a3%e6%9e%90\/sympy-%e3%81%a7%ef%bc%88%e3%81%82%e3%81%88%e3%81%a6%ef%bc%89%e6%95%b0%e5%80%a4%e8%a7%a3%e6%9e%90\/\">SymPy \u3066\u3099\uff08\u3042\u3048\u3066\uff09\u6570\u5024\u89e3\u6790<\/a>\n<ul>\n<li><a href=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/%e3%82%b3%e3%83%b3%e3%83%94%e3%83%a5%e3%83%bc%e3%82%bf%e6%bc%94%e7%bf%92\/python-%e3%81%a7%e6%95%b0%e5%80%a4%e8%a7%a3%e6%9e%90\/sympy-%e3%81%a7%ef%bc%88%e3%81%82%e3%81%88%e3%81%a6%ef%bc%89%e6%95%b0%e5%80%a4%e8%a7%a3%e6%9e%90\/#SymPy_nsolve\">SymPy \u306b\u3088\u308b\u65b9\u7a0b\u5f0f\u306e\u6570\u5024\u89e3\uff1ansolve()<\/a><\/li>\n<\/ul>\n<\/li>\n<\/ul>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[15]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"c1\"># 40\u00b0\u3092\u30e9\u30b8\u30a2\u30f3\u3078\u63db\u7b97\u3057\u3066\u5c0f\u6570\u70b9\u8868\u793a\u306b<\/span>\r\n<span class=\"n\">th1<\/span> <span class=\"o\">=<\/span> <span class=\"n\">N<\/span><span class=\"p\">(<\/span><span class=\"n\">rad<\/span><span class=\"p\">(<\/span><span class=\"mi\">40<\/span><span class=\"p\">))<\/span>\r\n<span class=\"n\">H<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.2<\/span>\r\n<span class=\"c1\"># theta = th1 \u306e\u4ed8\u8fd1\u3067\u65b9\u7a0b\u5f0f\u3092\u6570\u5024\u89e3\u3092\u63a2\u3059<\/span>\r\n<span class=\"n\">thetam<\/span> <span class=\"o\">=<\/span> <span class=\"n\">nsolve<\/span><span class=\"p\">(<\/span><span class=\"n\">dL<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">),<\/span> <span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">th1<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"nb\">print<\/span><span class=\"p\">(<\/span><span class=\"s1\">'H = <\/span><span class=\"si\">%.1f<\/span><span class=\"s1\"> \u306e\u3068\u304d'<\/span> <span class=\"o\">%<\/span> <span class=\"n\">H<\/span><span class=\"p\">)<\/span>\r\n<span class=\"c1\"># \u5f97\u3089\u308c\u305f\u6570\u5024\u89e3\u3092\u5ea6\u306b\u63db\u7b97<\/span>\r\n<span class=\"n\">thm<\/span> <span class=\"o\">=<\/span> <span class=\"n\">N<\/span><span class=\"p\">(<\/span><span class=\"n\">deg<\/span><span class=\"p\">(<\/span><span class=\"n\">thetam<\/span><span class=\"p\">))<\/span>\r\n<span class=\"nb\">print<\/span><span class=\"p\">(<\/span><span class=\"s1\">'  \u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b\u89d2\u5ea6\u306f <\/span><span class=\"si\">%.5f<\/span><span class=\"s1\">\u00b0'<\/span> <span class=\"o\">%<\/span> <span class=\"n\">thm<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">Lm<\/span> <span class=\"o\">=<\/span> <span class=\"n\">L<\/span><span class=\"p\">(<\/span><span class=\"n\">thetam<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nb\">print<\/span><span class=\"p\">(<\/span><span class=\"s1\">'  \u305d\u306e\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u306f <\/span><span class=\"si\">%.5f<\/span><span class=\"s1\">'<\/span> <span class=\"o\">%<\/span> <span class=\"n\">Lm<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nb\">print<\/span><span class=\"p\">()<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_subarea output_stream output_stdout output_text\">\n<pre>H = 0.2 \u306e\u3068\u304d\r\n  \u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b\u89d2\u5ea6\u306f 40.20297\u00b0\r\n  \u305d\u306e\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u306f 1.18322\r\n\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h3 id=\"\u25cb\u7df4\u7fd2\uff1a$H$-\u3092\u5909\u3048\u3066\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3092\u6c42\u3081\u308b\">\u25cb\u7df4\u7fd2\uff1a$H$ \u3092\u5909\u3048\u3066\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3092\u6c42\u3081\u308b<\/h3>\n<p>\u9ad8\u3055 $H$ \u304b\u3089\u306e\u659c\u65b9\u6295\u5c04\u306e\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u304c\u6700\u5927\u3068\u306a\u308b\u6253\u3061\u51fa\u3057\u89d2\u5ea6\u3068\uff0c\u305d\u306e\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3092\u6b21\u306e $H$ \u306b\u3064\u3044\u3066\u6c42\u3081\u3088\u3002<\/p>\n<p>$H = 0, 0.2, 0.4, 0.6, 0.8, 1.0$<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[16]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">Hs<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"mf\">0.2<\/span><span class=\"o\">*<\/span><span class=\"n\">i<\/span> <span class=\"k\">for<\/span> <span class=\"n\">i<\/span> <span class=\"ow\">in<\/span> <span class=\"nb\">range<\/span><span class=\"p\">(<\/span><span class=\"mi\">6<\/span><span class=\"p\">)]<\/span>\r\n<span class=\"n\">Hs<\/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 output_prompt\">Out[16]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\left[ 0.0, \\ 0.2, \\ 0.4, \\ 0.6, \\ 0.8, \\ 1.0\\right]$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[17]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"k\">for<\/span> <span class=\"n\">H<\/span> <span class=\"ow\">in<\/span> <span class=\"n\">Hs<\/span><span class=\"p\">:<\/span>\r\n    <span class=\"n\">thetam<\/span> <span class=\"o\">=<\/span> <span class=\"n\">nsolve<\/span><span class=\"p\">(<\/span><span class=\"n\">dL<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">),<\/span> <span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">th1<\/span><span class=\"p\">)<\/span>\r\n\r\n    <span class=\"nb\">print<\/span><span class=\"p\">(<\/span><span class=\"s1\">'H = <\/span><span class=\"si\">%.1f<\/span><span class=\"s1\"> \u306e\u3068\u304d'<\/span> <span class=\"o\">%<\/span> <span class=\"n\">H<\/span><span class=\"p\">)<\/span>\r\n    <span class=\"c1\"># \u5f97\u3089\u308c\u305f\u6570\u5024\u89e3\u3092\u5ea6\u306b\u63db\u7b97<\/span>\r\n    <span class=\"n\">thm<\/span> <span class=\"o\">=<\/span> <span class=\"n\">N<\/span><span class=\"p\">(<\/span><span class=\"n\">deg<\/span><span class=\"p\">(<\/span><span class=\"n\">thetam<\/span><span class=\"p\">))<\/span>\r\n    <span class=\"nb\">print<\/span><span class=\"p\">(<\/span><span class=\"s1\">'  \u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b\u89d2\u5ea6\u306f <\/span><span class=\"si\">%.5f<\/span><span class=\"s1\">\u00b0'<\/span> <span class=\"o\">%<\/span> <span class=\"n\">thm<\/span><span class=\"p\">)<\/span>\r\n    <span class=\"n\">Lm<\/span> <span class=\"o\">=<\/span> <span class=\"n\">L<\/span><span class=\"p\">(<\/span><span class=\"n\">thetam<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">)<\/span>\r\n    <span class=\"nb\">print<\/span><span class=\"p\">(<\/span><span class=\"s1\">'  \u305d\u306e\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u306f <\/span><span class=\"si\">%.5f<\/span><span class=\"s1\">'<\/span> <span class=\"o\">%<\/span> <span class=\"n\">Lm<\/span><span class=\"p\">)<\/span>\r\n    <span class=\"nb\">print<\/span><span class=\"p\">()<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_subarea output_stream output_stdout output_text\">\n<pre>H = 0.0 \u306e\u3068\u304d\r\n  \u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b\u89d2\u5ea6\u306f 45.00000\u00b0\r\n  \u305d\u306e\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u306f 1.00000\r\n\r\nH = 0.2 \u306e\u3068\u304d\r\n  \u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b\u89d2\u5ea6\u306f 40.20297\u00b0\r\n  \u305d\u306e\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u306f 1.18322\r\n\r\nH = 0.4 \u306e\u3068\u304d\r\n  \u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b\u89d2\u5ea6\u306f 36.69923\u00b0\r\n  \u305d\u306e\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u306f 1.34164\r\n\r\nH = 0.6 \u306e\u3068\u304d\r\n  \u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b\u89d2\u5ea6\u306f 33.98784\u00b0\r\n  \u305d\u306e\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u306f 1.48324\r\n\r\nH = 0.8 \u306e\u3068\u304d\r\n  \u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b\u89d2\u5ea6\u306f 31.80610\u00b0\r\n  \u305d\u306e\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u306f 1.61245\r\n\r\nH = 1.0 \u306e\u3068\u304d\r\n  \u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b\u89d2\u5ea6\u306f 30.00000\u00b0\r\n  \u305d\u306e\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u306f 1.73205\r\n\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h3 id=\"\u25cb\u7df4\u7fd2\uff1a\u9ad8\u3055-$H$-\u3068-\u6253\u3061\u51fa\u3057\u89d2\u5ea6-$\\theta_{\\rm-m}$-\u306e\u30b0\u30e9\u30d5\">\u25cb\u7df4\u7fd2\uff1a\u9ad8\u3055 $H$ \u3068 \u6253\u3061\u51fa\u3057\u89d2\u5ea6 $\\theta_{\\rm m}$ \u306e\u30b0\u30e9\u30d5<\/h3>\n<p>\u9ad8\u3055 $H$ \u3092\u6a2a\u8ef8\u306b\uff0c\u305d\u306e\u3068\u304d\u306e\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u304c\u6700\u5927\u3068\u306a\u308b\u6253\u3061\u51fa\u3057\u89d2\u5ea6 $\\theta_{\\rm m}$ \u3092\u7e26\u8ef8\u306b\u3057\u305f\u30b0\u30e9\u30d5\u3092\u63cf\u3051\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[18]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">Hs<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"mf\">0.2<\/span><span class=\"o\">*<\/span><span class=\"n\">i<\/span> <span class=\"k\">for<\/span> <span class=\"n\">i<\/span> <span class=\"ow\">in<\/span> <span class=\"nb\">range<\/span><span class=\"p\">(<\/span><span class=\"mi\">6<\/span><span class=\"p\">)]<\/span>\r\n<span class=\"n\">thms<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[]<\/span>\r\n<span class=\"n\">Lms<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[]<\/span>\r\n\r\n<span class=\"k\">for<\/span> <span class=\"n\">H<\/span> <span class=\"ow\">in<\/span> <span class=\"n\">Hs<\/span><span class=\"p\">:<\/span>\r\n    <span class=\"n\">thetam<\/span> <span class=\"o\">=<\/span> <span class=\"n\">nsolve<\/span><span class=\"p\">(<\/span><span class=\"n\">dL<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">),<\/span> <span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">th1<\/span><span class=\"p\">)<\/span>\r\n\r\n    <span class=\"n\">thm<\/span> <span class=\"o\">=<\/span> <span class=\"n\">N<\/span><span class=\"p\">(<\/span><span class=\"n\">deg<\/span><span class=\"p\">(<\/span><span class=\"n\">thetam<\/span><span class=\"p\">))<\/span>\r\n    <span class=\"n\">thms<\/span> <span class=\"o\">+=<\/span> <span class=\"p\">[<\/span><span class=\"n\">thm<\/span><span class=\"p\">]<\/span>\r\n    <span class=\"n\">Lm<\/span> <span class=\"o\">=<\/span> <span class=\"n\">L<\/span><span class=\"p\">(<\/span><span class=\"n\">thetam<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">)<\/span>\r\n    <span class=\"n\">Lms<\/span> <span class=\"o\">+=<\/span> <span class=\"p\">[<\/span><span class=\"n\">Lm<\/span><span class=\"p\">]<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[19]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">graphics<\/span><span class=\"p\">(<\/span>\r\n    <span class=\"n\">list_2d<\/span><span class=\"p\">(<\/span><span class=\"n\">Hs<\/span><span class=\"p\">,<\/span> <span class=\"n\">thms<\/span><span class=\"p\">,<\/span> <span class=\"n\">scatter<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">True<\/span><span class=\"p\">),<\/span> \r\n    <span class=\"n\">title<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u304c\u6700\u5927\u3068\u306a\u308b\u6253\u3061\u51fa\u3057\u89d2\u5ea6'<\/span><span class=\"p\">,<\/span> \r\n    <span class=\"n\">xlabel<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'\u898f\u683c\u5316\u3055\u308c\u305f\u9ad8\u3055 H'<\/span><span class=\"p\">,<\/span> \r\n    <span class=\"n\">ylabel<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'\u6253\u3061\u51fa\u3057\u89d2\u5ea6 (\u00b0)'<\/span><span class=\"p\">,<\/span> \r\n    <span class=\"n\">ylim<\/span> <span class=\"o\">=<\/span> <span class=\"p\">(<\/span><span class=\"mi\">29<\/span><span class=\"p\">,<\/span> <span class=\"mi\">46<\/span><span class=\"p\">)<\/span>\r\n<span class=\"p\">)<\/span><span class=\"p\">;<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_svg output_subarea \"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-10030\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/symtosha08.svg\" alt=\"\" width=\"640\" height=\"427\" \/><\/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=\"\u25cb\u7df4\u7fd2\uff1a\u9ad8\u3055-$H$-\u3068\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2-$L_{\\rm-m}$-\u306e\u30b0\u30e9\u30d5\">\u25cb\u7df4\u7fd2\uff1a\u9ad8\u3055 $H$ \u3068\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2 $L_{\\rm m}$ \u306e\u30b0\u30e9\u30d5<\/h3>\n<p>\u9ad8\u3055 $H$ \u3092\u6a2a\u8ef8\u306b\uff0c\u305d\u306e\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2 $L_{\\rm m} \\equiv L(\\theta_{\\rm m}, H)$ \u3092\u7e26\u8ef8\u306b\u3057\u305f\u30b0\u30e9\u30d5\u3092\u63cf\u3051\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[20]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">graphics<\/span><span class=\"p\">(<\/span>\r\n    <span class=\"n\">list_2d<\/span><span class=\"p\">(<\/span><span class=\"n\">Hs<\/span><span class=\"p\">,<\/span> <span class=\"n\">Lms<\/span><span class=\"p\">,<\/span> <span class=\"n\">scatter<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">True<\/span><span class=\"p\">),<\/span> \r\n    <span class=\"n\">title<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2'<\/span><span class=\"p\">,<\/span> \r\n    <span class=\"n\">xlabel<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'\u898f\u683c\u5316\u3055\u308c\u305f\u9ad8\u3055 H'<\/span><span class=\"p\">,<\/span> \r\n    <span class=\"n\">ylabel<\/span> <span class=\"o\">=<\/span> <span class=\"s1\">'\u898f\u683c\u5316\u3055\u308c\u305f\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2'<\/span>\r\n<span class=\"p\">)<\/span><span class=\"p\">;<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_svg output_subarea \"><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-10031\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/symtosha09.svg\" alt=\"\" width=\"640\" height=\"427\" \/><\/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=\"\u25cb\u7df4\u7fd2\uff1a\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3092\u6c42\u3081\u308b\">\u25cb\u7df4\u7fd2\uff1a\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3092\u6c42\u3081\u308b<\/h3>\n<p>$h = 5\\,\\mbox{(m)}$ \u306e\u9ad8\u3055\u304b\u3089\u901f\u3055 $v_0 = 10\\,\\mbox{(m\/s)}$ \u3067\u659c\u65b9\u6295\u5c04\u3057\u305f\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u306f\u4f55 $\\mbox{m}$ \u304b\u3002\u307e\u305f\u305d\u306e\u3068\u304d\u306e\u6253\u51fa\u3057\u89d2\u5ea6\u306f\u4f55\u5ea6\u304b\u3002<\/p>\n<p>\u30d2\u30f3\u30c8\uff1a<\/p>\n<p>$h = 5\\,\\mbox{(m)}$ \u306f\u898f\u683c\u5316\u3055\u308c\u305f\u9ad8\u3055 $H$ \u3067\u306f\u3044\u304f\u3089\u304b\uff0c\u307e\u305f\u898f\u683c\u5316\u3055\u308c\u305f\u6700\u5927\u5230\u9054\u8ddd\u96e2 $L_{\\rm m}$ \u3092\u6b21\u5143\u3092\u3082\u3063\u305f\u91cf\u306b\u76f4\u3059\u3068\u3044\u304f\u3089\u304b\uff0c\u3068\u3044\u3046\u3053\u3068\u3002<\/p>\n<p>\\begin{eqnarray}<br \/>\nH &amp;=&amp; \\frac{h}{\\ell} = \\frac{g h}{v_0^2} \\\\<br \/>\nx &amp;=&amp; \\ell X = \\frac{v_0^2}{g} X \\\\<br \/>\n\\therefore\\ \\ x_{\\rm m} &amp;=&amp; \\frac{v_0^2}{g} L_{\\rm m}<br \/>\n\\end{eqnarray}<\/p>\n<p>\u91cd\u529b\u52a0\u901f\u5ea6 $g = 9.80665 \\,(\\mbox{m}\/\\mbox{s}^2)$<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[21]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">h<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">5<\/span>\r\n<span class=\"n\">v0<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">10<\/span>\r\n<span class=\"n\">g<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">9.80665<\/span>\r\n\r\n<span class=\"c1\"># \u898f\u683c\u5316\u3055\u308c\u305f\u9ad8\u3055 H<\/span>\r\n<span class=\"n\">H<\/span> <span class=\"o\">=<\/span> <span class=\"n\">h<\/span> <span class=\"o\">*<\/span> <span class=\"n\">g<\/span><span class=\"o\">\/<\/span><span class=\"n\">v0<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span>\r\n\r\n<span class=\"n\">thetam<\/span> <span class=\"o\">=<\/span> <span class=\"n\">nsolve<\/span><span class=\"p\">(<\/span><span class=\"n\">dL<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">),<\/span> <span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">N<\/span><span class=\"p\">(<\/span><span class=\"n\">rad<\/span><span class=\"p\">(<\/span><span class=\"mi\">40<\/span><span class=\"p\">)))<\/span>\r\n\r\n<span class=\"n\">xm<\/span> <span class=\"o\">=<\/span> <span class=\"n\">v0<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"o\">\/<\/span><span class=\"n\">g<\/span> <span class=\"o\">*<\/span> <span class=\"n\">L<\/span><span class=\"p\">(<\/span><span class=\"n\">thetam<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nb\">print<\/span><span class=\"p\">(<\/span><span class=\"s1\">'\u9ad8\u3055 h = <\/span><span class=\"si\">%2d<\/span><span class=\"s1\"> m \u306e\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u306f'<\/span> <span class=\"o\">%<\/span> <span class=\"n\">h<\/span><span class=\"p\">,<\/span> <span class=\"n\">end<\/span><span class=\"o\">=<\/span><span class=\"s1\">' '<\/span><span class=\"p\">)<\/span>\r\n<span class=\"nb\">print<\/span><span class=\"p\">(<\/span><span class=\"s1\">'<\/span><span class=\"si\">%9.6f<\/span><span class=\"s1\"> m'<\/span> <span class=\"o\">%<\/span> <span class=\"n\">xm<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"nb\">print<\/span><span class=\"p\">(<\/span><span class=\"s1\">'\u305d\u306e\u3068\u304d\u306e\u6253\u3061\u51fa\u3057\u89d2\u5ea6\u306f <\/span><span class=\"si\">%9.6f<\/span><span class=\"s1\">\u00b0'<\/span> <span class=\"o\">%<\/span> <span class=\"n\">N<\/span><span class=\"p\">(<\/span><span class=\"n\">deg<\/span><span class=\"p\">(<\/span><span class=\"n\">thetam<\/span><span class=\"p\">)))<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_subarea output_stream output_stdout output_text\">\n<pre>\u9ad8\u3055 h =  5 m \u306e\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u306f 14.351088 m\r\n\u305d\u306e\u3068\u304d\u306e\u6253\u3061\u51fa\u3057\u89d2\u5ea6\u306f 35.395688\u00b0\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h3 id=\"\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b\u89d2\u5ea6\u3092\u89e3\u6790\u7684\u306b\u6c42\u3081\u308b\">\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b\u89d2\u5ea6\u3092\u89e3\u6790\u7684\u306b\u6c42\u3081\u308b<\/h3>\n<h4 id=\"\u89e3\u6cd51\uff1a\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b-$\\sin\\theta_{\\rm-m}$-\u3092\u89e3\u6790\u7684\u306b\u6c42\u3081\u308b\">\u89e3\u6cd51\uff1a\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b $\\sin\\theta_{\\rm m}$ \u3092\u89e3\u6790\u7684\u306b\u6c42\u3081\u308b<\/h4>\n<p>\u6c34\u5e73\u5230\u9054\u8ddd\u96e2 $L(\\theta)$ \u3092 $s \\equiv \\sin \\theta$ \u3092\u4f7f\u3063\u3066\u8868\u3057\u3066\u307f\u307e\u3059\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[22]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">var<\/span><span class=\"p\">(<\/span><span class=\"s1\">'H'<\/span><span class=\"p\">,<\/span> <span class=\"n\">positive<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">True<\/span><span class=\"p\">)<\/span> <span class=\"c1\"># \u6570\u5024\u3092\u5165\u308c\u3066\u305f\u306e\u3067\u521d\u671f\u5316<\/span>\r\n<span class=\"n\">display<\/span><span class=\"p\">(<\/span><span class=\"n\">L<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">))<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_latex output_subarea \">$\\displaystyle \\left(\\sqrt{2 H + \\sin^{2}{\\left(\\theta \\right)}} + \\sin{\\left(\\theta \\right)}\\right) \\cos{\\left(\\theta \\right)}$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[23]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">var<\/span><span class=\"p\">(<\/span><span class=\"s1\">'s'<\/span><span class=\"p\">,<\/span> <span class=\"n\">positive<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">True<\/span><span class=\"p\">)<\/span> <span class=\"c1\"># s &gt; 0 \u306e\u307f\u3092\u8003\u3048\u308b<\/span>\r\n\r\n<span class=\"k\">def<\/span> <span class=\"nf\">Ls<\/span><span class=\"p\">(<\/span><span class=\"n\">s<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"k\">return<\/span> <span class=\"n\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"mi\">1<\/span><span class=\"o\">-<\/span><span class=\"n\">s<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"p\">(<\/span><span class=\"n\">s<\/span> <span class=\"o\">+<\/span> <span class=\"n\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">H<\/span> <span class=\"o\">+<\/span> <span class=\"n\">s<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">))<\/span>\r\n<span class=\"n\">Ls<\/span><span class=\"p\">(<\/span><span class=\"n\">s<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[23]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\sqrt{1 -s^{2}} \\left(s + \\sqrt{2 H + s^{2}}\\right)$<\/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<p><code>dLs<\/code>$\\displaystyle \\equiv\\frac{d L(s)}{d s}$<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[24]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">dLs<\/span> <span class=\"o\">=<\/span> <span class=\"n\">diff<\/span><span class=\"p\">(<\/span><span class=\"n\">Ls<\/span><span class=\"p\">(<\/span><span class=\"n\">s<\/span><span class=\"p\">),<\/span> <span class=\"n\">s<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">dLs<\/span> <span class=\"o\">=<\/span> <span class=\"n\">dLs<\/span><span class=\"o\">.<\/span><span class=\"n\">simplify<\/span><span class=\"p\">()<\/span>\r\n<span class=\"n\">dLs<\/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 output_prompt\">Out[24]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\frac{\\left(s + \\sqrt{2 H + s^{2}}\\right) \\left(-s^{2} -s \\sqrt{2 H + s^{2}} + 1\\right)}{\\sqrt{1 -s^{2}} \\sqrt{2 H + s^{2}}}$<\/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<p><code>dLs = 0<\/code> \u3092 <code>s<\/code> \u306b\u3064\u3044\u3066\u89e3\u304f\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[25]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">sols<\/span> <span class=\"o\">=<\/span> <span class=\"n\">solve<\/span><span class=\"p\">(<\/span><span class=\"n\">dLs<\/span><span class=\"p\">,<\/span> <span class=\"n\">s<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">sols<\/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 output_prompt\">Out[25]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\left[ \\frac{\\sqrt{2}}{2 \\sqrt{H + 1}}\\right]$<\/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<p>\u3068\u3044\u3046\u3053\u3068\u3067\uff0c\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b $s_{\\rm m} = \\sin\\theta_{\\rm m} = $ <code>sols[0]<\/code> \u306f<\/p>\n<p>\\begin{eqnarray}<br \/>\ns_{\\rm m} = \\sin\\theta_{\\rm m} &amp;=&amp; \\frac{1}{\\sqrt{2 (1 + H)}} \\\\<br \/>\n\\therefore\\ \\ \\theta_{\\rm m} &amp;=&amp; \\arcsin\\left(\\frac{1}{\\sqrt{2 (1 + H)}}\\right)<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<h5 id=\"\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\">\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2<\/h5>\n<p>\u305d\u306e\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2 $L_{\\rm m}$ \u306f&#8230;<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[26]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">sm<\/span> <span class=\"o\">=<\/span> <span class=\"n\">sols<\/span><span class=\"p\">[<\/span><span class=\"mi\">0<\/span><span class=\"p\">]<\/span>\r\n<span class=\"n\">Ls<\/span><span class=\"p\">(<\/span><span class=\"n\">sm<\/span><span class=\"p\">)<\/span><span class=\"o\">.<\/span><span class=\"n\">factor<\/span><span class=\"p\">()<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[26]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\frac{\\sqrt{2 H + 1} \\left(\\sqrt{4 H^{2} + 4 H + 1} + 1\\right)}{2 \\left(H + 1\\right)}$<\/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<p>\u7c21\u5358\u5316\u3057\u3066\u304f\u308c\u306a\u3044\u306a\u3041&#8230;<\/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=\"\u89e3\u6cd52\uff1a\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b-$\\tan\\theta_{\\rm-m}$-\u3092\u89e3\u6790\u7684\u306b\u6c42\u3081\u308b\">\u89e3\u6cd52\uff1a\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b $\\tan\\theta_{\\rm m}$ \u3092\u89e3\u6790\u7684\u306b\u6c42\u3081\u308b<\/h4>\n<p>\u6c34\u5e73\u5230\u9054\u8ddd\u96e2 $L(\\theta)$ \u3092 $\\tan \\theta$ \u3092\u4f7f\u3063\u3066\u8868\u3057\u3066\u307f\u307e\u3059\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[27]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">var<\/span><span class=\"p\">(<\/span><span class=\"s1\">'H'<\/span><span class=\"p\">,<\/span> <span class=\"n\">positive<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">True<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">L<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">,<\/span> <span class=\"n\">H<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[27]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\left(\\sqrt{2 H + \\sin^{2}{\\left(\\theta \\right)}} + \\sin{\\left(\\theta \\right)}\\right) \\cos{\\left(\\theta \\right)}$<\/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<p>\u4eba\u529b\u3067\u6574\u7406\u3057\u3066\u307f\u3088\u3046\u3002<\/p>\n<p>\\begin{eqnarray}<br \/>\nL(\\theta) &amp;=&amp;\\left(\\sqrt{2 H + \\sin^{2}{\\left(\\theta \\right)}} + \\sin{\\left(\\theta \\right)}\\right) \\cos{\\left(\\theta \\right)} \\\\<br \/>\n&amp;=&amp;<br \/>\n\\left(\\sqrt{\\frac{2H}{\\cos^2\\theta}+ \\frac{\\sin^2\\theta}{\\cos^2\\theta}} + \\frac{\\sin\\theta}{\\cos\\theta}<br \/>\n\\right) \\cos^2\\theta \\\\<br \/>\n&amp;=&amp; \\left(\\sqrt{2H \\left(1+\\tan^2\\theta\\right)+ \\tan^2\\theta} + \\tan\\theta<br \/>\n\\right) \\frac{1}{1+\\tan^2\\theta} \\\\<br \/>\n&amp;=&amp; \\left(\\sqrt{2H \\left(1+t^2\\right)+ t^2} + t<br \/>\n\\right) \\frac{1}{1+t^2} \\\\<br \/>\n&amp;=&amp; \\frac{t + \\sqrt{2H + t^2 (1 + 2H)}}{1+t^2} \\\\<br \/>\nt &amp;\\equiv&amp; \\tan\\theta<br \/>\n\\end{eqnarray}<\/p>\n<p>\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3092\u5909\u6570 $t \\equiv \\tan\\theta$ \u306e\u95a2\u6570\u3068\u3057\u3066\u5b9a\u7fa9\u3057\u307e\u3059\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[28]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">var<\/span><span class=\"p\">(<\/span><span class=\"s1\">'t H'<\/span><span class=\"p\">,<\/span> <span class=\"n\">positive<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">True<\/span><span class=\"p\">)<\/span> <span class=\"c1\"># positive = True \u304c\u30ad\u30e2<\/span>\r\n<span class=\"k\">def<\/span> <span class=\"nf\">Lt<\/span><span class=\"p\">(<\/span><span class=\"n\">t<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"k\">return<\/span> <span class=\"p\">(<\/span><span class=\"n\">t<\/span><span class=\"o\">+<\/span><span class=\"n\">sqrt<\/span><span class=\"p\">((<\/span><span class=\"mi\">1<\/span><span class=\"o\">+<\/span><span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">H<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"n\">t<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">+<\/span> <span class=\"mi\">2<\/span><span class=\"o\">*<\/span><span class=\"n\">H<\/span><span class=\"p\">))<\/span><span class=\"o\">\/<\/span><span class=\"p\">(<\/span><span class=\"mi\">1<\/span> <span class=\"o\">+<\/span> <span class=\"n\">t<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">Lt<\/span><span class=\"p\">(<\/span><span class=\"n\">t<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[28]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\frac{t + \\sqrt{2 H + t^{2} \\left(2 H + 1\\right)}}{t^{2} + 1}$<\/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<p>$\\displaystyle \\frac{d L_{t}}{d t} = 0$ \u3068\u306a\u308b $t$ \u3092\u6c42\u3081\u307e\u3059\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[29]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">dLt<\/span> <span class=\"o\">=<\/span> <span class=\"n\">diff<\/span><span class=\"p\">(<\/span><span class=\"n\">Lt<\/span><span class=\"p\">(<\/span><span class=\"n\">t<\/span><span class=\"p\">),<\/span> <span class=\"n\">t<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">sols<\/span> <span class=\"o\">=<\/span> <span class=\"n\">solve<\/span><span class=\"p\">(<\/span><span class=\"n\">dLt<\/span><span class=\"p\">,<\/span> <span class=\"n\">t<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">sols<\/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 output_prompt\">Out[29]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\left[ \\frac{1}{\\sqrt{2 H + 1}}\\right]$<\/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<p>\u3068\u3044\u3046\u3053\u3068\u3067\uff0c\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b $t_{\\rm m} = \\tan\\theta_{\\rm m} = $ <code>sols[0]<\/code> \u306f<\/p>\n<p>\\begin{eqnarray}<br \/>\nt_{\\rm m} = \\tan\\theta_{\\rm m} &amp;=&amp; \\frac{1}{\\sqrt{1 + 2 H}} \\\\<br \/>\n\\therefore\\ \\ \\theta_{\\rm m} &amp;=&amp; \\arctan\\left(\\frac{1}{\\sqrt{1 + 2 H}}\\right)<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<h5 id=\"\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\">\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2<\/h5>\n<p>$\\displaystyle t = t_{\\rm m} \\equiv \\frac{1}{\\sqrt{1+2H}}$ \u306e\u3068\u304d\u6975\u5024\uff08\u6700\u5927\u5024\uff09\u3092\u6301\u3064\u3053\u3068\u304c\u308f\u304b\u308a\uff0c\u305d\u306e\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2 $L_{\\rm m} \\equiv L(t_{\\rm m}, H)$ \u306f&#8230;<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[30]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">tm<\/span> <span class=\"o\">=<\/span> <span class=\"n\">sols<\/span><span class=\"p\">[<\/span><span class=\"mi\">0<\/span><span class=\"p\">]<\/span>\r\n<span class=\"n\">Lt<\/span><span class=\"p\">(<\/span><span class=\"n\">tm<\/span><span class=\"p\">)<\/span><span class=\"o\">.<\/span><span class=\"n\">simplify<\/span><span class=\"p\">()<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[30]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\sqrt{2 H + 1}$<\/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=\"\u89e3\u6cd53\uff1a\u9670\u95a2\u6570\u5b9a\u7406\u3092\u4f7f\u3063\u3066\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3092\u6c42\u3081\u308b\">\u89e3\u6cd53\uff1a\u9670\u95a2\u6570\u5b9a\u7406\u3092\u4f7f\u3063\u3066\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3092\u6c42\u3081\u308b<\/h4>\n<p>\u3042\u3089\u305f\u3081\u3066\uff0c\uff08$H$ \u306f\u5909\u6570\u3067\u306f\u306a\u304f\u30d1\u30e9\u30e1\u30fc\u30bf\u3068\u3057\u3066\uff09<br \/>\n\u7121\u6b21\u5143\u5316\u3055\u308c\u305f\u89e3 $X(T, \\theta), Y(T, \\theta)$ \u306e\u5b9a\u7fa9\u3092\u66f8\u304d\u51fa\u3059\u3068&#8230;<\/p>\n<p>\\begin{eqnarray}<br \/>\nX(T, \\theta) &amp;\\equiv&amp; \\cos\\theta\\cdot T \\\\<br \/>\nY(T, \\theta)&amp;\\equiv&amp; H+ \\sin\\theta\\cdot T -\\frac{1}{2} T^2 \\\\<br \/>\n\\end{eqnarray}<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[31]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">var<\/span><span class=\"p\">(<\/span><span class=\"s1\">'H'<\/span><span class=\"p\">,<\/span> <span class=\"n\">positive<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">True<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"k\">def<\/span> <span class=\"nf\">X<\/span><span class=\"p\">(<\/span><span class=\"n\">T<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"k\">return<\/span> <span class=\"n\">cos<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">)<\/span> <span class=\"o\">*<\/span> <span class=\"n\">T<\/span>\r\n\r\n<span class=\"k\">def<\/span> <span class=\"nf\">Y<\/span><span class=\"p\">(<\/span><span class=\"n\">T<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"k\">return<\/span> <span class=\"n\">H<\/span> <span class=\"o\">+<\/span> <span class=\"n\">sin<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"n\">T<\/span> <span class=\"o\">-<\/span> <span class=\"n\">T<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"o\">\/<\/span><span class=\"mi\">2<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h5 id=\"\u6ede\u7a7a\u6642\u9593-$T_1$\">\u6ede\u7a7a\u6642\u9593 $T_1$<\/h5>\n<p>\u6ede\u7a7a\u6642\u9593 $T_1\\ (&gt; 0)$ \u306f\u4ee5\u4e0b\u306e\u5f0f\u3092\u6e80\u305f\u3059\u3002<\/p>\n<p>$$Y(T_1, \\theta) = 0$$<\/p>\n<p>\u3053\u308c\u304b\u3089\uff0c$T_1$ \u306f $\\theta$ \u306e\u9670\u95a2\u6570\u3068\u307f\u308b\u3053\u3068\u304c\u3067\u304d\u308b\u3002\uff08\u5fc5\u305a\u3057\u3082\u967d\u95a2\u6570\u3068\u3057\u3066\u89e3\u304f\u5fc5\u8981\u306f\u306a\u3044\u3002\uff09<\/p>\n<p>$$T_1 = T_1(\\theta)$$<\/p>\n<p>\u6c34\u5e73\u5230\u9054\u8ddd\u96e2 $L(\\theta,)$ \u306f\uff0c\u6ede\u7a7a\u6642\u9593 $T_1$ \u306e\u9593\u306b\u6c34\u5e73\u65b9\u5411\u306b\u9032\u3080\u8ddd\u96e2\u3067\u3042\u308b\u304b\u3089<\/p>\n<p>$$L(\\theta) \\equiv X(T_1(\\theta), \\theta)$$<\/p>\n<p>\u3067\u3042\u308b\u304b\u3089\uff0c$L(\\theta)$ \u304c\u6700\u5927\u3068\u306a\u308b\u89d2\u5ea6 $\\theta$ \u3092\u6c42\u3081\u308b\u554f\u984c\u306f\uff0c<br \/>\n\u4ee5\u4e0b\u306e\u9023\u7acb\u65b9\u7a0b\u5f0f\u3092 $\\theta$ \u306b\u3064\u3044\u3066\u89e3\u304f\u554f\u984c\u3067\u3042\u308b\u3002<\/p>\n<p>\\begin{eqnarray}<br \/>\nY(T_1(\\theta), \\theta) &amp;=&amp; 0 \\tag{1}\\\\<br \/>\n\\frac{d}{d\\theta} L(\\theta) &amp;=&amp; \\frac{d}{d\\theta} X(T_1(\\theta), \\theta) \\\\<br \/>\n&amp;=&amp;<br \/>\n\\frac{\\partial X}{\\partial T_1} \\frac{d T_1}{d\\theta} + \\frac{\\partial X}{\\partial \\theta} = 0<br \/>\n\\tag{2}<br \/>\n\\end{eqnarray}<\/p>\n<p>(2) \u5f0f\u3067\u5fc5\u8981\u3068\u306a\u308b $\\displaystyle \\frac{d T_1}{d\\theta}$ \u306b\u3064\u3044\u3066\u306f\uff0c\u4ee5\u4e0b\u306e\u3088\u3046\u306b (1) \u5f0f\u306e\u5fae\u5206\u3092\u3068\u3063\u3066\uff0c\u9670\u95a2\u6570\u5b9a\u7406\u304b\u3089\u6c42\u3081\u308b\u3053\u3068\u304c\u3067\u304d\u308b\u3002<\/p>\n<p>\\begin{eqnarray}<br \/>\nY(T_1(\\theta), \\theta) &amp;=&amp; 0 \\\\<br \/>\n\\therefore\\ \\ \\frac{dY}{d\\theta} &amp;=&amp; \\frac{\\partial Y}{\\partial T_1} \\frac{d T_1}{d\\theta} +<br \/>\n\\frac{\\partial Y}{\\partial \\theta} = 0 \\\\<br \/>\n\\therefore\\ \\ \\frac{d T_1}{d\\theta} &amp;=&amp; -\\frac{\\frac{\\partial Y}{\\partial \\theta}}{\\frac{\\partial Y}{\\partial T_1}}<br \/>\n\\end{eqnarray}<\/p>\n<p>\u3053\u308c\u3092 SymPy \u3067\u3084\u308b\u3068&#8230;<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[32]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"c1\"># T1 \u3092 theta \u306e\uff08\u9670\uff09\u95a2\u6570\u3068\u3057\u3066\u5ba3\u8a00<\/span>\r\n<span class=\"n\">T1<\/span> <span class=\"o\">=<\/span> <span class=\"n\">Function<\/span><span class=\"p\">(<\/span><span class=\"s1\">'T1'<\/span><span class=\"p\">)(<\/span><span class=\"n\">theta<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">dY<\/span> <span class=\"o\">=<\/span> <span class=\"n\">diff<\/span><span class=\"p\">(<\/span><span class=\"n\">Y<\/span><span class=\"p\">(<\/span><span class=\"n\">T1<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">),<\/span> <span class=\"n\">theta<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">sols<\/span> <span class=\"o\">=<\/span> <span class=\"n\">solve<\/span><span class=\"p\">(<\/span><span class=\"n\">dY<\/span><span class=\"p\">,<\/span> <span class=\"n\">diff<\/span><span class=\"p\">(<\/span><span class=\"n\">T1<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">))<\/span>\r\n<span class=\"n\">sols<\/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 output_prompt\">Out[32]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\left[ \\frac{T_{1}{\\left(\\theta \\right)} \\cos{\\left(\\theta \\right)}}{T_{1}{\\left(\\theta \\right)} -\\sin{\\left(\\theta \\right)}}\\right]$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[33]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"c1\"># \u89e3\u306f 1 \u500b\u3060\u3051\u306a\u306e\u3067 <\/span>\r\n<span class=\"n\">dT1<\/span> <span class=\"o\">=<\/span> <span class=\"n\">sols<\/span><span class=\"p\">[<\/span><span class=\"mi\">0<\/span><span class=\"p\">]<\/span>\r\n<span class=\"n\">Eq<\/span><span class=\"p\">(<\/span><span class=\"n\">Derivative<\/span><span class=\"p\">(<\/span><span class=\"n\">T1<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">),<\/span> <span class=\"n\">dT1<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[33]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\frac{d}{d \\theta} T_{1}{\\left(\\theta \\right)} = \\frac{T_{1}{\\left(\\theta \\right)} \\cos{\\left(\\theta \\right)}}{T_{1}{\\left(\\theta \\right)} -\\sin{\\left(\\theta \\right)}}$<\/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<h5 id=\"\u6c34\u5e73\u5230\u9054\u8ddd\u96e2-$L$\">\u6c34\u5e73\u5230\u9054\u8ddd\u96e2 $L$<\/h5>\n<p>\u6ede\u7a7a\u6642\u9593 $T_1(\\theta)$ \u306e\u9593\u306b\u6c34\u5e73\u65b9\u5411\u306b\u9032\u3080\u8ddd\u96e2\u306f<\/p>\n<p>$$L(\\theta) = X(T_1(\\theta), \\theta)$$<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[34]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"k\">def<\/span> <span class=\"nf\">L<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">):<\/span>\r\n    <span class=\"k\">return<\/span> <span class=\"n\">X<\/span><span class=\"p\">(<\/span><span class=\"n\">T1<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">L<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[34]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle T_{1}{\\left(\\theta \\right)} \\cos{\\left(\\theta \\right)}$<\/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<h5 id=\"$L$-\u304c\u6700\u5927\u3068\u306a\u308b\u89d2\u5ea6\u3092\u6c42\u3081\u308b\">$L$ \u304c\u6700\u5927\u3068\u306a\u308b\u89d2\u5ea6\u3092\u6c42\u3081\u308b<\/h5>\n<p>$\\displaystyle \\frac{dL}{d\\theta} = 0$ \u3068\u306a\u308b\u89d2\u5ea6\u3092\u6c42\u3081\u308b\u3002\u9670\u95a2\u6570\u5b9a\u7406\u306e\u7d50\u679c\u3092 <code>.subs()<\/code> \u3067\u4ee3\u5165\u3057\u3066&#8230;<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[35]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">dL<\/span> <span class=\"o\">=<\/span> <span class=\"n\">diff<\/span><span class=\"p\">(<\/span><span class=\"n\">L<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">),<\/span> <span class=\"n\">theta<\/span><span class=\"p\">)<\/span><span class=\"o\">.<\/span><span class=\"n\">subs<\/span><span class=\"p\">(<\/span><span class=\"n\">diff<\/span><span class=\"p\">(<\/span><span class=\"n\">T1<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">),<\/span> <span class=\"n\">dT1<\/span><span class=\"p\">)<\/span>\r\n<span class=\"n\">dL<\/span> <span class=\"o\">=<\/span> <span class=\"n\">simplify<\/span><span class=\"p\">(<\/span><span class=\"n\">dL<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">display<\/span><span class=\"p\">(<\/span><span class=\"n\">dL<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_latex output_subarea \">$\\displaystyle \\frac{\\left(-T_{1}{\\left(\\theta \\right)} \\sin{\\left(\\theta \\right)} + 1\\right) T_{1}{\\left(\\theta \\right)}}{T_{1}{\\left(\\theta \\right)} -\\sin{\\left(\\theta \\right)}}$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[36]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">Y<\/span><span class=\"p\">(<\/span><span class=\"n\">T1<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[36]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle H -\\frac{T_{1}^{2}{\\left(\\theta \\right)}}{2} + T_{1}{\\left(\\theta \\right)} \\sin{\\left(\\theta \\right)}$<\/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<p>$L$ \u304c\u6700\u5927\u3068\u306a\u308b\u89d2\u5ea6\u306f\uff0c\u9023\u7acb\u65b9\u7a0b\u5f0f<\/p>\n<p>\\begin{eqnarray}<br \/>\nY(T_1(\\theta), \\theta) &amp;=&amp; H + T_1(\\theta)\\,\\sin\\theta -\\frac{1}{2} T^2_1(\\theta)= 0 \\\\<br \/>\n\\frac{d}{d\\theta} L(\\theta) &amp;=&amp; \\frac{\\left(-T_{1}{\\left(\\theta \\right)} \\sin{\\left(\\theta \\right)} + 1\\right) T_{1}{\\left(\\theta \\right)}}{T_{1}{\\left(\\theta \\right)} -\\sin{\\left(\\theta \\right)}}=0<br \/>\n\\end{eqnarray}<\/p>\n<p>\u3092\uff0c\uff08$\\theta$ \u305d\u306e\u3082\u306e\u3067\u306f\u306a\u304f\uff09$\\sin\\theta$ \u3068 $T_1(\\theta)$ \u306b\u3064\u3044\u3066\u89e3\u3044\u3066\u6c42\u3081\u308b\u3053\u3068\u304c\u3067\u304d\u308b\u3002<\/p>\n<p>$\\displaystyle 0 &lt; \\theta \\leq \\frac{\\pi}{2}$ \u3067\u3042\u308b\u304b\u3089 $\\sin\\theta &gt; 0$ \u3067\u3042\u308b\u3002<\/p>\n<p>$ \\sin\\theta \\Rightarrow s$ \u3068\u7f6e\u304d\u63db\u3048\u3066\u5909\u6570 $s$ \u306b\u5bfe\u3057\u3066 $s &gt; 0$ \u3068\u3044\u3046\u6761\u4ef6\u3092\u3064\u3051\u3066\u89e3\u3044\u3066\u307f\u308b\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[37]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">var<\/span><span class=\"p\">(<\/span><span class=\"s1\">'s'<\/span><span class=\"p\">,<\/span> <span class=\"n\">positive<\/span> <span class=\"o\">=<\/span> <span class=\"kc\">True<\/span><span class=\"p\">)<\/span>\r\n\r\n<span class=\"n\">sols2<\/span> <span class=\"o\">=<\/span> <span class=\"n\">solve<\/span><span class=\"p\">([<\/span><span class=\"n\">Y<\/span><span class=\"p\">(<\/span><span class=\"n\">T1<\/span><span class=\"p\">,<\/span> <span class=\"n\">theta<\/span><span class=\"p\">)<\/span><span class=\"o\">.<\/span><span class=\"n\">subs<\/span><span class=\"p\">(<\/span><span class=\"n\">sin<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">),<\/span> <span class=\"n\">s<\/span><span class=\"p\">),<\/span> \r\n               <span class=\"n\">dL<\/span><span class=\"o\">.<\/span><span class=\"n\">subs<\/span><span class=\"p\">(<\/span><span class=\"n\">sin<\/span><span class=\"p\">(<\/span><span class=\"n\">theta<\/span><span class=\"p\">),<\/span> <span class=\"n\">s<\/span><span class=\"p\">)],<\/span> <span class=\"p\">[<\/span><span class=\"n\">s<\/span><span class=\"p\">,<\/span> <span class=\"n\">T1<\/span><span class=\"p\">])<\/span>\r\n<span class=\"n\">sols2<\/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 output_prompt\">Out[37]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\left[ \\left( \\frac{\\sqrt{2}}{2 \\sqrt{H + 1}}, \\ \\sqrt{2 H + 2}\\right)\\right]$<\/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<p>\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2\u3068\u306a\u308b $\\sin\\theta_{\\rm m}$ \u3068\u305d\u306e\u3068\u304d\u306e\u6ede\u7a7a\u6642\u9593 $T_1$ \u306f\uff0c<br \/>\n$s_m = \\sin\\theta_m = $ <code>sols2[0][0]<\/code>\uff0c$T_1 = $ <code>sols2[0][1]<\/code> \u3067\u3042\u308a\uff0c<\/p>\n<p>\\begin{eqnarray}<br \/>\ns_m = \\sin\\theta_{\\rm m} &amp;=&amp; \\frac{1}{\\sqrt{2 (1+H)}} \\\\<br \/>\n\\therefore\\ \\ \\theta_{\\rm m} &amp;=&amp; \\arcsin \\left( \\frac{1}{\\sqrt{2 (1+H)}}\\right) \\\\<br \/>\nT_1 &amp;=&amp; \\sqrt{2 (1 + H)}<br \/>\n\\end{eqnarray}<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[38]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">thetam<\/span> <span class=\"o\">=<\/span> <span class=\"n\">asin<\/span><span class=\"p\">(<\/span><span class=\"n\">sols2<\/span><span class=\"p\">[<\/span><span class=\"mi\">0<\/span><span class=\"p\">][<\/span><span class=\"mi\">0<\/span><span class=\"p\">])<\/span>\r\n<span class=\"n\">display<\/span><span class=\"p\">(<\/span><span class=\"n\">thetam<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_latex output_subarea \">$\\displaystyle \\operatorname{asin}{\\left(\\frac{\\sqrt{2}}{2 \\sqrt{H + 1}} \\right)}$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[39]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">T1m<\/span> <span class=\"o\">=<\/span> <span class=\"n\">sols2<\/span><span class=\"p\">[<\/span><span class=\"mi\">0<\/span><span class=\"p\">][<\/span><span class=\"mi\">1<\/span><span class=\"p\">]<\/span>\r\n<span class=\"n\">display<\/span><span class=\"p\">(<\/span><span class=\"n\">T1m<\/span><span class=\"p\">)<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt\"><\/div>\n<div class=\"output_latex output_subarea \">$\\displaystyle \\sqrt{2 H + 2}$<\/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<p>\u305d\u306e\u3068\u304d\u306e\u6700\u5927\u6c34\u5e73\u5230\u9054\u8ddd\u96e2 $L_{\\rm m} = X(T_1(\\theta_{\\rm m}), \\theta_{\\rm m})$ \u306f&#8230;<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[40]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-ipython3\">\n<pre><span class=\"n\">X<\/span><span class=\"p\">(<\/span><span class=\"n\">T1m<\/span><span class=\"p\">,<\/span> <span class=\"n\">thetam<\/span><span class=\"p\">)<\/span><span class=\"o\">.<\/span><span class=\"n\">simplify<\/span><span class=\"p\">()<\/span>\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[40]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">$\\displaystyle \\sqrt{2 H + 1}$<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"excerpt":{"rendered":"<p>\u5730\u9762\u304b\u3089\u9ad8\u3055 $h$ \u306e\u5730\u70b9\u304b\u3089\u7a7a\u6c17\u62b5\u6297\u306a\u3057\u306e\u659c\u65b9\u6295\u5c04\u3092\u884c\u3046\u3068\uff0c\u521d\u901f\u5ea6\u306e\u5927\u304d\u3055\u3092\u4e00\u5b9a\u3068\u3057\u305f\u5834\u5408\uff0c\u6c34\u5e73\u65b9\u5411\u306e\u5230\u9054\u8ddd\u96e2\u304c\u6700\u5927\u3068\u306a\u308b\u306e\u306f\u6253\u3061\u4e0a\u3052\u89d2\u5ea6\uff08\u4ef0\u89d2\uff09\u304c\u4f55\u5ea6\u306e\u3068\u304d\u304b\u3092 Python \u306e SymPy \u3092\u4f7f\u3063\u3066\uff08\u304b\u3064 SciPy \u3084 NumPy \u3092\u4f7f\u308f\u305a\u306b\uff09\u6c42\u3081\u308b\u3002<\/p>\n<p style=\"text-align: right;\">update: 2025.01.31<\/p><p><a class=\"more-link btn\" href=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/%e3%82%b3%e3%83%b3%e3%83%94%e3%83%a5%e3%83%bc%e3%82%bf%e6%bc%94%e7%bf%92\/%e9%ab%98%e3%81%95-h-%e3%81%8b%e3%82%89%e3%81%ae%e6%96%9c%e6%96%b9%e6%8a%95%e5%b0%84%e3%81%ae%e6%9c%80%e5%a4%a7%e5%88%b0%e9%81%94%e8%b7%9d%e9%9b%a2%e3%82%92%e6%b1%82%e3%82%81%e3%82%8b%e6%ba%96\/sympy-%e3%81%a7%e9%ab%98%e3%81%95-h-%e3%81%8b%e3%82%89%e3%81%ae%e6%96%9c%e6%96%b9%e6%8a%95%e5%b0%84%e3%81%ae%e6%9c%80%e5%a4%a7%e5%88%b0%e9%81%94%e8%b7%9d%e9%9b%a2%e3%82%92%e6%b1%82%e3%82%81%e3%82%8b\/\">\u7d9a\u304d\u3092\u8aad\u3080<\/a><\/p>\n","protected":false},"author":33,"featured_media":0,"parent":5052,"menu_order":10,"comment_status":"closed","ping_status":"closed","template":"","meta":{"inline_featured_image":false,"footnotes":""},"class_list":["post-7371","page","type-page","status-publish","hentry","nodate","item-wrap"],"aioseo_notices":[],"_links":{"self":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/7371","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/types\/page"}],"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=7371"}],"version-history":[{"count":24,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/7371\/revisions"}],"predecessor-version":[{"id":10166,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/7371\/revisions\/10166"}],"up":[{"embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/5052"}],"wp:attachment":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/media?parent=7371"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}