Anastasia Kvasha_HW8

# integral approximations for average degree <k>

f = lambda x: 0.75 * (x**(1 - 2.5))
intervalBegin = 1
intervalEnd = 10000
iterations = 100000
step = (intervalEnd - intervalBegin) / iterations
integral = 0
for i in range(1, iterations):
    integral += step * f(intervalBegin + (i - 1) * step)
print("Integral is equal to: ", integral)


# numerical calculations for average degree <k>

summa = 0
for i in range(1, 10001):
    calc = 0.75 * (i**(1 - 2.5))
    summa += calc
print(summa)



# integral approximations for second moment <k2>

f = lambda x: 0.75 * (x**(2 - 2.5))
intervalBegin = 1
intervalEnd = 10000
iterations = 100000
step = (intervalEnd - intervalBegin) / iterations
integral = 0
for i in range(1, iterations):
    integral += step * f(intervalBegin + (i - 1) * step)
print("Integral is equal to: ", integral)


# numerical calculations for second moment <k2>

summa = 0
for i in range(1, 10001):
    calc = 0.75 * (i**(2 - 2.5))
    summa += calc
print(summa)