Страница 4 из 16
Попробуем воспроизвести ее в виде программы, чтобы примерно оценить объем вычислений.
Первым шагом необходимо вычислить √12. Возникает резонный вопрос — как это сделать? Оказывается, уже в Вавилоне был известен метод вычисления квадратного корня, который сейчас так и называется «вавилонским». Суть его в вычислении √S по простой формуле:
Здесь x0 — любое приближенное значение, например для √12 можно взять 3.
Запишем формулу в виде программы:
from decimal import Decimal
print ("Квадратный корень:")
number = Decimal(12)
result = Decimal(3)
for p in range(1, 9):
result = (result + number / result)/Decimal(2)
difference = result**2 - number
print (p, result, difference)
sqrt12 = result
Результаты весьма интересны: