xxxxxxxxxx
import math
def R(n):
h = (math.pi / 2 - 0) / (2 ** n)
if n == 0:
return math.pi / 4 * (math.sin(0) + math.sin(3.14 / 2))
else:
sigma = 0
for k in range(1, 2 ** (n -1) + 1):
sigma += math.sin((2 * k - 1) * h) # a tidak ditulis karena a = 0
return R(n - 1) / 2 + h * sigma
for n in range(0, 3 + 1):
print("R(" + str(n) + ", 0) = " + str(R(n)))
print("")
print("Estimasi Error = |R(3,0) – R(2,0)| = " + str(R(3) - R(2)))
xxxxxxxxxx
R(0, 0) = 0.7853979143726935
R(1, 0) = 0.9480593244561425
R(2, 0) = 0.9871157387165868
R(3, 0) = 0.9967851407580752
Estimasi Error = |R(3,0) – R(2,0)| = 0.009669402041488473
>>>