Разработка и анализ математических моделей с использованием MATLAB и MAPLE. Коробейников А.Г. - 38 стр.

UptoLike

Составители: 

37
> k:=1.45:
> #Решение обыкновенного дифференциального уравнения методом
интегрирующего множителя
> P:=t->k/(m[0]-beta*t):
> Q:=t->beta*c/(m[0]-beta*t)-g:
> my_sys_ode:=diff(v(t),t)+P(t)*v(t)=Q(t):
> #Расчет интегрирующего множителя
> ro:=simplify(exp(int(P(t),t))):
> #Умножаем правую часть на интегрирующий множитель
> Q(t):=simplify(Q(t)*ro):
> #Интегрируем правую часть с добавлением константы
> ro_v:=simplify(int(Q(t),t) + С1):
> #Получаем решение - делим на интегрирующий множитель
> V:=ro_v/ro:
> `С1`:=C1:
> rrr:=eval(V,t=0)=0:
> C1:=solve(rrr,C1):
> #Высота
> h:=int(V,t=0..70);
:= h 41436.56270
>
plot(V,t = 0..70);