Ramanujan's constant: Difference between revisions

</pre>

=={{header|Python}}==

{{libheader|mpmath}}

<lang Python>from mpmath import mp

heegner = [19,43,67,163]

mp.dps = 50

x = mp.exp(mp.pi*mp.sqrt(163))

print("calculated Ramanujan's constant: {}".format(x))

print("Heegner numbers yielding 'almost' integers:")

for i in heegner:

print(" for {}: {} ~ {} error: {}".format(str(i),mp.exp(mp.pi*mp.sqrt(i)),round(mp.exp(mp.pi*mp.sqrt(i))),(mp.pi*mp.sqrt(i)) - round(mp.pi*mp.sqrt(i))))

</lang>

{{out}}

<pre>

calculated Ramanujan's constant: 262537412640768743.99999999999925007259719818568888

Heegner numbers yielding 'almost' integers:

for 19: 885479.77768015431949753789348171962682071428650187 ~ 885480 error: 0.30611510123230903757863689092534707729405221250933

for 43: 884736743.99977746603490666193746207858537684739915 ~ 884736744 error: -0.39919930568613989412676260444831671571796782935998

for 67: 147197952743.99999866245422450682926131257862850819 ~ 147197952744 error: -0.28495586484466040200154673774982799575003729030943

for 163: 262537412640768743.99999999999925007259719818568888 ~ 262537412640768736 error: 0.10916999113251975535008362290414005390053481224586