Section 13.7

Approximate Integration of Differential Equations

[29/62%]
Show
LaTeX
Problem 5342

Find the power-series expansion of the general integral of the equation y′−xy=1y' - xy = 1.

?
Problem 5343

Find the power-series expansion of the general integral of the equation (1+x)y′−ny=0(1+x)y' - ny = 0.

?
Problem 5344

Find the power-series expansion of the general integral of the equation y′′+y=0y'' + y = 0.

?
Problem 5345

Find the power-series expansion of the general integral of the equation y′′−xy=0y'' - xy = 0.

?
Problem 5346

Prove that the corresponding particular integral of the differential equation

x(x−1)y′′+((α+β+1)x−γ)y′+αβy=0 x(x-1)y'' + \left((\alpha +\beta +1)x-\gamma \right)y' + \alpha \beta y = 0

can be represented by the indicated power series:

y=1+α⋅β1⋅γx+α(α+1)1⋅2⋅β(β+1)γ(γ+1)x2+α(α+1)(α+2)1⋅2⋅3⋅β(β+1)(β+2)γ(γ+1)(γ+2)x3+… y = 1 + \frac{\alpha \cdot \beta }{1\cdot \gamma }x + \frac{\alpha (\alpha +1)}{1\cdot 2}\cdot \frac{\beta (\beta +1)}{\gamma (\gamma +1)}x^{2} + \frac{\alpha (\alpha +1)(\alpha +2)}{1\cdot 2\cdot 3}\cdot \frac{\beta (\beta +1)(\beta +2)}{\gamma (\gamma +1)(\gamma +2)}x^{3} + \ldots

