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

# Пределы изменения переменной величины
# В данной задаче переменной величиной является время
t = np.arange(0, 10, 0.01)

# Создание дифференциальной функции
def invest_function(s, t):
    dsdt = - k * s * t
    return dsdt

# Определение начальных условий и параметров
s_0 = 1000
k = 0.08

# Решение дифференциального уравнения функцией odeint
solve_invest = odeint(invest_function, s_0, t)

# Построение решения в виде графика функции
plt.plot(t, solve_invest[:,0], label='Кривая падения инвестиций')

plt.xlabel('Время, годы')
plt.ylabel('Объем инвестиций')
plt.title('Инвестициив в производство')
plt.legend()

plt.show()