Задача 176. Филиппов А.Ф. Сборник задач по дифференциальным уравнениям

За время Δt (где Δt очень мало и выражено в долях года) из каждого грамма радия распадается 0,00044 Δt грамма и образуется 0,00043 Δt грамма радона. Из каждого грамма радона за время Δt распадается 70 Δt грамма. В начале опыта имелось некоторое количество x0 чистого радия. Когда количество образовавшегося и еще не распавшегося радона будет наибольшим?

Решение
Пусть x(t) - количество радия в момент времени t.
За время Δt из каждого грамма радия распадается 0,00044 Δt грамма, соответственно за время Δt количество радия уменьшится на Δt0,00044x(t). Таким образом:
x(t+Δt)x(t)=Δt0,00044x(t)

Поделив обе части уравнения на Δt и перейдя к пределу при Δt0 получим дифференциальное уравнение:
x=0,00044x
Разделим переменные:
dxx=0,00044 dt
Переменные разделены. Интегрируем обе части уравнения:
dxx=0,00044 dt
ln|x|=0,00044t+lnC
x=Ce0,00044t
В момент времени t=0 было x0 чистого радия, следовательно x(0)=x0. Получаем C=x0.
Таким образом, количество радия изменяется со временем по закону:
x=x0e0,00044t

Пусть y(t) - количество радона в момент времени t.
За время Δt из каждого грамма радия образуется 0,00043 Δt грамма радона и из каждого грамма радона за время Δt распадается 70 Δt грамма радона. Таким образом:
y(t+Δt)y(t)=0,00043Δtx(t)70Δty(t)
Поделив обе части уравнения на Δt и перейдя к пределу при Δt0, и подставив x(t) получим дифференциальное уравнение:
y=0,00043x0e0,00044t70y
y+70y=0,00043x0e0,00044t

Данное уравнение является линейным дифференциальным уравнением.
Найдем решение однородного уравнения:
y+70y=0
Разделим переменные:
dyy=70dt
Переменные разделены. Интегрируем обе части уравнения:
dyy=70dt
ln|y|=70t+lnC
y=Ce70t
Таким образом, решение однородного уравнения: y=Ce70t.

Считая постоянную C функцией от t, подставим решение однородного уравнения в исходное уравнение.
Так как y=Ce70t70Ce70t, то:
Ce70t70Ce70t+70Ce70t=0,00043x0e0,00044t
Ce70t=0,00043x0e0,00044t
C=0,00043x0e(700,00044)t
C=0,00043x0e(700,00044)tdt
C=0,00043x0700,00044e(700,00044)t+C1

Подставим полученную функцию в решение однородного уравнения:
y=Ce70t=(0,00043x0700,00044e(700,00044)t+C1)e70t
y=0,00043x0700,00044e0,00044t+C1e70t

В момент времени t=0 был только чистый радий и не было радона, следовательно y(0)=0. Получаем:
0=0,00043x0700,00044+C1
C1=0,00043x0700,00044

Таким образом, количество радона изменяется со временем по закону:
y=0,00043x0700,00044(e0,00044te70t)

Чтобы найти наибольшее количество образовавшегося и еще не распавшегося радона, найдем производную y:
y=0,00043x0700,00044(0,00044e0,00044t+70e70t)

Приравняв производную к нулю, получим:
0,00044e0,00044t+70e70t=0
e(700,00044)t=700,00044
t=ln70ln(0,00044)700,000440,171 года62,45 дня

Таким образом, количество образовавшегося и еще не распавшегося радона будет наибольшим через 62,45 дня.