(Gauss's hypergeometric series).

?
Problem 5347

Prove that the corresponding particular integral of the differential equation

x2y′′+xy′+(x2−n2)y=0 x^{2}y'' + xy' + (x^{2}-n^{2})y = 0

can be represented by the indicated power series:

y=∑ν=0∞(−1)ν(x2)n+2νν! Γ(n+ν+1) y = \sum _{\nu =0}^{\infty } \frac{(-1)^{\nu }\left(\frac{x}{2}\right)^{n+2\nu }}{\nu !\, \Gamma (n+\nu +1)}

(the Bessel function Jn(x)J_{n}(x)).

?
Problem 5348

Prove that the corresponding particular integral of the differential equation

x(x4−1)u′′−(3x4−1)u′+x3u=0 x(x^{4}-1)u'' - (3x^{4}-1)u' + x^{3}u = 0

can be represented by the indicated power series:

u=1+(14)2x4+(1⋅54⋅8)2x8+(1⋅5⋅94⋅8⋅12)2x12+…;∣x∣<1. u = 1 + \left(\frac{1}{4}\right)^{2}x^{4} + \left(\frac{1\cdot 5}{4\cdot 8}\right)^{2}x^{8} + \left(\frac{1\cdot 5\cdot 9}{4\cdot 8\cdot 12}\right)^{2}x^{12} + \ldots ; \qquad \left|x\right| < 1.
?
Problem 5349

Prove that the corresponding particular integral of the differential equation

x2(x4−1)v′′′+3x(3x4−1)v′′+(19x4−1)v′+8x3v=0 x^{2}(x^{4}-1)v''' + 3x(3x^{4}-1)v'' + (19x^{4}-1)v' + 8x^{3}v = 0

can be represented by the indicated power series:

v=1+(12)3x4+(1⋅32⋅4)3x8+(1⋅3⋅52⋅4⋅6)3x12+…;∣x∣<1. v = 1 + \left(\frac{1}{2}\right)^{3}x^{4} + \left(\frac{1\cdot 3}{2\cdot 4}\right)^{3}x^{8} + \left(\frac{1\cdot 3\cdot 5}{2\cdot 4\cdot 6}\right)^{3}x^{12} + \ldots ; \qquad \left|x\right| < 1.
?
Problem 5350

Prove Gauss's equality:

(1+(14)2x4+(1⋅54⋅8)2x8+(1⋅5⋅94⋅8⋅12)2x12+…)2=1+(12)3x4+(1⋅32⋅4)3x8+(1⋅3⋅52⋅4⋅6)3x12+…;∣x∣<1. \left(1 + \left(\frac{1}{4}\right)^{2}x^{4} + \left(\frac{1\cdot 5}{4\cdot 8}\right)^{2}x^{8} + \left(\frac{1\cdot 5\cdot 9}{4\cdot 8\cdot 12}\right)^{2}x^{12} + \ldots \right)^{2} = 1 + \left(\frac{1}{2}\right)^{3}x^{4} + \left(\frac{1\cdot 3}{2\cdot 4}\right)^{3}x^{8} + \left(\frac{1\cdot 3\cdot 5}{2\cdot 4\cdot 6}\right)^{3}x^{12} + \ldots ; \qquad \left|x\right|<1.
?
Problem 5351

Prove that, applying Picard's method to the equation y′=yy' = y, with a=0a=0, b=1b=1, we obtain:

y=1+x+x22!+…+xn−1(n−1)!+xnn!+… y = 1 + x + \frac{x^{2}}{2!} + \ldots + \frac{x^{n-1}}{(n-1)!} + \frac{x^{n}}{n!} + \ldots
?
Remark.
?

Picard's method of finding the integral of the equation y′=φ(x,y)y' = \varphi (x,y), which equals bb at x=ax=a, consists in successively constructing a sequence of functions y0,y1,y2,…y_{0}, y_{1}, y_{2}, \ldots by the formula yn=b+∫axφ(x,yn−1) dxy_{n} = b + \displaystyle \int_{a}^{x} \varphi (x,y_{n-1})\, dx; y0=by_{0}=b. If yn→yy_{n} \to y as n→∞n \to \infty, then yy is the desired integral.

Problem 5352

Prove that in the preceding example (problem 5351), yn→yy_{n} \to y as n→∞n \to \infty for an arbitrary choice of the function y0y_{0}.

?
Problem 5353

Prove that, replacing the equation

y′′+p(x)y′+q(x)y=0 y'' + p(x)y' + q(x)y = 0

by the system of equations

y′=u;u′=−p(x)u−q(x)y y' = u; \qquad u' = -p(x)u - q(x)y

and applying Picard's method to solve these equations, we obtain a convergent process if the radius vector from the origin to the point xx does not pass, in the complex plane of the variable xx, through any singular point of the functions p(x)p(x) and q(x)q(x).

?
Problem 5354

Solve the equation y′−y=axyy' - y = axy by the method of successive approximations, setting y0=0y_{0}=0 and computing the quantities yny_{n} by the formula yn′−yn=axyn−1y_{n}' - y_{n} = axy_{n-1} (at x=0x=0 we must have y=Cy=C).

?
Problem 5355

In the same way (see problem 5354), solve the equation y′′+a2y=σxyy'' + a^{2}y = \sigma xy, where σ\sigma is a small parameter. Here y=0y=0, y′=ay'=a at x=0x=0. It is useful to apply the theorem: if y′′+a2y=f(x)y''+a^{2}y=f(x), then

y=C1sin⁡ax+C2cos⁡ax+1a∫0xf(t)sin⁡a(x−t) dt. y = C_{1}\sin ax + C_{2}\cos ax + \frac{1}{a}\int _{0}^{x} f(t)\sin a(x-t)\, dt.
?
Problem 5356

Using the substitution x=n+1nx = n + \dfrac {1}{\sqrt{n}}, show that the integral of the equation y′′+xy=0y'' + xy = 0 on the interval (n,π+n)(n, \pi +n), where nn is large, has approximately n\sqrt{n} roots.

?
Problem 5357

Prove that the integral of Bessel's equation

x2y′′+xy′+(x2−n2)y=0 x^{2}y'' + xy' + (x^{2}-n^{2})y = 0

on an interval (a,a+b)(a, a+b) far from the origin has approximately bπ\dfrac {b}{\pi } roots.

?
Problem 5358

The integral of the equation y′′+y=σy′2y'' + y = \sigma y'^{2} has initial conditions: x=0x=0, y=1y=1, y′=0y'=0. Find the first three terms of its expansion in powers of the small parameter σ\sigma.

?
Problem 5359

The Chebyshev--Hermite polynomials Hn(x)=ex22n!dne−x22dxnH_{n}(x) = \dfrac {e^{\frac{x^{2}}{2}}}{n!}\dfrac {d^{n}e^{-\frac{x^{2}}{2}}}{dx^{n}} satisfy the equation

Hn′′(x)−xHn′(x)+nHn(x)=0. H_{n}''(x) - xH_{n}'(x) + nH_{n}(x) = 0.

By the substitution Hn(x)=f(x)⋅yH_{n}(x) = f(x)\cdot y, transform this equation to the form

y′′+(−x24+n+12)y=0. y'' + \left(-\frac{x^{2}}{4} + n + \frac{1}{2}\right)y = 0.
?
Problem 5360

Find the asymptotic expression, as n→∞n \to \infty, for the Chebyshev--Hermite polynomials Hn(x)=ex22n!dne−x22dxnH_{n}(x) = \dfrac {e^{\frac{x^{2}}{2}}}{n!}\dfrac {d^{n}e^{-\frac{x^{2}}{2}}}{dx^{n}} (see problem 5359).

?
Problem 5361

Find, for 0<x<π0<x<\pi, the integral of the equation

y′′+4y=πx−x2, y'' + 4y = \pi x - x^{2},

by expanding it in a series of the form ∑ansin⁡nx\sum a_{n}\sin nx. Boundary conditions: y=0y=0 at x=0x=0 and x=πx=\pi.

?
Problem 5362

Find the integral of the equation yIV+y=π−∣x∣y^{\mathrm{IV}} + y = \pi - \left|x\right| for −π<x<π-\pi <x<\pi and boundary conditions y′=y′′′=0y' = y''' = 0 at x=±πx = \pm \pi.

?
Problem 5363

Construct a table of values of yy at intervals of 0,050{,}05, given y′+2xy=0y'+2xy=0 and y=1y=1 at x=0x=0. Here the values y−2,y−1,y0,y1,y2y_{-2}, y_{-1}, y_{0}, y_{1}, y_{2} are to be computed by means of a series expansion, in the Adams formula terms up through the third differences are to be retained, and the computations are to be carried out to the sixth decimal digit.

?
Remark.
?

Numerical integration of the differential equation y′=f(x,y)y'=f(x,y) consists in this: given the initial condition y=y0y=y_{0} at x=x0x=x_{0}, one constructs a table of values of yy for values of xx equal to the numbers x0,x0+h,x0+2h,…x_{0}, x_{0}+h, x_{0}+2h, \ldots, where hh is a chosen number, usually fairly small.

The Adams--Störmer formula: Δyn=un+12Δun−1+512Δ2un−2+38Δ3un−3+251720Δ4un−4+13Δ5un−5+…\Delta y_{n} = u_{n} + \dfrac {1}{2}\Delta u_{n-1} + \dfrac {5}{12}\Delta^{2}u_{n-2} + \dfrac {3}{8}\Delta^{3}u_{n-3} + \dfrac {251}{720}\Delta^{4}u_{n-4} + \dfrac {1}{3}\Delta^{5}u_{n-5} + \ldots, where u=hf(x,y)u = hf(x,y); this formula makes it possible to compute yn+1y_{n+1}.

Problem 5364

Find by the Adams--Störmer method, at x=2x=2, the integral of the equation y′=x+yy' = x+y that equals unity at x=1x=1; here use the formula for Δy\Delta y taken up through the fourth differences inclusive, and carry out the computations to the fifth digit.

?
Problem 5365

Find a table of values of the integral of the equation

y′=x+yforx=0; 0,1; 0,2; … 1, y' = \sqrt{x} + \sqrt{y} \quad \text{for} \quad x = 0;\ 0{,}1;\ 0{,}2;\ \ldots \ 1,

given y=0y=0 at x=0x=0.

?
Remark.
?

The expansion of yy near the start has the form y=23x32+4723x74+17x2+14923x94+…y = \dfrac {2}{3}x^{\frac32} + \dfrac {4}{7}\sqrt{\dfrac {2}{3}}x^{\frac74} + \dfrac {1}{7}x^{2} + \dfrac {1}{49}\sqrt{\dfrac {2}{3}}x^{\frac94} + \ldots It is advantageous to introduce a new variable by setting x=ξ4x = \xi^{4}.

Problem 5366

Construct a table of values of the functions sin⁡x\sin x and cos⁡x\cos x at intervals of 0,10{,}1 from 00 to 11, using the equations (sin⁡x)′=cos⁡x(\sin x)' = \cos x; (cos⁡x)′=−sin⁡x(\cos x)' = -\sin x and applying the Adams--Störmer method.

?
Problem 5367

Construct a table of values of the Bessel functions J0(x)J_{0}(x) and J0′(x)J_{0}'(x) for values of xx at intervals of one tenth from zero to one, replacing the equation for J0(x)J_{0}(x), which has the form

xJ0′′(x)+J0′(x)+xJ0(x)=0, xJ_{0}''(x) + J_{0}'(x) + xJ_{0}(x) = 0,

by the system: xv′+v+xu=0xv' + v + xu = 0; u′=vu' = v.

?
Remark.
?

For second-order equations of the form y′′=f(x,y,y′)y'' = f(x,y,y') one may use the formulas:

Δ2yn−1=h2fn+112h2Δ2fn−2+112h2Δ3fn−3+19240h2Δ4fn−4+340h2Δ5fn−5,Δyn=hyn′+h22fn+h26Δfn−1+h28Δ2fn−2+19h2180Δ3fn−3+3h232Δ4fn−4,Δyn′=hfn+h2Δfn−1+512hΔ2fn−2+38hΔ3fn−3+251720hΔ4fn−4. \begin{align} \Delta ^{2}y_{n-1} & = h^{2}f_{n} + \frac{1}{12}h^{2}\Delta ^{2}f_{n-2} + \frac{1}{12}h^{2}\Delta ^{3}f_{n-3} + \frac{19}{240}h^{2}\Delta ^{4}f_{n-4} + \frac{3}{40}h^{2}\Delta ^{5}f_{n-5}, \\ \Delta y_{n} & = hy_{n}' + \frac{h^{2}}{2}f_{n} + \frac{h^{2}}{6}\Delta f_{n-1} + \frac{h^{2}}{8}\Delta ^{2}f_{n-2} + \frac{19h^{2}}{180}\Delta ^{3}f_{n-3} + \frac{3h^{2}}{32}\Delta ^{4}f_{n-4}, \\ \Delta y_{n}' & = hf_{n} + \frac{h}{2}\Delta f_{n-1} + \frac{5}{12}h\Delta ^{2}f_{n-2} + \frac{3}{8}h\Delta ^{3}f_{n-3} + \frac{251}{720}h\Delta ^{4}f_{n-4}. \end{align}
Problem 5368

Construct a table of values, at intervals of 0,10{,}1 on the interval (0,1)(0,1), of the function y=ch⁡xy = \operatorname {ch}x, satisfying the equation y′′=yy'' = y and the initial conditions: at x=1x=1, y=1y=1, y′=0y'=0.

?
Problem 5369

Construct a table of values of the function J0(x)J_{0}(x) for values of xx at intervals of 0,10{,}1, knowing that J0(x)J_{0}(x) satisfies the differential equation xJ0′′(x)+J0′(x)+xJ0(x)=0xJ_{0}''(x) + J_{0}'(x) + xJ_{0}(x) = 0, with initial values at x=0x=0: J0(0)=1J_{0}(0)=1, J0′(0)=0J_{0}'(0)=0.

?
Problem 5370

Knowing the values: cos⁡1=0,54030\cos 1 = 0{,}54030; cos⁡1,1=0,45360\cos 1{,}1 = 0{,}45360; cos⁡1,2=0,36236\cos 1{,}2 = 0{,}36236; cos⁡1,3=0,26750\cos 1{,}3 = 0{,}26750; cos⁡1,4=0,16997\cos 1{,}4 = 0{,}16997, using the equation y′′=−yy'' = -y, find the values of cos⁡x\cos x at x=1,5x=1{,}5; 1,61{,}6; 1,71{,}7; 1,81{,}8; 1,91{,}9; 2,02{,}0.

?