Tak sme sa priatelia s Elí zamýšľali nad tým, kedy jej funkcia dosiahne pí 3.14 Napadlo ma nasimulovať to a keďže sa pýtala na optimálnu lambdu, napadli ma RSA funkcie. Takto je totiž ich matematika definovaná, aby nikdy nedosiahla optimálnu lambdu pre dešifrovanie vzorca:

import math
import matplotlib.pyplot as plt
import sys
n=1
k=1
lamb=3*1000000000 #warning here try negative value
ypsilons=[]
func=1/math.factorial(n)*(1/(n+lamb)-4/(2*n+1))*((2*n+1)**2/(4*(n+lamb))-n)
for k in range(k-1,10000):
    if abs(func)>0:
        print(func)
        ypsilons.append(func) #check func where is optimal lambda
    func=1/math.factorial(n)*(1/(n+lamb)-4/(2*n+1))*((2*n+1)**2/(4*(n+lamb))-n)
    n+=1 #due to factiorial definition n>=0
    lamb=10*5000000000000.1 #check negative lambda
print("lambda over sys.maxint is "+str(sys.maxsize<lamb))
plt.plot(ypsilons)
plt.ylabel('So called pi')
plt.show()
#result: there is no optimal lambda to achieve pi=3.14 or similar value
#lesna.vevericka@gmail.com

Ak chceme reálnu funkciu, tak za optimálnu lambdu môžeme prehlásiť február 85, čo je môj dátum narodenia, obrázok:

import math
import matplotlib.pyplot as plt
import sys
n=1
k=1
lamb=10 #warning here try negative value
ypsilons=[]
func=1/math.factorial(n)*(1/(n+lamb)-4/(2*n+1))*((2*n+1)**2/(4*(n+lamb))-n)
for k in range(k-1,10000):
    if abs(func)>0:
        print(func)
        ypsilons.append(func) #check func where is optimal lambda
    func+=1/math.factorial(n)*(1/(n+lamb)-4/(2*n+1))*((2*n+1)**2/(4*(n+lamb))-n)
    n+=1 #due to factiorial definition n>=0
    lamb+=285 #check negative lambda
print("lambda over sys.maxint is "+str(sys.maxsize<lamb))
#result: optimal lambda
print("optimal lambda is "+str(lamb))
plt.plot(ypsilons)
plt.ylabel('So called pi')
plt.show()
#lesna.vevericka@gmail.com

---