{"id":4950,"date":"2023-12-18T10:39:50","date_gmt":"2023-12-18T01:39:50","guid":{"rendered":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/?page_id=4950"},"modified":"2023-12-18T10:41:01","modified_gmt":"2023-12-18T01:41:01","slug":"maxima-%e3%81%a7%e9%81%8b%e5%8b%95%e6%96%b9%e7%a8%8b%e5%bc%8f%e3%82%92%e8%a7%a3%e3%81%84%e3%81%a6%e6%8c%af%e5%b9%85%e3%81%a8%e5%91%a8%e6%9c%9f%e3%81%ae%e9%96%a2%e4%bf%82%e3%82%92%e8%aa%bf%e3%81%b9","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\/%e5%8f%82%e8%80%83%ef%bc%9amaxima-%e7%b7%a8\/maxima-%e3%81%a7%e9%81%8b%e5%8b%95%e6%96%b9%e7%a8%8b%e5%bc%8f%e3%82%92%e8%a7%a3%e3%81%84%e3%81%a6%e6%8c%af%e5%b9%85%e3%81%a8%e5%91%a8%e6%9c%9f%e3%81%ae%e9%96%a2%e4%bf%82%e3%82%92%e8%aa%bf%e3%81%b9\/","title":{"rendered":"Maxima \u3067\u904b\u52d5\u65b9\u7a0b\u5f0f\u3092\u89e3\u3044\u3066\u632f\u5e45\u3068\u5468\u671f\u306e\u95a2\u4fc2\u3092\u8abf\u3079\u308b"},"content":{"rendered":"<p>Maxima \u3092\u4f7f\u3063\u3066\u904b\u52d5\u65b9\u7a0b\u5f0f\u3092\u6570\u5024\u7684\u306b\u89e3\u304d\uff0c\u5358\u632f\u308a\u5b50\u306e\u632f\u5e45\u3068\u5468\u671f\u306e\u95a2\u4fc2\u3092\u8abf\u3079\u308b\u3002\u5c0e\u51fa\u306b\u3064\u3044\u3066\u306f\uff0c\u4ee5\u4e0b\u3092\u53c2\u7167\u3002<\/p>\n<ul>\n<li>\u300c<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\/%e5%8d%98%e6%8c%af%e3%82%8a%e5%ad%90%e3%81%ae%e9%81%8b%e5%8b%95%e6%96%b9%e7%a8%8b%e5%bc%8f%e3%82%92%e6%95%b0%e5%80%a4%e7%9a%84%e3%81%ab%e8%a7%a3%e3%81%8f%e6%ba%96%e5%82%99\/\">\u5358\u632f\u308a\u5b50\u306e\u904b\u52d5\u65b9\u7a0b\u5f0f\u3092\u6570\u5024\u7684\u306b\u89e3\u304f\u6e96\u5099<\/a>\u300d<\/li>\n<\/ul>\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=\"\u4e00\u69d8\u91cd\u529b\u5834\u4e2d\u306e\u5358\u632f\u308a\u5b50\">\u4e00\u69d8\u91cd\u529b\u5834\u4e2d\u306e\u5358\u632f\u308a\u5b50<\/h3>\n<p>\u5358\u632f\u308a\u5b50\u306e\u3072\u3082\u306e\u9577\u3055\u3092 $\\ell$\uff0c\u91cd\u529b\u52a0\u901f\u5ea6\u306e\u5927\u304d\u3055\u3092 $g$\uff0c\u925b\u76f4\u4e0b\u5411\u304d\u304b\u3089\u306e\u632f\u308c\u89d2\u3092 $\\theta$ \u3068\u3059\u308b\u3068\uff0c\u5358\u632f\u308a\u5b50\u306e\u904b\u52d5\u65b9\u7a0b\u5f0f\u306f<\/p>\n<p>$$\\frac{d^2\\theta}{dt^2} = &#8211; \\frac{g}{\\ell} \\sin\\theta$$<\/p>\n<p>\u3068\u306a\u308b\u3002<\/p>\n<h3 id=\"\u632f\u5e45\u304c\u5c0f\u3055\u3044\u5834\u5408\u306b\u306f\u5358\u632f\u52d5\u3068\u306a\u308b\u3053\u3068\">\u632f\u5e45\u304c\u5c0f\u3055\u3044\u5834\u5408\u306b\u306f\u5358\u632f\u52d5\u3068\u306a\u308b\u3053\u3068<\/h3>\n<p>$|\\theta| \\ll 1$ \u306e\u5834\u5408\u306b\u306f\uff0c$\\sin\\theta \\simeq \\theta$ \u3068\u8fd1\u4f3c\u3059\u308b\u3053\u3068\u3067\uff0c\u904b\u52d5\u65b9\u7a0b\u5f0f\u306f<\/p>\n<p>$$\\frac{d^2\\theta}{dt^2} = &#8211; \\frac{g}{\\ell} \\theta$$<\/p>\n<p>\u3068\u306a\u308a\uff0c\u3053\u308c\u306f\u5358\u632f\u52d5\u3068\u306a\u308b\u3002\u5358\u632f\u52d5\u306e\u5468\u671f $\\tau_0$ \u306f<\/p>\n<p>$$\\tau_0 = 2 \\pi \\sqrt{\\frac{\\ell}{g}}$$<\/p>\n<p>\u3068\u306a\u308a\uff0c\u5358\u632f\u308a\u5b50\u306e\u3072\u3082\u306e\u9577\u3055\u3068\u91cd\u529b\u52a0\u901f\u5ea6\u3060\u3051\u3067\u6c7a\u307e\u308a\uff0c\u521d\u671f\u6761\u4ef6\u3067\u4e0e\u3048\u3089\u308c\u308b\u632f\u5e45 $\\theta_0$ \u306b\u4f9d\u5b58\u3057\u306a\u3044\u3002\u3053\u308c\u3092\u300c\u632f\u308a\u5b50\u306e\u7b49\u6642\u6027\u300d\u3068\u8a00\u3063\u305f\u308a\u3059\u308b\u3002<\/p>\n<h3 id=\"\u7121\u6b21\u5143\u5316\">\u7121\u6b21\u5143\u5316<\/h3>\n<p>\u5358\u632f\u52d5\u306e\u5468\u671f $\\tau_0$ \u3067\u7121\u6b21\u5143\u5316\u3057\u305f\u6642\u9593 $T$ \u3092<\/p>\n<p>$$T \\equiv \\frac{t}{\\tau_0}$$<\/p>\n<p>\u3068\u5b9a\u7fa9\u3059\u308b\u3002\u3053\u306e\u7121\u6b21\u5143\u5316\u3055\u308c\u305f\u6642\u9593\u5ea7\u6a19\u3067\u66f8\u3044\u305f\u904b\u52d5\u65b9\u7a0b\u5f0f\u306f<\/p>\n<p>$$\\frac{d^2 \\theta}{dT^2} = &#8211; 4 \\pi^2 \\sin\\theta$$<\/p>\n<p>\u3068\u306a\u308b\u3002<\/p>\n<h3 id=\"\u554f\u984c\uff1a\u632f\u5e45\u304c\u5927\u304d\u3044\u5834\u5408\u306e\u5468\u671f\u306f\uff1f\">\u554f\u984c\uff1a\u632f\u5e45\u304c\u5927\u304d\u3044\u5834\u5408\u306e\u5468\u671f\u306f\uff1f<\/h3>\n<p>\u632f\u5e45\u304c\u5927\u304d\u3044\uff0c\u3064\u307e\u308a $\\sin\\theta \\simeq \\theta$ \u3068\u8fd1\u4f3c\u3067\u304d\u306a\u3044\u5834\u5408\u306f\uff0c\u632f\u308a\u5b50\u306e\u5468\u671f\u306f\u4e00\u822c\u306b\u632f\u5e45\u306b\u4f9d\u5b58\u3059\u308b\u306e\u3067\u306f\u306a\u3044\u304b\u3068\u8003\u3048\u3089\u308c\u308b\u3002<\/p>\n<p>\u305d\u3053\u3067\uff0c$T = 0$ \u3067\u306e\u632f\u5e45 $\\theta_0$ \u3092\u305f\u3068\u3048\u3070 10\u00b0 \u304b\u3089 90\u00b0 \u307e\u3067 10\u00b0 \u304d\u3056\u307f\u3067\u5927\u304d\u304f\u3057\u3066\u3044\u3063\u305f\u5834\u5408\uff0c\u3075\u308a\u3053\u306e\u5468\u671f\u306f\u3069\u3046\u306a\u308b\u304b\uff1f\u3000\u3068\u3044\u3046\u306e\u304c\u554f\u984c\u3002<\/p>\n<h3 id=\"\u30eb\u30f3\u30b2\u30fb\u30af\u30c3\u30bf\u6cd5\u3067\u89e3\u304f\u305f\u3081\u306e\u6e96\u5099\">\u30eb\u30f3\u30b2\u30fb\u30af\u30c3\u30bf\u6cd5\u3067\u89e3\u304f\u305f\u3081\u306e\u6e96\u5099<\/h3>\n<p>2\u968e\u5e38\u5fae\u5206\u65b9\u7a0b\u5f0f\u3092\u30eb\u30f3\u30b2\u30fb\u30af\u30c3\u30bf\u6cd5\u3067\u6570\u5024\u7684\u306b\u89e3\u304f\u305f\u3081\u306b\u306f\uff0c\u9023\u7acb1\u968e\u5fae\u5206\u65b9\u7a0b\u5f0f\u7cfb\u306b\u306a\u304a\u3059\u3002<\/p>\n<p>\\begin{eqnarray}<br \/>\n\\frac{d\\theta}{dT} &amp;=&amp; F_1(V) = V \\\\<br \/>\n\\frac{dV}{dT} &amp;=&amp; F_2(\\theta) = &#8211; 4 \\pi^2 \\sin\\theta<br \/>\n\\end{eqnarray}<\/p>\n<p>\u3053\u308c\u3092\u521d\u671f\u6761\u4ef6 $T = 0$ \u3067<\/p>\n<p>\\begin{eqnarray}<br \/>\n\\theta(0) &amp;=&amp; \\theta_0 \\\\<br \/>\nV(0) &amp;=&amp; 0<br \/>\n\\end{eqnarray}<\/p>\n<p>\u3068\u3057\uff0c$\\theta_0 = 10^{\\circ}, \\ 60^{\\circ}, \\ 80^{\\circ}$ \u306e\u5834\u5408\u306b $T_0 = 0$ \u304b\u3089 $ T_1 = 2$ \u307e\u3067\u89e3\u304f\u3002<\/p>\n<p>\uff08\u304f\u3069\u3044\u3088\u3046\u3067\u3059\u304c\uff0c\u898f\u683c\u5316\u3055\u308c\u305f\u6642\u9593\u3067 $T_1 = 2$ \u3068\u3044\u3046\u306e\u306f\uff0c\u5358\u632f\u52d5\u306e\u5834\u5408\u306e\u5468\u671f $\\tau_0$ \u306e 2 \u500d\u306e\u6642\u9593\u307e\u3067\uff0c\u3068\u3044\u3046\u610f\u5473\u3067\u3059\u3088\u3002\uff09<\/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[1]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"cm\">\/* \u9023\u7acb\u5fae\u5206\u65b9\u7a0b\u5f0f\u306e\u53f3\u8fba *\/<\/span>\r\n<span class=\"nf\">F1<\/span><span class=\"p\">(<\/span><span class=\"nv\">V<\/span><span class=\"p\">)<\/span><span class=\"o\">:=<\/span> <span class=\"nv\">V<\/span>;\r\n<span class=\"nf\">F2<\/span><span class=\"p\">(<\/span><span class=\"nv\">theta<\/span><span class=\"p\">)<\/span><span class=\"o\">:=<\/span> <span class=\"o\">-<\/span> <span class=\"mi\">4<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">%pi<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">sin<\/span><span class=\"p\">(<\/span><span class=\"nv\">theta<\/span><span class=\"p\">)<\/span>; \r\n\r\n<span class=\"cm\">\/* \u5358\u632f\u52d5\u306e\u5834\u5408 <\/span>\r\n<span class=\"cm\">F2(theta):= - 4 * %pi**2 * theta; *\/<\/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[1]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{1}$}F_{1}\\left(V\\right):=V\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[1]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{2}$}F_{2}\\left(\\vartheta\\right):=\\left(-4\\right)\\,\\pi^2\\,\\sin \\vartheta\\]<\/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>\u307e\u305a\uff0c$\\theta_0 = 80^{\\circ}$ \u306e\u5834\u5408\u306b\uff0c\u523b\u307f\u5e45 $h$ \u3068\u6570\u5024\u8a08\u7b97\u306e\u7cbe\u5ea6\u306b\u3064\u3044\u3066\u78ba\u8a8d\u3057\u3066\u304a\u304f\u3002<\/p>\n<p>Maxima \u306e Runge-Kutta \u6cd5\u306e\u95a2\u6570 <code>rk()<\/code> \u3067\u306f\uff0c\u6570\u5024\u8a08\u7b97\u306e\u7cbe\u5ea6\u306b\u95a2\u4fc2\u3059\u308b\u306e\u306f\u523b\u307f\u5e45 $h$ \u3060\u3051\u3067\u3042\u308b\u304b\u3089\u306d\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-maxima\">\n<pre><span class=\"cm\">\/* \u5ea6\u3092\u30e9\u30b8\u30a2\u30f3\u306b\u5909\u63db\u3059\u308b\u95a2\u6570\u3092\u6e96\u5099\u3057\u3066\u304a\u304f *\/<\/span>\r\n<span class=\"nf\">radians<\/span><span class=\"p\">(<\/span><span class=\"nv\">deg<\/span><span class=\"p\">)<\/span><span class=\"o\">:=<\/span> <span class=\"nv\">deg<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">float<\/span><span class=\"p\">(<\/span><span class=\"nv\">%pi<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span>180$\r\n\r\n<span class=\"cm\">\/* \u521d\u671f\u6761\u4ef6\u3068\u8a08\u7b97\u7bc4\u56f2 *\/<\/span>\r\n<span class=\"nv\">T0<\/span><span class=\"o\">:<\/span> 0$\r\n<span class=\"nv\">th<\/span><span class=\"o\">:<\/span> 80$\r\n<span class=\"nv\">theta0<\/span><span class=\"o\">:<\/span> <span class=\"nf\">radians<\/span><span class=\"p\">(<\/span><span class=\"nv\">th<\/span><span class=\"p\">)<\/span>$ <span class=\"cm\">\/* \u30e9\u30b8\u30a2\u30f3\u306b\u306a\u304a\u3059 *\/<\/span>\r\n<span class=\"nv\">V0<\/span><span class=\"o\">:<\/span> 0$\r\n<span class=\"nv\">T1<\/span><span class=\"o\">:<\/span> 2$\r\n\r\n<span class=\"cm\">\/* N \u3092\u5909\u3048\u3066\u6570\u5024\u8a08\u7b97... *\/<\/span>\r\n<span class=\"k\">for<\/span> <span class=\"nv\">N<\/span> <span class=\"nv\">in<\/span> <span class=\"p\">[<\/span><span class=\"mi\">200<\/span>, <span class=\"mi\">2000<\/span>, 20000<span class=\"p\">]<\/span> <span class=\"k\">do<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"cm\">\/* \u523b\u307f\u5e45 *\/<\/span>\r\n  <span class=\"nv\">h<\/span><span class=\"o\">:<\/span> <span class=\"nf\">float<\/span><span class=\"p\">((<\/span><span class=\"nv\">T1<\/span> <span class=\"o\">-<\/span> <span class=\"nv\">T0<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"nv\">N<\/span><span class=\"p\">)<\/span>, \r\n  \r\n  <span class=\"cm\">\/* Runge-Kutta \u6cd5\u3067\u6570\u5024\u8a08\u7b97 *\/<\/span>\r\n  <span class=\"nv\">ansN<\/span><span class=\"o\">:<\/span> <span class=\"nf\">rk<\/span><span class=\"p\">([<\/span><span class=\"nf\">F1<\/span><span class=\"p\">(<\/span><span class=\"nv\">V<\/span><span class=\"p\">)<\/span>, <span class=\"nf\">F2<\/span><span class=\"p\">(<\/span><span class=\"nv\">theta<\/span><span class=\"p\">)]<\/span>, \r\n           <span class=\"p\">[<\/span><span class=\"nv\">theta<\/span>, <span class=\"nv\">V<\/span><span class=\"p\">]<\/span>, <span class=\"p\">[<\/span><span class=\"nv\">theta0<\/span>, <span class=\"nv\">V0<\/span><span class=\"p\">]<\/span>, \r\n           <span class=\"p\">[<\/span><span class=\"nv\">T<\/span>, <span class=\"nv\">T0<\/span>, <span class=\"nv\">T1<\/span>, <span class=\"nv\">h<\/span><span class=\"p\">])<\/span>,\r\n           \r\n  <span class=\"cm\">\/* \u6700\u5f8c\u306e\u5024\u3092\u8868\u793a\u3057\uff0c\u8aa4\u5dee\u3092\u78ba\u8a8d\u3002 *\/<\/span>\r\n  <span class=\"nf\">printf<\/span><span class=\"p\">(<\/span><span class=\"no\">true<\/span>, <span class=\"s\">\"N = ~6d, h = ~,8f, V(T=2) = ~15,13f~%\"<\/span>, \r\n                <span class=\"nv\">N<\/span>, <span class=\"nv\">h<\/span>, <span class=\"nv\">ansN<\/span><span class=\"p\">[<\/span><span class=\"nf\">length<\/span><span class=\"p\">(<\/span><span class=\"nv\">ansN<\/span><span class=\"p\">)][<\/span>3<span class=\"p\">])<\/span>\r\n<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>N =    200, h = 0.01000000, V(T=2) = 8.0634656639326\r\nN =   2000, h = 0.00100000, V(T=2) = 8.0634655783821\r\nN =  20000, h = 0.00010000, V(T=2) = 8.0634655783315\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<p>\u4e0a\u306e\u7d50\u679c\u304b\u3089 $N=2000$ \u3067\u5c0f\u6570\u70b9\u4ee5\u4e0b10\u6841\u7a0b\u5ea6\u306e\u7cbe\u5ea6\u304c\u3042\u308b\u3068\u308f\u304b\u308b\u3002\u3053\u3053\u3067\u306f $N=20000$ \u3067\u8a08\u7b97\u3092\u884c\u3046\u3053\u3068\u306b\u3057\u3066\u307f\u308b\u3002\u30b0\u30e9\u30d5\u3092\u63cf\u304f\u3060\u3051\u306a\u3089 $N$ \u3092\u305d\u3093\u306a\u306b\u5927\u304d\u304f\u3068\u308b\u5fc5\u8981\u306f\u3042\u308a\u307e\u305b\u3093\u304c\uff0c\u3042\u3068\u3067\u7406\u7531\u304c\u308f\u304b\u308a\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[3]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"cm\">\/* \u521d\u671f\u6761\u4ef6\u3068\u8a08\u7b97\u7bc4\u56f2 *\/<\/span>\r\n<span class=\"nv\">T0<\/span><span class=\"o\">:<\/span> 0$\r\n<span class=\"nv\">V0<\/span><span class=\"o\">:<\/span> 0$\r\n<span class=\"nv\">T1<\/span><span class=\"o\">:<\/span> 2$\r\n<span class=\"cm\">\/* N \u3092\u5927\u304d\u3081\u306b\u3068\u308b\u7406\u7531\u306f\uff0c\u3042\u3068\u3067\u308f\u304b\u308b *\/<\/span>\r\n<span class=\"nv\">N<\/span><span class=\"o\">:<\/span> 20000$\r\n<span class=\"cm\">\/* \u523b\u307f\u5e45 *\/<\/span>\r\n<span class=\"nv\">h<\/span><span class=\"o\">:<\/span> <span class=\"nf\">float<\/span><span class=\"p\">((<\/span><span class=\"nv\">T1<\/span> <span class=\"o\">-<\/span> <span class=\"nv\">T0<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"nv\">N<\/span><span class=\"p\">)<\/span>$ \r\n\r\n<span class=\"cm\">\/* \u632f\u5e45\u306e\u521d\u671f\u5024\u3092\u5909\u3048\u3066\u6570\u5024\u8a08\u7b97... *\/<\/span>\r\n<span class=\"k\">for<\/span> <span class=\"nv\">th<\/span> <span class=\"nv\">in<\/span> <span class=\"p\">[<\/span><span class=\"mi\">10<\/span>, <span class=\"mi\">60<\/span>, 80<span class=\"p\">]<\/span> <span class=\"k\">do<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"nv\">theta0<\/span><span class=\"o\">:<\/span> <span class=\"nf\">radians<\/span><span class=\"p\">(<\/span><span class=\"nv\">th<\/span><span class=\"p\">)<\/span>,\r\n  <span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">]<\/span><span class=\"o\">:<\/span> <span class=\"nf\">rk<\/span><span class=\"p\">([<\/span><span class=\"nf\">F1<\/span><span class=\"p\">(<\/span><span class=\"nv\">V<\/span><span class=\"p\">)<\/span>, <span class=\"nf\">F2<\/span><span class=\"p\">(<\/span><span class=\"nv\">theta<\/span><span class=\"p\">)]<\/span>, \r\n                 <span class=\"p\">[<\/span><span class=\"nv\">theta<\/span>, <span class=\"nv\">V<\/span><span class=\"p\">]<\/span>, <span class=\"p\">[<\/span><span class=\"nv\">theta0<\/span>, <span class=\"nv\">V0<\/span><span class=\"p\">]<\/span>, \r\n                 <span class=\"p\">[<\/span><span class=\"nv\">T<\/span>, <span class=\"nv\">T0<\/span>, <span class=\"nv\">T1<\/span>, <span class=\"nv\">h<\/span><span class=\"p\">])<\/span>\r\n<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<p>\u6570\u5024\u8a08\u7b97\u3067\u304d\u305f\u3089\uff0c\u7d50\u679c\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[4]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"cm\">\/* \u3042\u3089\u304b\u3058\u3081\u30b0\u30e9\u30d5\u306b\u3059\u308b points \u306e\u30ea\u30b9\u30c8\u3092\u4f5c\u6210\u3057\u3066\u304a\u304f *\/<\/span>\r\n\r\n<span class=\"nv\">i<\/span><span class=\"o\">:<\/span> 0$\r\n<span class=\"k\">for<\/span> <span class=\"nv\">th<\/span> <span class=\"nv\">in<\/span> <span class=\"p\">[<\/span><span class=\"mi\">10<\/span>, <span class=\"mi\">60<\/span>, 80<span class=\"p\">]<\/span> <span class=\"k\">do<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"cm\">\/* \u7dda\u306e\u8272 *\/<\/span>\r\n  <span class=\"nv\">i<\/span><span class=\"o\">:<\/span> <span class=\"nv\">i<\/span> <span class=\"o\">+<\/span> <span class=\"mi\">1<\/span>,\r\n  <span class=\"nv\">col<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">]<\/span><span class=\"o\">:<\/span> <span class=\"nv\">i<\/span>,\r\n\r\n  <span class=\"cm\">\/* \u51e1\u4f8b\u306b\u632f\u5e45\u306e\u521d\u671f\u5024\u3092\u3044\u308c\u308b *\/<\/span>\r\n  <span class=\"nv\">legend<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">]<\/span><span class=\"o\">:<\/span> <span class=\"nf\">concat<\/span><span class=\"p\">(<\/span><span class=\"s\">\"\u03b8_0 = \"<\/span>, <span class=\"nv\">th<\/span>, <span class=\"s\">\"\u00b0\"<\/span><span class=\"p\">)<\/span>,\r\n\r\n  <span class=\"cm\">\/* \u305f\u304f\u3055\u3093\u70b9\u3092\u63cf\u304f\u3068\u6642\u9593\u304c\u304b\u304b\u308b\u306e\u3067<\/span>\r\n<span class=\"cm\">     \u30d7\u30ed\u30c3\u30c8\u3059\u308b\u30c7\u30fc\u30bf\u3092\u5c11\u3057\u9593\u5f15\u304f *\/<\/span>\r\n  <span class=\"nv\">pts<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">]<\/span><span class=\"o\">:<\/span> <span class=\"nf\">makelist<\/span><span class=\"p\">([<\/span><span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">][<\/span><span class=\"nv\">i<\/span><span class=\"p\">][<\/span>1<span class=\"p\">]<\/span>, <span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">][<\/span><span class=\"nv\">i<\/span><span class=\"p\">][<\/span>2<span class=\"p\">]]<\/span>, <span class=\"nv\">i<\/span>, <span class=\"mi\">1<\/span>, <span class=\"nv\">N<\/span><span class=\"o\">+<\/span><span class=\"mi\">1<\/span>, <span class=\"mi\">100<\/span><span class=\"p\">)<\/span>\r\n<span class=\"p\">)<\/span>$\r\n\r\n<span class=\"nv\">tmp<\/span><span class=\"o\">:<\/span> <span class=\"nf\">makelist<\/span><span class=\"p\">([<\/span><span class=\"nv\">key<\/span> <span class=\"o\">=<\/span><span class=\"nv\">legend<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">]<\/span>, \r\n               <span class=\"nv\">color<\/span> <span class=\"o\">=<\/span> <span class=\"nv\">col<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">]<\/span>, \r\n               <span class=\"nf\">points<\/span><span class=\"p\">(<\/span><span class=\"nv\">pts<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">])]<\/span>, <span class=\"nv\">th<\/span>, <span class=\"p\">[<\/span><span class=\"mi\">80<\/span>, <span class=\"mi\">60<\/span>, 10<span class=\"p\">])<\/span>$\r\n<span class=\"nf\">draw2d<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"nv\">font<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"Arial\"<\/span>, <span class=\"nv\">font_size<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">16<\/span>, \r\n  <span class=\"nv\">title<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u5358\u632f\u308a\u5b50\u306e\u632f\u5e45\u3068\u5468\u671f\"<\/span>,\r\n  <span class=\"nv\">user_preamble<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"set key samplen 1\"<\/span>,\r\n  <span class=\"cm\">\/* \u8868\u793a\u7bc4\u56f2 *\/<\/span>\r\n  <span class=\"nv\">xrange<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"mi\">0<\/span>, 2<span class=\"p\">]<\/span>, <span class=\"nv\">yrange<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"o\">-<\/span><span class=\"mf\">1.5<\/span>, 1<span class=\"o\">.<\/span>5<span class=\"p\">]<\/span>, \r\n  <span class=\"cm\">\/* \u30b0\u30ea\u30c3\u30c9  *\/<\/span>\r\n  <span class=\"nv\">grid<\/span> <span class=\"o\">=<\/span> <span class=\"no\">true<\/span>, \r\n  <span class=\"nv\">xaxis<\/span> <span class=\"o\">=<\/span> <span class=\"no\">true<\/span>, \r\n  <span class=\"cm\">\/* \u5ea7\u6a19\u8ef8\u306e\u30e9\u30d9\u30eb *\/<\/span>\r\n  <span class=\"nv\">xlabel<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u898f\u683c\u5316\u3055\u308c\u305f\u6642\u9593 T\"<\/span>, <span class=\"nv\">ylabel<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u632f\u5e45 \u03b8 (rad)\"<\/span>, \r\n      \r\n  <span class=\"nv\">point_type<\/span> <span class=\"o\">=<\/span> <span class=\"nv\">dot<\/span>, \r\n  <span class=\"nv\">points_joined<\/span> <span class=\"o\">=<\/span> <span class=\"no\">true<\/span>, \r\n  <span class=\"nv\">line_width<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">1.5<\/span>,\r\n  <span class=\"nv\">tmp<\/span>  \r\n<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 \">\n<p><!--?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?--><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-5009\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/mfuri01A.svg\" alt=\"\" width=\"640\" height=\"533\" \/><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<p>\u7e26\u8ef8\u306e\u632f\u5e45\u3092\u30e9\u30b8\u30a2\u30f3\u304b\u3089\u5ea6\u306b\u5909\u3048\u305f\u30b0\u30e9\u30d5\u3092\u63cf\u304d\u307e\u3059\u3002<\/p>\n<p>\u306a\u304a\uff0c\u4ee5\u4e0b\u3067\u306f\u76f4\u63a5\u4f7f\u3044\u307e\u305b\u3093\u304c\uff0c\u53c2\u8003\u307e\u3067\u306b Maxima \u3067\u30ea\u30b9\u30c8\u5185\u306e\u5168\u3066\u306e\u8981\u7d20\u30922\u500d\u3059\u308b\u4f8b\u306f\u4ee5\u4e0b\u306e\u901a\u308a\u3067\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[5]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">a<\/span><span class=\"o\">:<\/span> <span class=\"p\">[<\/span><span class=\"mi\">1<\/span>, <span class=\"mi\">2<\/span>, 3<span class=\"p\">]<\/span>;\r\n<span class=\"nv\">a<\/span><span class=\"o\">*<\/span><span class=\"mi\">2<\/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\">\\[\\tag{${\\it \\%o}_{20}$}\\left[ 1 , 2 , 3 \\right] \\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[5]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{21}$}\\left[ 2 , 4 , 6 \\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[6]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"cm\">\/* \u30e9\u30b8\u30a2\u30f3\u304b\u3089\u5ea6\u3078\u5909\u63db\u3059\u308b\u95a2\u6570\u3092\u6e96\u5099\u3057\u3066\u304a\u304f *\/<\/span>\r\n<span class=\"nf\">degrees<\/span><span class=\"p\">(<\/span><span class=\"nv\">rad<\/span><span class=\"p\">)<\/span><span class=\"o\">:=<\/span> <span class=\"nv\">rad<\/span> <span class=\"o\">*<\/span> 180<span class=\"o\">\/<\/span><span class=\"nf\">float<\/span><span class=\"p\">(<\/span><span class=\"nv\">%pi<\/span><span class=\"p\">)<\/span>$\r\n\r\n<span class=\"cm\">\/* \u3042\u3089\u304b\u3058\u3081\u30b0\u30e9\u30d5\u306b\u3059\u308b points \u306e\u30ea\u30b9\u30c8\u3092\u4f5c\u6210\u3057\u3066\u304a\u304f *\/<\/span>\r\n\r\n<span class=\"nv\">i<\/span><span class=\"o\">:<\/span> 0$\r\n<span class=\"k\">for<\/span> <span class=\"nv\">th<\/span> <span class=\"nv\">in<\/span> <span class=\"p\">[<\/span><span class=\"mi\">10<\/span>, <span class=\"mi\">60<\/span>, 80<span class=\"p\">]<\/span> <span class=\"k\">do<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"cm\">\/* \u7dda\u306e\u8272 *\/<\/span>\r\n  <span class=\"nv\">i<\/span><span class=\"o\">:<\/span> <span class=\"nv\">i<\/span> <span class=\"o\">+<\/span> <span class=\"mi\">1<\/span>,\r\n  <span class=\"nv\">col<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">]<\/span><span class=\"o\">:<\/span> <span class=\"nv\">i<\/span>,\r\n  <span class=\"cm\">\/* \u51e1\u4f8b *\/<\/span>\r\n  <span class=\"nv\">legend<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">]<\/span><span class=\"o\">:<\/span> <span class=\"nf\">concat<\/span><span class=\"p\">(<\/span><span class=\"s\">\"\u03b8_0 = \"<\/span>, <span class=\"nv\">th<\/span>, <span class=\"s\">\"\u00b0\"<\/span><span class=\"p\">)<\/span>,\r\n  <span class=\"cm\">\/* \u30d7\u30ed\u30c3\u30c8\u3059\u308b\u30c7\u30fc\u30bf\u3092\u5c11\u3057\u9593\u5f15\u304d\uff0c\u7e26\u8ef8\u3092\u30e9\u30b8\u30a2\u30f3\u304b\u3089\u5ea6\u306b\u3002 *\/<\/span>\r\n  <span class=\"nv\">pts<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">]<\/span><span class=\"o\">:<\/span> <span class=\"nf\">makelist<\/span><span class=\"p\">([<\/span><span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">][<\/span><span class=\"nv\">i<\/span><span class=\"p\">][<\/span>1<span class=\"p\">]<\/span>, \r\n                     <span class=\"nf\">degrees<\/span><span class=\"p\">(<\/span><span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">][<\/span><span class=\"nv\">i<\/span><span class=\"p\">][<\/span>2<span class=\"p\">])]<\/span>, <span class=\"nv\">i<\/span>, <span class=\"mi\">1<\/span>, <span class=\"nv\">N<\/span><span class=\"o\">+<\/span><span class=\"mi\">1<\/span>, <span class=\"mi\">100<\/span><span class=\"p\">)<\/span> \r\n<span class=\"p\">)<\/span>$\r\n\r\n<span class=\"nv\">tmp<\/span><span class=\"o\">:<\/span> <span class=\"nf\">makelist<\/span><span class=\"p\">([<\/span><span class=\"nv\">key<\/span> <span class=\"o\">=<\/span><span class=\"nv\">legend<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">]<\/span>, \r\n               <span class=\"nv\">color<\/span> <span class=\"o\">=<\/span> <span class=\"nv\">col<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">]<\/span>, \r\n               <span class=\"nf\">points<\/span><span class=\"p\">(<\/span><span class=\"nv\">pts<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">])]<\/span>, <span class=\"nv\">th<\/span>, <span class=\"p\">[<\/span><span class=\"mi\">80<\/span>,<span class=\"mi\">60<\/span>,10<span class=\"p\">])<\/span>$\r\n<span class=\"nf\">draw2d<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"nv\">font<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"Arial\"<\/span>, <span class=\"nv\">font_size<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">16<\/span>, \r\n  <span class=\"nv\">title<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u5358\u632f\u308a\u5b50\u306e\u632f\u5e45\u3068\u5468\u671f\"<\/span>,\r\n  <span class=\"nv\">user_preamble<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"set key samplen 1\"<\/span>,\r\n  <span class=\"cm\">\/* \u8868\u793a\u7bc4\u56f2 *\/<\/span>\r\n  <span class=\"nv\">xrange<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"mi\">0<\/span>, 2<span class=\"p\">]<\/span>, <span class=\"nv\">yrange<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"o\">-<\/span><span class=\"mi\">90<\/span>, 90<span class=\"p\">]<\/span>, \r\n  <span class=\"cm\">\/* \u30b0\u30ea\u30c3\u30c9  *\/<\/span>\r\n  <span class=\"nv\">grid<\/span> <span class=\"o\">=<\/span> <span class=\"no\">true<\/span>, \r\n  <span class=\"nv\">xaxis<\/span> <span class=\"o\">=<\/span> <span class=\"no\">true<\/span>, \r\n  <span class=\"cm\">\/* \u5ea7\u6a19\u8ef8\u306e\u30e9\u30d9\u30eb *\/<\/span>\r\n  <span class=\"nv\">xlabel<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u898f\u683c\u5316\u3055\u308c\u305f\u6642\u9593 T\"<\/span>, <span class=\"nv\">ylabel<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u632f\u5e45 \u03b8 (\u00b0)\"<\/span>, \r\n      \r\n  <span class=\"nv\">point_type<\/span> <span class=\"o\">=<\/span> <span class=\"nv\">dot<\/span>, \r\n  <span class=\"nv\">points_joined<\/span> <span class=\"o\">=<\/span> <span class=\"no\">true<\/span>, \r\n  <span class=\"nv\">line_width<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">1.5<\/span>,\r\n  <span class=\"nv\">tmp<\/span>  \r\n<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 \">\n<p><!--?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?--><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-5010\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/mfuri02A.svg\" alt=\"\" width=\"640\" height=\"533\" \/><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h3 id=\"\u5468\u671f\u306f\u901f\u5ea6\u306e\u5411\u304d\u304c\u5909\u308f\u308b\u6642\u523b\u304b\u3089...\">\u5468\u671f\u306f\u901f\u5ea6\u306e\u5411\u304d\u304c\u5909\u308f\u308b\u6642\u523b\u304b\u3089&#8230;<\/h3>\n<p>\u521d\u901f\u5ea6 $V(0) = 0$ \u3067\u306f\u3058\u307e\u308b\u632f\u308a\u5b50\u306e\u904b\u52d5\u306f\uff0c\u534a\u5468\u671f\u3054\u3068\u306b $V = 0$ \u3068\u306a\u308b\u306f\u305a\u3067\u3042\u308b\u3002\u5b9f\u969b\u306e\u6570\u5024\u8a08\u7b97\u3067\u306f\uff0c\u523b\u307f\u5e45 $h$ \u3054\u3068\u306e\u5024\u3092\u6c42\u3081\u3066\u3044\u308b\u306e\u3067\uff0c\u53b3\u5bc6\u306b $V=0$ \u3068\u306a\u308b\u77ac\u9593\u304c\u5f97\u3089\u308c\u308b\u3068\u306f\u9650\u3089\u306a\u3044\u3002\u3057\u305f\u304c\u3063\u3066\uff0c\u901f\u5ea6\u306e\u5411\u304d\u304c\u5909\u308f\u308b\u6642\u523b\uff0c\u3064\u307e\u308a $V_i \\times V_{i+1} \\leq 0$ \u3068\u306a\u308b\u6642\u523b\u3092\u63a2\u3059\u3053\u3068\u306b\u3059\u308b\u3002<\/p>\n<p>\u305f\u3060\u30b0\u30e9\u30d5\u3092\u63cf\u304f\u305f\u3081\u3060\u3051\u306a\u3089\uff0c<code>N: 20000<\/code> \u306a\u3069\u3068\u3044\u3046\u5927\u304d\u306a\u5024\u3092\u3068\u308b\u5fc5\u8981\u306f\u306a\u3044\u306e\u3060\u304c\uff0c\u3053\u306e\u3088\u3046\u306b\u901f\u5ea6\u306e\u5411\u304d\u304c\u5909\u308f\u308b\u77ac\u9593\u3092\u3068\u3089\u3048\u305f\u3044\u3068\u304d\u306b\u306f\uff0c\u306a\u308b\u3079\u304f\u523b\u307f\u5e45\u3092\u5c0f\u3055\u304f\u3059\u308b\u3068\u826f\u3044\u306f\u305a\u3002<\/p>\n<p>\u3053\u306e\u3088\u3046\u306a\u7406\u7531\u3067 <code>N<\/code> \u306e\u5024\u3092\u5927\u304d\u304f\u3057\u3066\u3044\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[7]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"k\">for<\/span> <span class=\"nv\">th<\/span> <span class=\"nv\">in<\/span> <span class=\"p\">[<\/span><span class=\"mi\">10<\/span>, <span class=\"mi\">60<\/span>, 80<span class=\"p\">]<\/span> <span class=\"k\">do<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"nf\">printf<\/span><span class=\"p\">(<\/span><span class=\"no\">true<\/span>, <span class=\"s\">\"\u03b8_0 = ~2d\u00b0 \u306e\u3068\u304d\uff0c\u901f\u5ea6\u306e\u5411\u304d\u304c\u5909\u308f\u308b\u6642\u523b~%\"<\/span>, <span class=\"nv\">th<\/span><span class=\"p\">)<\/span>,\r\n  <span class=\"k\">for<\/span> <span class=\"nv\">i<\/span><span class=\"o\">:<\/span><span class=\"mi\">2<\/span> <span class=\"k\">thru<\/span> <span class=\"nv\">N<\/span> <span class=\"k\">do<\/span><span class=\"p\">(<\/span>\r\n    <span class=\"cm\">\/* if then \u306e\u5f8c\u306b\u8907\u6570\u884c\u306e\u5b9f\u884c\u6587\u3092\u66f8\u304f\u5834\u5408\u306f () \u3067\u56f2\u3080 *\/<\/span>\r\n    <span class=\"k\">if<\/span> <span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">][<\/span><span class=\"nv\">i<\/span><span class=\"p\">][<\/span>3<span class=\"p\">]<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">][<\/span><span class=\"nv\">i<\/span><span class=\"o\">+<\/span>1<span class=\"p\">][<\/span>3<span class=\"p\">]<\/span> <span class=\"o\">&lt;=<\/span> <span class=\"mi\">0<\/span> <span class=\"k\">then<\/span> <span class=\"p\">(<\/span>\r\n      <span class=\"nv\">t<\/span><span class=\"o\">:<\/span> <span class=\"p\">(<\/span><span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">][<\/span><span class=\"nv\">i<\/span><span class=\"p\">][<\/span>1<span class=\"p\">]<\/span> <span class=\"o\">+<\/span> <span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">][<\/span><span class=\"nv\">i<\/span><span class=\"o\">+<\/span>1<span class=\"p\">][<\/span>1<span class=\"p\">])<\/span><span class=\"o\">\/<\/span><span class=\"mi\">2<\/span>,\r\n      <span class=\"nf\">printf<\/span><span class=\"p\">(<\/span><span class=\"no\">true<\/span>, <span class=\"s\">\"    T = ~,5f\"<\/span>, <span class=\"nv\">t<\/span><span class=\"p\">)<\/span>\r\n    <span class=\"p\">)<\/span>\r\n  <span class=\"p\">)<\/span>, \r\n  <span class=\"nf\">printf<\/span><span class=\"p\">(<\/span><span class=\"no\">true<\/span>, <span class=\"s\">\"~%\"<\/span><span class=\"p\">)<\/span> <span class=\"cm\">\/* \u6539\u884c *\/<\/span>\r\n<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>\u03b8_0 = 10\u00b0 \u306e\u3068\u304d\uff0c\u901f\u5ea6\u306e\u5411\u304d\u304c\u5909\u308f\u308b\u6642\u523b\r\n    T = 0.50095    T = 1.00195    T = 1.50285\r\n\u03b8_0 = 60\u00b0 \u306e\u3068\u304d\uff0c\u901f\u5ea6\u306e\u5411\u304d\u304c\u5909\u308f\u308b\u6642\u523b\r\n    T = 0.53655    T = 1.07315    T = 1.60975\r\n\u03b8_0 = 80\u00b0 \u306e\u3068\u304d\uff0c\u901f\u5ea6\u306e\u5411\u304d\u304c\u5909\u308f\u308b\u6642\u523b\r\n    T = 0.56875    T = 1.13745    T = 1.70625\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<p>\u4ee5\u4e0a\u306e\u7d50\u679c\u304b\u3089\uff0c\u4f8b\u3048\u3070\u632f\u5e45\u306e\u521d\u671f\u6761\u4ef6\u3092 $\\theta_0 = 60^{\\circ}$ \u3068\u3057\u305f\u6642\u306e\u5468\u671f\u306f $T = 1.07315$\uff0c\u3059\u306a\u308f\u3061\u5358\u632f\u52d5\u306e\u5834\u5408\u306e\u5468\u671f<\/p>\n<p>$$T_0 = 2\\pi \\sqrt{\\frac{\\ell}{g}}$$<\/p>\n<p>\u306e $1.07315$ \u500d\u3060\u3068\u3044\u3046\u3053\u3068\u306b\u306a\u308b\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h3 id=\"\u554f\u984c\">\u554f\u984c<\/h3>\n<p>\u3067\u306f\u672c\u984c\u3002$\\theta_0 = 10^{\\circ}$ \u304b\u3089 $10^{\\circ}$ \u304d\u3056\u307f\u3067 $90^{\\circ}$ \u307e\u3067\uff0c $T = 0$ \u304b\u3089 $ T = 2$ \u307e\u3067\u89e3\u304d\uff0c\u632f\u5e45 $\\theta_0$ \u306b\u5bfe\u3059\u308b\u5468\u671f $T$ \u3092\u6c42\u3081\u3066\u307f\u3088\u3002<\/p>\n<p>\u30d2\u30f3\u30c8\uff1a<\/p>\n<div class=\"highlight\">\n<pre><span class=\"k\">for<\/span> <span class=\"nv\">th<\/span> <span class=\"nv\">in<\/span> <span class=\"p\">[<\/span><span class=\"mi\">10<\/span>, <span class=\"mi\">60<\/span>, 80<span class=\"p\">]<\/span>\r\n<\/pre>\n<\/div>\n<p>\u306a\u3069\u3068\u306a\u3063\u3066\u3044\u308b\u3068\u3053\u308d\u3092<\/p>\n<div class=\"highlight\">\n<pre><span class=\"k\">for<\/span> <span class=\"nv\">th<\/span><span class=\"o\">:<\/span><span class=\"mi\">10<\/span> <span class=\"k\">thru<\/span> <span class=\"mi\">90<\/span> <span class=\"k\">step<\/span> <span class=\"mi\">10<\/span>\r\n<\/pre>\n<\/div>\n<p>\u3068\u306a\u3069\u3059\u308c\u3070\u3044\u3044\u304b\u3082\u3002\u307e\u305f\uff0c\u4e0a\u8a18\u306e\u4f8b\u3067\u306f\u534a\u5468\u671f\u3054\u3068\u306b <code>T =<\/code> \u3068\u8868\u793a\u3057\u3066\u3044\u308b\u304c\uff0c\u5468\u671f\u306f2\u56de\u76ee\u306b\u901f\u5ea6\u306e\u5411\u304d\u304c\u5909\u308f\u308b\u6642\u9593\u3060\u304b\u3089\uff0c\u305d\u306e\u6642\u9593\u3060\u3051\u3092\u8868\u793a\u3059\u308b\u3088\u3046\u306b\u3057\u3066\u304f\u3060\u3055\u3044\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing text_cell rendered\">\n<div class=\"prompt input_prompt\"><\/div>\n<div class=\"inner_cell\">\n<div class=\"text_cell_render border-box-sizing rendered_html\">\n<h4 id=\"\u89e3\u7b54\u4f8b\">\u89e3\u7b54\u4f8b<\/h4>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[8]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"cm\">\/* \u632f\u5e45\u3068\u5468\u671f\u306e\u30b0\u30e9\u30d5\u3092\u63cf\u304f\u305f\u3081\u306b x, y \u306b\u8ffd\u52a0\u3057\u3066\u3044\u304f *\/<\/span>\r\n<span class=\"nv\">x<\/span><span class=\"o\">:<\/span> <span class=\"p\">[]<\/span>$ <span class=\"cm\">\/* \u632f\u5e45 *\/<\/span>\r\n<span class=\"nv\">y<\/span><span class=\"o\">:<\/span> <span class=\"p\">[]<\/span>$ <span class=\"cm\">\/* \u5468\u671f *\/<\/span>\r\n\r\n<span class=\"cm\">\/* \u53c2\u8003\u3068\u3057\u3066\uff0c\u307b\u307c\u5358\u632f\u52d5\u3068\u306a\u308b\u3067\u3042\u308d\u3046 \u03b8 = 1\u00b0\u306e\u5834\u5408 *\/<\/span>\r\n<span class=\"nv\">th<\/span><span class=\"o\">:<\/span> 1$\r\n<span class=\"nv\">theta0<\/span><span class=\"o\">:<\/span> <span class=\"nf\">radians<\/span><span class=\"p\">(<\/span><span class=\"nv\">th<\/span><span class=\"p\">)<\/span>$\r\n<span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">]<\/span><span class=\"o\">:<\/span> <span class=\"nf\">rk<\/span><span class=\"p\">([<\/span><span class=\"nf\">F1<\/span><span class=\"p\">(<\/span><span class=\"nv\">V<\/span><span class=\"p\">)<\/span>, <span class=\"nf\">F2<\/span><span class=\"p\">(<\/span><span class=\"nv\">theta<\/span><span class=\"p\">)]<\/span>, \r\n               <span class=\"p\">[<\/span><span class=\"nv\">theta<\/span>, <span class=\"nv\">V<\/span><span class=\"p\">]<\/span>, <span class=\"p\">[<\/span><span class=\"nv\">theta0<\/span>, <span class=\"nv\">V0<\/span><span class=\"p\">]<\/span>, \r\n               <span class=\"p\">[<\/span><span class=\"nv\">T<\/span>, <span class=\"nv\">T0<\/span>, <span class=\"nv\">T1<\/span>, <span class=\"nv\">h<\/span><span class=\"p\">])<\/span>$\r\n\r\n<span class=\"nf\">printf<\/span><span class=\"p\">(<\/span><span class=\"no\">true<\/span>, <span class=\"s\">\"\u03b8_0 = ~2d\u00b0 \u306e\u3068\u304d  \"<\/span>, <span class=\"nv\">th<\/span><span class=\"p\">)<\/span>$\r\n<span class=\"nv\">count<\/span><span class=\"o\">:<\/span> 0$\r\n<span class=\"k\">for<\/span> <span class=\"nv\">i<\/span><span class=\"o\">:<\/span><span class=\"mi\">2<\/span> <span class=\"k\">thru<\/span> <span class=\"nv\">N<\/span> <span class=\"k\">do<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"cm\">\/* if then \u306e\u5f8c\u306b\u8907\u6570\u306e\u5b9f\u884c\u6587\u3092\u66f8\u304f\u5834\u5408\u306f () \u3067\u56f2\u3080 *\/<\/span>\r\n  <span class=\"k\">if<\/span> <span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">][<\/span><span class=\"nv\">i<\/span><span class=\"p\">][<\/span>3<span class=\"p\">]<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">][<\/span><span class=\"nv\">i<\/span><span class=\"o\">+<\/span>1<span class=\"p\">][<\/span>3<span class=\"p\">]<\/span> <span class=\"o\">&lt;=<\/span> <span class=\"mi\">0<\/span> <span class=\"k\">then<\/span> <span class=\"p\">(<\/span>\r\n    <span class=\"nv\">count<\/span><span class=\"o\">:<\/span> <span class=\"nv\">count<\/span> <span class=\"o\">+<\/span> <span class=\"mi\">1<\/span>,\r\n    <span class=\"cm\">\/* \u5468\u671f\u306f2\u56de\u76ee\u306b\u901f\u5ea6\u306e\u5411\u304d\u304c\u5909\u308f\u308b\u6642\u9593\u3060\u304b\u3089 *\/<\/span>\r\n    <span class=\"k\">if<\/span> <span class=\"nv\">count<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">2<\/span> <span class=\"k\">then<\/span> <span class=\"p\">(<\/span>\r\n      <span class=\"nv\">t<\/span><span class=\"o\">:<\/span> <span class=\"p\">(<\/span><span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">][<\/span><span class=\"nv\">i<\/span><span class=\"p\">][<\/span>1<span class=\"p\">]<\/span> <span class=\"o\">+<\/span> <span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">][<\/span><span class=\"nv\">i<\/span><span class=\"o\">+<\/span>1<span class=\"p\">][<\/span>1<span class=\"p\">])<\/span><span class=\"o\">\/<\/span><span class=\"mi\">2<\/span>,\r\n      <span class=\"nv\">x<\/span><span class=\"o\">:<\/span> <span class=\"nf\">append<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span>, <span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">])<\/span>, \r\n      <span class=\"nv\">y<\/span><span class=\"o\">:<\/span> <span class=\"nf\">append<\/span><span class=\"p\">(<\/span><span class=\"nv\">y<\/span>, <span class=\"p\">[<\/span><span class=\"nv\">t<\/span><span class=\"p\">])<\/span>,\r\n      <span class=\"nf\">printf<\/span><span class=\"p\">(<\/span><span class=\"no\">true<\/span>, <span class=\"s\">\"\u5468\u671f T = ~,5f~%\"<\/span>, <span class=\"nv\">t<\/span><span class=\"p\">)<\/span>, \r\n      <span class=\"nf\">return<\/span><span class=\"p\">()<\/span> <span class=\"cm\">\/* \u5468\u671f\u3092\u8868\u793a\u3057\u305f\u3089 for \u30eb\u30fc\u30d7\u304b\u3089\u629c\u3051\u308b *\/<\/span>\r\n    <span class=\"p\">)<\/span>\r\n  <span class=\"p\">)<\/span>\r\n<span class=\"p\">)<\/span>$\r\n\r\n<span class=\"cm\">\/* \u632f\u5e45\u306e\u521d\u671f\u5024\u3092\u5909\u3048\u3066\u6570\u5024\u8a08\u7b97... *\/<\/span>\r\n<span class=\"k\">for<\/span> <span class=\"nv\">th<\/span><span class=\"o\">:<\/span><span class=\"mi\">10<\/span> <span class=\"k\">thru<\/span> <span class=\"mi\">90<\/span> <span class=\"k\">step<\/span> <span class=\"mi\">10<\/span> <span class=\"k\">do<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"nv\">theta0<\/span><span class=\"o\">:<\/span> <span class=\"nf\">radians<\/span><span class=\"p\">(<\/span><span class=\"nv\">th<\/span><span class=\"p\">)<\/span>,\r\n  <span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">]<\/span><span class=\"o\">:<\/span> <span class=\"nf\">rk<\/span><span class=\"p\">([<\/span><span class=\"nf\">F1<\/span><span class=\"p\">(<\/span><span class=\"nv\">V<\/span><span class=\"p\">)<\/span>, <span class=\"nf\">F2<\/span><span class=\"p\">(<\/span><span class=\"nv\">theta<\/span><span class=\"p\">)]<\/span>, \r\n                 <span class=\"p\">[<\/span><span class=\"nv\">theta<\/span>, <span class=\"nv\">V<\/span><span class=\"p\">]<\/span>, <span class=\"p\">[<\/span><span class=\"nv\">theta0<\/span>, <span class=\"nv\">V0<\/span><span class=\"p\">]<\/span>, \r\n                 <span class=\"p\">[<\/span><span class=\"nv\">T<\/span>, <span class=\"nv\">T0<\/span>, <span class=\"nv\">T1<\/span>, <span class=\"nv\">h<\/span><span class=\"p\">])<\/span>\r\n<span class=\"p\">)<\/span>$\r\n\r\n<span class=\"k\">for<\/span> <span class=\"nv\">th<\/span><span class=\"o\">:<\/span><span class=\"mi\">10<\/span> <span class=\"k\">thru<\/span> <span class=\"mi\">90<\/span> <span class=\"k\">step<\/span> <span class=\"mi\">10<\/span> <span class=\"k\">do<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"nf\">printf<\/span><span class=\"p\">(<\/span><span class=\"no\">true<\/span>, <span class=\"s\">\"\u03b8_0 = ~2d\u00b0 \u306e\u3068\u304d  \"<\/span>, <span class=\"nv\">th<\/span><span class=\"p\">)<\/span>, \r\n  <span class=\"nv\">count<\/span><span class=\"o\">:<\/span> <span class=\"mi\">0<\/span>,\r\n  <span class=\"k\">for<\/span> <span class=\"nv\">i<\/span><span class=\"o\">:<\/span><span class=\"mi\">2<\/span> <span class=\"k\">thru<\/span> <span class=\"nv\">N<\/span> <span class=\"k\">do<\/span><span class=\"p\">(<\/span>\r\n    <span class=\"cm\">\/* if then \u306e\u5f8c\u306b\u8907\u6570\u306e\u5b9f\u884c\u6587\u3092\u66f8\u304f\u5834\u5408\u306f () \u3067\u56f2\u3080 *\/<\/span>\r\n    <span class=\"k\">if<\/span> <span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">][<\/span><span class=\"nv\">i<\/span><span class=\"p\">][<\/span>3<span class=\"p\">]<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">][<\/span><span class=\"nv\">i<\/span><span class=\"o\">+<\/span>1<span class=\"p\">][<\/span>3<span class=\"p\">]<\/span> <span class=\"o\">&lt;=<\/span> <span class=\"mi\">0<\/span> <span class=\"k\">then<\/span> <span class=\"p\">(<\/span>\r\n      <span class=\"nv\">count<\/span><span class=\"o\">:<\/span> <span class=\"nv\">count<\/span> <span class=\"o\">+<\/span> <span class=\"mi\">1<\/span>,\r\n      <span class=\"cm\">\/* \u5468\u671f\u306f2\u56de\u76ee\u306b\u901f\u5ea6\u306e\u5411\u304d\u304c\u5909\u308f\u308b\u6642\u9593\u3060\u304b\u3089 *\/<\/span>\r\n      <span class=\"k\">if<\/span> <span class=\"nv\">count<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">2<\/span> <span class=\"k\">then<\/span> <span class=\"p\">(<\/span>\r\n        <span class=\"nv\">t<\/span><span class=\"o\">:<\/span> <span class=\"p\">(<\/span><span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">][<\/span><span class=\"nv\">i<\/span><span class=\"p\">][<\/span>1<span class=\"p\">]<\/span> <span class=\"o\">+<\/span> <span class=\"nv\">ansdeg<\/span><span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">][<\/span><span class=\"nv\">i<\/span><span class=\"o\">+<\/span>1<span class=\"p\">][<\/span>1<span class=\"p\">])<\/span><span class=\"o\">\/<\/span><span class=\"mi\">2<\/span>,\r\n        <span class=\"nv\">x<\/span><span class=\"o\">:<\/span> <span class=\"nf\">append<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span>, <span class=\"p\">[<\/span><span class=\"nv\">th<\/span><span class=\"p\">])<\/span>, \r\n        <span class=\"nv\">y<\/span><span class=\"o\">:<\/span> <span class=\"nf\">append<\/span><span class=\"p\">(<\/span><span class=\"nv\">y<\/span>, <span class=\"p\">[<\/span><span class=\"nv\">t<\/span><span class=\"p\">])<\/span>,\r\n        <span class=\"nf\">printf<\/span><span class=\"p\">(<\/span><span class=\"no\">true<\/span>, <span class=\"s\">\"\u5468\u671f T = ~,5f~%\"<\/span>, <span class=\"nv\">t<\/span><span class=\"p\">)<\/span>, \r\n        <span class=\"nf\">return<\/span><span class=\"p\">()<\/span> <span class=\"cm\">\/* \u5468\u671f\u3092\u8868\u793a\u3057\u305f\u3089 for \u30eb\u30fc\u30d7\u304b\u3089\u629c\u3051\u308b *\/<\/span>\r\n      <span class=\"p\">)<\/span>\r\n    <span class=\"p\">)<\/span>\r\n  <span class=\"p\">)<\/span>\r\n<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>\u03b8_0 =  1\u00b0 \u306e\u3068\u304d  \u5468\u671f T = 1.00005\r\n\u03b8_0 = 10\u00b0 \u306e\u3068\u304d  \u5468\u671f T = 1.00195\r\n\u03b8_0 = 20\u00b0 \u306e\u3068\u304d  \u5468\u671f T = 1.00765\r\n\u03b8_0 = 30\u00b0 \u306e\u3068\u304d  \u5468\u671f T = 1.01745\r\n\u03b8_0 = 40\u00b0 \u306e\u3068\u304d  \u5468\u671f T = 1.03135\r\n\u03b8_0 = 50\u00b0 \u306e\u3068\u304d  \u5468\u671f T = 1.04975\r\n\u03b8_0 = 60\u00b0 \u306e\u3068\u304d  \u5468\u671f T = 1.07315\r\n\u03b8_0 = 70\u00b0 \u306e\u3068\u304d  \u5468\u671f T = 1.10215\r\n\u03b8_0 = 80\u00b0 \u306e\u3068\u304d  \u5468\u671f T = 1.13745\r\n\u03b8_0 = 90\u00b0 \u306e\u3068\u304d  \u5468\u671f T = 1.18035\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-maxima\">\n<pre><span class=\"cm\">\/* \u30aa\u30d7\u30b7\u30e7\u30f3\u3092\u8a2d\u5b9a\u3057\u3066\u30b0\u30e9\u30d5\u3092\u63cf\u304f *\/<\/span>\r\n<span class=\"nf\">draw2d<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"nv\">font<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"Arial\"<\/span>, <span class=\"nv\">font_size<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">16<\/span>, \r\n  <span class=\"nv\">title<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u5358\u632f\u308a\u5b50\u306e\u632f\u5e45\u3068\u5468\u671f\u306e\u95a2\u4fc2\"<\/span>,\r\n  <span class=\"cm\">\/* \u8868\u793a\u7bc4\u56f2 *\/<\/span>\r\n  <span class=\"nv\">xrange<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"o\">-<\/span><span class=\"mi\">2<\/span>, 92<span class=\"p\">]<\/span>, <span class=\"nv\">yrange<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"mf\">0.95<\/span>, 1<span class=\"o\">.<\/span>2<span class=\"p\">]<\/span>, \r\n  <span class=\"cm\">\/* \u30b0\u30ea\u30c3\u30c9  *\/<\/span>\r\n  <span class=\"nv\">grid<\/span> <span class=\"o\">=<\/span> <span class=\"no\">true<\/span>, \r\n  <span class=\"cm\">\/* \u5ea7\u6a19\u8ef8\u306e\u30e9\u30d9\u30eb *\/<\/span>\r\n  <span class=\"nv\">xlabel<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u632f\u5e45 \u03b8 (\u00b0)\"<\/span>, <span class=\"nv\">ylabel<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u898f\u683c\u5316\u3055\u308c\u305f\u5468\u671f T\"<\/span>, \r\n  \r\n  <span class=\"cm\">\/* \u70b9\u3092\u5857\u308a\u3064\u3076\u3057\u305f\u4e38\u306b *\/<\/span>\r\n  <span class=\"nv\">point_type<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">7<\/span>,\r\n  <span class=\"cm\">\/* \u70b9\u306e\u5927\u304d\u3055 *\/<\/span>\r\n  <span class=\"nv\">point_size<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.5<\/span>,\r\n  <span class=\"nf\">points<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span>, <span class=\"nv\">y<\/span><span class=\"p\">)<\/span>\r\n<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 \">\n<p><!--?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?--><\/p>\n<p><img loading=\"lazy\" decoding=\"async\" class=\"aligncenter size-large wp-image-5011\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/mfuri03A.svg\" alt=\"\" width=\"640\" height=\"533\" \/><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"excerpt":{"rendered":"<p>Maxima \u3092\u4f7f\u3063\u3066\u904b\u52d5\u65b9\u7a0b\u5f0f\u3092\u6570\u5024\u7684\u306b\u89e3\u304d\uff0c\u5358\u632f\u308a\u5b50\u306e\u632f\u5e45\u3068\u5468\u671f\u306e\u95a2\u4fc2\u3092\u8abf\u3079\u308b\u3002\u5c0e\u51fa\u306b\u3064\u3044\u3066\u306f\uff0c\u4ee5\u4e0b\u3092\u53c2\u7167\u3002<\/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\/%e5%8f%82%e8%80%83%ef%bc%9amaxima-%e7%b7%a8\/maxima-%e3%81%a7%e9%81%8b%e5%8b%95%e6%96%b9%e7%a8%8b%e5%bc%8f%e3%82%92%e8%a7%a3%e3%81%84%e3%81%a6%e6%8c%af%e5%b9%85%e3%81%a8%e5%91%a8%e6%9c%9f%e3%81%ae%e9%96%a2%e4%bf%82%e3%82%92%e8%aa%bf%e3%81%b9\/\">\u7d9a\u304d\u3092\u8aad\u3080<\/a><\/p>\n<ul>\n<li>\u300c\u5358\u632f\u308a\u5b50\u306e\u904b\u52d5\u65b9\u7a0b\u5f0f\u3092\u6570\u5024\u7684\u306b\u89e3\u304f\u6e96\u5099\u300d<\/li>\n<\/ul>\n","protected":false},"author":33,"featured_media":0,"parent":7137,"menu_order":40,"comment_status":"closed","ping_status":"closed","template":"","meta":{"inline_featured_image":false,"footnotes":""},"class_list":["post-4950","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\/4950","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=4950"}],"version-history":[{"count":13,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/4950\/revisions"}],"predecessor-version":[{"id":5090,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/4950\/revisions\/5090"}],"up":[{"embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/7137"}],"wp:attachment":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/media?parent=4950"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}