{"id":4425,"date":"2022-12-08T15:24:18","date_gmt":"2022-12-08T06:24:18","guid":{"rendered":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/?page_id=4425"},"modified":"2022-12-13T22:24:19","modified_gmt":"2022-12-13T13:24:19","slug":"%e8%a3%9c%e8%b6%b3%ef%bc%9a%e8%a7%92%e5%be%84%e8%b7%9d%e9%9b%a2%e3%81%8c%e6%a5%b5%e5%a4%a7%e3%81%a8%e3%81%aa%e3%82%8b%e8%b5%a4%e6%96%b9%e5%81%8f%e7%a7%bb%e3%82%92%e6%b1%82%e3%82%81%e3%82%8b","status":"publish","type":"page","link":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/%e4%b8%80%e8%88%ac%e7%9b%b8%e5%af%be%e8%ab%96%e7%9a%84%e5%ae%87%e5%ae%99%e8%ab%96\/%e8%a7%92%e5%be%84%e8%b7%9d%e9%9b%a2\/%e8%a3%9c%e8%b6%b3%ef%bc%9a%e8%a7%92%e5%be%84%e8%b7%9d%e9%9b%a2%e3%81%8c%e6%a5%b5%e5%a4%a7%e3%81%a8%e3%81%aa%e3%82%8b%e8%b5%a4%e6%96%b9%e5%81%8f%e7%a7%bb%e3%82%92%e6%b1%82%e3%82%81%e3%82%8b\/","title":{"rendered":"\u88dc\u8db3\uff1aMaxima \u3067\u89d2\u5f84\u8ddd\u96e2\u304c\u6975\u5927\u3068\u306a\u308b\u8d64\u65b9\u504f\u79fb\u3092\u6c42\u3081\u308b"},"content":{"rendered":"<p>\u30b0\u30e9\u30d5\u304b\u3089\u308f\u304b\u308b\u3088\u3046\u306b\uff0c\u89d2\u5f84\u8ddd\u96e2 $d_A(z)$ \u306f\u6975\u5927\u5024\u3092\u6301\u3063\u3066\u3044\u308b\u3088\u3046\u3060\u3002$d_A(z)$ \u304c\u6975\u5927\u3068\u306a\u308b $z$ \u3092\u6c42\u3081\u3066\u307f\u308b\u3002<br \/>\n<!--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<p>\u30b0\u30e9\u30d5\u304b\u3089\u308f\u304b\u308b\u3088\u3046\u306b\uff0c\u89d2\u5f84\u8ddd\u96e2 $d_A(z)$ \u306f\u6975\u5927\u5024\u3092\u6301\u3063\u3066\u3044\u308b\u3088\u3046\u3060\u3002$d_A(z)$ \u304c\u6975\u5927\u3068\u306a\u308b $z$ \u3092\u6c42\u3081\u3066\u307f\u308b\u3002<\/p>\n<p><img decoding=\"async\" src=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-content\/uploads\/sites\/76\/dA1.svg\" alt=\"dA\" width=\"100%\" align=\"middle\" \/><\/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=\"$\\Omega_{\\Lambda}-=-0$-\u306e\u5834\u5408\u306e\u89d2\u5f84\u8ddd\u96e2\">$\\Omega_{\\Lambda} = 0$ \u306e\u5834\u5408\u306e\u89d2\u5f84\u8ddd\u96e2<\/h3>\n<p>$H_0$ \u3092\u7701\u7565\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[1]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">dA<\/span><span class=\"p\">(<\/span><span class=\"nv\">z<\/span>, <span class=\"nv\">Omega<\/span><span class=\"p\">)<\/span><span class=\"o\">:=<\/span> 2<span class=\"o\">\/<\/span><span class=\"p\">(<\/span><span class=\"nv\">Omega<\/span><span class=\"o\">**<\/span>2<span class=\"o\">*<\/span><span class=\"p\">(<\/span>1<span class=\"o\">+<\/span><span class=\"nv\">z<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span><span class=\"p\">)<\/span> <span class=\"o\">*<\/span> \r\n  <span class=\"p\">(<\/span>2<span class=\"o\">-<\/span><span class=\"nv\">Omega<\/span><span class=\"o\">+<\/span> <span class=\"nv\">Omega<\/span><span class=\"o\">*<\/span><span class=\"nv\">z<\/span> <span class=\"o\">-<\/span><span class=\"p\">(<\/span>2<span class=\"o\">-<\/span><span class=\"nv\">Omega<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">(<\/span>1<span class=\"o\">+<\/span><span class=\"nv\">Omega<\/span><span class=\"o\">*<\/span><span class=\"nv\">z<\/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}$}{\\it dA}\\left(z , \\Omega\\right):=\\frac{2}{\\Omega^2\\,\\left(1+z\\right)^2}\\,\\left(2-\\Omega+\\Omega\\,z+\\left(-\\left(2-\\Omega\\right)\\right)\\,\\sqrt{1+\\Omega\\,z}\\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=\"$\\Omega_{\\rm-m}-=-1$-\u306e\u5834\u5408\u306e-$d_A(z)$-\u304c\u6975\u5927\u3068\u306a\u308b\u8d64\u65b9\u504f\u79fb\">$\\Omega_{\\rm m} = 1$ \u306e\u5834\u5408\u306e $d_A(z)$ \u304c\u6975\u5927\u3068\u306a\u308b\u8d64\u65b9\u504f\u79fb<\/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[2]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">dA<\/span><span class=\"p\">(<\/span><span class=\"nv\">z<\/span>, <span class=\"mi\">1<\/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[2]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{2}$}\\frac{2\\,\\left(-\\sqrt{z+1}+z+1\\right)}{\\left(z+1\\right)^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[3]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">ddA<\/span><span class=\"o\">:<\/span> <span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nf\">dA<\/span><span class=\"p\">(<\/span><span class=\"nv\">z<\/span>, <span class=\"mi\">1<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">z<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">ratsimp<\/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}$}-\\frac{2\\,\\sqrt{z+1}-3}{\\sqrt{z+1}\\,\\left(z^2+2\\,z+1\\right)}\\]<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[4]:<\/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\">ddA<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>, <span class=\"nv\">z<\/span><span class=\"p\">)<\/span>;\r\n<\/pre>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"output_wrapper\">\n<div class=\"output\">\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[4]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{4}$}\\left[ z=\\frac{5}{4} \\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=\"$\\Omega_{\\rm-m}-=-0.1$-\u306e\u5834\u5408\u306b\u306f\u6570\u5024\u7684\u306b\u6c42\u3081\u308b\">$\\Omega_{\\rm m} = 0.1$ \u306e\u5834\u5408\u306b\u306f\u6570\u5024\u7684\u306b\u6c42\u3081\u308b<\/h4>\n<p>$\\Omega_{\\Lambda} = 0$ \u306e\u5834\u5408\u3067\u3042\u3063\u3066\u3082\uff0c$d_A(z)$ \u304c\u6975\u5927\u3068\u306a\u308b $z$ \u306f3\u6b21\u65b9\u7a0b\u5f0f\u3092\u89e3\u3044\u3066\u6c42\u3081\u308b\u5fc5\u8981\u304c\u3042\u308b\u3002\u305d\u3053\u3067\uff0c\u89e3\u306e\u516c\u5f0f\u304b\u3089\u53b3\u5bc6\u89e3\u3092\u6c42\u3081\u308b\u3088\u308a\u3082\uff0c\u4ee5\u4e0b\u306e\u3088\u3046\u306b\u6570\u5024\u89e3\u3092\u6c42\u3081\u3066\u307f\u308b\u3053\u3068\u306b\u3059\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[5]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">ddA<\/span><span class=\"o\">:<\/span> <span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nf\">dA<\/span><span class=\"p\">(<\/span><span class=\"nv\">z<\/span>, 1<span class=\"o\">\/<\/span><span class=\"mi\">10<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">z<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">ratsimp<\/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}$}-\\frac{\\sqrt{z+10}\\,\\left(2\\,10^{\\frac{3}{2}}\\,z+74\\,10^{\\frac{3}{2}}\\right)-570\\,z-7410}{\\sqrt{z+10}\\,\\left(\\sqrt{10}\\,z^3+3\\,\\sqrt{10}\\,z^2+3\\,\\sqrt{10}\\,z+\\sqrt{10}\\right)}\\]<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[6]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"cm\">\/* \u89d2\u5f84\u8ddd\u96e2\u306e\u30b0\u30e9\u30d5\u304b\u3089\u6975\u5927\u306f 2 \u304b\u3089 5 \u307e\u3067\u306e\u9593\u306b\u3042\u308a\u305d\u3046 *\/<\/span>\r\n<span class=\"nf\">find_root<\/span><span class=\"p\">(<\/span><span class=\"nv\">ddA<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>, <span class=\"nv\">z<\/span>, <span class=\"mi\">2<\/span>, <span class=\"mi\">5<\/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}$}2.835430200415257\\]<\/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=\"3\u6b21\u65b9\u7a0b\u5f0f\u3092\u89e3\u304f\u4f8b\">3\u6b21\u65b9\u7a0b\u5f0f\u3092\u89e3\u304f\u4f8b<\/h4>\n<p>$\\Omega_{\\Lambda} = 0$ \u306e\u5834\u5408\uff0c$d_A(z)$ \u304c\u6975\u5927\u3068\u306a\u308b $z$ \u30923\u6b21\u65b9\u7a0b\u5f0f\u3092\u89e3\u3044\u3066\u6c42\u3081\u3066\u307f\u308b\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[7]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">ddA<\/span><span class=\"o\">:<\/span> <span class=\"nf\">diff<\/span><span class=\"p\">(<\/span><span class=\"nf\">dA<\/span><span class=\"p\">(<\/span><span class=\"nv\">z<\/span>, <span class=\"nv\">Omega<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">z<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">ratsimp<\/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[7]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{7}$}-\\frac{\\sqrt{\\Omega\\,z+1}\\,\\left(2\\,\\Omega\\,z-6\\,\\Omega+8\\right)+\\left(3\\,\\Omega^2-6\\,\\Omega\\right)\\,z-\\Omega^2+6\\,\\Omega-8}{\\sqrt{\\Omega\\,z+1}\\,\\left(\\Omega^2\\,z^3+3\\,\\Omega^2\\,z^2+3\\,\\Omega^2\\,z+\\Omega^2\\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[8]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"cm\">\/* \u5206\u5b50\u304c\u30bc\u30ed\u306b\u306a\u308c\u3070\u3044\u3044\u306e\u3060\u304b\u3089... *\/<\/span>\r\n<span class=\"nv\">tmp<\/span><span class=\"o\">:<\/span> <span class=\"nf\">num<\/span><span class=\"p\">(<\/span><span class=\"nv\">ddA<\/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}$}-\\sqrt{\\Omega\\,z+1}\\,\\left(2\\,\\Omega\\,z-6\\,\\Omega+8\\right)-\\left(3\\,\\Omega^2-6\\,\\Omega\\right)\\,z+\\Omega^2-6\\,\\Omega+8\\]<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[9]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"cm\">\/* \u30eb\u30fc\u30c8\u306e\u9805\u3092\u4e8c\u4e57\u3057\u3066 3\u6b21\u65b9\u7a0b\u5f0f\u306b *\/<\/span>\r\n<span class=\"nv\">eq<\/span><span class=\"o\">:<\/span> <span class=\"p\">(<\/span><span class=\"nv\">tmp<\/span> <span class=\"o\">-<\/span> <span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"p\">(<\/span>3<span class=\"o\">*<\/span><span class=\"nv\">Omega<\/span><span class=\"o\">**<\/span>2<span class=\"o\">-<\/span>6<span class=\"o\">*<\/span><span class=\"nv\">Omega<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nv\">z<\/span><span class=\"o\">+<\/span> <span class=\"nv\">Omega<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">-<\/span> 6<span class=\"o\">*<\/span><span class=\"nv\">Omega<\/span> <span class=\"o\">+<\/span> <span class=\"mi\">8<\/span><span class=\"p\">))<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> \r\n         <span class=\"o\">-<\/span> <span class=\"p\">(<\/span><span class=\"o\">-<\/span><span class=\"p\">(<\/span>3<span class=\"o\">*<\/span><span class=\"nv\">Omega<\/span><span class=\"o\">**<\/span>2<span class=\"o\">-<\/span>6<span class=\"o\">*<\/span><span class=\"nv\">Omega<\/span><span class=\"p\">)<\/span><span class=\"o\">*<\/span><span class=\"nv\">z<\/span><span class=\"o\">+<\/span> <span class=\"nv\">Omega<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">-<\/span> 6<span class=\"o\">*<\/span><span class=\"nv\">Omega<\/span> <span class=\"o\">+<\/span> <span class=\"mi\">8<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>, <span class=\"nv\">ratsimp<\/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}_{9}$}4\\,\\Omega^3\\,z^3+\\left(12\\,\\Omega^3-9\\,\\Omega^4\\right)\\,z^2+\\left(6\\,\\Omega^4-12\\,\\Omega^3\\right)\\,z-\\Omega^4+12\\,\\Omega^3-16\\,\\Omega^2=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[10]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nv\">sols<\/span><span class=\"o\">:<\/span> <span class=\"nf\">solve<\/span><span class=\"p\">(<\/span><span class=\"nv\">eq<\/span>, <span class=\"nv\">z<\/span><span class=\"p\">)<\/span>$\r\n<span class=\"nv\">sols<\/span><span class=\"p\">[<\/span>3<span class=\"p\">]<\/span>, <span class=\"nv\">ratsimp<\/span>;\r\n<span class=\"nf\">ev<\/span><span class=\"p\">(<\/span><span class=\"nv\">sols<\/span><span class=\"p\">[<\/span>3<span class=\"p\">]<\/span>, <span class=\"nv\">Omega<\/span><span class=\"o\">=<\/span><span class=\"mi\">1<\/span><span class=\"p\">)<\/span>;\r\n<span class=\"nf\">ev<\/span><span class=\"p\">(<\/span><span class=\"nv\">sols<\/span><span class=\"p\">[<\/span>3<span class=\"p\">]<\/span>, <span class=\"nv\">Omega<\/span><span class=\"o\">=<\/span><span class=\"mf\">0.1<\/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[10]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{11}$}z=\\frac{\\left(27\\,\\Omega^4-144\\,\\Omega^3+\\sqrt{1-\\Omega}\\,\\left(64\\,\\Omega^2-192\\,\\Omega+128\\right)+272\\,\\Omega^2-256\\,\\Omega+128\\right)^{\\frac{2}{3}}+\\Omega^{\\frac{1}{3}}\\,\\left(3\\,\\Omega-4\\right)\\,\\left(27\\,\\Omega^4-144\\,\\Omega^3+\\sqrt{1-\\Omega}\\,\\left(64\\,\\Omega^2-192\\,\\Omega+128\\right)+272\\,\\Omega^2-256\\,\\Omega+128\\right)^{\\frac{1}{3}}+\\Omega^{\\frac{2}{3}}\\,\\left(9\\,\\Omega^2-32\\,\\Omega+32\\right)}{4\\,\\Omega^{\\frac{1}{3}}\\,\\left(27\\,\\Omega^4-144\\,\\Omega^3+\\sqrt{1-\\Omega}\\,\\left(64\\,\\Omega^2-192\\,\\Omega+128\\right)+272\\,\\Omega^2-256\\,\\Omega+128\\right)^{\\frac{1}{3}}}\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[10]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{12}$}z=\\frac{5}{4}\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[10]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{13}$}z=2.835430200415417\\]<\/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=\"$\\Omega_{\\rm-m}-+-\\Omega_{\\Lambda}-=-1$-\u306e\u5834\u5408\u306e\u89d2\u5f84\u8ddd\u96e2\">$\\Omega_{\\rm m} + \\Omega_{\\Lambda} = 1$ \u306e\u5834\u5408\u306e\u89d2\u5f84\u8ddd\u96e2<\/h3>\n<p>$H_0$ \u3092\u7701\u7565\u3057\u3066&#8230;<\/p>\n<p>\\begin{eqnarray}<br \/>\nd_A &amp;=&amp; \\frac{1}{1+z} \\int_0^z \\frac{dz}{\\sqrt{(1-\\Omega_{\\rm m}) + \\Omega_{\\rm m} (1+z)^3} }<br \/>\n\\end{eqnarray}<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[11]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">f<\/span><span class=\"p\">(<\/span><span class=\"nv\">z<\/span>, <span class=\"nv\">Omega<\/span><span class=\"p\">)<\/span><span class=\"o\">:=<\/span> 1<span class=\"o\">\/<\/span><span class=\"p\">(<\/span><span class=\"nf\">sqrt<\/span><span class=\"p\">((<\/span>1<span class=\"o\">-<\/span><span class=\"nv\">Omega<\/span><span class=\"p\">)<\/span><span class=\"o\">+<\/span><span class=\"nv\">Omega<\/span><span class=\"o\">*<\/span><span class=\"p\">(<\/span>1<span class=\"o\">+<\/span><span class=\"nv\">z<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">3<\/span><span class=\"p\">))<\/span>;\r\n<span class=\"nv\">rombergtol<\/span><span class=\"o\">:<\/span> 1<span class=\"nv\">E<\/span><span class=\"o\">-<\/span>12$\r\n\r\n<span class=\"cm\">\/* dA \u306e\u5fae\u5206 *\/<\/span>\r\n<span class=\"nf\">ddA<\/span><span class=\"p\">(<\/span><span class=\"nv\">z<\/span>, <span class=\"nv\">Omega<\/span><span class=\"p\">)<\/span><span class=\"o\">:=<\/span> <span class=\"o\">-<\/span> 1<span class=\"o\">\/<\/span><span class=\"p\">(<\/span>1<span class=\"o\">+<\/span><span class=\"nv\">z<\/span><span class=\"p\">)<\/span><span class=\"o\">**<\/span><span class=\"mi\">2<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">romberg<\/span><span class=\"p\">(<\/span><span class=\"nf\">f<\/span><span class=\"p\">(<\/span><span class=\"nv\">x<\/span>, <span class=\"nv\">Omega<\/span><span class=\"p\">)<\/span>, <span class=\"nv\">x<\/span>, <span class=\"mi\">0<\/span>, <span class=\"nv\">z<\/span><span class=\"p\">)<\/span> <span class=\"o\">+<\/span> \r\n                  1<span class=\"o\">\/<\/span><span class=\"p\">(<\/span>1<span class=\"o\">+<\/span><span class=\"nv\">z<\/span><span class=\"p\">)<\/span> <span class=\"o\">*<\/span> <span class=\"nf\">f<\/span><span class=\"p\">(<\/span><span class=\"nv\">z<\/span>, <span class=\"nv\">Omega<\/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[11]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{14}$}f\\left(z , \\Omega\\right):=\\frac{1}{\\sqrt{1-\\Omega+\\Omega\\,\\left(1+z\\right)^3}}\\]<\/div>\n<\/div>\n<div class=\"output_area\">\n<div class=\"prompt output_prompt\">Out[11]:<\/div>\n<div class=\"output_latex output_subarea output_execute_result\">\\[\\tag{${\\it \\%o}_{16}$}{\\it ddA}\\left(z , \\Omega\\right):=\\frac{-1}{\\left(1+z\\right)^2}\\,{\\it romberg}\\left(f\\left(x , \\Omega\\right) , x , 0 , z\\right)+\\frac{1}{1+z}\\,f\\left(z , \\Omega\\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=\"$\\Omega_{\\rm-m}-=-0.1$-\u306e\u5834\u5408\u306b\u6570\u5024\u7684\u306b\u6c42\u3081\u308b\">$\\Omega_{\\rm m} = 0.1$ \u306e\u5834\u5408\u306b\u6570\u5024\u7684\u306b\u6c42\u3081\u308b<\/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[12]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"cm\">\/* \u89d2\u5f84\u8ddd\u96e2\u306e\u30b0\u30e9\u30d5\u304b\u3089\u6975\u5927\u306f 1 \u304b\u3089 3 \u307e\u3067\u306e\u9593\u306b\u3042\u308a\u305d\u3046 *\/<\/span>\r\n<span class=\"nf\">find_root<\/span><span class=\"p\">(<\/span><span class=\"nf\">ddA<\/span><span class=\"p\">(<\/span><span class=\"nv\">z<\/span>, 1<span class=\"o\">\/<\/span><span class=\"mi\">10<\/span><span class=\"p\">)<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>, <span class=\"nv\">z<\/span>, <span class=\"mi\">1<\/span>, <span class=\"mi\">3<\/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}_{17}$}2.0793683111316\\]<\/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=\"Planck-2018-results-$\\Omega_{\\rm-m}-=-0.315$-\u306e\u5834\u5408\">Planck 2018 results $\\Omega_{\\rm m} = 0.315$ \u306e\u5834\u5408<\/h4>\n<p><a href=\"https:\/\/en.wikipedia.org\/wiki\/Angular_diameter_distance\">Wikipedia \u82f1\u8a9e\u7248\u306e Angular diameter distance<\/a> \u306b\u3088\u308b\u3068\uff0c$\\Lambda \\mbox{CDM}$ \u30e2\u30c7\u30eb\u3067\u306f $z$ \u304c\u7d04 $1.5$ \u3050\u3089\u3044\u306e\u3068\u3053\u308d\u3067 $d_A(z)$ \u306e\u6975\u5927\u304c\u3042\u308b\u3088\u3046\u306a\u8a18\u8ff0\u304c\u3042\u308b\u3002\u8abf\u3079\u3066\u307f\u308b\u3068\uff0c\u4ee5\u4e0b\u306e\u3088\u3046\u306b $z = 1.59$ \u3042\u305f\u308a\u304c\u6975\u5927\u5024\u3092\u4e0e\u3048\u308b\u3053\u3068\u304c\u308f\u304b\u308b\u3002<\/p>\n<\/div>\n<\/div>\n<\/div>\n<div class=\"cell border-box-sizing code_cell rendered\">\n<div class=\"input\">\n<div class=\"prompt input_prompt\">In\u00a0[13]:<\/div>\n<div class=\"inner_cell\">\n<div class=\"input_area\">\n<div class=\" highlight hl-maxima\">\n<pre><span class=\"nf\">find_root<\/span><span class=\"p\">(<\/span><span class=\"nf\">ddA<\/span><span class=\"p\">(<\/span><span class=\"nv\">z<\/span>, 315<span class=\"o\">\/<\/span><span class=\"mi\">1000<\/span><span class=\"p\">)<\/span> <span class=\"o\">=<\/span> <span class=\"mi\">0<\/span>, <span class=\"nv\">z<\/span>, <span class=\"mi\">1<\/span>, <span class=\"mi\">5<\/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}_{18}$}1.587997010828216\\]<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n","protected":false},"excerpt":{"rendered":"<p>\u30b0\u30e9\u30d5\u304b\u3089\u308f\u304b\u308b\u3088\u3046\u306b\uff0c\u89d2\u5f84\u8ddd\u96e2 $d_A(z)$ \u306f\u6975\u5927\u5024\u3092\u6301\u3063\u3066\u3044\u308b\u3088\u3046\u3060\u3002$d_A(z)$ \u304c\u6975\u5927\u3068\u306a\u308b $z$ \u3092\u6c42\u3081\u3066\u307f\u308b\u3002<\/p><p><a class=\"more-link btn\" href=\"https:\/\/home.hirosaki-u.ac.jp\/relativity\/%e4%b8%80%e8%88%ac%e7%9b%b8%e5%af%be%e8%ab%96%e7%9a%84%e5%ae%87%e5%ae%99%e8%ab%96\/%e8%a7%92%e5%be%84%e8%b7%9d%e9%9b%a2\/%e8%a3%9c%e8%b6%b3%ef%bc%9a%e8%a7%92%e5%be%84%e8%b7%9d%e9%9b%a2%e3%81%8c%e6%a5%b5%e5%a4%a7%e3%81%a8%e3%81%aa%e3%82%8b%e8%b5%a4%e6%96%b9%e5%81%8f%e7%a7%bb%e3%82%92%e6%b1%82%e3%82%81%e3%82%8b\/\">\u7d9a\u304d\u3092\u8aad\u3080<\/a><\/p>\n","protected":false},"author":33,"featured_media":0,"parent":1551,"menu_order":4,"comment_status":"closed","ping_status":"closed","template":"","meta":{"inline_featured_image":false,"footnotes":""},"class_list":["post-4425","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\/4425","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=4425"}],"version-history":[{"count":7,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/4425\/revisions"}],"predecessor-version":[{"id":4480,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/4425\/revisions\/4480"}],"up":[{"embeddable":true,"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/pages\/1551"}],"wp:attachment":[{"href":"https:\/\/home.hirosaki-u.ac.jp\/relativity\/wp-json\/wp\/v2\/media?parent=4425"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}