Для знаходження розв'язку лінійного неоднорідного диференціального рівняння (ЛНДР) методом операційного числення спершу знайдемо загальний розв'язок відповідного однорідного рівняння:
y′′−9y=0.
Характеристичне рівняння цього однорідного рівняння має вигляд:
r2−9=0.
Розв'язавши його, отримаємо два корені: r1=3 і r2=−3. Загальний розв'язок однорідного рівняння виглядає так:
yh(t)=c1e3t+c2e−3t,
де c1 і c2 - це довільні константи.
Тепер перейдемо до неоднорідного рівняння:
y′′−9y=10sin(t).
Для знаходження частинного розв'язку неоднорідного рівняння використовуємо метод варіації сталих. Припустимо, що частинний розв'язок має вигляд:
yp(t)=Asin(t)+Bcos(t).
Тепер знайдемо похідні від yp(t) та підставимо їх в неоднорідне рівняння:
yp′′(t)=−Asin(t)−Bcos(t),
yp′′(t)−9yp(t)=−Asin(t)−Bcos(t)−9(Asin(t)+Bcos(t))=(−8A−9B)sin(t)+(−8B+9A)cos(t).
Підставимо це у неоднорідне рівняння:
(−8A−9B)sin(t)+(−8B+9A)cos(t)=10sin(t). Рівняння має виконуватися для будь-якого t, тому коефіцієнти при синусі та косинусі мають дорівнювати:
−8A−9B=10,−8B+9A=0. Розв'язавши цю систему рівнянь, отримаємо A=1790 і B=−1780.
Таким чином, частинний розв'язок має вигляд:
yp(t)=1790sin(t)−1780cos(t).
Тепер можемо знайти загальний розв'язок неоднорідного рівняння, який задовольняє початковим умовам y(0)=1 і y′(0)=2:
y(t)=yh(t)+yp(t)=c1e3t+c2e−3t+1790sin(t)−1780cos(t).
Зараз можна використовувати початкові умови для знаходження значень констант c1 і c2:
y(0)=c1+c2+1790⋅0−1780⋅1=1,
y′(0)=3c1−3c2+
0
0