Problem 46117. Test approximations of the prime counting function

Created by ChrisR

Solve Later

{"parts":[{"partUri":"/matlab/document.xml","contentType":"application/vnd.mathworks.matlab.code.document+xml","content":"<?xml version=\"1.0\" encoding=\"UTF-8\"?><w:document xmlns:w=\"http://schemas.openxmlformats.org/wordprocessingml/2006/main\"><w:body><w:p><w:pPr><w:pStyle w:val=\"text\"/><w:jc w:val=\"left\"/></w:pPr><w:hyperlink w:docLocation=\"https://www.mathworks.com/matlabcentral/cody/problems/241\"><w:r><w:t>Cody Problem 241</w:t></w:r></w:hyperlink><w:r><w:t>, which is based on </w:t></w:r><w:hyperlink w:docLocation=\"http://projecteuler.net/problem=7\"><w:r><w:t>Project Euler Problem 7</w:t></w:r></w:hyperlink><w:r><w:t>, asks us to identify the </w:t></w:r><w:customXml w:element=\"equation\"><w:customXmlPr><w:attr w:name=\"displayStyle\" w:val=\"false\"/><w:attr w:name=\"altTextString\" w:val=\"N\"/></w:customXmlPr><w:r><w:t>N</w:t></w:r></w:customXml><w:r><w:t>th prime number. That is, the problem seeks the inverse of the prime counting function </w:t></w:r><w:customXml w:element=\"equation\"><w:customXmlPr><w:attr w:name=\"displayStyle\" w:val=\"false\"/><w:attr w:name=\"altTextString\" w:val=\"pi(n)\"/></w:customXmlPr><w:r><w:t>\\pi(n)</w:t></w:r></w:customXml><w:r><w:t>, which provides the number of primes less than or equal to </w:t></w:r><w:customXml w:element=\"equation\"><w:customXmlPr><w:attr w:name=\"displayStyle\" w:val=\"false\"/><w:attr w:name=\"altTextString\" w:val=\"n\"/></w:customXmlPr><w:r><w:t>n</w:t></w:r></w:customXml><w:r><w:t>. The </w:t></w:r><w:hyperlink w:docLocation=\"https://mathworld.wolfram.com/PrimeNumberTheorem.html\"><w:r><w:t>Prime Number Theorem</w:t></w:r></w:hyperlink><w:r><w:t> gives approximate forms of </w:t></w:r><w:customXml w:element=\"equation\"><w:customXmlPr><w:attr w:name=\"displayStyle\" w:val=\"false\"/><w:attr w:name=\"altTextString\" w:val=\"pi(n)\"/></w:customXmlPr><w:r><w:t>\\pi(n)</w:t></w:r></w:customXml><w:r><w:t> for large </w:t></w:r><w:customXml w:element=\"equation\"><w:customXmlPr><w:attr w:name=\"displayStyle\" w:val=\"false\"/><w:attr w:name=\"altTextString\" w:val=\"n\"/></w:customXmlPr><w:r><w:t> n</w:t></w:r></w:customXml><w:r><w:t>. Two such approximations are </w:t></w:r><w:customXml w:element=\"equation\"><w:customXmlPr><w:attr w:name=\"displayStyle\" w:val=\"false\"/><w:attr w:name=\"altTextString\" w:val=\"n/ln(n)\"/></w:customXmlPr><w:r><w:t>n/\\ln(n)</w:t></w:r></w:customXml><w:r><w:t> and the </w:t></w:r><w:hyperlink w:docLocation=\"https://en.wikipedia.org/wiki/Logarithmic_integral_function#Offset_logarithmic_integral\"><w:r><w:t>offset logarithmic integral</w:t></w:r></w:hyperlink><w:r><w:t> </w:t></w:r><w:customXml w:element=\"equation\"><w:customXmlPr><w:attr w:name=\"displayStyle\" w:val=\"false\"/><w:attr w:name=\"altTextString\" w:val=\"Li(n) = li(n) - li(2)\"/></w:customXmlPr><w:r><w:t>{\\rm Li(n)} ={\\rm li(n)} - {\\rm li(2)}</w:t></w:r></w:customXml><w:r><w:t>, where </w:t></w:r><w:customXml w:element=\"equation\"><w:customXmlPr><w:attr w:name=\"displayStyle\" w:val=\"false\"/><w:attr w:name=\"altTextString\" w:val=\"li(x) = integral of 1/ln(t) from t = 0 to t = infinity\"/></w:customXmlPr><w:r><w:t>{\\rm li}(x) = \\int_0^\\infty dt/\\ln(t)</w:t></w:r></w:customXml><w:r><w:t> (See </w:t></w:r><w:hyperlink w:docLocation=\"https://www.mathworks.com/matlabcentral/cody/problems/46066-evaluate-the-logarithmic-integral\"><w:r><w:t>Cody Problem 46066</w:t></w:r></w:hyperlink><w:r><w:t>).</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/><w:jc w:val=\"left\"/></w:pPr><w:r><w:t>Test these approximations by computing two ratios: </w:t></w:r><w:customXml w:element=\"equation\"><w:customXmlPr><w:attr w:name=\"displayStyle\" w:val=\"false\"/><w:attr w:name=\"altTextString\" w:val=\"r1 = [n/ln(n)]/pi(n)\"/></w:customXmlPr><w:r><w:t>r_1 =[n/\\ln(n)]/ \\pi(n)</w:t></w:r></w:customXml><w:r><w:t> and </w:t></w:r><w:customXml w:element=\"equation\"><w:customXmlPr><w:attr w:name=\"displayStyle\" w:val=\"false\"/><w:attr w:name=\"altTextString\" w:val=\"r2 = Li(n)/pi(n)\"/></w:customXmlPr><w:r><w:t>r_2 = {\\rm Li}(n)/\\pi(n)</w:t></w:r></w:customXml><w:r><w:t>. Do not round the approximations to integers. For </w:t></w:r><w:customXml w:element=\"equation\"><w:customXmlPr><w:attr w:name=\"displayStyle\" w:val=\"false\"/><w:attr w:name=\"altTextString\" w:val=\"n = 100\"/></w:customXmlPr><w:r><w:t>n = 100</w:t></w:r></w:customXml><w:r><w:t>, you will find that the first approximation is about 13% low and the second is about 16% high. However, for </w:t></w:r><w:customXml w:element=\"equation\"><w:customXmlPr><w:attr w:name=\"displayStyle\" w:val=\"false\"/><w:attr w:name=\"altTextString\" w:val=\"n = 10^8\"/></w:customXmlPr><w:r><w:t>n = 10^8</w:t></w:r></w:customXml><w:r><w:t>, the first approximation is 6% low and the second is only 0.01% high. </w:t></w:r></w:p></w:body></w:document>","relationship":null}],"relationships":[{"relationshipType":"http://schemas.mathworks.com/matlab/code/2013/relationships/document","target":"/matlab/document.xml","relationshipId":"rId1"}]}

