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

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

def sopr_func(v, t):
    dVdt = F / m - (gamma * v**2) / m
    return dVdt

# Определение начальных условий и параметров
F = 8                      # Cила, H
m = 15                     # Масса, кг
gamma = 0.06                 # Коэффициент
v_0 = 0                    # Начальная скорость

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

# Построение решения в виде графика функции
plt.plot(t, solve[:,0], color='c', label='Изменение скорости')

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

plt.show()