# внутренний цикл
def func(a, b, f):
    i = a
    temp = []
    while i <= b:
        # результат каждого вычисления записываем в массив temp
        temp.append(a * f + b * i * f)
        i = i + 1

    return temp

# функция
def psi(a, b):
    f = 0.001
    j = 1
    p = [] #
    while f <= 1.0 and j <= 1000:

        # присваиваем сумму результатов внутреннего цикла в  p[j]
        p.append(sum(func(a, b, f)))

        f += 0.001
        j += 1

    return p

print(psi(2, 5))