<div style = "text-align: start; line-height: 20.4333px; min-height: 0px; white-space: normal; color: rgb(0, 0, 0); font-family: Menlo, Monaco, Consolas, monospace; font-style: normal; font-size: 14px; font-weight: 400; text-decoration: rgb(0, 0, 0); white-space: normal; "><div style="block-size: 162.817px; display: block; min-width: 0px; padding-block-start: 0px; padding-top: 0px; perspective-origin: 407px 81.4083px; transform-origin: 407px 81.4083px; vertical-align: baseline; "><div style="block-size: 90px; font-family: Helvetica, Arial, sans-serif; line-height: 21px; margin-block-end: 9px; margin-block-start: 2px; margin-bottom: 9px; margin-inline-end: 10px; margin-inline-start: 4px; margin-left: 4px; margin-right: 10px; margin-top: 2px; perspective-origin: 384px 45px; text-align: left; transform-origin: 384px 45px; white-space: pre-wrap; margin-left: 4px; margin-top: 2px; margin-bottom: 9px; margin-right: 10px; "><a target='_blank' href = "https://www.mathworks.com/matlabcentral/cody/problems/241">Cody Problem 241</a><span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 61.4583px 7.79167px; transform-origin: 61.4583px 7.79167px; unicode-bidi: normal; ">, which is based on <a target='_blank' href = "http://projecteuler.net/problem=7">Project Euler Problem 7</a><span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 73.1083px 7.79167px; transform-origin: 73.1083px 7.79167px; unicode-bidi: normal; ">, asks us to identify the N<span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 91px 7.79167px; transform-origin: 91px 7.79167px; unicode-bidi: normal; ">th prime number. That is, the problem seeks the inverse of the prime counting function <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAD4AAAAlCAYAAAAa0FocAAACX0lEQVRoge2YXZGDMBSFPw84wAAGUFAFdYADHNRCNSChHmoBDVjoPoSzvaVJubCkywNnhtmZpHv/cnJuEjhw4MCBPKiAYiNb5Ya2suIMXDe0VwG38e9usXXSQgn07DT5EyG4XLSsMttfhRIYgDqznwvQZfaxCFfCPsyNglDg0xd8zaIEHuRfbeFKoPy/49uBVHy30EkM5FHyrD5LQgtqHV+sldSE6p8X+KzH31t7JdB88DPFjZUsKwgVeyz4yoidFh/trsA9Ya+djA8JXzG/i1pbMQbREyonh7fxU4C9GUvRSsXzBFBMfOn/b4REL2Z+TrW9BX/BjUAreAqFbUXNONY6bT2cfqX+D560tj1ZyQzMF/KEr0C/KHjdQ1qxxoz145h3r3kTP/NM/Mz7KUy2POcBaYtncd5gqac9JQZ4hWNJ4iryMNqf0nTgfRFS+FPiorStsPbZxWljSeJKLBZwbeY8TFvTTX4hStsgFJxXNDp8iYtJKTa1H+ZiWCVuEK+wBMMjLtMA5lbJtqzYKqmTeA8lsjfX9t4git7NWBcZu8wY96qrEhsic1ZrZGcuoS5hK4mSZ4LT6ov6olvDaxFiUNCfRMYmFltRqb2YVjn89iy4nlp6T1eWyZw+zx7qZgIVK1L2LNP02PBp60gv3MJm++gQMW6Z4G0r1m6KnraNxbTDHp09jxktC2kOofpNIsiC54lqqWj0pFugLkGpFSrG+QafqPas7N85UOO7XPwVDTt5hLDI/R6md71dvrR2+LVhKe6sPKl9A7rnbx1gDptZ0LLd+/eJHbyvHThw4MC/4wdzuO6gMG8oWgAAAABJRU5ErkJggg==" alt="pi(n)" style="width: 31px; height: 18.5px;" width="31" height="18.5"><span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 185.15px 7.79167px; transform-origin: 185.15px 7.79167px; unicode-bidi: normal; ">, which provides the number of primes less than or equal to n<span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 17.8833px 7.79167px; transform-origin: 17.8833px 7.79167px; unicode-bidi: normal; ">. The <a target='_blank' href = "https://mathworld.wolfram.com/PrimeNumberTheorem.html">Prime Number Theorem</a><span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 87.9px 7.79167px; transform-origin: 87.9px 7.79167px; unicode-bidi: normal; "> gives approximate forms of <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAD4AAAAlCAYAAAAa0FocAAACX0lEQVRoge2YXZGDMBSFPw84wAAGUFAFdYADHNRCNSChHmoBDVjoPoSzvaVJubCkywNnhtmZpHv/cnJuEjhw4MCBPKiAYiNb5Ya2suIMXDe0VwG38e9usXXSQgn07DT5EyG4XLSsMttfhRIYgDqznwvQZfaxCFfCPsyNglDg0xd8zaIEHuRfbeFKoPy/49uBVHy30EkM5FHyrD5LQgtqHV+sldSE6p8X+KzH31t7JdB88DPFjZUsKwgVeyz4yoidFh/trsA9Ya+djA8JXzG/i1pbMQbREyonh7fxU4C9GUvRSsXzBFBMfOn/b4REL2Z+TrW9BX/BjUAreAqFbUXNONY6bT2cfqX+D560tj1ZyQzMF/KEr0C/KHjdQ1qxxoz145h3r3kTP/NM/Mz7KUy2POcBaYtncd5gqac9JQZ4hWNJ4iryMNqf0nTgfRFS+FPiorStsPbZxWljSeJKLBZwbeY8TFvTTX4hStsgFJxXNDp8iYtJKTa1H+ZiWCVuEK+wBMMjLtMA5lbJtqzYKqmTeA8lsjfX9t4git7NWBcZu8wY96qrEhsic1ZrZGcuoS5hK4mSZ4LT6ov6olvDaxFiUNCfRMYmFltRqb2YVjn89iy4nlp6T1eWyZw+zx7qZgIVK1L2LNP02PBp60gv3MJm++gQMW6Z4G0r1m6KnraNxbTDHp09jxktC2kOofpNIsiC54lqqWj0pFugLkGpFSrG+QafqPas7N85UOO7XPwVDTt5hLDI/R6md71dvrR2+LVhKe6sPKl9A7rnbx1gDptZ0LLd+/eJHbyvHThw4MC/4wdzuO6gMG8oWgAAAABJRU5ErkJggg==" alt="pi(n)" style="width: 31px; height: 18.5px;" width="31" height="18.5"><span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 29.5583px 7.79167px; transform-origin: 29.5583px 7.79167px; unicode-bidi: normal; "> for large n<span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 96.8583px 7.79167px; transform-origin: 96.8583px 7.79167px; unicode-bidi: normal; ">. Two such approximations are <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAGgAAAAlCAYAAACu2qwTAAADf0lEQVRoge1a0ZGDIBTcHuyABmzAClJBOkgHdpAWrCElXA9pwRpsIfeBO3lHEB6oaObYmcxczggP9rHsQ4GKioqKih3QHh1AxTIaABOA7uhAFNgyRgM79tPjCmA8OggFBthYt0IL4AdfoB4PAHflb1sANwDP+b5S2JocwsAm52lJMgBeiAdoYAkZ59+/UI6g+859tbDjOqXc3ZAmbyS0FEEX2P1x78nbOwmy8YQlKQUlCRoB9AX6oVG6FOhLDa4Gk3hfKYKuKLN6iAEnM0t3WBeTilIEjbCTVgot7LhOU26MyHNGMYKaud3O+d7DyqlmxXKy9nBuIUxYkRQtrEZKnZSDT9HPDnYCcuTDR1ADO/kPcf2Gt0y9xEdja3uky2839yfbNnMcvaJPwCpKkszd5pvkAJlVdDjympakAfmZskSQgZVNSYQsBDtxLSatHHMMA6zRkXNAUnvn/xPihPOe5MRlENw0Wf13sKQwCG3BucaxhCROkjDgc6BcYbHNf4QdswYNPokf5r/dpImNmQQl70MsEB/zzbKwkpOiaZiykwstQT77LrM6JDmaVUbI2oxyJmNjnxpHyGRPSl43AK6cnACAtKMdH85G0FW0SWWR80C51LTH+JNqLxnAiM+BUzY0tpdysObs6WwEDXgnqJu8wHuf1hTkWQRJZ+QLOjQhLlKPdnw4G0HSKLkTK+PRJCV/n2TvGcDk6UQGoLGkOUc7LkoQNEFHEOslqstSf9qkTDYJrityQYeiCSD3aMdFCYJ+oDMysj1f1tP9akuK5PpLWkTfgBgAN/1Qw3forWsIJQjiuGOmR5YfLqT9piuLTfxjoa1oAL4VIt1dNwcUekKYe7TjogRBtLshqZEE+FYIzRXdbYt4go5IOGOUAfhcxU1cN3PnSwRQq9eeDMuYcgiKKYLEhHA5IAt0H5E0V0+8H8qF+kw+/5P22heAvB5zcQPWnz4bvC0tM9M9F5SO84m/ktLhr+PynTS4MYf2Vmmvfe24scY2/h6J8tbNN4U8OSvnmLaufRhlRCzuh4mxdL0N3B9KKkr4Utw8JF7K+Eb0oVGOUg8HP8BD1W9Ej22MTQxb1IfZeKDsg6+tsUXtFoKBv8Ysgi2Odo5Gg31fjQqZq93xLS8lxmCwz0uGe71vp8YTB218O4Cb/la44OD3D7Y62qnYCS3WPfepqKioqPg/+AVFtJvsMIygjwAAAABJRU5ErkJggg==" alt="n/ln(n)" style="width: 52px; height: 18.5px;" width="52" height="18.5"><span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 27.225px 7.79167px; transform-origin: 27.225px 7.79167px; unicode-bidi: normal; "> and the <a target='_blank' href = "https://en.wikipedia.org/wiki/Logarithmic_integral_function#Offset_logarithmic_integral">offset logarithmic integral</a><span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 1.94167px 7.79167px; transform-origin: 1.94167px 7.79167px; unicode-bidi: normal; "> <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAPwAAAAlCAYAAABmmoSAAAAFV0lEQVR4nO2c4bHqIBCFTw/pIA3YgBVYgR3YgR3YgjVYgj3YgjXYgvcH2ckmJrCQBcG730xm3nvxEXLCgWUhAQzDMAzDMAzDMAzDMLbRAdgplqdZ1i9ieufHNF6hA3CH7g1dhsP4JIfeRwC3oWyjUo2PAM6BI6XMU0SlcghDXABcM5Rbih7jQ36t/GYP95x6YZk59T4AeKAt0+/g2usDTuclcmlMz5e8dkBYu00a71gBb3a8hn/fR5Z3ZGVIR9cHnOC5uGUuPxdXTJ/LkuF7dv4hLDe33iesG6cmerh6PjFquFTvHBp3cP54LxzkPR+bNeY39RZccI1Yw1/gesKc9HAitjr/oQapYfgSegP5OxVNuIYahpdofB/Kew5/ng+4Eg9u1phfbIs5jsMRCjlIyNgoIoUz5L1zbVDj8IX0J4TDzZJ67+HqKw2Bv43P8ICuxqfhN8eF/3tjdXnC76HNGmsZXsoVZUYbYHwQc5FbIGR4KSX1BlyDbSV/EjK8FInGT/hHZj7FCHXOmzQuafhvGPCGsg1eCw3Df0PvE1ydW0jgaRheovEe4UiTIoA3wiH7Jo1LGp5uShKOdHAi7md/P0MWZqVcsyY0DB977wdMG+4BYyZZ2jZ2aCeq0jC8ROM9wqP2HnLDb9JYcw5/g1+8O1w4skaHMRPJb/4I1/B5XZ/C+pKQLTRAjmQOf4Vfz5DegNPlilHfG5yuPMSMTeq+0UZYL5nDa2gsgRte0q6TNd5i+B7TxhLqLd/wh9fdUCZfuqCMJtWNCyMN1WOWCzk9wvsVJEdKwmzN8BdMzeiLACQa0RIt1/uFsc49pisGkjCSMtC1s9ZmtTWWQJGCVLdkjTVGeAoxfIanuY7EeNzUV3w2MooApA0w9aHwemw5UpZRfCN8B/+yHRCnNzBdgpqHp3x+Kem8qO6142uzOTT2QW1aGokma6wV0kvCozdku/hC85kz4ur8QFpvqDXCp+gaCuklIb9Ub8D//GLml8AYodWeuAu1WW2N16COI2ZQStb4Pxi+lRGH07Lh6fnETGU0OtfYPE0thqct1DGJ5RSNAZjha+W/GV5j+hSbba/B8LRLNdZ7ZngPWpnUkvyC4WNGLI0RPnZr+LcNv4NsD/0SKRoDKGd4Sux9w/Cx8yOixiy99HyM3oCu4VuJqLYaPlbj+f9NNTtQOGlHrwEulRNalpOEXTkMn7JmWWuWXnIeiAtztQ3fQkS11fBURuxUYqvZgQ0apxh+6aUUieGl2XJNw1MvnGK6lrP0QNzqhKbhXyvl1IaG4WNXgDrIzB6KGpI1jjU89U7zCkkMT0sJITQNn5oU+TYahpfqDegZvhf+rgY0DB+jcQfXOYTW7a+B3yRrHPs+PN/iOp+XSgxPo23oOqEGxnfihYxMHztoDXpXeovhpXoD/ufHd+JJXuxISiYVpoOO4aUak9kpn7R20GYfn35JGu8w3hA/1irCf8N3uHWYjrihcCX0al+H6V76+c4veh+Yzi/txONlvdDGaEPQ+wTzHAC/x0PgPEfyKuX8+fHOvMf0Qw1LO/E4d9QfztOWcF+b1dSYm11ySD6kEaWx5Jt2kuUP3xx3zWQUYq+J55sHr13Pd61WXtUkfLp3gvNzQnqvtQV6xjHXKvmxjVQkbVZbY/o2nkbOpwWNP7gj/5dlaS90a2/J5aCE3oAbdVJWQ34B09hDie/NXVB/aFmKEnofEP480y9jGgegd65zVL7FzybnJqfe9Aptaysh2pjGAfZwoZCmQJSMNLN/kkPvnN+7bxHTOMAOuvPstcSK4dDWO+ZzY/8F09gwDMMwDMMwDMMwDMMwjCl/+3GdfYQc30QAAAAASUVORK5CYII=" alt="Li(n) = li(n) - li(2)" style="width: 126px; height: 18.5px;" width="126" height="18.5"><span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 24.8917px 7.79167px; transform-origin: 24.8917px 7.79167px; unicode-bidi: normal; ">, where <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAPQAAAA2CAYAAADwD7RBAAAHYElEQVR4nO2db5X7LBCFr4c4qIEaiIJVUAd1EAe1UA2VEA+1UA210P2Q3JMpCwEmJCXdec7Jed9f6J8UuMMwDCxgGIZhGIZhGIZhGIZhGIZhGIZhGIZhGIZhGMa/pB0vwzB2zBnAE8BrvO54F3Y73nsBeADoADQbP6NhGAn8YBDqEYNIr+O/nwAO433+u8ck/CtM1IaxKWcMIuzH//cJ8DGWSW4YRNuP5TdR1mAarc09N4yNuGMYRcl1vCc5YBDmj3O/wSBkitqFo3pX5EkNw5ilwyA4jsh0m91RtR3vnTyfQdE+4B/ZXwAuhZ53Cx4Y6sH4IFf4O5uWHn9Ho2+jwTTfde9xXkwo9Bv+ivaMyQj4hBsyBDXyg+G31z7nP+Hdq0rlgKG9ixmsI4ZGD7lgDK7kUFrMpF/pc2sh1x2ma805doNphO8wBcHkPPs03i/NEcPzlx75r8gTSqw/r4FWzITTJLWoDxgsu1zuuHle14pyN/gS4or13DkGdb41oHNB3uhJA+BebKuTuMd2eaLcaHAYP1f2o5IGl95JzDNL7c9rcMLf+IaGI8JTpGSkYH0VIF23FAt0KvFQEVrswwXT0CM/At1iaJt+/K/73qNTfkB56BW8Cn9+rjcR68+lOaKsgezgD2RmMVcBXNu8Iv7QtKZbuMRXbGeBt4QjzN48EHoWj8Kfe0O+t7eloGkkS+JbisyiVAV0KN+gIRgQ2lvHj8G22Jv3wbXtOfFdkNdeDaakmRy2EjS9gVyPJFYPzPxT94ESFcDReZFlycRdq9070l3cE1wP962JE87lc8R5hm6A2ErQveI7UuqBhkzt6ZaoAE2DLcVdr907exW0DLz52oJxmFxx3qEbIFL6c4P3vtogb3MLjVjO8+XUQ48FgbZYBRzGh5lbBmCEMRdWpM9tOWK+gimAb1mbpoHaatqyBNn5Zaqp5IgpyMf+xffFjDAFowmwhfrzAYPxuYlymf/Oy13v90FxpgxgmnpYNFj5KqDB+4+PWTw3GWKOBsM8gvMudoZGlMsKCAmWrolmvVFa5CVXyYguG3FxlHMluETFufIN70tWsh1+MI0ybMfHeC8lyHWBvh5C/fWMKXjHXWncfdZiiqinGFXqIoa2HrgcqYoRhQTd4n3ZKiRoCit1PttgEoK0jrR29/G7eud+6Nk1DS/d2yVXyZgBf2+NgqZbfcf7qCE7qq/zyeWsnNHmAf0cMtZfpbDcwUIKfq7f0RikklsP7J+q5JjUCgiVL/lyKawThgrNCXTR2uVCN2jpVdLdp6BrW46j0fXlF/CZQ9MtlufMB9kntLGRJf05NZEq1/Dm1gOnHLsTtPz8O/IDAayob8Dnun4aGfDyTS84QodyGPjenN+Um+rpUpugtfWg9tY+LWg5X86Nkn+ToDWNviZyOcr3TLKjzu34yp0LpqR6zlGboLX1sFtBa+dZwHQKx96R4qll80nMnZYd1Td6y/loKiU2jtQmaE098DtU069PC1oG3nKjejnRdUltUW7ZkWrIfpMGJuT+MtIbCg7NueMhNKmeLrUJWlMPzIT8yByar9HMe+Tm/Rd0ebsaK1ZblFs+Tw0b+aXX5HN/peB9bSbLU+tIm+rpsoWgebxTDE09yOdQ9a8SgtZEm7kH9ARdNHRJJLC2KLcMPvk4oPy69xwy/8A3DYrlCcjfkypQbaqnyxaCZuQ/hqYe5PuKrUPnlAPpP1Ai3SvfPPqMecEsWnyvDP5+d/7I0yxuGH7vHUOdrZ3uKgU7dxoKy7lPncjTR4Fpi+ccd5QJCG4h6FTBaeoBmObdKkoIWh4nO/eaDu/7c4m7Hp2yaZxZSt+AL0uMHszduffE+ptSZCBHBum4o44jOBNN7ng3vjxJhdtu3YQUlyWpni5bCLpJeA2QXw/yfaqAmLS2vsSBNlIueWJ+Diytuu+znpFyl0fk+/YEBSQbkSJ3Ow2t/ppzbTn3o6GhcXFTc31ZXbIs5QCADmVO/Vjan1NTnfna2DPn1gMQPs01+iaOku51HssvM+U+eIZVCJ7GEDrk/SdS7r62lEWvAQpEupy07q5bx0679igtPakb3kXbYspHDi1ZpbYlsCzVE4j3Z2CK1ciLUxkE3j/3/DQOcyLNrQe+p4oNOnQRt9gT3aOeBIwSUNDs1DJpIzTKlBjRaoArHXvcBntD2VRdTqmqiQu1WP+8rzPS5yJ7gYJmQ0rXcM5t/IY62PNxUqUFyB1sVUHXZg14MmINa7UlcQX6nwS9NNXz0xQ5rRNTILjKNu1Q3tIwGPNtYgamwAn5L4LmIfp7h/EfbXswL6Lq9jyh7Hw65eTR2jngbxCJAr04r0sR9N6p0sVUwiBiLgzoVS1m4y/yT7tKUTNq7Qqdy3huQzPCX+NBCDmUSvU0jI8gR1Yp3tB6OtdE3fml/FM3e6bF90TqjX8K1255DBPXPn3uFg2A68Zx9465aIZRAScMo2uHeHSXf4uKyyLMId5zVNgw/jVM6O8QzswyDMMwDMMwDMMwDMMwDMMwDMMwDMPYml8KELv7nlwFLgAAAABJRU5ErkJggg==" alt="li(x) = integral of 1/ln(t) from t = 0 to t = infinity" style="width: 122px; height: 27px;" width="122" height="27"><span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 18.6667px 7.79167px; transform-origin: 18.6667px 7.79167px; unicode-bidi: normal; "> (See <a target='_blank' href = "https://www.mathworks.com/matlabcentral/cody/problems/46066-evaluate-the-logarithmic-integral">Cody Problem 46066</a><span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 4.275px 7.79167px; transform-origin: 4.275px 7.79167px; unicode-bidi: normal; ">).</div><div style="block-size: 63.8167px; font-family: Helvetica, Arial, sans-serif; line-height: 21px; margin-block-end: 9px; margin-block-start: 2px; margin-bottom: 9px; margin-inline-end: 10px; margin-inline-start: 4px; margin-left: 4px; margin-right: 10px; margin-top: 2px; perspective-origin: 384px 31.9083px; text-align: left; transform-origin: 384px 31.9083px; white-space: pre-wrap; margin-left: 4px; margin-top: 2px; margin-bottom: 9px; margin-right: 10px; "><span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 160.642px 7.79167px; transform-origin: 160.642px 7.79167px; unicode-bidi: normal; ">Test these approximations by computing two ratios: <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAQAAAAAoCAYAAAAR+iSJAAAGTElEQVR4nO2d8ZGyMBDFXw90YAM2YAVWYAd2YAe2YA2WYA/XgjXYgt8f+D4XJGQ3gRB0fzPM3A0HJJvNyyZZOMBxHMdxHMdxHOenaADsAkezQFk2hZ+5JuZqq03m9c6K2QF4Bo5d4bKcANwKPzOFDaYTKovozdVWN7S2L81SdnQEdKoauAM4LF2ICFu0HWaqEbMBcEVevXMEYPO6vnTnqdGOP0ktArBFWw6NQzQA9gAuAB6va0swtdNK/pDuvDkCcHw9uyS12vEnqUUALmgVPMYRrfPI8LeEAGwwr9g0aCOglPvnCMAdrU1LUbMdf5JaBOABm3JfUFYA/jD/PHmP1nmtI2OqADDqKhn+12zHn6QGAdijFQALJ5QTgAPa8pVwqBuAs/GaVAE4Qxd1TUXtdvxJahCAy+uwUFIA7ii3Sn6AfVROFYDSi6612/EnWVoAGqQ5cCkBWMKRHrB1lBT7MeoqFSavwY4/ydICcEA7MljRCMAWXQfg7ye0HUDDNbF8OVifmSIAKVFXDmuw4wcbtIY94a1cDdpV0xO+Y5HBKgDMRjuiW3/ZuSx2SZ2rhQRg97rf/XWOo1x/5+CpfO4Tto6yRSsuUmAatEKnFR7WTTtaWgWgQWsXrQhy2/WkOEL3XIMd/9PfaqITbV8/WxxIC4Um97BWViMA7EDsVE+8lZUJGLJjafeVc5JQQgKwxXuriW3HHYZmoLxjHYe2iYWRQ1uTnFsz1JbnYs67U/4dsQqAJeo64bP8Y8fQmsJa7PgBC3ND61h3dLefppxfSIfOOazzYUsEsEe37g3azn4WP/O8JgrISUKJTQFuI+e30LUh561aB2L9OWCwo+3QtV1s4NB2GGIVgKuiDMA70eqGd+f7e/1O+z7E76HknrXY8QOOeCd0M4z2+AyBc5kqArCWySIAbEh2qhu6SSRSMDXkJKFoBSC0vchrx7bB+Axt56K/XF/XyL1omccfux8XRrURpqWMjLpiA8UJXdswkmJ9KKKatl6LHTvQUCxIyQWTklgEgNEPt3P6nccSGVlSf4eoTQCkvxzxHrH699KuvFuE1NK5tFGXvB+FX9qKvqAR8LXYsYMc7e7QzVO5QLimxUGLAFCZ7/jMtJJhtWYakpuEUpsA9P2l3zG47qCt81wC8Ad71MWQXM7vOSXQ9Iu12HHwQdoQgllOOaPaEmgFQHbwocaUyqzBmvrbpzYBkP4y5HByVNPwhD7q1JYxJfWX18gRl/Nw7frNWuzYQa4ixwx2RndxIkUAat4FAFqDs35DBr2NnOuTkvrbp4QAsM4ax5X+0i+PnLdq2meuRcCUqIuhvuyM7KTa8q3FjoMPsqhHjgDUvgsgV9X7xm/EOc2oPkUSSgkBoG1io03MX854h7QaeD9thKTtXNbUX+YLSBvIttZOO9Zix48HWTtUDRHAHLsAstGHOsxQ9BNSaN4reW/2RQkB0K4ix/yFc2jeJzZ6ccTU+p6mM7KNLP4hp3Uss4wEea8DxttzLXb8eJA1jTBHAJZCIwCygw+paT9MPCA8wqem/vYpIQDAe9879jchf5Gr2ju8k6rGnJL771o0AmCJuhp07Suvk3kwzJTVLJKvwY4fD7LOHb5VAKjMoa0XholnxDu4NgklRmy0GBOAWEQj4SgSalN5ryF/kSPmBrov1jxgn3rGBECb+ivD/v7oDwynVGv6yRrsCKC7DWF9weJbBYDbfyEF76cIh+qvTUKJ0U8LvaDrpHIu2Z97NuiOYg+MtzPLHHK2mL/I85Z5sHXqOVYH7lBpkAPgUL2PvfPaRcU12BFA96UHK98qALRHyMn4IlBMkaf4/hxfBhk6tgivpRwj18dCyVC55fNC8OUxzcr1BfZ965gAWKMutmfIJrRhyhuINdsxm28VgKlISUKpAevbc6mkfqJrTAAYWtfwjbza7ZiNC0CYxRplIkp8Zy5VIMcEgOm0tVCzHbNxAQhT+vtzc3DEfHU4Iz03YkwAaoy6arVjNi4AYaxJKLVyxPQONrZtqiEkAEv90w8NNdoxC7nyXMN8S0t/xVwe1kWeEKW/Pzc3e0xnmw30I7S1rc6o+1+tLWXHyRn6XNJUFZubEv8cdIrUX8feVt8SdTkr54x1RUXfwhXfE3U5juM4juM4juM4juM4juM4juMY+Afhj6HXww+3bQAAAABJRU5ErkJggg==" alt="r1 = [n/ln(n)]/pi(n)" style="width: 128px; height: 20px;" width="128" height="20"><span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 15.5583px 7.79167px; transform-origin: 15.5583px 7.79167px; unicode-bidi: normal; "> and <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAMkAAAAoCAYAAABU3t4iAAAFz0lEQVR4nO1cbZHjMAx9HMKgBEqgCIKgDMqgDJZCMRTCciiFYAiF3g/nXdScY0uu40339GYys9NsHdnW06dTwOFwOBwOh8PhcOwEx58WwPH/4jBdtbCFMp8APDcYV/vsWjgA6CqO52iAI4Bv1N24r+mqidt0tcYNwLnieFxv94oKXDNXycacAVygV/gtCEJ8oa5SjwD6iuNpUJsgxAHAACdKFieEDRgRwghej+medQHPYgyNFT9Mz95yo+4IpH0XXKeW+EKQfyscEYjioZcCUrmfKF80K0keCB5rS9Qi4h1tQ60eQe6tFXhrIv4aMCF94n1reZ6u3ObSMrewYlcEQpaiQ1iblqHJgO0NCBDm9hNh5MehJkm0aKUEQPAmT5TH9mcEeVuhpQEBgodsOb+PRGuSMCyrWfLN4Y5QIChBi7BQYkDb0O6IsB81y8y/Dq1Jckd7y3VBGTHphVoRmgq7RUUrhRFvEPOAoERXzAvVISz6Fb+jMlA7J7kjnQw+YdsQ7sGyStUj7IGmesU5WpXPms8cMetG7orpzhV2UrJKKXOmg5BDk0t9o8BwXaYvSuXppgfKkmnNhhXJ+O5ltXrvkuSAoPRyXdZIwmflwpcTwvrLMUmsI4LiyoqcZh9K9muAjljs9zyV15pCcowcbvh3Dbjv18XnI/I6we8UGX0KzQ4l40UKUDNWXU6u9LJWYWp5EoYKKZIwH9FWU77EmOzdDJgblrK3kwP3UQvOJ6c4NJwPzIo7Ts+ihaaM/GzN+w3Qey45f87rNv19wOva5dabuleUl3CCdLu0Kj1snWUNankSq0w1w60cSaybQYNETz4svmsxVlRiLW7I9xCWMjH3kfLcoZ+zhcjMl56YQyspL9daUynrYTNeUSFaN5NaYs8koZGihZQWWMqtGU8byhCa/sHyYCblPYr7lnW1kEQ2b1mmlmSQUVAO2jA4KcSAdFzXYU4kS/KCn8ReSSKN1BAZk1ZbK7OFJOx4W0BrLPMNyqg1sBaSSC+79LDAnM9ZChtmktBN5hI+GZMO4jsta+vvYK8koYKthQxUeu2xCksFpyRyoL7IfWeOog1jLCSRRY2lrsk91eSopdW/v0KkqgOxOFl6IEsi9KnVLYkcSaj4mnWRRiq2eal7a7JpFLDkGIr0epwbE39L95wJfw6ySBIj/jVxL4aixF0qTsqinBH3MtxgS8nx06pbJ/w7vxxJ+CxNCJBSgl7c1xoFbdhTcgxFJskEq0uWihrL3trnrRkJejCtNyzpz7yUz1JKd0K6KWQhyadVt75hJwmtdG5dJAlihOL+sFrVIT1vWl4NOR8K+SR4SHA5r+VnHfJKy3nl9lCWmmPycO0Y5uWU/74ylkqI0uMTJMknnK4sIQlDyqV7zpEEmPsFKUgjFdtgWZoH8i8nUd6cl7UeQznhtanH78nw64aZIFrjkAp7JAlipONcZQM8V/qOFUaSkBMsTb61bnMPkDlUrjvbIV1/15CEeUnKWlLxYmSS8XgP3Wu12vNiF+h7KcvwWBLgtLhnyU1GpMkkvWyMTAz1H5ibrynjUHRerDTxJkiyT/AiPV4rctzM78i1PAYhFViSh2OszT93dD1npOT9HCEp2whdqMVuvgby5EWMgLL6NMLWQE0RetlgTcmlee4VBQZd9jtK8CmNx3dzH9kwW/ufNYW7Yd1iy/HWvBr3SKN42vcztMdQiA7pl8w4D+vJjJyRPU/jrhkZGiztc1u+2wNgPgXrSKPVG3Es0WtCiT29ynrFe29TanFB49cW/MV6GxjqbbleFsUf0f49jhQeqPMjFmto8UMcL2AFYbnhPK7iiOOC7ax3j/ierP3v3got9IJbKbE8sLs56EGW8THLfv5aZBoX1M/hrL/npamQ/QQO2OaH5HIl86og29e63/6CvQ496hoTS1OVfYe9GjMm4rVQe61VD0xVgPa68I4ZrX8NxeH4OLCM63A4HA6Hw+FwOBwOh8PhcDh+Hf4AMC7RQCoQmHQAAAAASUVORK5CYII=" alt="r2 = Li(n)/pi(n)" style="width: 100.5px; height: 20px;" width="100.5" height="20"><span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 57.95px 7.79167px; transform-origin: 57.95px 7.79167px; unicode-bidi: normal; ">. Do not round the approximations to integers. For <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAGcAAAAkCAYAAACUjSQ7AAACLklEQVRoge2ZfZGDMBDFnwccxEANoAAFOKgDHNQCGiqhHmoBDVjo/UH2sqX52BCuR4b9zWRu5rIvJVk2eQCgKIqiKIqiKEeltU1KA6ADMNjWfUl7KloADwAvAFehZgAwAxitvgPwtC210CXa02DgkkJNkpzRxo6r/zdYFv2F8CKXaE+DwZII+itNTsdiG09/b/smT3+J9rS0kCdnsnGPQH/DxrrtqD0t0uTwO3+IxD3hKmAP7amRJmdkcbFzgcdddtCKuNiB+eANlr2yZjsoTQ43D7GFGzzjlWiDXPHpaHrb18E5DMldscbA+fySlvNs4kOaHD5P6QKTKyvRJqG9cIarlgnOp285yPiilDTps4nkOqTJMZE47v7uO2iTkNO4Y5kMt3t8cjl38V6Vk7U3e8hNzpwx3jo5W7RRDN4vniqGoFKkqqqNLcmJ3f094snJ1UbhggmfE7jnDHZApMnh52vuuVGijUKLH3qA2mvv/y+qdWuAy/rsGZhPLFauPmpzazcWF/tNX1yJVnThvjKjwbY80dbm1jphHFUJP4NLtEF4Jn3lSBabLHRO9dTm1oCDvVuLvevhLq61gz9QvljfJic5VAEh19UjfOeXaD/gmfS9rOMPTAZLIntP3NGRbjlE7JsMVceW7zkp7RvcQvsOJ95fq1ujr5A0hydkBoO2+/XXzBnpxS3R/tLC7eshrrY/16kdgdRZltpaDNxL31zXWKJVFEVRFEVRFEVRlL/lB90inPH0D2r4AAAAAElFTkSuQmCC" alt="n = 100" style="width: 51.5px; height: 18px;" width="51.5" height="18"><span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 257.5px 7.79167px; transform-origin: 257.5px 7.79167px; unicode-bidi: normal; ">, you will find that the first approximation is about 13% low and the second is about 16% high. However, for <img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAGIAAAAmCAYAAAA/bE50AAACk0lEQVRoge2abZGDMBCGHw9xgIEaqIIqwME5wAEWqgEJeKiFasBC70eyw14INCnctRz7zDA3kyYk7JvsR+bAMAzDMIx3cgLO4a/xBi7ADajxQnRAD1TvXNTRqIAH01PQh8f9+YoOSoMXIjb4XPuhcZS7CYd3M2eWjSkGv0TtV2AonPPf4vCGGoCvzDFnvEu5hbFtGN+SFlNc0xDGgndTj4I5/zUiwIN8o9Shb+zbxbADaTFk3AMfqO+MohyWCr97a7x7yBXipPqm0k9xQbeZ8V9q/MDUVR2aM/lC9CwbWlxQ6l0nvPEb9R5zTYpcIfRpaBf63ZmK5UJ7r9paxpNhxR35QrSqX73Qr1P9JFZIfIjHSd+mdNGOqXoV/mP2WiHmCqHdydIObpgKNpe+ytxLJwzCQPFpkl3IYqtocaVHTOfga561GyBXCJ1d5QohBr6Q3vnSvnTCgHG3i9+T46YDT2riHLQB1jxrg12uEHrOXCE61X7F21HGOnwc6Sjgxpg7V/zc/fpDniqrODHetax51qaApUI8q4T1+2IjS7osRWCJvXDRQuUGUdC58R7vTF4RYuk7L8wLsQr94h6vqOaqftsjf+WaVqMrz3vidwlixSnYh/CbWdOmxZoO1LE/flbyL7G3rElvyFwhNrvC0IZOHTOJD69c5e4ta6oz+8nJmbv8ewkdiFPqyqRx3Mhhb1mTY3TDS98rfTaND1rd1MLitLVjf/cmJZd+4nbmMied2GxqB3lpagfoo+pCn5KC7lMoEQLGmirVVzbupomLVjd1/LXbKq4QPwTHz0u6jue1kFTE2hPIRiy9XciiZvTDcwuSfwkpqhA/BLlDi58u/PZMEKmSZVzD/tyyYRiGYRiGYRiGsSXfQQpnPGbnc4wAAAAASUVORK5CYII=" alt="n = 10^8" style="width: 49px; height: 19px;" width="49" height="19"><span style="block-size: auto; display: inline; margin-block-end: 0px; margin-block-start: 0px; margin-bottom: 0px; margin-inline-end: 0px; margin-inline-start: 0px; margin-left: 0px; margin-right: 0px; margin-top: 0px; perspective-origin: 216.658px 7.79167px; transform-origin: 216.658px 7.79167px; unicode-bidi: normal; ">, the first approximation is 6% low and the second is only 0.01% high. </div></div></div>

Solve

Solution Stats

27 Solutions
12 Solvers

Last Solution submitted on Jun 23, 2021

Last 200 Solutions

Problem Comments

3 Comments

3 Comments

ChrisR on 3 Jan 2021

Test added to discourage lookup table solutions.

Tim on 1 Feb 2021

In the problem description, r1 and r2 appear to be upside down.

ChrisR on 26 Feb 2021

Corrected. Thanks, Tim.

Problem Recent Solvers12

Problem Tags

special functions

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Problem 46117. Test approximations of the prime counting function

Solution Stats

Last 200 Solutions

Problem Comments

Problem Recent Solvers12

Suggested Problems

More from this Author131

Problem Tags

Community Treasure Hunt

Problem 46117. Test approximations of the prime counting function

Solution Stats

Last 200 Solutions

Problem Comments

Problem Recent Solvers12

Suggested Problems

More from this Author131

Problem Tags

Community Treasure Hunt

Players