{"id":2395,"date":"2023-04-13T10:36:30","date_gmt":"2023-04-13T01:36:30","guid":{"rendered":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/?page_id=2395"},"modified":"2023-07-13T16:13:45","modified_gmt":"2023-07-13T07:13:45","slug":"maxima-jupyter-%e3%81%a6%e3%82%99%e5%b8%b8%e5%be%ae%e5%88%86%e6%96%b9%e7%a8%8b%e5%bc%8f","status":"publish","type":"page","link":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/%e7%90%86%e5%b7%a5%e7%b3%bb%e3%81%ae%e6%95%b0%e5%ad%a6c\/maxima-%e3%81%a7%e7%90%86%e5%b7%a5%e7%b3%bb%e3%81%ae%e6%95%b0%e5%ad%a6c\/maxima-jupyter-%e3%81%a6%e3%82%99%e5%b8%b8%e5%be%ae%e5%88%86%e6%96%b9%e7%a8%8b%e5%bc%8f\/","title":{"rendered":"Maxima \u3066\u3099\u5e38\u5fae\u5206\u65b9\u7a0b\u5f0f"},"content":{"rendered":"<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=\"\u5909\u6570\u5206\u96e2\u6cd5\">\u5909\u6570\u5206\u96e2\u6cd5<\/h3>\n<h4 id=\"\u4f8b\u984c\">\u4f8b\u984c<\/h4>\n<p>$$ \\frac{dy}{dx} = &#8211; 2 x\\, y $$<\/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<p>\u89e3\u304f\u3079\u304d\u65b9\u7a0b\u5f0f\u3092\uff0c\u5909\u6570 <code>eq<\/code> \u306b\u4ee3\u5165\u3057\u307e\u3059\u3002<code>diff()<\/code> \u306fMaxima \u3067\u5fae\u5206\u3092\u884c\u3046\u95a2\u6570\u3067\u3059\u304c\uff0c\u5fae\u5206\u65b9\u7a0b\u5f0f\u3092\u5b9a\u7fa9\u3059\u308b\u3068\u304d\u306f <code>'<\/code> \u3067\u59cb\u307e\u308b\u306e\u304c\u5927\u4e8b\u3067\u3059\u3002<\/p>\n<p><code>'diff(y, x)<\/code> \u3067\u306f\u306a\u304f\uff0c<code>diff(y, x)<\/code> \u3068 <code>'<\/code> \u3092\u6700\u521d\u306b\u3064\u3051\u308b\u306e\u3092\u5fd8\u308c\u308b\u3068\u3069\u3046\u306a\u308b\u304b\u3068\u3044\u3046\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[1]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/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[1]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{1}$}0\\]<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[2]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"cm\">\/* \u89e3\u304f\u3079\u304d\u65b9\u7a0b\u5f0f\u3092\uff0c\u5909\u6570 eq \u306b\u4ee3\u5165\u3002*\/<\/span>\r\n\r\n<span class=\"nv\">eq<\/span><span class=\"o\">:<\/span> <span class=\"o\">'<\/span><span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">=<\/span> <span class=\"o\">-<\/span>2<span class=\"o\">*<\/span><span class=\"nv\">x<\/span><span class=\"o\">*<\/span><span class=\"nv\">y<\/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[2]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{2}$}\\frac{d}{d\\,x}\\,y=-2\\,x\\,y\\]<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[3]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"cm\">\/* \u5fae\u5206\u65b9\u7a0b\u5f0f\u3092\u89e3\u304f\u95a2\u6570\u306f ode2() \u3002*\/<\/span>\r\n<span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq<\/span>, <span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/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\">\\[\\tag{${\\it \\%o}_{3}$}y={\\it \\%c}\\,e^ {- x^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>\u4e0a\u8a18\u306e\u7b54\u3048\u306e <code>%c<\/code> \u306f\u7a4d\u5206\u5b9a\u6570 $C$ \u3092\u8868\u3057\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=\"\u30de\u30eb\u30b5\u30b9\u306e\u4eba\u53e3\u30e2\u30c7\u30eb\">\u30de\u30eb\u30b5\u30b9\u306e\u4eba\u53e3\u30e2\u30c7\u30eb<\/h3>\n<p>$$\\frac{dN}{dt} = \\gamma \\, N$$<\/p>\n<p>\u3092\u521d\u671f\u6761\u4ef6 $t = t_0$ \u3067 $N(t_0) = N_0$ \u3068\u3057\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[4]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">eq2<\/span><span class=\"o\">:<\/span> <span class=\"o\">'<\/span><span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nv\">N<\/span>, <span class=\"nv\">t<\/span><span class=\"p\">)<\/span> <span class=\"o\">=<\/span> <span class=\"nv\">gamma<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">N<\/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\">\\[\\tag{${\\it \\%o}_{4}$}\\frac{d}{d\\,t}\\,N=N\\,\\gamma\\]<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[5]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq2<\/span>, <span class=\"nv\">N<\/span>, <span class=\"nv\">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[5]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{5}$}N={\\it \\%c}\\,e^{t\\,\\gamma}\\]<\/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\">\/* 1\u968e\u5fae\u5206\u65b9\u7a0b\u5f0f\u306e\u521d\u671f\u6761\u4ef6 *\/<\/span>\r\n\r\n<span class=\"nf\">ic1<\/span><span class=\"p\">(<\/span><span class=\"nv\">%<\/span>, <span class=\"nv\">t<\/span> <span class=\"o\">=<\/span> <span class=\"nv\">t0<\/span>, <span class=\"nv\">N<\/span> <span class=\"o\">=<\/span> <span class=\"nv\">N0<\/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[6]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{6}$}N=N_{0}\\,e^{t\\,\\gamma-t_{0}\\,\\gamma}\\]<\/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=\"\u53c2\u8003\uff1a\u7c73\u56fd\u306e\u4eba\u53e3\u3068\u30de\u30eb\u30b5\u30b9\u30e2\u30c7\u30eb\">\u53c2\u8003\uff1a\u7c73\u56fd\u306e\u4eba\u53e3\u3068\u30de\u30eb\u30b5\u30b9\u30e2\u30c7\u30eb<\/h4>\n<ul>\n<li>\u53c2\u8003\uff1a<a href=\"https:\/\/www.nippyo.co.jp\/shop\/book\/1240.html\">\u300c\u5fae\u5206\u65b9\u7a0b\u5f0f\u3067\u6570\u5b66\u30e2\u30c7\u30eb\u3092\u4f5c\u308d\u3046\u300d \u30d0\u30fc\u30b8\u30a7\u30b9\uff0c\u30dc\u30ea\u30fc\u8457\uff0c\u65e5\u672c\u8a55\u8ad6\u793e<\/a><\/li>\n<\/ul>\n<p>\u4e0a\u8a18\u306b\u3088\u308c\u3070\uff0c1790\u5e74\u304b\u30891930\u5e74\u306e\u30a2\u30e1\u30ea\u30ab\u306e\u4eba\u53e3\u306f\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u306a\u3063\u3066\u3044\u308b\uff08\u4eba\u53e3\u306e\u5358\u4f4d\u306f\u767e\u4e07\u4eba\uff09\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=\"nv\">usa<\/span><span class=\"o\">:<\/span> <span class=\"p\">[<\/span>\r\n<span class=\"cm\">\/* \u897f\u66a6, \u4eba\u53e3\uff08\u767e\u4e07\u4eba\uff09*\/<\/span>\r\n  <span class=\"p\">[<\/span><span class=\"mi\">1790<\/span>,  3<span class=\"o\">.<\/span>9<span class=\"p\">]<\/span>, \r\n  <span class=\"p\">[<\/span><span class=\"mi\">1800<\/span>,  5<span class=\"o\">.<\/span>3<span class=\"p\">]<\/span>,\r\n  <span class=\"p\">[<\/span><span class=\"mi\">1810<\/span>,  7<span class=\"o\">.<\/span>2<span class=\"p\">]<\/span>,\r\n  <span class=\"p\">[<\/span><span class=\"mi\">1820<\/span>,  9<span class=\"o\">.<\/span>6<span class=\"p\">]<\/span>,\r\n  <span class=\"p\">[<\/span><span class=\"mi\">1830<\/span>, 12<span class=\"o\">.<\/span>9<span class=\"p\">]<\/span>,\r\n  <span class=\"p\">[<\/span><span class=\"mi\">1840<\/span>, 17<span class=\"o\">.<\/span>1<span class=\"p\">]<\/span>,\r\n  <span class=\"p\">[<\/span><span class=\"mi\">1850<\/span>, 23<span class=\"o\">.<\/span>2<span class=\"p\">]<\/span>,\r\n  <span class=\"p\">[<\/span><span class=\"mi\">1860<\/span>, 31<span class=\"o\">.<\/span>4<span class=\"p\">]<\/span>,\r\n  <span class=\"p\">[<\/span><span class=\"mi\">1870<\/span>, 38<span class=\"o\">.<\/span>6<span class=\"p\">]<\/span>,\r\n  <span class=\"p\">[<\/span><span class=\"mi\">1880<\/span>, 50<span class=\"o\">.<\/span>2<span class=\"p\">]<\/span>,\r\n  <span class=\"p\">[<\/span><span class=\"mi\">1890<\/span>, 62<span class=\"o\">.<\/span>9<span class=\"p\">]<\/span>,\r\n  <span class=\"p\">[<\/span><span class=\"mi\">1900<\/span>, 76<span class=\"o\">.<\/span>0<span class=\"p\">]<\/span>,\r\n  <span class=\"p\">[<\/span><span class=\"mi\">1910<\/span>, 92<span class=\"o\">.<\/span>0<span class=\"p\">]<\/span>,\r\n  <span class=\"p\">[<\/span><span class=\"mi\">1920<\/span>,106<span class=\"o\">.<\/span>5<span class=\"p\">]<\/span>,\r\n  <span class=\"p\">[<\/span><span class=\"mi\">1930<\/span>,123<span class=\"o\">.<\/span>2<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>$t_0 = 1790$ \uff08\u5e74\uff09\u3068\u3059\u308b\u3068\uff0c$N_0 = 3.9$ \u3002\u3057\u305f\u304c\u3063\u3066\uff0c\u30de\u30eb\u30b5\u30b9\u306e\u4eba\u53e3\u30e2\u30c7\u30eb\u3092 $N_m(t)$ \u3068\u3059\u308b\u3068<\/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-maxima\">\n<pre><span class=\"nv\">t0<\/span><span class=\"o\">:<\/span> <span class=\"nv\">usa<\/span><span class=\"p\">[<\/span>1<span class=\"p\">][<\/span>1<span class=\"p\">]<\/span>; <span class=\"cm\">\/* 1\u884c\u76ee\u306e1\u5217\u76ee *\/<\/span>\r\n<span class=\"nv\">N0<\/span><span class=\"o\">:<\/span> <span class=\"nv\">usa<\/span><span class=\"p\">[<\/span>1<span class=\"p\">][<\/span>2<span class=\"p\">]<\/span>; <span class=\"cm\">\/* 1\u884c\u76ee\u306e2\u5217\u76ee *\/<\/span>\r\n\r\n<span class=\"nf\">kill<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span>, <span class=\"nv\">gamma<\/span><span class=\"p\">)<\/span>$\r\n\r\n<span class=\"nf\">Nm<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)<\/span> <span class=\"o\">:=<\/span> <span class=\"nv\">N0<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">exp<\/span><span class=\"p\">(<\/span><span class=\"nv\">gamma1<\/span><span class=\"o\">*<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"o\">-<\/span><span class=\"nv\">t0<\/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[8]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{8}$}1790\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[8]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{9}$}3.9\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[8]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{11}$}{\\it Nm}\\left(t\\right):=N_{0}\\,\\exp \\left(\\gamma_{1}\\,\\left(t-t_{0}\\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>\u6b8b\u308a\u306e\u30d1\u30e9\u30e1\u30fc\u30bf\u30fc $\\gamma_1$ \u306f\uff0c\u5225\u306e\u6642\u523b $t_1$ \u306b\u304a\u3051\u308b $N_1 = N_m(t_1)$ \u304b\u3089\u6c42\u3081\u308b\u3002<\/p>\n<p>\\begin{eqnarray}<br \/>\nN_1 &amp;=&amp; N_0 e^{\\gamma_1 (t_1 &#8211; t_0)} \\\\<br \/>\n\\log \\frac{N_1}{N_0} &amp;=&amp; \\gamma_1 (t_1 &#8211; t_0) \\\\<br \/>\n\\therefore\\ \\ \\gamma_1 &amp;=&amp; \\frac{1}{t_1 &#8211; t_0} \\log \\frac{N_1}{N_0}<br \/>\n\\end{eqnarray}<\/p>\n<p>\u305f\u3068\u3048\u3070\uff0c$t_1 = 1830$ \u3068\u3059\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[9]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">t1<\/span><span class=\"o\">:<\/span> <span class=\"nv\">usa<\/span><span class=\"p\">[<\/span>5<span class=\"p\">][<\/span>1<span class=\"p\">]<\/span>; <span class=\"cm\">\/* 5\u884c\u76ee\u306e1\u5217\u76ee *\/<\/span>\r\n<span class=\"nv\">N1<\/span><span class=\"o\">:<\/span> <span class=\"nv\">usa<\/span><span class=\"p\">[<\/span>5<span class=\"p\">][<\/span>2<span class=\"p\">]<\/span>; <span class=\"cm\">\/* 5\u884c\u76ee\u306e2\u5217\u76ee *\/<\/span>\r\n\r\n<span class=\"nv\">gamma1<\/span><span class=\"o\">:<\/span> 1<span class=\"o\">\/<\/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=\"nf\">log<\/span><span class=\"p\">(<\/span><span class=\"nv\">N1<\/span><span class=\"o\">\/<\/span><span class=\"nv\">N0<\/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[9]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{12}$}1830\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[9]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{13}$}12.9\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[9]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{14}$}0.02990626895580065\\]<\/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>\u3053\u306e\u3088\u3046\u306b\u3057\u3066\u6c42\u3081\u3089\u308c\u305f $N(t)$ \u3092\uff0c\u4eba\u53e3\u30c7\u30fc\u30bf <code>usa<\/code> \u3068\u5171\u30b0\u30e9\u30d5\u306b\u3057\u3066\u307f\u307e\u3059\u3002<\/p>\n<p>Maxima \u3067\u30b0\u30e9\u30d5\u3092\u4f5c\u6210\u3059\u308b\u95a2\u6570\u3068\u3057\u3066 <code>plot2d()<\/code> \u306a\u3069\u306e\u3044\u308f\u3086\u308b plot \u7cfb\u3068\uff0c<code>draw2d()<\/code> \u306a\u3069\u306e\u3044\u308f\u3086\u308b draw \u7cfb\u304c\u3042\u308a\u307e\u3059\u3002\u30d1\u30e9\u30e1\u30fc\u30bf\u8a2d\u5b9a\u7b49\u306e\u66f8\u5f0f\u304c\u7570\u306a\u308a\u307e\u3059\u306e\u3067\uff0c\u305d\u308c\u305e\u308c\u306e\u4f8b\u3092\u3042\u3052\u3066\u304a\u304d\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-maxima\">\n<pre><span class=\"cm\">\/* plot2d() \u306b\u3088\u308b\u30b0\u30e9\u30d5 *\/<\/span>\r\n\r\n<span class=\"nf\">plot2d<\/span><span class=\"p\">([[<\/span><span class=\"nv\">discrete<\/span>, <span class=\"nv\">usa<\/span><span class=\"p\">]<\/span>, <span class=\"nf\">Nm<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)]<\/span>, <span class=\"p\">[<\/span><span class=\"nv\">t<\/span>, <span class=\"mi\">1790<\/span>, 1940<span class=\"p\">]<\/span>, \r\n       <span class=\"p\">[<\/span><span class=\"nv\">style<\/span>, <span class=\"p\">[<\/span><span class=\"nv\">points<\/span>, 1<span class=\"o\">.<\/span>5<span class=\"p\">]<\/span>, <span class=\"nv\">lines<\/span><span class=\"p\">]<\/span>, \r\n       <span class=\"p\">[<\/span><span class=\"nv\">legend<\/span>, <span class=\"s\">\"\u30a2\u30e1\u30ea\u30ab\u306e\u4eba\u53e3\"<\/span>, <span class=\"s\">\"\u30de\u30eb\u30b5\u30b9\u30e2\u30c7\u30eb\"<\/span><span class=\"p\">]<\/span>, \r\n       <span class=\"p\">[<\/span><span class=\"nv\">xlabel<\/span>, <span class=\"s\">\"\u5e74\"<\/span><span class=\"p\">]<\/span>, <span class=\"p\">[<\/span><span class=\"nv\">ylabel<\/span>, <span class=\"s\">\"\u4eba\u53e3\uff08\u5358\u4f4d\uff1a\u767e\u4e07\u4eba\uff09\"<\/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 \">\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-6120\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/mmathc101.svg\" alt=\"\" width=\"640\" height=\"481\" \/><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[11]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"cm\">\/* draw2d() \u306b\u3088\u308b\u30b0\u30e9\u30d5 *\/<\/span>\r\n<span class=\"nf\">draw2d<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"nv\">color<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"red\"<\/span>,\r\n  <span class=\"nv\">key<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u30de\u30eb\u30b5\u30b9\u30e2\u30c7\u30eb\"<\/span>, \r\n  <span class=\"nf\">explicit<\/span><span class=\"p\">(<\/span><span class=\"nf\">Nm<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">t<\/span>, <span class=\"mi\">1790<\/span>, <span class=\"mi\">1940<\/span><span class=\"p\">)<\/span>, \r\n  \r\n  <span class=\"nv\">color<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"blue\"<\/span>, \r\n  <span class=\"nv\">point_type<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">7<\/span>,\r\n  <span class=\"nv\">point_size<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.7<\/span>,\r\n  <span class=\"nv\">key<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u30a2\u30e1\u30ea\u30ab\u306e\u4eba\u53e3\"<\/span>, \r\n  <span class=\"nf\">points<\/span><span class=\"p\">(<\/span><span class=\"nv\">usa<\/span><span class=\"p\">)<\/span>, \r\n  \r\n  <span class=\"nv\">xlabel<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u5e74\"<\/span>, <span class=\"nv\">ylabel<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u4eba\u53e3\uff08\u5358\u4f4d\uff1a\u767e\u4e07\u4eba\uff09\"<\/span>,\r\n  <span class=\"nv\">yrange<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"o\">-<\/span><span class=\"mi\">20<\/span>, 410<span class=\"p\">]<\/span>, <span class=\"nv\">xaxis<\/span> <span class=\"o\">=<\/span> <span class=\"no\">true<\/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-6121\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/mmathc102.svg\" alt=\"\" width=\"640\" height=\"481\" \/><\/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=\"\u30f4\u30a7\u30a2\u30d5\u30eb\u30b9\u30c8\u306b\u3088\u308b\u4fee\u6b63\u4eba\u53e3\u30e2\u30c7\u30eb\">\u30f4\u30a7\u30a2\u30d5\u30eb\u30b9\u30c8\u306b\u3088\u308b\u4fee\u6b63\u4eba\u53e3\u30e2\u30c7\u30eb<\/h3>\n<p>$$<br \/>\n\\frac{dN}{dt} = \\gamma N \\left(1 &#8211; \\frac{N}{N_{\\rm max}}\\right)<br \/>\n$$<\/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-maxima\">\n<pre><span class=\"cm\">\/* \u5ff5\u306e\u305f\u3081\uff0c\u5909\u6570\u3092\u521d\u671f\u5316 *\/<\/span>\r\n<span class=\"nf\">kill<\/span><span class=\"p\">(<\/span><span class=\"nv\">N<\/span>, <span class=\"nv\">t<\/span>, <span class=\"nv\">Nmax<\/span>, <span class=\"nv\">N0<\/span>, <span class=\"nv\">t0<\/span><span class=\"p\">)<\/span>$\r\n\r\n<span class=\"nv\">eq3<\/span><span class=\"o\">:<\/span> <span class=\"o\">'<\/span><span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nv\">N<\/span>, <span class=\"nv\">t<\/span><span class=\"p\">)<\/span> <span class=\"o\">=<\/span> <span class=\"nv\">gamma2<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">N<\/span> <span class=\"o\">*<\/span> <span class=\"p\">(<\/span><span class=\"mi\">1<\/span> <span class=\"o\">-<\/span> <span class=\"nv\">N<\/span><span class=\"o\">\/<\/span><span class=\"nv\">Nmax<\/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[12]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{20}$}\\frac{d}{d\\,t}\\,N=N\\,\\left(1-\\frac{N}{{\\it Nmax}}\\right)\\,\\gamma_{2}\\]<\/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[13]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">ans<\/span><span class=\"o\">:<\/span> <span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq3<\/span>, <span class=\"nv\">N<\/span>, <span class=\"nv\">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[13]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{21}$}-\\frac{\\log \\left(N-{\\it Nmax}\\right)-\\log N}{\\gamma_{2}}=t+{\\it \\%c}\\]<\/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>\u3069\u3046\u3084\u3089\uff0cMaxima \u306f\u3061\u3083\u3093\u3068 $N= \\dots$ \u306e\u3088\u3046\u306b\u89e3\u3044\u3066\u306f\u304f\u308c\u306a\u3044\u3088\u3046\u3067\u3059\u3002\u3042\u3068\u306f\u4eba\u9593\u306e\u624b\u3067\u6574\u7406\u3057\u3066\u3044\u304f\u3053\u3068\u306b\u306a\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<p>\u4e0a\u5f0f <code>ans<\/code> \u306e\u4e21\u8fba\u306b $\\gamma_2$ \u3092\u304b\u3051\u305f\u5f0f\u3092 <code>ans2<\/code> \u3068\u3057\uff0c<\/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-maxima\">\n<pre><span class=\"nv\">ans2<\/span><span class=\"o\">:<\/span> <span class=\"nv\">gamma2<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">ans<\/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\">\\[\\tag{${\\it \\%o}_{22}$}\\log N-\\log \\left(N-{\\it Nmax}\\right)=\\gamma_{2}\\,\\left(t+{\\it \\%c}\\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>\u4e21\u8fba\u306e exponential \u3092\u3068\u3063\u305f\u5f0f\u3092 <code>ans3<\/code> \u3068\u3059\u308b\u3002<code>lhs()<\/code> \u306f\u5de6\u8fba\uff08left-hand-side\uff09\uff0c<code>rhs()<\/code> \u306f\u53f3\u8fba\uff08right-hand-side\uff09\u3092\u3068\u308b\u95a2\u6570\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[15]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">ans3<\/span><span class=\"o\">:<\/span> <span class=\"nf\">exp<\/span><span class=\"p\">(<\/span><span class=\"nf\">lhs<\/span><span class=\"p\">(<\/span><span class=\"nv\">ans2<\/span><span class=\"p\">))<\/span> <span class=\"o\">=<\/span> <span class=\"nf\">exp<\/span><span class=\"p\">(<\/span><span class=\"nf\">rhs<\/span><span class=\"p\">(<\/span><span class=\"nv\">ans2<\/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[15]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{23}$}\\frac{N}{N-{\\it Nmax}}=e^{\\gamma_{2}\\,\\left(t+{\\it \\%c}\\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>\u3053\u306e <code>ans3<\/code> \u5f0f\u3092 <code>N<\/code> \u306b\u3064\u3044\u3066\u89e3\u304f\u3002\u65b9\u7a0b\u5f0f\u3092\u89e3\u304f\u306e\u306f <code>solve()<\/code> \u95a2\u6570\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[16]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">sol<\/span><span class=\"o\">:<\/span> <span class=\"nf\">solve<\/span><span class=\"p\">(<\/span><span class=\"nv\">ans3<\/span>, <span class=\"nv\">N<\/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[16]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{24}$}\\left[ N=\\frac{{\\it Nmax}\\,e^{\\gamma_{2}\\,t+{\\it \\%c}\\,\\gamma_{2}}}{e^{\\gamma_{2}\\,t+{\\it \\%c}\\,\\gamma_{2}}-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>\u6700\u5f8c\u306b\u521d\u671f\u6761\u4ef6\uff0c$t = t_0$ \u3067 $N = N_0$ \u3092\u5165\u308c\u308b\u3002\uff08\u4e00\u822c\u306b <code>solve()<\/code> \u306e\u89e3\u306f\u8907\u6570\u304b\u3082\u77e5\u308c\u306a\u3044\u306e\u3067\u30ea\u30b9\u30c8\u5f62\u5f0f <code>[... ]<\/code> \u3067\u51fa\u529b\u3055\u308c\u308b\u30021\u500b\u76ee\u3092\u3068\u308b\u306e\u304c <code>sol[1]<\/code>\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[17]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">sol1<\/span><span class=\"o\">:<\/span> <span class=\"nf\">ic1<\/span><span class=\"p\">(<\/span><span class=\"nv\">sol<\/span><span class=\"p\">[<\/span>1<span class=\"p\">]<\/span>, <span class=\"nv\">t<\/span> <span class=\"o\">=<\/span> <span class=\"nv\">t0<\/span>, <span class=\"nv\">N<\/span> <span class=\"o\">=<\/span> <span class=\"nv\">N0<\/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[17]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{25}$}N=\\frac{N_{0}\\,{\\it Nmax}\\,e^{\\gamma_{2}\\,t}}{\\left({\\it Nmax}-N_{0}\\right)\\,e^{\\gamma_{2}\\,t_{0}}+N_{0}\\,e^{\\gamma_{2}\\,t}}\\]<\/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[18]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">define<\/span><span class=\"p\">(<\/span><span class=\"nf\">Nv<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)<\/span>, <span class=\"nf\">rhs<\/span><span class=\"p\">(<\/span><span class=\"nv\">sol1<\/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[18]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{26}$}{\\it Nv}\\left(t\\right):=\\frac{N_{0}\\,{\\it Nmax}\\,e^{\\gamma_{2}\\,t}}{\\left({\\it Nmax}-N_{0}\\right)\\,e^{\\gamma_{2}\\,t_{0}}+N_{0}\\,e^{\\gamma_{2}\\,t}}\\]<\/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=\"\u53c2\u8003\uff1a\u7c73\u56fd\u306e\u4eba\u53e3\u3068\u30f4\u30a7\u30a2\u30d5\u30eb\u30b9\u30c8\u30e2\u30c7\u30eb\">\u53c2\u8003\uff1a\u7c73\u56fd\u306e\u4eba\u53e3\u3068\u30f4\u30a7\u30a2\u30d5\u30eb\u30b9\u30c8\u30e2\u30c7\u30eb<\/h4>\n<p>\u521d\u671f\u6761\u4ef6\u3092 $t_0 = 1790$\uff08\u5e74\uff09\u306e\u3068\u304d $N(t_0) = N_0 = 3.9$\uff08\u767e\u4e07\u4eba\uff09\u3068\u3057\u307e\u3059\u3002<\/p>\n<p>$t_1$\uff08\u5e74\uff09\u3068 $t_2$\uff08\u5e74\uff09\u306e\u5024\u3092\u4f7f\u3063\u3066 $\\gamma$ \u3068 $N_{\\rm max}$ \u3092\u6c42\u3081\u307e\u3059\u3002\u6307\u6570\u95a2\u6570\u3092\u542b\u3080\u9023\u7acb\u65b9\u7a0b\u5f0f\u306f\u306a\u304b\u306a\u304b\u89e3\u3044\u3066\u304f\u308c\u306a\u3044\u306e\u3067\uff0c\u7c21\u5358\u306a\u4ee3\u6570\u65b9\u7a0b\u5f0f\u306e\u5f62\u306b\u3057\u3066\u89e3\u304d\u307e\u3059\u3002<\/p>\n<p>$Nv(t)$ \u306e\u5206\u5b50\u5206\u6bcd\u3092 $e^{\\gamma t}$ \u3067\u5272\u308a\uff0c\u3055\u3089\u306b<\/p>\n<p>\\begin{eqnarray}<br \/>\nn_0 &amp;\\equiv&amp; \\frac{N_0}{N_{\\rm max}}<br \/>\n\\end{eqnarray}<\/p>\n<p>\u3068\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[19]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">t0<\/span><span class=\"o\">:<\/span> <span class=\"nv\">usa<\/span><span class=\"p\">[<\/span>1<span class=\"p\">][<\/span>1<span class=\"p\">]<\/span>; <span class=\"cm\">\/* 1\u884c\u76ee\u306e1\u5217\u76ee *\/<\/span>\r\n<span class=\"nv\">N0<\/span><span class=\"o\">:<\/span> <span class=\"nv\">usa<\/span><span class=\"p\">[<\/span>1<span class=\"p\">][<\/span>2<span class=\"p\">]<\/span>; <span class=\"cm\">\/* 1\u884c\u76ee\u306e2\u5217\u76ee *\/<\/span>\r\n\r\n<span class=\"cm\">\/* n0 = N0\/Nmax \u3068\u304a\u3044\u305f *\/<\/span>\r\n\r\n<span class=\"nf\">Nv1<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)<\/span><span class=\"o\">:=<\/span> <span class=\"nv\">N0<\/span><span class=\"o\">\/<\/span><span class=\"p\">(<\/span><span class=\"nv\">n0<\/span> <span class=\"o\">+<\/span> <span class=\"p\">(<\/span>1<span class=\"o\">-<\/span><span class=\"nv\">n0<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">exp<\/span><span class=\"p\">(<\/span><span class=\"nv\">gamma2<\/span><span class=\"o\">*<\/span><span class=\"p\">(<\/span><span class=\"nv\">t0<\/span><span class=\"o\">-<\/span><span class=\"nv\">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[19]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{27}$}1790\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[19]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{28}$}3.9\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[19]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{29}$}{\\it Nv}_{1}\\left(t\\right):=\\frac{N_{0}}{n_{0}+\\left(1-n_{0}\\right)\\,\\exp \\left(\\gamma_{2}\\,\\left(t_{0}-t\\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 = 1850$\uff08\u5e74\uff09\u3068 $t_2 = 1910$\uff08\u5e74\uff09\u306e\u5024\u3092\u5165\u308c\u3066\u9023\u7acb\u65b9\u7a0b\u5f0f\u306e\u5f62\u306b\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[20]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">t1<\/span><span class=\"o\">:<\/span> <span class=\"nv\">usa<\/span><span class=\"p\">[<\/span>7<span class=\"p\">][<\/span>1<span class=\"p\">]<\/span>;\r\n<span class=\"nv\">t2<\/span><span class=\"o\">:<\/span> <span class=\"nv\">usa<\/span><span class=\"p\">[<\/span>13<span class=\"p\">][<\/span>1<span class=\"p\">]<\/span>;\r\n\r\n<span class=\"nv\">eq1<\/span><span class=\"o\">:<\/span> <span class=\"nv\">usa<\/span><span class=\"p\">[<\/span>7<span class=\"p\">][<\/span>2<span class=\"p\">]<\/span> <span class=\"o\">=<\/span> <span class=\"nf\">Nv1<\/span><span class=\"p\">(<\/span><span class=\"nv\">t1<\/span><span class=\"p\">)<\/span>;\r\n<span class=\"nv\">eq2<\/span><span class=\"o\">:<\/span> <span class=\"nv\">usa<\/span><span class=\"p\">[<\/span>13<span class=\"p\">][<\/span>2<span class=\"p\">]<\/span> <span class=\"o\">=<\/span> <span class=\"nf\">Nv1<\/span><span class=\"p\">(<\/span><span class=\"nv\">t2<\/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[20]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{30}$}1850\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[20]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{31}$}1910\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[20]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{32}$}23.2=\\frac{3.9}{n_{0}+e^ {- 60\\,\\gamma_{2} }\\,\\left(1-n_{0}\\right)}\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[20]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{33}$}92.0=\\frac{3.9}{n_{0}+e^ {- 120\\,\\gamma_{2} }\\,\\left(1-n_{0}\\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\u306a\u9023\u7acb\u65b9\u7a0b\u5f0f\u306b\u3059\u308b\u305f\u3081\uff0c\u3055\u3089\u306b $T \\equiv e^{-60\\, \\gamma_2}$ \u3068\u3057\uff0c$n_0$ \u3068 $T$ \u306b\u95a2\u3059\u308b\u30b7\u30f3\u30d7\u30eb\u306a\u9023\u7acb\u65b9\u7a0b\u5f0f\u306e\u5f62\u306b\u3057\u3066\uff0c<code>solve()<\/code> \u3067\u89e3\u304d\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[21]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">kill<\/span><span class=\"p\">(<\/span><span class=\"nv\">T<\/span>, <span class=\"nv\">n0<\/span><span class=\"p\">)<\/span>$\r\n\r\n<span class=\"nv\">eq1<\/span><span class=\"o\">:<\/span> <span class=\"nf\">subst<\/span><span class=\"p\">(<\/span><span class=\"nf\">exp<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span>60<span class=\"o\">*<\/span><span class=\"nv\">gamma2<\/span><span class=\"p\">)<\/span><span class=\"o\">=<\/span><span class=\"nv\">T<\/span>, <span class=\"nv\">eq1<\/span><span class=\"p\">)<\/span>;\r\n<span class=\"nv\">eq2<\/span><span class=\"o\">:<\/span> <span class=\"nf\">subst<\/span><span class=\"p\">(<\/span><span class=\"nf\">exp<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span>120<span class=\"o\">*<\/span><span class=\"nv\">gamma2<\/span><span class=\"p\">)<\/span><span class=\"o\">=<\/span><span class=\"nv\">T<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span>, <span class=\"nv\">eq2<\/span><span class=\"p\">)<\/span>;\r\n\r\n<span class=\"nv\">ans<\/span><span class=\"o\">:<\/span> <span class=\"nf\">solve<\/span><span class=\"p\">([<\/span><span class=\"nv\">eq1<\/span>, <span class=\"nv\">eq2<\/span><span class=\"p\">]<\/span>, <span class=\"p\">[<\/span><span class=\"nv\">n0<\/span>, <span class=\"nv\">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\"><\/div>\n<div class=\"output_subarea output_stream output_stdout output_text\">\n<pre>rat: replaced 23.2 by 116\/5 = 23.2\r\n\r\nrat: replaced -3.9 by -39\/10 = -3.9\r\n\r\nrat: replaced 92.0 by 92\/1 = 92.0\r\n\r\nrat: replaced -3.9 by -39\/10 = -3.9\r\n<\/pre>\n<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[21]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{35}$}23.2=\\frac{3.9}{n_{0}+T\\,\\left(1-n_{0}\\right)}\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[21]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{36}$}92.0=\\frac{3.9}{n_{0}+T^2\\,\\left(1-n_{0}\\right)}\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[21]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{37}$}\\left[ \\left[ n_{0}=\\frac{87477}{4371112} , T=\\frac{3354}{22195} \\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>$n_0, \\ T$ \u304b\u3089\u3082\u3068\u306e $N_{\\rm max}, \\ \\gamma_2$ \u306e\u5024\u306b\u306a\u304a\u3059\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[22]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">Nmax<\/span><span class=\"o\">:<\/span> <span class=\"nf\">float<\/span><span class=\"p\">(<\/span><span class=\"nv\">N0<\/span><span class=\"o\">\/<\/span><span class=\"nf\">ev<\/span><span class=\"p\">(<\/span><span class=\"nv\">n0<\/span>, <span class=\"nv\">ans<\/span><span class=\"p\">[<\/span>1<span class=\"p\">][<\/span>1<span class=\"p\">]))<\/span>;\r\n<span class=\"nv\">gamma2<\/span><span class=\"o\">:<\/span> <span class=\"nf\">float<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span>1<span class=\"o\">\/<\/span>60<span class=\"o\">*<\/span><span class=\"nf\">log<\/span><span class=\"p\">(<\/span><span class=\"nf\">ev<\/span><span class=\"p\">(<\/span><span class=\"nv\">T<\/span>, <span class=\"nv\">ans<\/span><span class=\"p\">[<\/span>1<span class=\"p\">][<\/span>2<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[22]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{38}$}194.8779313419527\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[22]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{39}$}0.03149522291474969\\]<\/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-maxima\">\n<pre><span class=\"cm\">\/* plot2d() \u306b\u3088\u308b\u30b0\u30e9\u30d5 *\/<\/span>\r\n\r\n<span class=\"nf\">plot2d<\/span><span class=\"p\">([<\/span><span class=\"nf\">Nm<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)<\/span>, <span class=\"p\">[<\/span><span class=\"nv\">discrete<\/span>, <span class=\"nv\">usa<\/span><span class=\"p\">]<\/span>, <span class=\"nf\">Nv<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)]<\/span>, <span class=\"p\">[<\/span><span class=\"nv\">t<\/span>, <span class=\"mi\">1790<\/span>, 1940<span class=\"p\">]<\/span>, \r\n       <span class=\"p\">[<\/span><span class=\"nv\">point_type<\/span>,<span class=\"nv\">bullet<\/span><span class=\"p\">]<\/span>,\r\n       <span class=\"p\">[<\/span><span class=\"nv\">style<\/span>, <span class=\"nv\">lines<\/span>, <span class=\"p\">[<\/span><span class=\"nv\">points<\/span>, 1<span class=\"o\">.<\/span>5<span class=\"p\">]<\/span>, <span class=\"nv\">lines<\/span><span class=\"p\">]<\/span>, \r\n       <span class=\"p\">[<\/span><span class=\"nv\">color<\/span>, <span class=\"nv\">red<\/span>, <span class=\"nv\">blue<\/span>, <span class=\"nv\">green<\/span><span class=\"p\">]<\/span>,\r\n       <span class=\"p\">[<\/span><span class=\"nv\">legend<\/span>, <span class=\"s\">\"\u30de\u30eb\u30b5\u30b9\u30e2\u30c7\u30eb\"<\/span>, <span class=\"s\">\"\u30a2\u30e1\u30ea\u30ab\u306e\u4eba\u53e3\"<\/span>, <span class=\"s\">\"\u30f4\u30a7\u30a2\u30d5\u30eb\u30b9\u30c8\u30e2\u30c7\u30eb\"<\/span><span class=\"p\">]<\/span>, \r\n       <span class=\"p\">[<\/span><span class=\"nv\">xlabel<\/span>, <span class=\"s\">\"\u5e74\"<\/span><span class=\"p\">]<\/span>, <span class=\"p\">[<\/span><span class=\"nv\">ylabel<\/span>, <span class=\"s\">\"\u4eba\u53e3\uff08\u5358\u4f4d\uff1a\u767e\u4e07\u4eba\uff09\"<\/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 \">\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-6122\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/mmathc103.svg\" alt=\"\" width=\"640\" height=\"481\" \/><\/p>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[24]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"cm\">\/* draw2d() \u306b\u3088\u308b\u30b0\u30e9\u30d5 *\/<\/span>\r\n<span class=\"nf\">draw2d<\/span><span class=\"p\">(<\/span>\r\n  <span class=\"nv\">color<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"red\"<\/span>,\r\n  <span class=\"nv\">key<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u30de\u30eb\u30b5\u30b9\u30e2\u30c7\u30eb\"<\/span>, \r\n  <span class=\"nf\">explicit<\/span><span class=\"p\">(<\/span><span class=\"nf\">Nm<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">t<\/span>, <span class=\"mi\">1790<\/span>, <span class=\"mi\">1940<\/span><span class=\"p\">)<\/span>, \r\n  \r\n  <span class=\"nv\">color<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"blue\"<\/span>, \r\n  <span class=\"nv\">point_type<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">7<\/span>,\r\n  <span class=\"nv\">point_size<\/span> <span class=\"o\">=<\/span> <span class=\"mf\">0.7<\/span>,\r\n  <span class=\"nv\">key<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u30a2\u30e1\u30ea\u30ab\u306e\u4eba\u53e3\"<\/span>, \r\n  <span class=\"nf\">points<\/span><span class=\"p\">(<\/span><span class=\"nv\">usa<\/span><span class=\"p\">)<\/span>, \r\n  \r\n  <span class=\"nv\">color<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"green\"<\/span>,\r\n  <span class=\"nv\">key<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u30f4\u30a7\u30a2\u30d5\u30eb\u30b9\u30c8\u30e2\u30c7\u30eb\"<\/span>,\r\n  <span class=\"nf\">explicit<\/span><span class=\"p\">(<\/span><span class=\"nf\">Nv<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">t<\/span>, <span class=\"mi\">1790<\/span>, <span class=\"mi\">1940<\/span><span class=\"p\">)<\/span>, \r\n  \r\n  <span class=\"nv\">xlabel<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u5e74\"<\/span>, <span class=\"nv\">ylabel<\/span> <span class=\"o\">=<\/span> <span class=\"s\">\"\u4eba\u53e3\uff08\u5358\u4f4d\uff1a\u767e\u4e07\u4eba\uff09\"<\/span>,\r\n  <span class=\"nv\">yrange<\/span> <span class=\"o\">=<\/span> <span class=\"p\">[<\/span><span class=\"o\">-<\/span><span class=\"mi\">20<\/span>, 410<span class=\"p\">]<\/span>, <span class=\"nv\">xaxis<\/span> <span class=\"o\">=<\/span> <span class=\"no\">true<\/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-6123\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/mmathc104.svg\" alt=\"\" width=\"640\" height=\"481\" \/><\/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=\"1\u968e\u7dda\u5f62\u5fae\u5206\u65b9\u7a0b\u5f0f\u3068\u7a4d\u5206\u56e0\u5b50\u6cd5\">1\u968e\u7dda\u5f62\u5fae\u5206\u65b9\u7a0b\u5f0f\u3068\u7a4d\u5206\u56e0\u5b50\u6cd5<\/h3>\n<p>Maxima \u3067\u306f\uff0c1\u968e\u7dda\u5f62\u5fae\u5206\u65b9\u7a0b\u5f0f\u306f\uff0c\u7279\u306b\u7a4d\u5206\u56e0\u5b50\u6cd5\u3092 coding \u3057\u306a\u304f\u3066\u3082\uff0c<code>ode2()<\/code> \u95a2\u6570\u3067\u89e3\u304f\u3053\u3068\u304c\u3067\u304d\u307e\u3059\u3002<\/p>\n<h4 id=\"\u4f8b\u984c\">\u4f8b\u984c<\/h4>\n<p>$$ \\frac{dy}{dx} + \\frac{y}{x} = \\frac{\\sin x}{x} $$<\/p>\n<p>Maxima \u3067\u306e\u4e09\u89d2\u95a2\u6570 $\\sin x, \\cos x, \\tan x$ \u306f <code>sin(x)<\/code>, <code>cos(x)<\/code>, <code>tan(x)<\/code> \u3068\u66f8\u304d\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[25]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">eq<\/span><span class=\"o\">:<\/span> <span class=\"o\">'<\/span><span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"nv\">y<\/span><span class=\"o\">\/<\/span><span class=\"nv\">x<\/span> <span class=\"o\">=<\/span> <span class=\"nf\">sin<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"nv\">x<\/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\">\\[\\tag{${\\it \\%o}_{44}$}\\frac{d}{d\\,x}\\,y+\\frac{y}{x}=\\frac{\\sin x}{x}\\]<\/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[26]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq<\/span>, <span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/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\">\\[\\tag{${\\it \\%o}_{45}$}y=\\frac{{\\it \\%c}-\\cos x}{x}\\]<\/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=\"\u53c2\u8003\uff1a\u3042\u3048\u3066\u7a4d\u5206\u56e0\u5b50\u6cd5\u306e\u516c\u5f0f\u3067...\">\u53c2\u8003\uff1a\u3042\u3048\u3066\u7a4d\u5206\u56e0\u5b50\u6cd5\u306e\u516c\u5f0f\u3067&#8230;<\/h4>\n<p>\u3053\u306e\u4f8b\u984c\u3092\u3042\u3048\u3066\u7a4d\u5206\u56e0\u5b50\u6cd5\u306e\u516c\u5f0f\u306b\u5f93\u3063\u3066\u3084\u3063\u3066\u307f\u308b\u3002<\/p>\n<p>\u4e00\u822c\u7684\u306a1\u968e\u7dda\u5f62\u5fae\u5206\u65b9\u7a0b\u5f0f\uff08\u975e\u540c\u6b21\u65b9\u7a0b\u5f0f\u306e\u683c\u597d\u3057\u3066\u3044\u308b\u3082\u306e\uff09\u3092<\/p>\n<p>$$\\frac{dy}{dx} + P(x) y = Q(x)$$<\/p>\n<p>\u3068\u66f8\u304f\u3068\uff0c\u7a4d\u5206\u56e0\u5b50 $g(x)$ \u306f<\/p>\n<p>$$g(x) = \\exp \\left\\{\\int^x P(x&#8217;)\\,dx&#8217; \\right\\}$$<\/p>\n<p>\u89e3\u306f<\/p>\n<p>$$y = \\frac{1}{g(x)} \\left\\{\\int^x g(x&#8217;) Q(x&#8217;) \\, dx&#8217; + C \\right\\}$$<\/p>\n<p>\u3053\u306e\u4f8b\u984c\u3067\u306f<\/p>\n<p>$$P(x) = \\frac{1}{x}, \\ Q(x) = \\frac{\\sin x}{x}$$<\/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-maxima\">\n<pre><span class=\"nf\">P<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">:=<\/span> 1<span class=\"o\">\/<\/span><span class=\"nv\">x<\/span>;\r\n<span class=\"nf\">Q<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">:=<\/span> <span class=\"nf\">sin<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"nv\">x<\/span>;\r\n\r\n<span class=\"nf\">define<\/span><span class=\"p\">(<\/span><span class=\"nf\">g<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span>, <span class=\"nf\">exp<\/span><span class=\"p\">(<\/span><span class=\"nf\">integrate<\/span><span class=\"p\">(<\/span><span class=\"nf\">P<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)))<\/span>;\r\n\r\n<span class=\"nv\">y<\/span> <span class=\"o\">=<\/span> 1<span class=\"o\">\/<\/span><span class=\"nf\">g<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">*<\/span> <span class=\"p\">(<\/span><span class=\"nf\">integrate<\/span><span class=\"p\">(<\/span><span class=\"nf\">g<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">Q<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"nv\">C<\/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\">\\[\\tag{${\\it \\%o}_{46}$}P\\left(x\\right):=\\frac{1}{x}\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[27]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{47}$}Q\\left(x\\right):=\\frac{\\sin x}{x}\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[27]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{48}$}g\\left(x\\right):=x\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[27]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{49}$}y=\\frac{C-\\cos x}{x}\\]<\/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\u3082\u7c21\u5358\u306a\u5b9a\u6570\u4fc2\u65702\u968e\u5fae\u5206\u65b9\u7a0b\u5f0f\">\u6700\u3082\u7c21\u5358\u306a\u5b9a\u6570\u4fc2\u65702\u968e\u5fae\u5206\u65b9\u7a0b\u5f0f<\/h3>\n<p>\u6700\u3082\u7c21\u5358\u306a\u5b9a\u6570\u4fc2\u65702\u968e\u7dda\u5f62\u5fae\u5206\u65b9\u7a0b\u5f0f $y^{&#8221;} + y = 0$ \u3068 $y^{&#8221;} &#8211; y = 0$ \u3092\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[28]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">eq1<\/span><span class=\"o\">:<\/span> <span class=\"o\">'<\/span><span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span>, <span class=\"mi\">2<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"nv\">y<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/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\">\\[\\tag{${\\it \\%o}_{50}$}\\frac{d^2}{d\\,x^2}\\,y+y=0\\]<\/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[29]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq1<\/span>, <span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/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[29]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{51}$}y={\\it \\%k}_{1}\\,\\sin x+{\\it \\%k}_{2}\\,\\cos x\\]<\/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[30]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">eq2<\/span><span class=\"o\">:<\/span> <span class=\"o\">'<\/span><span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span>, <span class=\"mi\">2<\/span><span class=\"p\">)<\/span> <span class=\"o\">-<\/span> <span class=\"nv\">y<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/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\">\\[\\tag{${\\it \\%o}_{52}$}\\frac{d^2}{d\\,x^2}\\,y-y=0\\]<\/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[31]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq2<\/span>, <span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/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[31]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{53}$}y={\\it \\%k}_{1}\\,e^{x}+{\\it \\%k}_{2}\\,e^ {- x }\\]<\/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=\"\u53c2\u8003\uff1a\u810a\u9ac4\u53cd\u5c04\u306b\u3088\u3089\u305a\u306b\u89e3\u304f\">\u53c2\u8003\uff1a\u810a\u9ac4\u53cd\u5c04\u306b\u3088\u3089\u305a\u306b\u89e3\u304f<\/h4>\n<p>$y^{&#8221;} + y = 0$ \u3092\u810a\u9ac4\u53cd\u5c04\u306b\u3088\u3089\u305a\u306b\u89e3\u304f\u3002\u4e21\u8fba\u306b $2 y&#8217;$ \u3092\u304b\u3051\u3066\u6574\u7406\u3059\u308b\u3068<\/p>\n<p>$$\\left(\\left(y&#8217;\\right)^2 + y^2\\right)&#8217; = 0$$<\/p>\n<p>\u5fae\u5206\u3057\u3066\u30bc\u30ed\u3068\u3044\u3046\u3053\u3068\u306f\uff0c\u304b\u3063\u3053\u306e\u4e2d\u8eab\u306f\u5b9a\u6570\u3067\u3042\u308b\u3068\u3044\u3046\u3053\u3068\u306a\u306e\u3067\uff0c<\/p>\n<p>$$ \\left(y&#8217;\\right)^2 + y^2 = \\mbox{const.} \\equiv a^2$$<\/p>\n<p>\u3068\u304a\u3051\u308b\u3002\u3064\u307e\u308a<\/p>\n<p>$$\\frac{dy}{dx} = \\pm \\sqrt{a^2 &#8211; y^2}$$<\/p>\n<p>\u3092\u89e3\u3051\u3070\u3088\u3044\u3002\u307e\u305a\u306f $+$ \u306e\u5f0f\u304b\u3089&#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-maxima\">\n<pre><span class=\"nf\">assume<\/span><span class=\"p\">(<\/span><span class=\"nv\">a<\/span> <span class=\"o\">&gt;<\/span> <span class=\"mi\">0<\/span><span class=\"p\">)<\/span>$\r\n<span class=\"nv\">eq1<\/span><span class=\"o\">:<\/span> <span class=\"o\">'<\/span><span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">=<\/span> <span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"nv\">a<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">-<\/span> <span class=\"nv\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/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[32]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{55}$}\\frac{d}{d\\,x}\\,y=\\sqrt{a^2-y^2}\\]<\/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-maxima\">\n<pre><span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq1<\/span>, <span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/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\">\\[\\tag{${\\it \\%o}_{56}$}\\arcsin \\left(\\frac{y}{a}\\right)=x+{\\it \\%c}\\]<\/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>$y = $ \u306e\u5f62\u306b\u89e3\u304d\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[34]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">solve<\/span><span class=\"p\">(<\/span><span class=\"nv\">%<\/span>, <span class=\"nv\">y<\/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\">\\[\\tag{${\\it \\%o}_{57}$}\\left[ y=a\\,\\sin \\left(x+{\\it \\%c}\\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>$-$ \u306e\u5f0f\u306b\u3064\u3044\u3066\u3082\u540c\u69d8\u306b&#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-maxima\">\n<pre><span class=\"nv\">eq2<\/span><span class=\"o\">:<\/span> <span class=\"o\">'<\/span><span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">=<\/span> <span class=\"o\">-<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"nv\">a<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">-<\/span> <span class=\"nv\">y<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>;\r\n<span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq2<\/span>, <span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span>;\r\n<span class=\"nf\">solve<\/span><span class=\"p\">(<\/span><span class=\"nv\">%<\/span>, <span class=\"nv\">y<\/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[35]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{58}$}\\frac{d}{d\\,x}\\,y=-\\sqrt{a^2-y^2}\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[35]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{59}$}-\\arcsin \\left(\\frac{y}{a}\\right)=x+{\\it \\%c}\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[35]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{60}$}\\left[ y=-a\\,\\sin \\left(x+{\\it \\%c}\\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>\u3068\u3044\u3046\u3053\u3068\u3067 2 \u3064\u306e\u89e3\u3092\u307e\u3068\u3081\u308b\u3068\uff0c\u7a4d\u5206\u5b9a\u6570 $C_1, C_2$ \u3092\u4f7f\u3063\u3066<\/p>\n<p>\\begin{eqnarray}<br \/>\ny &amp;=&amp; C_1 \\sin(x + C_2) \\\\<br \/>\n&amp;=&amp; C_1 \\left(\\sin x \\cos C_2 + \\cos x \\sin C_2 \\right) \\\\<br \/>\n&amp;=&amp; (C_1 \\cos C_2) \\sin x + (C_1 \\sin C_2) \\cos x \\\\<br \/>\n&amp;\\equiv&amp; A \\cos x + B \\sin x<br \/>\n\\end{eqnarray}<\/p>\n<p>\u306e\u3088\u3046\u306b\u66f8\u3051\u308b\u3068\u3044\u3046\u3053\u3068\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=\"\u6700\u3082\u7c21\u5358\u306a\u5b9a\u6570\u4fc2\u65702\u968e\u5fae\u5206\u65b9\u7a0b\u5f0f\uff1a\u7d9a\u304d\">\u6700\u3082\u7c21\u5358\u306a\u5b9a\u6570\u4fc2\u65702\u968e\u5fae\u5206\u65b9\u7a0b\u5f0f\uff1a\u7d9a\u304d<\/h3>\n<p>$y^{&#8221;} + K y = 0$ \u3042\u308b\u3044\u306f\u79fb\u9805\u3057\u3066 $y^{&#8221;} = \u2013 K y$ \u3092\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[36]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">eq<\/span><span class=\"o\">:<\/span> <span class=\"o\">'<\/span><span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span>, <span class=\"mi\">2<\/span><span class=\"p\">)<\/span> <span class=\"o\">=<\/span> <span class=\"o\">-<\/span> <span class=\"nv\">K<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">y<\/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\">\\[\\tag{${\\it \\%o}_{61}$}\\frac{d^2}{d\\,x^2}\\,y=-K\\,y\\]<\/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=\"$K-&gt;-0$-\u306e\u5834\u5408\">$K &gt; 0$ \u306e\u5834\u5408<\/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[37]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq<\/span>, <span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span>;\r\n\r\n<span class=\"cm\">\/* positive; \u3068\u56de\u7b54\u3059\u308b *\/<\/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>\u001bXIs K positive, negative or zero?\r\n\u001b\\Is K positive, negative or zero?\r\npositive;\r\n\r\n<\/pre>\n<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[37]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{62}$}y={\\it \\%k}_{1}\\,\\sin \\left(\\sqrt{K}\\,x\\right)+{\\it \\%k}_{2}\\,\\cos \\left(\\sqrt{K}\\,x\\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<h4 id=\"$K-&lt;-0$-\u306e\u5834\u5408\">$K &lt; 0$ \u306e\u5834\u5408<\/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[38]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">ans3<\/span><span class=\"o\">:<\/span> <span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq<\/span>, <span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span>;\r\n\r\n<span class=\"cm\">\/* negative; \u3068\u56de\u7b54\u3059\u308b *\/<\/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>\u001bXIs K positive, negative or zero?\r\n\u001b\\Is K positive, negative or zero?\r\nnegative;\r\n\r\n<\/pre>\n<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[38]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{63}$}y={\\it \\%k}_{1}\\,e^{i\\,\\sqrt{K}\\,x}+{\\it \\%k}_{2}\\,e^ {- i\\,\\sqrt{K}\\,x }\\]<\/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\u8a18\u306e\u89e3\u306f\u9593\u9055\u3044\u3067\u306f\u306a\u3044\u304b\u3082\u77e5\u308c\u306a\u3044\u304c\uff0c\u3084\u3084\u3053\u3057\u3044\u3002<\/p>\n<p>$K \\equiv -|K|$ \u3068\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[39]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">ev<\/span><span class=\"p\">(<\/span><span class=\"nv\">%<\/span>, <span class=\"nv\">K<\/span><span class=\"o\">=-<\/span><span class=\"nf\">abs<\/span><span class=\"p\">(<\/span><span class=\"nv\">K<\/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[39]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{64}$}y={\\it \\%k}_{2}\\,e^{\\sqrt{\\left| K\\right| }\\,x}+{\\it \\%k}_{1}\\,e^ {- \\sqrt{\\left| K\\right| }\\,x }\\]<\/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=\"$K-=-0$-\u306e\u5834\u5408\">$K = 0$ \u306e\u5834\u5408<\/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[40]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq<\/span>, <span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span>;\r\n\r\n<span class=\"cm\">\/* zero; \u3068\u56de\u7b54\u3059\u308b *\/<\/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>\u001bXIs K positive, negative or zero?\r\n\u001b\\Is K positive, negative or zero?\r\nzero; \r\n\r\n<\/pre>\n<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[40]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{65}$}y={\\it \\%k}_{2}\\,x+{\\it \\%k}_{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<h3 id=\"\u4eba\u985e\u306e\u81f3\u5b9d\uff1a\u30aa\u30a4\u30e9\u30fc\u306e\u516c\u5f0f\">\u4eba\u985e\u306e\u81f3\u5b9d\uff1a\u30aa\u30a4\u30e9\u30fc\u306e\u516c\u5f0f<\/h3>\n<p>Maxima \u304c\u30aa\u30a4\u30e9\u30fc\u306e\u516c\u5f0f\u3092\u77e5\u3063\u3066\u3044\u308b\u304b\u78ba\u8a8d\u3059\u308b\u3002<\/p>\n<p>Maxima \u3067\u306f\u865a\u6570\u5358\u4f4d $i$ \u306f <code>%i<\/code> \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[41]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">exp<\/span><span class=\"p\">(<\/span><span class=\"nv\">%i<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">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[41]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{66}$}e^{i\\,\\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><code>demoivre()<\/code> \u95a2\u6570\u306f\u8907\u7d20\u6307\u6570\u95a2\u6570\u3092\u4e09\u89d2\u95a2\u6570\u3067\u8868\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[42]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">exp<\/span><span class=\"p\">(<\/span><span class=\"nv\">%i<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">theta<\/span><span class=\"p\">)<\/span> <span class=\"o\">=<\/span> <span class=\"nf\">demoivre<\/span><span class=\"p\">(<\/span><span class=\"nf\">exp<\/span><span class=\"p\">(<\/span><span class=\"nv\">%i<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">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[42]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{67}$}e^{i\\,\\vartheta}=i\\,\\sin \\vartheta+\\cos \\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<h4 id=\"\u30aa\u30a4\u30e9\u30fc\u306e\u7b49\u5f0f\">\u30aa\u30a4\u30e9\u30fc\u306e\u7b49\u5f0f<\/h4>\n<p>$\\theta = \\pi$ \u306e\u3068\u304d\uff0c$e^{i \\pi} + 1 = 0$\u3002$\\pi$ \u306f Maxima \u3067\u306f <code>%pi<\/code> \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[43]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">exp<\/span><span class=\"p\">(<\/span><span class=\"nv\">%i<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">%pi<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"mi\">1<\/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[43]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{68}$}0\\]<\/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=\"\u30aa\u30a4\u30e9\u30fc\u306e\u516c\u5f0f\u304b\u3089\u307f\u305f\u4e09\u89d2\u95a2\u6570\u3068\u53cc\u66f2\u7dda\u95a2\u6570\u306e\u95a2\u4fc2\">\u30aa\u30a4\u30e9\u30fc\u306e\u516c\u5f0f\u304b\u3089\u307f\u305f\u4e09\u89d2\u95a2\u6570\u3068\u53cc\u66f2\u7dda\u95a2\u6570\u306e\u95a2\u4fc2<\/h4>\n<p>\u4e09\u89d2\u95a2\u6570 $\\cos x, \\sin x$ \u3092 <code>exponentialize()<\/code> \u3067\uff08\u8907\u7d20\uff09\u6307\u6570\u95a2\u6570\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[44]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">cos<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">=<\/span> <span class=\"nf\">exponentialize<\/span><span class=\"p\">(<\/span><span class=\"nf\">cos<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">))<\/span>;\r\n<span class=\"nf\">sin<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">=<\/span> <span class=\"nf\">exponentialize<\/span><span class=\"p\">(<\/span><span class=\"nf\">sin<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/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[44]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{69}$}\\cos x=\\frac{e^{i\\,x}+e^ {- i\\,x }}{2}\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[44]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{70}$}\\sin x=-\\frac{i\\,\\left(e^{i\\,x}-e^ {- i\\,x }\\right)}{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>\u53cc\u66f2\u7dda\u95a2\u6570 $\\cosh x, \\sinh x$ \u3082\u6307\u6570\u95a2\u6570\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[45]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">cosh<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">=<\/span> <span class=\"nf\">exponentialize<\/span><span class=\"p\">(<\/span><span class=\"nf\">cosh<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">))<\/span>;\r\n<span class=\"nf\">sinh<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">=<\/span> <span class=\"nf\">exponentialize<\/span><span class=\"p\">(<\/span><span class=\"nf\">sinh<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/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[45]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{71}$}\\cosh x=\\frac{e^{x}+e^ {- x }}{2}\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[45]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{72}$}\\sinh x=\\frac{e^{x}-e^ {- x }}{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>\u3055\u3066\uff0c\u4e09\u89d2\u95a2\u6570\u3084\u53cc\u66f2\u7dda\u95a2\u6570\u306e\u5909\u6570\uff08\u5f15\u6570\uff09\u304c\u865a\u6570\u3067\u3082\u3044\u3044\u306e\u3060\u3068\u62e1\u5f35\u3059\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[46]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">cosh<\/span><span class=\"p\">(<\/span><span class=\"nv\">%i<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">theta<\/span><span class=\"p\">)<\/span>;\r\n<span class=\"nf\">sinh<\/span><span class=\"p\">(<\/span><span class=\"nv\">%i<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">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[46]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{73}$}\\cos \\vartheta\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[46]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{74}$}i\\,\\sin \\vartheta\\]<\/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[47]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">cos<\/span><span class=\"p\">(<\/span><span class=\"nv\">%i<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">x<\/span><span class=\"p\">)<\/span>;\r\n<span class=\"nf\">sin<\/span><span class=\"p\">(<\/span><span class=\"nv\">%i<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">x<\/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[47]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{75}$}\\cosh x\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[47]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{76}$}i\\,\\sinh x\\]<\/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=\"\u4f8b\u984c\">\u4f8b\u984c<\/h4>\n<p>\u5fae\u5206\u65b9\u7a0b\u5f0f $y^{&#8221;} + K y = 0$ \u306e $K &gt; 0$ \u306e\u5834\u5408\u306e\u89e3\u306f<\/p>\n<p>\\begin{eqnarray}<br \/>\ny &amp;=&amp; A \\cos\\left( \\sqrt{K} x\\right) + \\frac{B}{\\sqrt{K}} \\sin\\left( \\sqrt{K} x\\right)<br \/>\n\\end{eqnarray}<\/p>\n<p>\u306e\u3088\u3046\u306b\u66f8\u3051\u308b\u3002\u3053\u306e\u3053\u3068\u3092\u4f7f\u3063\u3066\uff0c$K&lt;0$ \u306e\u5834\u5408\u3068 $K = 0$ \u306e\u5834\u5408\u306e\u89e3\u3092\u4e0a\u8a18\u306e\u89e3\u304b\u3089\u76f4\u63a5\u6c42\u3081\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<h5 id=\"$K-&gt;-0$-\u306e\u5834\u5408\">$K &gt; 0$ \u306e\u5834\u5408<\/h5>\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[48]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">y<\/span> <span class=\"o\">=<\/span> <span class=\"nv\">y1<\/span><span class=\"o\">:<\/span> <span class=\"nv\">A<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">cos<\/span><span class=\"p\">(<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"nv\">K<\/span><span class=\"p\">)<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"nv\">B<\/span><span class=\"o\">\/<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"nv\">K<\/span><span class=\"p\">)<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">sin<\/span><span class=\"p\">(<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"nv\">K<\/span><span class=\"p\">)<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">x<\/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[48]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{77}$}y=\\frac{B\\,\\sin \\left(\\sqrt{K}\\,x\\right)}{\\sqrt{K}}+A\\,\\cos \\left(\\sqrt{K}\\,x\\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=\"$K-&lt;-0$-\u306e\u5834\u5408\">$K &lt; 0$ \u306e\u5834\u5408<\/h5>\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[49]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">y<\/span> <span class=\"o\">=<\/span> <span class=\"nf\">ev<\/span><span class=\"p\">(<\/span><span class=\"nv\">y1<\/span>, <span class=\"nv\">K<\/span><span class=\"o\">=-<\/span><span class=\"nf\">abs<\/span><span class=\"p\">(<\/span><span class=\"nv\">K<\/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[49]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{78}$}y=\\frac{B\\,\\sinh \\left(\\sqrt{\\left| K\\right| }\\,x\\right)}{\\sqrt{\\left| K\\right| }}+A\\,\\cosh \\left(\\sqrt{\\left| K\\right| }\\,x\\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=\"$K-=-0$-\u306e\u5834\u5408\">$K = 0$ \u306e\u5834\u5408<\/h5>\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[50]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">y<\/span> <span class=\"o\">=<\/span> <span class=\"nf\">limit<\/span><span class=\"p\">(<\/span><span class=\"nv\">y1<\/span>, <span class=\"nv\">K<\/span>, <span class=\"mi\">0<\/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[50]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{79}$}y=B\\,x+A\\]<\/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=\"\u5b9a\u6570\u4fc2\u65702\u968e\u7dda\u5f62\u540c\u6b21\u65b9\u7a0b\u5f0f\">\u5b9a\u6570\u4fc2\u65702\u968e\u7dda\u5f62\u540c\u6b21\u65b9\u7a0b\u5f0f<\/h3>\n<p>\u5b9a\u6570\u4fc2\u65702\u968e\u7dda\u5f62\u5fae\u5206\u65b9\u7a0b\u5f0f\uff08\u540c\u6b21\u65b9\u7a0b\u5f0f\uff09\u306f\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u66f8\u3051\u308b\u3002<\/p>\n<p>$$ \\frac{d^2 y}{dx^2} + 2 b \\frac{dy}{dx} + cy = 0$$<\/p>\n<p>$b, c$ \u306f\u5b9a\u6570\u3002\u4e00\u822c\u89e3 $y$ \u306f\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u5834\u5408\u5206\u3051\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[51]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">eq<\/span><span class=\"o\">:<\/span> <span class=\"o\">'<\/span><span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span>, <span class=\"mi\">2<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"mi\">2<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">b<\/span> <span class=\"o\">*<\/span> <span class=\"o\">'<\/span><span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"nv\">c<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">y<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/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[51]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{80}$}\\frac{d^2}{d\\,x^2}\\,y+2\\,b\\,\\left(\\frac{d}{d\\,x}\\,y\\right)+c\\,y=0\\]<\/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=\"$c---b^2-&gt;-0$-\u306e\u5834\u5408\">$c &#8211; b^2 &gt; 0$ \u306e\u5834\u5408<\/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[52]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq<\/span>, <span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">factor<\/span>;\r\n\r\n<span class=\"cm\">\/* c - b**2 \u306f positive; \u3068\u7b54\u3048\u308b *\/<\/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>\u001bX        2\r\nIs c - b  positive, negative or zero?\r\n\u001b\\        2\r\nIs c - b  positive, negative or zero?\r\npositive;\r\n\r\n<\/pre>\n<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[52]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{81}$}y=e^ {- b\\,x }\\,\\left({\\it \\%k}_{1}\\,\\sin \\left(\\sqrt{c-b^2}\\,x\\right)+{\\it \\%k}_{2}\\,\\cos \\left(\\sqrt{c-b^2}\\,x\\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<h4 id=\"$c---b^2-&lt;-0$-\u306e\u5834\u5408\">$c &#8211; b^2 &lt; 0$ \u306e\u5834\u5408<\/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[53]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq<\/span>, <span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">factor<\/span>;\r\n\r\n<span class=\"cm\">\/* c - b**2 \u306f negative; \u3068\u7b54\u3048\u308b *\/<\/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>\u001bX        2\r\nIs c - b  positive, negative or zero?\r\n\u001b\\        2\r\nIs c - b  positive, negative or zero?\r\nnegative;\r\n\r\n<\/pre>\n<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[53]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{82}$}y=\\left({\\it \\%k}_{1}\\,e^{2\\,\\sqrt{b^2-c}\\,x}+{\\it \\%k}_{2}\\right)\\,e^{-\\sqrt{b^2-c}\\,x-b\\,x}\\]<\/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=\"$c---b^2-=-0$-\u306e\u5834\u5408\">$c &#8211; b^2 = 0$ \u306e\u5834\u5408<\/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[54]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nf\">ev<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq<\/span>, <span class=\"nv\">c<\/span><span class=\"o\">=<\/span><span class=\"nv\">b<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/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>\u001bXIs b zero or nonzero?\r\n\u001b\\Is b zero or nonzero?\r\nnonzero;\r\n\r\n<\/pre>\n<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[54]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{83}$}y=\\left({\\it \\%k}_{2}\\,x+{\\it \\%k}_{1}\\right)\\,e^ {- b\\,x }\\]<\/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=\"\u3072\u3068\u3064\u3060\u3051\u306e\u89e3\u3067\u3059\u307e\u3059\u306b\u306f\u2026\">\u3072\u3068\u3064\u3060\u3051\u306e\u89e3\u3067\u3059\u307e\u3059\u306b\u306f\u2026<\/h4>\n<p>\u4ee5\u4e0b\u306e\u8868\u73fe\u3072\u3068\u3064\u30673\u3064\u306e\u5834\u5408\u306e\u5168\u3066\u306b\u5bfe\u5fdc\u3057\u305f\u89e3\u306b\u306a\u308b\u3053\u3068\u3092\u793a\u3057\u3066\u304a\u304d\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[55]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">y<\/span><span class=\"p\">(<\/span><span class=\"nv\">mu<\/span><span class=\"p\">)<\/span><span class=\"o\">:=<\/span> <span class=\"nf\">exp<\/span><span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"nv\">b<\/span><span class=\"o\">*<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">*<\/span> <span class=\"p\">(<\/span><span class=\"nv\">A<\/span><span class=\"o\">*<\/span><span class=\"nf\">cos<\/span><span class=\"p\">(<\/span><span class=\"nv\">mu<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"nv\">B<\/span><span class=\"o\">\/<\/span><span class=\"nv\">mu<\/span><span class=\"o\">*<\/span><span class=\"nf\">sin<\/span><span class=\"p\">(<\/span><span class=\"nv\">mu<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">x<\/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[55]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{84}$}y\\left(\\mu\\right):=\\exp \\left(\\left(-b\\right)\\,x\\right)\\,\\left(A\\,\\cos \\left(\\mu\\,x\\right)+\\frac{B}{\\mu}\\,\\sin \\left(\\mu\\,x\\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>$c-b^2 &gt; 0$ \u306e\u5834\u5408\u306f $\\mu \\equiv \\sqrt{c-b^2}$ \u3068\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[56]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">y<\/span> <span class=\"o\">=<\/span> <span class=\"nf\">y<\/span><span class=\"p\">(<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"nv\">c<\/span><span class=\"o\">-<\/span><span class=\"nv\">b<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/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[56]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{85}$}y=e^ {- b\\,x }\\,\\left(\\frac{B\\,\\sin \\left(\\sqrt{c-b^2}\\,x\\right)}{\\sqrt{c-b^2}}+A\\,\\cos \\left(\\sqrt{c-b^2}\\,x\\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>$c-b^2 &lt; 0$ \u306e\u5834\u5408\u306f $\\mu \\equiv \\sqrt{c-b^2} = i \\sqrt{b^2-c}$ \u3068\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[57]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">y<\/span> <span class=\"o\">=<\/span> <span class=\"nf\">y<\/span><span class=\"p\">(<\/span><span class=\"nv\">%i<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span><span class=\"nv\">b<\/span><span class=\"o\">**<\/span>2<span class=\"o\">-<\/span><span class=\"nv\">c<\/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[57]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{86}$}y=e^ {- b\\,x }\\,\\left(\\frac{B\\,\\sinh \\left(\\sqrt{b^2-c}\\,x\\right)}{\\sqrt{b^2-c}}+A\\,\\cosh \\left(\\sqrt{b^2-c}\\,x\\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\u5f8c\u306b $c-b^2 = 0$ \u306e\u5834\u5408\u306f $\\mu \\rightarrow 0$ \u306e\u6975\u9650\u3092\u3068\u3063\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[58]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">y<\/span> <span class=\"o\">=<\/span> <span class=\"nf\">limit<\/span><span class=\"p\">(<\/span><span class=\"nf\">y<\/span><span class=\"p\">(<\/span><span class=\"nv\">mu<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">mu<\/span>, <span class=\"mi\">0<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">factor<\/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>\u001bXIs b x an integer?\r\n\u001b\\Is b x an integer?\r\nno;\r\n\r\n<\/pre>\n<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[58]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{87}$}y=\\left(B\\,x+A\\right)\\,e^ {- b\\,x }\\]<\/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=\"\u5b9a\u6570\u4fc2\u65702\u968e\u7dda\u5f62\u975e\u540c\u6b21\u65b9\u7a0b\u5f0f\">\u5b9a\u6570\u4fc2\u65702\u968e\u7dda\u5f62\u975e\u540c\u6b21\u65b9\u7a0b\u5f0f<\/h3>\n<p>\u4eba\u529b\u3067\u89e3\u304f\u969b\u306b\u306f\uff0c\u540c\u6b21\u65b9\u7a0b\u5f0f\u3068\u975e\u540c\u6b21\u65b9\u7a0b\u5f0f\u3068\u3067\u306f\uff0c\u89e3\u304f\u624b\u9593\u304c\u305a\u3044\u3076\u3093\u9055\u3063\u305f\u304c\uff0cMaxima \u3067\u306f\u3069\u3061\u3089\u3082\u540c\u3058\u3002<code>ode2()<\/code> \u3092\u4f7f\u3046\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=\"\u4f8b\u984c\">\u4f8b\u984c<\/h4>\n<p>\u975e\u540c\u6b21\u65b9\u7a0b\u5f0f $y^{&#8221;} + a^2 y = \\sin b x$ \u3092\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[59]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"cm\">\/* a, b \u306f\u6b63\uff08\u30bc\u30ed\u3067\u306a\u3044\uff09\u3068\u4eee\u5b9a\u3059\u308b\u3002*\/<\/span>\r\n<span class=\"nf\">assume<\/span><span class=\"p\">(<\/span><span class=\"nv\">a<\/span><span class=\"o\">&gt;<\/span><span class=\"mi\">0<\/span>, <span class=\"nv\">b<\/span><span class=\"o\">&gt;<\/span><span class=\"mi\">0<\/span><span class=\"p\">)<\/span>$\r\n\r\n<span class=\"nv\">eq<\/span><span class=\"o\">:<\/span> <span class=\"o\">'<\/span><span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span>, <span class=\"mi\">2<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"nv\">a<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">y<\/span> <span class=\"o\">=<\/span> <span class=\"nf\">sin<\/span><span class=\"p\">(<\/span><span class=\"nv\">b<\/span><span class=\"o\">*<\/span><span class=\"nv\">x<\/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[59]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{89}$}\\frac{d^2}{d\\,x^2}\\,y+a^2\\,y=\\sin \\left(b\\,x\\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[60]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq<\/span>, <span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/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[60]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{90}$}y=-\\frac{\\sin \\left(b\\,x\\right)}{b^2-a^2}+{\\it \\%k}_{1}\\,\\sin \\left(a\\,x\\right)+{\\it \\%k}_{2}\\,\\cos \\left(a\\,x\\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>$a = b$ \u306e\u5834\u5408\u306f\u4e0a\u8a18\u306e\u89e3\u306e\u5206\u6bcd\u304c\u30bc\u30ed\u306b\u306a\u3063\u3066\u3057\u307e\u3046\u306e\u3067\uff0c\u5225\u9014\u8a08\u7b97\u3059\u308b\u5fc5\u8981\u304c\u3042\u308b\u3002<\/p>\n<p>\u3053\u306e\u72b6\u6cc1\u306f\uff0c\u529b\u5b66\u3067\u306f\u56fa\u6709\u632f\u52d5\u6570 $a$<br \/>\n\u306b\u7b49\u3057\u3044\u632f\u52d5\u6570\u306e\u5916\u529b\u304c\u52a0\u3048\u3089\u308c\u305f\u6642\u306b\u8d77\u3053\u308b\u300c\u5171\u9cf4\uff08\u5171\u632f\uff09\u300d\u3068\u547c\u3070\u308c\u308b\u73fe\u8c61\u3067\u3042\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[61]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">eq1<\/span><span class=\"o\">:<\/span> <span class=\"nf\">ev<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq<\/span>, <span class=\"nv\">b<\/span> <span class=\"o\">=<\/span> <span class=\"nv\">a<\/span><span class=\"p\">)<\/span>; <span class=\"cm\">\/* ev(eq) \u306f b = a \u306e\u5834\u5408\u3092\u8a55\u4fa1\u3057\u3066\u304f\u308c\u308b\u3002*\/<\/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[61]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{91}$}\\frac{d^2}{d\\,x^2}\\,y+a^2\\,y=\\sin \\left(a\\,x\\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[62]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq1<\/span>, <span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/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[62]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{92}$}y={\\it \\%k}_{1}\\,\\sin \\left(a\\,x\\right)-\\frac{x\\,\\cos \\left(a\\,x\\right)}{2\\,a}+{\\it \\%k}_{2}\\,\\cos \\left(a\\,x\\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>\u7a4d\u5206\u5b9a\u6570 $\\% k_1, \\% k_2 $\u304c\u3064\u304f\u306e\u306f\u540c\u6b21\u65b9\u7a0b\u5f0f\u306e\u4e00\u822c\u89e3\u3002\u975e\u540c\u6b21\u65b9\u7a0b\u5f0f\u306e\u7279\u89e3\u306f\uff0c\u632f\u5e45\u304c $x$ \u306b\u6bd4\u4f8b\u3057\u3066\u5358\u8abf\u5897\u52a0\u3057\u3066\u3044\u304f\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=\"\u3042\u3048\u3066\u975e\u540c\u6b21\u65b9\u7a0b\u5f0f\u306e\u30ed\u30f3\u30b9\u30ad\u30a2\u30f3\u3092\u4f7f\u3063\u3066\u7279\u6b8a\u89e3\u3092\u6c42\u3081\u308b\">\u3042\u3048\u3066\u30ed\u30f3\u30b9\u30ad\u30a2\u30f3\u3092\u4f7f\u3063\u305f\u516c\u5f0f\u3067\u7279\u6b8a\u89e3\u3092\u6c42\u3081\u308b<\/h4>\n<p>\u975e\u540c\u6b21\u65b9\u7a0b\u5f0f\u3082 <code>ode2()<\/code> \u3067\u89e3\u304c\u3082\u3068\u307e\u3063\u305f\u308f\u3051\u3060\u304c\uff0c\u305d\u3053\u3092\u3042\u3048\u3066\u30ed\u30f3\u30b9\u30ad\u30a2\u30f3\u3092\u4f7f\u3063\u305f\u516c\u5f0f\u3067\u7279\u6b8a\u89e3\u3092\u6c42\u3081\u3066\u307f\u308b\u3002<\/p>\n<p>\u307e\u305a\uff0c\u540c\u6b21\u65b9\u7a0b\u5f0f\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[63]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">eq0<\/span><span class=\"o\">:<\/span> <span class=\"o\">'<\/span><span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/span>, <span class=\"mi\">2<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> <span class=\"nv\">b<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">*<\/span> <span class=\"nv\">y<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/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[63]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{93}$}\\frac{d^2}{d\\,x^2}\\,y+b^2\\,y=0\\]<\/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[64]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">ode2<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq0<\/span>, <span class=\"nv\">y<\/span>, <span class=\"nv\">x<\/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[64]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{94}$}y={\\it \\%k}_{1}\\,\\sin \\left(b\\,x\\right)+{\\it \\%k}_{2}\\,\\cos \\left(b\\,x\\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>\u4e0a\u8a18\u306e\u3088\u3046\u306b\u540c\u6b21\u65b9\u7a0b\u5f0f\u306e1\u6b21\u72ec\u7acb\u306a\u57fa\u672c\u89e3\u306f\u305d\u308c\u305e\u308c&#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[65]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">y1<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">:=<\/span> <span class=\"nf\">sin<\/span><span class=\"p\">(<\/span><span class=\"nv\">a<\/span><span class=\"o\">*<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span>;\r\n<span class=\"nf\">y2<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">:=<\/span> <span class=\"nf\">cos<\/span><span class=\"p\">(<\/span><span class=\"nv\">a<\/span><span class=\"o\">*<\/span><span class=\"nv\">x<\/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[65]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{95}$}y_{1}\\left(x\\right):=\\sin \\left(a\\,x\\right)\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[65]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{96}$}y_{2}\\left(x\\right):=\\cos \\left(a\\,x\\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>\u30ed\u30f3\u30b9\u30ad\u30a2\u30f3\u306f $W(x) \\equiv y_1 y_2&#8242; &#8211; y_1&#8242; y_2 = \\cdots$<\/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[66]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">W<\/span><span class=\"o\">:<\/span> <span class=\"nf\">y1<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nf\">y2<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">-<\/span> <span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nf\">y1<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">y2<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span>;\r\n<span class=\"nv\">W<\/span><span class=\"o\">:<\/span> <span class=\"nf\">trigsimp<\/span><span class=\"p\">(<\/span><span class=\"nv\">W<\/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[66]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{97}$}-a\\,\\sin ^2\\left(a\\,x\\right)-a\\,\\cos ^2\\left(a\\,x\\right)\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[66]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{98}$}-a\\]<\/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>\u975e\u540c\u6b21\u65b9\u7a0b\u5f0f\u306e\u7279\u6b8a\u89e3 $y_s(x)$ \u306f\uff0c\u975e\u540c\u6b21\u9805\u3092 $R(x) = \\sin b x$ \u3068\u3057\u3066<\/p>\n<p>$$y_s(x) = y_2(x) \\int^x \\frac{R(x\u2019) y_1 (x\u2019)}{W(x\u2019)} dx\u2019 \u2013 y_1(x) \\int^x \\frac{R(x\u2019) y_2(x\u2019)}{W(x\u2019)} dx\u2019$$<\/p>\n<p>$b \\neq a$ \u3092\u4eee\u5b9a\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[67]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">R<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">:=<\/span> <span class=\"nf\">sin<\/span><span class=\"p\">(<\/span><span class=\"nv\">b<\/span><span class=\"o\">*<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span>;\r\n\r\n<span class=\"nv\">ys<\/span><span class=\"o\">:<\/span> <span class=\"nf\">y2<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">integrate<\/span><span class=\"p\">(<\/span><span class=\"nf\">R<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">y1<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"nv\">W<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">-<\/span> <span class=\"nf\">y1<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">integrate<\/span><span class=\"p\">(<\/span><span class=\"nf\">R<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">y2<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"nv\">W<\/span>, <span class=\"nv\">x<\/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[67]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{99}$}R\\left(x\\right):=\\sin \\left(b\\,x\\right)\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[67]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{100}$}\\frac{\\sin \\left(a\\,x\\right)\\,\\left(-\\frac{\\cos \\left(\\left(b+a\\right)\\,x\\right)}{2\\,\\left(b+a\\right)}-\\frac{\\cos \\left(\\left(b-a\\right)\\,x\\right)}{2\\,\\left(b-a\\right)}\\right)}{a}-\\frac{\\cos \\left(a\\,x\\right)\\,\\left(\\frac{\\sin \\left(\\left(a-b\\right)\\,x\\right)}{2\\,\\left(a-b\\right)}-\\frac{\\sin \\left(\\left(b+a\\right)\\,x\\right)}{2\\,\\left(b+a\\right)}\\right)}{a}\\]<\/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>expand(ys)<\/code> \u3067\u304b\u3063\u3053\u3092\u306f\u305a\u3057\u3066\u5c55\u958b\u3057\u3066&#8230;<br \/>\n<code>trigexpand(%)<\/code> \u3067\u4e09\u89d2\u95a2\u6570\u306e\u52a0\u6cd5\u5b9a\u7406\u3092\u4f7f\u3044&#8230;<br \/>\n\u6700\u5f8c\u306b <code>trigsimp(%)<\/code> \u3067\u7c21\u5358\u5316\u3059\u308b\u3068&#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<p>\u6700\u5f8c\u306b <code>trigsimp(%)<\/code> \u3067\u7c21\u5358\u5316\u3059\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[68]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">expand<\/span><span class=\"p\">(<\/span><span class=\"nv\">ys<\/span><span class=\"p\">)<\/span>$\r\n<span class=\"nf\">trigexpand<\/span><span class=\"p\">(<\/span><span class=\"nv\">%<\/span><span class=\"p\">)<\/span>$\r\n<span class=\"nf\">trigsimp<\/span><span class=\"p\">(<\/span><span class=\"nv\">%<\/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[68]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{103}$}-\\frac{\\sin \\left(b\\,x\\right)}{b^2-a^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>\u7121\u4e8b\uff0c\u975e\u540c\u6b21\u65b9\u7a0b\u5f0f\u306e\u7279\u6b8a\u89e3<br \/>\n$y_s$ \u304c\u51fa\u307e\u3057\u305f\u3002\u305f\u3060\u3057\uff0c\u7a4d\u5206\u306e\u4e0b\u9650\u306e\u4efb\u610f\u6027\u304c\u3042\u308b\u306e\u3067\uff0c\u4ee5\u4e0b\u306e\u3088\u3046\u306b<\/p>\n<p>$$y_s(x) = y_2(x) \\int_0^x \\frac{R(t) y_1 (t)}{W(t)} dt \u2013 y_1(x) \\int_0^x \\frac{R(t) y_2(t)}{W(t)} dt$$<\/p>\n<p>\u3068\u3059\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[69]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">assume<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"o\">&gt;<\/span><span class=\"mi\">0<\/span><span class=\"p\">)<\/span>$\r\n<span class=\"nv\">ys<\/span><span class=\"o\">:<\/span> <span class=\"nf\">y2<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">integrate<\/span><span class=\"p\">(<\/span><span class=\"nf\">R<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">y1<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"nv\">W<\/span>, <span class=\"nv\">t<\/span>, <span class=\"mi\">0<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">-<\/span> <span class=\"nf\">y1<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">integrate<\/span><span class=\"p\">(<\/span><span class=\"nf\">R<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">y2<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"nv\">W<\/span>, <span class=\"nv\">t<\/span>, <span class=\"mi\">0<\/span>, <span class=\"nv\">x<\/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[69]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{105}$}\\frac{\\cos \\left(a\\,x\\right)\\,\\left(\\left(b-a\\right)\\,\\sin \\left(\\left(b+a\\right)\\,x\\right)+\\left(b+a\\right)\\,\\sin \\left(\\left(a-b\\right)\\,x\\right)\\right)}{a\\,\\left(2\\,b^2-2\\,a^2\\right)}+\\frac{\\sin \\left(a\\,x\\right)\\,\\left(\\frac{b}{b^2-a^2}-\\frac{\\left(b-a\\right)\\,\\cos \\left(\\left(b+a\\right)\\,x\\right)+\\left(b+a\\right)\\,\\cos \\left(\\left(b-a\\right)\\,x\\right)}{2\\,b^2-2\\,a^2}\\right)}{a}\\]<\/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>\u306a\u3093\u3060\u304b\u3093\u3060\u3067\u7c21\u5358\u5316\u3057\u3066\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[70]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">expand<\/span><span class=\"p\">(<\/span><span class=\"nv\">ys<\/span><span class=\"p\">)<\/span>$\r\n<span class=\"nf\">trigexpand<\/span><span class=\"p\">(<\/span><span class=\"nv\">%<\/span><span class=\"p\">)<\/span>$\r\n<span class=\"nf\">trigsimp<\/span><span class=\"p\">(<\/span><span class=\"nv\">%<\/span><span class=\"p\">)<\/span>$\r\n<span class=\"nf\">expand<\/span><span class=\"p\">(<\/span><span class=\"nv\">%<\/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[70]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{109}$}\\frac{b\\,\\sin \\left(a\\,x\\right)}{a\\,b^2-a^3}-\\frac{a\\,\\sin \\left(b\\,x\\right)}{a\\,b^2-a^3}\\]<\/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>&#8230; \u3068\u306a\u308a\uff0c$\\displaystyle y_s = \\frac{\\sin (b x)}{a^2-b^2}$ \u306e\u307b\u304b\u306b $\\sin (a x)$ \u306b\u6bd4\u4f8b\u3059\u308b\u9805\u3082\u73fe\u308c\u307e\u3059\u3002\u3057\u304b\u3057\uff0c$\\sin (a x)$ \u306b\u6bd4\u4f8b\u3059\u308b\u9805\u306f\uff0c\u540c\u6b21\u65b9\u7a0b\u5f0f\u306e\u4e00\u822c\u89e3\u306b\u7d44\u307f\u8fbc\u307e\u308c\u308b\u306e\u3067\uff0c$\\sin (b x)$ \u306b\u6bd4\u4f8b\u3059\u308b\u9805\u304c\u7d50\u5c40\u975e\u540c\u6b21\u65b9\u7a0b\u5f0f\u306e\u7279\u6b8a\u89e3\u3068\u306a\u308a\u307e\u3059\u3002<\/p>\n<p>\u3053\u306e\u3088\u3046\u306b\uff0c\u30ed\u30f3\u30b9\u30ad\u30a2\u30f3\u3092\u4f7f\u3063\u305f\u516c\u5f0f\u3067\u7279\u6b8a\u89e3\u3092\u6c42\u3081\u308b\u3068\uff0c\u4e0d\u5b9a\u7a4d\u5206\u306e\u4e0d\u5b9a\u6027\uff08\u7a4d\u5206\u306e\u4e0b\u9650\u306e\u4e0d\u5b9a\u6027\uff09\u306b\u3088\u308a\uff0c\u540c\u6b21\u65b9\u7a0b\u5f0f\u306e\u4e00\u822c\u89e3\u306b\u7d44\u307f\u8fbc\u307e\u308c\u308b\u9805\u304c\u3042\u3089\u308f\u308c\u308b\u3053\u3068\u3082\u3042\u308a\u307e\u3059\u306e\u3067\uff0c\u305d\u306e\u3078\u3093\u306f\u5404\u81ea\u51e6\u7406\u3057\u3066\u3044\u305f\u3060\u304f\u3068\u3044\u3046\u3053\u3068\u3067&#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<p>$b = a$ \u306e\u5834\u5408\u306f\uff0c\u975e\u540c\u6b21\u9805 $R(x) = \\sin (a x)$ \u3068\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[71]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">R<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">:=<\/span> <span class=\"nf\">sin<\/span><span class=\"p\">(<\/span><span class=\"nv\">a<\/span><span class=\"o\">*<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span>$\r\n\r\n<span class=\"nv\">ys<\/span><span class=\"o\">:<\/span> <span class=\"nf\">y2<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">integrate<\/span><span class=\"p\">(<\/span><span class=\"nf\">R<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">y1<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"nv\">W<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">-<\/span> <span class=\"nf\">y1<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">integrate<\/span><span class=\"p\">(<\/span><span class=\"nf\">R<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">y2<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"nv\">W<\/span>, <span class=\"nv\">x<\/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[71]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{111}$}-\\frac{\\cos \\left(a\\,x\\right)\\,\\left(a\\,x-\\frac{\\sin \\left(2\\,a\\,x\\right)}{2}\\right)}{2\\,a^2}-\\frac{\\cos ^2\\left(a\\,x\\right)\\,\\sin \\left(a\\,x\\right)}{2\\,a^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>expand(ys)<\/code> \u3067\u304b\u3063\u3053\u3092\u306f\u305a\u3057\u3066\u5c55\u958b\u3057\u3066&#8230;<br \/>\n<code>trigexpand(%)<\/code> \u3067\u4e09\u89d2\u95a2\u6570\u306e\u52a0\u6cd5\u5b9a\u7406\u3092\u4f7f\u3044&#8230;<br \/>\n\u6700\u5f8c\u306b <code>trigsimp(%)<\/code> \u3067\u7c21\u5358\u5316\u3059\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[72]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">expand<\/span><span class=\"p\">(<\/span><span class=\"nv\">ys<\/span><span class=\"p\">)<\/span>$\r\n<span class=\"nf\">trigexpand<\/span><span class=\"p\">(<\/span><span class=\"nv\">%<\/span><span class=\"p\">)<\/span>$\r\n<span class=\"nf\">trigsimp<\/span><span class=\"p\">(<\/span><span class=\"nv\">%<\/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[72]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{114}$}-\\frac{x\\,\\cos \\left(a\\,x\\right)}{2\\,a}\\]<\/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>\u7121\u4e8b\uff0c\u975e\u540c\u6b21\u65b9\u7a0b\u5f0f\u306e\u7279\u6b8a\u89e3<br \/>\n$y_s$ \u304c\u51fa\u307e\u3057\u305f\u3002\u305f\u3060\u3057\uff0c\u7a4d\u5206\u306e\u4e0b\u9650\u306e\u4efb\u610f\u6027\u304c\u3042\u308b\u306e\u3067\uff0c\u4ee5\u4e0b\u306e\u3088\u3046\u306b<\/p>\n<p>$$y_s(x) = y_2(x) \\int_0^x \\frac{R(t) y_1 (t)}{W(t)} dt \u2013 y_1(x) \\int_0^x \\frac{R(t) y_2(t)}{W(t)} dt$$<\/p>\n<p>\u3068\u3059\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[73]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">assume<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span> <span class=\"o\">&gt;<\/span> <span class=\"mi\">0<\/span><span class=\"p\">)<\/span>$\r\n\r\n<span class=\"nv\">ys<\/span><span class=\"o\">:<\/span> <span class=\"nf\">y2<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">integrate<\/span><span class=\"p\">(<\/span><span class=\"nf\">R<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">y1<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"nv\">W<\/span>, <span class=\"nv\">t<\/span>, <span class=\"mi\">0<\/span>, <span class=\"nv\">x<\/span><span class=\"p\">)<\/span> <span class=\"o\">-<\/span> <span class=\"nf\">y1<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">integrate<\/span><span class=\"p\">(<\/span><span class=\"nf\">R<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">y2<\/span><span class=\"p\">(<\/span><span class=\"nv\">t<\/span><span class=\"p\">)<\/span><span class=\"o\">\/<\/span><span class=\"nv\">W<\/span>, <span class=\"nv\">t<\/span>, <span class=\"mi\">0<\/span>, <span class=\"nv\">x<\/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[73]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{116}$}\\frac{\\cos \\left(a\\,x\\right)\\,\\left(\\sin \\left(2\\,a\\,x\\right)-2\\,a\\,x\\right)}{4\\,a^2}+\\frac{\\left(\\frac{1}{2\\,a}-\\frac{\\cos ^2\\left(a\\,x\\right)}{2\\,a}\\right)\\,\\sin \\left(a\\,x\\right)}{a}\\]<\/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>\u306a\u3093\u3060\u304b\u3093\u3060\u3067\u7c21\u5358\u5316\u3057\u3066\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[74]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">expand<\/span><span class=\"p\">(<\/span><span class=\"nv\">ys<\/span><span class=\"p\">)<\/span>$\r\n<span class=\"nf\">trigexpand<\/span><span class=\"p\">(<\/span><span class=\"nv\">%<\/span><span class=\"p\">)<\/span>$\r\n<span class=\"nf\">trigsimp<\/span><span class=\"p\">(<\/span><span class=\"nv\">%<\/span><span class=\"p\">)<\/span>$\r\n<span class=\"nf\">expand<\/span><span class=\"p\">(<\/span><span class=\"nv\">%<\/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[74]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{120}$}\\frac{\\sin \\left(a\\,x\\right)}{2\\,a^2}-\\frac{x\\,\\cos \\left(a\\,x\\right)}{2\\,a}\\]<\/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>&#8230; \u3068\u306a\u308a\uff0c\u3084\u306f\u308a\u975e\u540c\u6b21\u65b9\u7a0b\u5f0f\u306e\u7279\u6b8a\u89e3\u306e\u9805\u306e\u307b\u304b\u306b\uff0c$\\sin (a x)$ \u306b\u6bd4\u4f8b\u3059\u308b\u9805\uff0c\u3064\u307e\u308a\u540c\u6b21\u65b9\u7a0b\u5f0f\u306e\u4e00\u822c\u89e3\u306b\u7d44\u307f\u8fbc\u307e\u308c\u308b\u9805\u3082\u73fe\u308c\u307e\u3059\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"excerpt":{"rendered":"","protected":false},"author":33,"featured_media":0,"parent":6118,"menu_order":10,"comment_status":"closed","ping_status":"closed","template":"","meta":{"inline_featured_image":false,"footnotes":""},"class_list":["post-2395","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\/2395","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=2395"}],"version-history":[{"count":11,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/2395\/revisions"}],"predecessor-version":[{"id":6665,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/2395\/revisions\/6665"}],"up":[{"embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/6118"}],"wp:attachment":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/media?parent=2395"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}