Одномерная краевая задача теплопроводности

Рассмотрим одномерную задачу теплопроводности на отрезке от α\alpha до β\beta. Уравнение теплопроводности (1.8) в этом случае принимает вид

T(x,t)t=a22T(x,t)x2+f(x,t),\frac{\partial T(x, t)}{\partial t} = a^2 \cdot \frac{\partial^2 T(x, t)}{\partial x^2} + f(x, t),
(2.23)

Краевые условия Дирихле (1.13) в одномерном случае имеют вид

{T(x,0)=T0(x),T(α,t)=Φα(t),T(β,t)=Φβ(t)\begin{cases} T(x, 0) = T_0(x),\\ T(\alpha, t) = \Phi_\alpha(t),\\ T(\beta, t) = \Phi_\beta(t) \end{cases}
(2.24)

Для избавления от неоднородности введем T(x,t)=T^(x,t)+U(x,t)T(x,t) = \widehat{T}(x,t) + U(x,t), выберем функцию U(x,t)=a1(t)+a2(t)xU(x,t) = a_1(t) + a_2(t) \cdot x, и вычислим коэффициенты, в соответствии с (2.2)

{Φα(t)a1(t)a2(t)α=0Φβ(t)a1(t)a2(t)β=0\begin{cases} \Phi_\alpha(t) - a_1(t) - a_2(t) \cdot \alpha = 0\\ \Phi_\beta(t) - a_1(t) - a_2(t) \cdot \beta = 0 \end{cases}
(2.25)
{a1(t)=Φα(t)βΦβ(t)αβαa2(t)=Φβ(t)Φα(t)βα\begin{cases} a_1(t) = \frac{\displaystyle \Phi_\alpha(t) \cdot \beta - \Phi_\beta(t) \cdot \alpha}{\displaystyle \beta - \alpha}\\ a_2(t) = \frac{\displaystyle \Phi_\beta(t) - \Phi_\alpha(t)}{\displaystyle \beta - \alpha} \end{cases}
(2.26)

Тогда уравнение (2.23) и краевые условия (2.24) примут вид

{T^(x,t)t=a22T^(x,t)x2+f^(x,t),T^(x,0)=T0(x)a1(0)a2(0)x,T^(α,t)=0,T^(β,t)=0,\begin{cases} \frac{\displaystyle \partial \widehat{T}(x, t)}{\displaystyle \partial t} = a^2 \cdot \frac{\displaystyle \partial^2 \widehat{T}(x, t)}{\displaystyle \partial x^2} + \widehat{f}(x, t),\\ \widehat{T}(x, 0) = T_0(x) - a_1(0) - a_2(0) \cdot x,\\ \widehat{T}(\alpha, t) = 0, \widehat{T}(\beta, t) = 0, \end{cases}
(2.27)

где f^(x,t)=f(x,t)da1(t)dtda2(t)dtx\widehat{f}(x, t) = f(x, t) - \frac{\displaystyle da_1(t)}{\displaystyle dt} - \frac{\displaystyle da_2(t)}{\displaystyle dt} \cdot x.

Краевые условия Неймана (1.14) в одномерном случае имеют вид

{T(x,0)=T0(x),T(x,t)xx=α=Φα(t),T(x,t)xx=β=Φβ(t)\begin{cases} T(x, 0) = T_0(x),\\ \frac{\displaystyle \partial T(x, t)}{\displaystyle \partial x} \bigg|_{x=\alpha} = \Phi_\alpha(t),\\ \frac{\displaystyle \partial T(x, t)}{\displaystyle \partial x} \bigg|_{x=\beta} = \Phi_\beta(t) \end{cases}
(2.28)

Для избавления от неоднородности введем T(x,t)=T^(x,t)+U(x,t)T(x,t) = \widehat{T}(x,t) + U(x,t), выберем функцию U(x,t)=a1(t)x+a2(t)x2U(x,t) = a_1(t) \cdot x + a_2(t) \cdot x^2, и вычислим коэффициенты, в соответствии с (2.3)

{Φα(t)a1(t)2a2(t)α=0Φβ(t)a1(t)2a2(t)β=0\begin{cases} \Phi_\alpha(t) - a_1(t) - 2 \cdot a_2(t) \cdot \alpha = 0\\ \Phi_\beta(t) - a_1(t) - 2 \cdot a_2(t) \cdot \beta = 0 \end{cases}
(2.29)
{a1(t)=Φα(t)βΦβ(t)αβαa2(t)=Φβ(t)Φα(t)2(βα)\begin{cases} a_1(t) = \frac{\displaystyle \Phi_\alpha(t) \cdot \beta - \Phi_\beta(t) \cdot \alpha}{\displaystyle \beta - \alpha}\\ a_2(t) = \frac{\displaystyle \Phi_\beta(t) - \Phi_\alpha(t)}{\displaystyle 2 \cdot (\beta - \alpha)} \end{cases}
(2.30)

Тогда уравнение (2.23) и краевые условия (2.28) примут вид

{T^(x,t)t=a22T^(x,t)x2+f^(x,t),T^(x,0)=T0(x)a1(0)xa2(0)x2,T^(x,t)xx=α=0,T^(x,t)xx=β=0,\begin{cases} \frac{\displaystyle \partial \widehat{T}(x, t)}{\displaystyle \partial t} = a^2 \cdot \frac{\displaystyle \partial^2 \widehat{T}(x, t)}{\displaystyle \partial x^2} + \widehat{f}(x, t),\\ \widehat{T}(x, 0) = T_0(x) - a_1(0) \cdot x - a_2(0) \cdot x^2,\\ \frac{\displaystyle \partial \widehat{T}(x, t)}{\displaystyle \partial x} \bigg|_{x=\alpha} = 0, \frac{\displaystyle \partial \widehat{T}(x, t)}{\displaystyle \partial x} \bigg|_{x=\beta} = 0, \end{cases}
(2.31)

где f^(x,t)=f(x,t)da1(t)dtxda2(t)dtx2+2a2a2(t)\widehat{f}(x, t) = f(x, t) - \frac{\displaystyle da_1(t)}{\displaystyle dt} \cdot x - \frac{\displaystyle da_2(t)}{\displaystyle dt} \cdot x^2 + 2 \cdot a^2 \cdot a_2(t).

Уравнение (2.8) примет вид

d2Ψ(x)dx2+γ2Ψ(x)=0.\frac{\displaystyle d^2 \Psi(x)}{\displaystyle dx^2} + \gamma^2 \cdot \Psi(x) = 0.
(2.32)

Общее решение этого уравнения имеет вид

Ψ(x)=c1cos[γ(xα)]+c2sin[γ(xα)],\Psi(x) = c_1 \cdot \cos[\gamma \cdot (x - \alpha)] + c_2 \cdot \sin[\gamma \cdot (x - \alpha)],
(2.33)

где c1c_1 и c2c_2 — коэффициенты, определяемые из граничных условий, а вес, обеспечивающий ортогональность собственных функций, ρ(x)=1\rho(x) = 1.

Подставим общее решение в граничные условия (2.27) и получим

{c1cos(γ0)+c2sin(γ0)=0,c1cos[γ(βα)]+c2sin[γ(βα)]=0.\begin{cases} c_1 \cdot \cos(\gamma \cdot 0) + c_2 \cdot \sin(\gamma \cdot 0) = 0,\\ c_1 \cdot \cos[\gamma \cdot (\beta - \alpha)] + c_2 \cdot \sin[\gamma \cdot (\beta - \alpha)] = 0. \end{cases}
(2.34)

Понятно, что из первого уравнения c1=0c_1 = 0 и решение имеет смысл только тогда, когда c20c_2 \neq 0, что возможно только тогда, когда sin[γ(βα)]=0\sin[\gamma \cdot (\beta - \alpha)] = 0

γ(βα)=πn,n(1..).\gamma \cdot (\beta - \alpha) = \pi \cdot n, n \in (1..\infty).
(2.35)

Нулевое и отрицательные значения nn не годятся: при n=0n = 0 имеем γ=0\gamma = 0, и Ψ(x)=c2sin[0(xα)]0\Psi(x) = c_2 \cdot \sin[0 \cdot (x - \alpha)] \equiv 0 — тождественный нуль, не являющийся собственной функцией; отрицательным же nn отвечает то же собственное значение γ2\gamma^2 (ведь γn=γn\gamma_{-n} = -\gamma_n), а sin[γn(xα)]=sin[γn(xα)]\sin[-\gamma_n \cdot (x - \alpha)] = -\sin[\gamma_n \cdot (x - \alpha)] линейно зависимо с решением для положительного nn и нового вклада не даёт. Поэтому n(1..)n \in (1..\infty).

Таким образом, собственные значения и собственные функции имеют вид

{γn=πnβα,n(1..),Ψn(x)=c2nsin[πnβα(xα)],\begin{cases} \gamma_n = \frac{\displaystyle \pi \cdot n}{\displaystyle \beta - \alpha}, n \in (1..\infty),\\ \Psi_n(x) = c_{2n} \cdot \sin \left[ \frac{\displaystyle \pi \cdot n}{\displaystyle \beta - \alpha} \cdot (x - \alpha) \right], \end{cases}
(2.36)

где c2nc_{2n} — произвольный постоянный множитель: собственная функция определена с точностью до него. Положим c2n=1c_{2n} = 1 — этот множитель всё равно вбирается в коэффициенты разложения T0nT_{0n} и fnf_n, как и в общем решении.

Для разложения функций в ряд Фурье по Ψn(x)\Psi_n(x) необходимо вычислить норму Ψn(x)2\|\Psi_n(x)\|^2, вес для декартовых координат ρ(x)=1\rho(x) = 1.

Ψn(x)2=αβΨn(x)2dx=αβsin[πnβα(xα)]2dx=βα2.\|\Psi_n(x)\|^2 = \int_\alpha^\beta \Psi_n(x)^2 \,dx = \int_\alpha^\beta \sin \left[ \frac{\displaystyle \pi \cdot n}{\displaystyle \beta - \alpha} \cdot (x - \alpha) \right]^2 \,dx = \frac{\displaystyle \beta - \alpha}{\displaystyle 2}.
(2.37)

Теперь подставим решение в граничные условия (2.31) и получим

{c1sin(γn0)+c2cos(γn0)=0,c1sin[γn(βα)]+c2cos[γn(βα)]=0.\begin{cases} - c_1 \cdot \sin(\gamma_n \cdot 0) + c_2 \cdot \cos(\gamma_n \cdot 0) = 0,\\ - c_1 \cdot \sin[\gamma_n \cdot (\beta - \alpha)] + c_2 \cdot \cos[\gamma_n \cdot (\beta - \alpha)] = 0. \end{cases}
(2.38)

Понятно, что из первого уравнения c2=0c_2 = 0 и решение имеет смысл только тогда, когда sin[γn(βα)]=0\sin[\gamma_n \cdot (\beta - \alpha)] = 0

γn(βα)=πn,n(0..).\gamma_n \cdot (\beta - \alpha) = \pi \cdot n, n \in (0..\infty).
(2.39)

Здесь, в отличие от задачи Дирихле, мода n=0n = 0 сохраняется: при γ0=0\gamma_0 = 0 собственная функция Ψ0(x)=c10cos[0(xα)]=c10\Psi_0(x) = c_{10} \cdot \cos[0 \cdot (x - \alpha)] = c_{10} — ненулевая постоянная, то есть полноценная собственная функция (постоянная мода). Поэтому индексация начинается с нуля.

Таким образом, собственные значения и собственные функции имеют вид

{γn=πnβα,n(0..),Ψn(x)=c1ncos[πnβα(xα)],\begin{cases} \gamma_n = \frac{\displaystyle \pi \cdot n}{\displaystyle \beta - \alpha}, n \in (0..\infty),\\ \Psi_n(x) = c_{1n} \cdot \cos \left[ \frac{\displaystyle \pi \cdot n}{\displaystyle \beta - \alpha} \cdot (x - \alpha) \right], \end{cases}
(2.40)

где c1nc_{1n} — произвольный постоянный множитель: собственная функция определена с точностью до него. Положим c1n=1c_{1n} = 1 — этот множитель всё равно вбирается в коэффициенты разложения T0nT_{0n} и fnf_n, как и в общем решении.

Норму собственных функций вычислим так же, как для случая с синусом, понизив степень косинуса по формуле половинного угла. Для n1n \ge 1

Ψn(x)2=αβΨn(x)2dx=αβcos[γn(xα)]2dx=αβ1+cos[2γn(xα)]2dx=βα2,n(1..).\|\Psi_n(x)\|^2 = \int_\alpha^\beta \Psi_n(x)^2 \,dx = \int_\alpha^\beta \cos \left[ \gamma_n \cdot (x - \alpha) \right]^2 \,dx = \int_\alpha^\beta \frac{\displaystyle 1 + \cos \left[ 2 \cdot \gamma_n \cdot (x - \alpha) \right]}{\displaystyle 2} \,dx = \frac{\displaystyle \beta - \alpha}{\displaystyle 2}, \quad n \in (1..\infty).
(2.41)

поскольку интеграл от косинуса по целому числу полупериодов равен нулю: sin[2γn(βα)]=sin(2πn)=0\sin \left[ 2 \cdot \gamma_n \cdot (\beta - \alpha) \right] = \sin(2 \cdot \pi \cdot n) = 0. Особый случай n=0n = 0 (при γ0=0\gamma_0 = 0 собственная функция Ψ0(x)=1\Psi_0(x) = 1) вычислим отдельно

Ψ0(x)2=αβΨ0(x)2dx=αβ1dx=βα.\|\Psi_0(x)\|^2 = \int_\alpha^\beta \Psi_0(x)^2 \,dx = \int_\alpha^\beta 1 \,dx = \beta - \alpha.
(2.42)

Финальное решение одномерной краевой задачи Дирихле с неоднородными граничными условиями примет вид

a1(t)=Φα(t)βΦβ(t)αβα,a2(t)=Φβ(t)Φα(t)βα,U(x,t)=a1(t)+a2(t)x,T^(x,0)=T0(x)U(x,0),f^(x,t)=f(x,t)da1(t)dtda2(t)dtx,γn=πnβα,Ψn(x)=sin[γn(xα)],Ψn(x)2=βα2,n(1..),T0n=αβT^(ω,0)Ψn(ω)dω,fn(τ)=αβf^(ω,τ)Ψn(ω)dω,T^(x,t)=2βαn=1Ψn(x)ea2γn2t[T0n+0tea2γn2τfn(τ)dτ],T(x,t)=U(x,t)+T^(x,t).\begin{aligned} & a_1(t) = \frac{\displaystyle \Phi_\alpha(t) \cdot \beta - \Phi_\beta(t) \cdot \alpha}{\displaystyle \beta - \alpha}, \quad a_2(t) = \frac{\displaystyle \Phi_\beta(t) - \Phi_\alpha(t)}{\displaystyle \beta - \alpha},\\ & U(x, t) = a_1(t) + a_2(t) \cdot x,\\ & \widehat{T}(x, 0) = T_0(x) - U(x, 0),\\ & \widehat{f}(x, t) = f(x, t) - \frac{\displaystyle da_1(t)}{\displaystyle dt} - \frac{\displaystyle da_2(t)}{\displaystyle dt} \cdot x,\\ & \gamma_n = \frac{\displaystyle \pi \cdot n}{\displaystyle \beta - \alpha}, \quad \Psi_n(x) = \sin \left[ \gamma_n \cdot (x - \alpha) \right], \quad \|\Psi_n(x)\|^2 = \frac{\displaystyle \beta - \alpha}{\displaystyle 2}, \quad n \in (1..\infty),\\ & T_{0n} = \int_\alpha^\beta \widehat{T}(\omega, 0) \cdot \Psi_n(\omega) \,d\omega, \quad f_n(\tau) = \int_\alpha^\beta \widehat{f}(\omega, \tau) \cdot \Psi_n(\omega) \,d\omega,\\ & \widehat{T}(x, t) = \frac{\displaystyle 2}{\displaystyle \beta - \alpha} \cdot \sum_{n=1}^{\infty} \Psi_n(x) \cdot e^{- a^2 \cdot \gamma_n^2 \cdot t} \cdot \left[ T_{0n} + \int_0^t e^{a^2 \cdot \gamma_n^2 \cdot \tau} \cdot f_n(\tau) \,d\tau \right],\\ & T(x, t) = U(x, t) + \widehat{T}(x, t). \end{aligned}
(2.43)

В задаче Неймана, в отличие от Дирихле, спектр содержит нулевую моду (γ0=0\gamma_0 = 0, Ψ0(x)=1\Psi_0(x) = 1) — случай γn=0\gamma_n = 0, который в общем решении (2.18) был оставлен для отдельного рассмотрения; именно она даёт первое (незатухающее) слагаемое. Финальное решение одномерной краевой задачи Неймана с неоднородными граничными условиями примет вид

a1(t)=Φα(t)βΦβ(t)αβα,a2(t)=Φβ(t)Φα(t)2(βα),U(x,t)=a1(t)x+a2(t)x2,T^(x,0)=T0(x)U(x,0),f^(x,t)=f(x,t)da1(t)dtxda2(t)dtx2+2a2a2(t),γn=πnβα,Ψn(x)=cos[γn(xα)],n(0..),Ψn(x)2=βα2(n1),Ψ0(x)2=βα,T0n=αβT^(ω,0)Ψn(ω)dω,fn(τ)=αβf^(ω,τ)Ψn(ω)dω,T^(x,t)=1βα[T00+0tf0(τ)dτ]+2βαn=1Ψn(x)ea2γn2t[T0n+0tea2γn2τfn(τ)dτ],T(x,t)=U(x,t)+T^(x,t).\begin{aligned} & a_1(t) = \frac{\displaystyle \Phi_\alpha(t) \cdot \beta - \Phi_\beta(t) \cdot \alpha}{\displaystyle \beta - \alpha}, \quad a_2(t) = \frac{\displaystyle \Phi_\beta(t) - \Phi_\alpha(t)}{\displaystyle 2 \cdot (\beta - \alpha)},\\ & U(x, t) = a_1(t) \cdot x + a_2(t) \cdot x^2,\\ & \widehat{T}(x, 0) = T_0(x) - U(x, 0),\\ & \widehat{f}(x, t) = f(x, t) - \frac{\displaystyle da_1(t)}{\displaystyle dt} \cdot x - \frac{\displaystyle da_2(t)}{\displaystyle dt} \cdot x^2 + 2 \cdot a^2 \cdot a_2(t),\\ & \gamma_n = \frac{\displaystyle \pi \cdot n}{\displaystyle \beta - \alpha}, \quad \Psi_n(x) = \cos \left[ \gamma_n \cdot (x - \alpha) \right], \quad n \in (0..\infty),\\ & \|\Psi_n(x)\|^2 = \frac{\displaystyle \beta - \alpha}{\displaystyle 2} \quad (n \ge 1), \qquad \|\Psi_0(x)\|^2 = \beta - \alpha,\\ & T_{0n} = \int_\alpha^\beta \widehat{T}(\omega, 0) \cdot \Psi_n(\omega) \,d\omega, \quad f_n(\tau) = \int_\alpha^\beta \widehat{f}(\omega, \tau) \cdot \Psi_n(\omega) \,d\omega,\\ & \widehat{T}(x, t) = \frac{\displaystyle 1}{\displaystyle \beta - \alpha} \cdot \left[ T_{00} + \int_0^t f_0(\tau) \,d\tau \right] + \frac{\displaystyle 2}{\displaystyle \beta - \alpha} \cdot \sum_{n=1}^{\infty} \Psi_n(x) \cdot e^{- a^2 \cdot \gamma_n^2 \cdot t} \cdot \left[ T_{0n} + \int_0^t e^{a^2 \cdot \gamma_n^2 \cdot \tau} \cdot f_n(\tau) \,d\tau \right],\\ & T(x, t) = U(x, t) + \widehat{T}(x, t). \end{aligned}
(2.44)

Чтобы получить стационарные решения, устремим в найденных решениях время к бесконечности (tt \to \infty), как это сделано для общего случая в (2.18). Слагаемое с начальным условием занулится, поскольку ea2γn2t0e^{- a^2 \cdot \gamma_n^2 \cdot t} \to 0; источники и граничные условия перестают зависеть от времени, а интеграл по времени для γn0\gamma_n \neq 0 даёт множитель 1a2γn2\frac{1}{a^2 \cdot \gamma_n^2}:

0tea2γn2(tτ)dτ=1a2γn2ea2γn2(tτ)0t=1a2γn2(1ea2γn2t)1a2γn2приt,γn0.\int_0^t e^{- a^2 \cdot \gamma_n^2 \cdot (t - \tau)} \,d\tau = \frac{\displaystyle 1}{\displaystyle a^2 \cdot \gamma_n^2} \cdot e^{- a^2 \cdot \gamma_n^2 \cdot (t - \tau)} \bigg|_0^t = \frac{\displaystyle 1}{\displaystyle a^2 \cdot \gamma_n^2} \cdot (1 - e^{- a^2 \cdot \gamma_n^2 \cdot t}) \rightarrow \frac{\displaystyle 1}{\displaystyle a^2 \cdot \gamma_n^2} \quad \text{при} \quad t \rightarrow \infty, \quad \gamma_n \neq 0.

Стационарное решение одномерной краевой задачи Дирихле примет вид

a1=ΦαβΦβαβα,a2=ΦβΦαβα,U(x)=a1+a2x,γn=πnβα,Ψn(x)=sin[γn(xα)],n(1..),fn=αβf(ω)Ψn(ω)dω,T(x)=U(x)+2βαn=1Ψn(x)a2γn2fn.\begin{aligned} & a_1 = \frac{\displaystyle \Phi_\alpha \cdot \beta - \Phi_\beta \cdot \alpha}{\displaystyle \beta - \alpha}, \quad a_2 = \frac{\displaystyle \Phi_\beta - \Phi_\alpha}{\displaystyle \beta - \alpha}, \quad U(x) = a_1 + a_2 \cdot x,\\ & \gamma_n = \frac{\displaystyle \pi \cdot n}{\displaystyle \beta - \alpha}, \quad \Psi_n(x) = \sin \left[ \gamma_n \cdot (x - \alpha) \right], \quad n \in (1..\infty),\\ & f_n = \int_\alpha^\beta f(\omega) \cdot \Psi_n(\omega) \,d\omega,\\ & T(x) = U(x) + \frac{\displaystyle 2}{\displaystyle \beta - \alpha} \cdot \sum_{n=1}^{\infty} \frac{\displaystyle \Psi_n(x)}{\displaystyle a^2 \cdot \gamma_n^2} \cdot f_n. \end{aligned}
(2.45)

Стационарное решение одномерной краевой задачи Неймана не такое простое. В отличие от задачи Дирихле, чисто неймановская стационарная задача разрешима не всегда. Приведённая функция T^\widehat{T} удовлетворяет однородным условиям Неймана T^(α)=T^(β)=0\widehat{T}'(\alpha) = \widehat{T}'(\beta) = 0 и стационарному уравнению a2T^(x)+f^(x)=0a^2 \cdot \widehat{T}''(x) + \widehat{f}(x) = 0. Проинтегрируем его по отрезку [α,β][\alpha, \beta]:

0=αβ(a2T^(x)+f^(x))dx=a2[T^(β)T^(α)]+αβf^(x)dx=αβf^(x)dx.0 = \int_\alpha^\beta \left( a^2 \cdot \widehat{T}''(x) + \widehat{f}(x) \right) dx = a^2 \cdot \left[ \widehat{T}'(\beta) - \widehat{T}'(\alpha) \right] + \int_\alpha^\beta \widehat{f}(x) \,dx = \int_\alpha^\beta \widehat{f}(x) \,dx.

Граничный член обратился в нуль в силу однородных условий Неймана, и осталось условие разрешимости (2.22): αβf^(ω)dω=0\int_\alpha^\beta \widehat{f}(\omega) \,d\omega = 0 — суммарный приведённый источник должен обращаться в нуль (иначе тепло накапливается, среднее неограниченно растёт и стационара нет). При выполнении этого условия решение определяется лишь с точностью до произвольной аддитивной постоянной CC: ей отвечает нулевая мода (γ0=0\gamma_0 = 0, Ψ0(x)=1\Psi_0(x) = 1), амплитуда которой стационарным уравнением не фиксируется. С этой оговоркой решение примет вид

a1=ΦαβΦβαβα,a2=ΦβΦα2(βα),U(x)=a1x+a2x2,f^(x)=f(x)+2a2a2,γn=πnβα,Ψn(x)=cos[γn(xα)],n(1..),fn=αβf^(ω)Ψn(ω)dω,T(x)=U(x)+2βαn=1Ψn(x)a2γn2fn+C,C=const.\begin{aligned} & a_1 = \frac{\displaystyle \Phi_\alpha \cdot \beta - \Phi_\beta \cdot \alpha}{\displaystyle \beta - \alpha}, \quad a_2 = \frac{\displaystyle \Phi_\beta - \Phi_\alpha}{\displaystyle 2 \cdot (\beta - \alpha)}, \quad U(x) = a_1 \cdot x + a_2 \cdot x^2,\\ & \widehat{f}(x) = f(x) + 2 \cdot a^2 \cdot a_2,\\ & \gamma_n = \frac{\displaystyle \pi \cdot n}{\displaystyle \beta - \alpha}, \quad \Psi_n(x) = \cos \left[ \gamma_n \cdot (x - \alpha) \right], \quad n \in (1..\infty),\\ & f_n = \int_\alpha^\beta \widehat{f}(\omega) \cdot \Psi_n(\omega) \,d\omega,\\ & T(x) = U(x) + \frac{\displaystyle 2}{\displaystyle \beta - \alpha} \cdot \sum_{n=1}^{\infty} \frac{\displaystyle \Psi_n(x)}{\displaystyle a^2 \cdot \gamma_n^2} \cdot f_n + C, \quad C = \text{const}. \end{aligned}
(2.46)