Когда возникает ошибка «Недопустимый индекс в скалярной переменной»?

Question

import numpy as np
import scipy.integrate as integrate
import matplotlib.pyplot as plt

t = [0]
RM = 0.4
RM1 = [0]
tr = 0.4
for i in range(0,100):
    # Time module
    k = i / 10
    t.append(k)
    # RM1 module
    y = (RM1[i]/tr) + RM
    RM1.append(y)
    RM1 = integrate.simps(RM1,t)
plt.plot(t, RM1, '-b')
plt.title(r'RM1')
plt.xlabel(r't (in s)')
plt.grid()

Для приведенного выше кода, почему отображается ошибка «Недопустимый индекс для скалярной переменной»? Я пытаюсь получить результаты для сравнения. RM1 = интеграл (Rm1 / tr + RM). Но так как позже мне нужно будет написать RM как еще одно разностное выражение (здесь я считал его постоянным), чтобы запустить их одновременно, я использую эту форму. Также были бы полезны любые предложения по запуску нескольких зависимых дифференциальных уравнений.