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

# Пределы изменения переменной величины
t = np.linspace(0, 600, 1000)  # Секунды

def loco_func(v, t):
    dvdt = (b - k * v) / m
    return dvdt

# Определение начальных условий и параметров
b = 8000
k = 100
m = 15 * 10**3
v_0 = 0

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

# Построение решения в виде графика функции
plt.plot(t, solve[:,0])

plt.xlabel('Время, t, с')
plt.ylabel('Скорость, v, м/с')
plt.grid()

plt.show()