Comprehensive Guide to Mathematical Analysis II: Integral Calculus, Series, and Differential Equations

Fundamental Definitions of Intervallic Subdivisions and Step Functions

In the study of integral calculus, we begin by examining the segment [a,b][a,b] as a closed and bounded interval in R\mathbb{R} where a<ba < b. A subdivision of this interval, denoted as σ={x0,x1,…,xn}\sigma = \{x_0, x_1, \dots, x_n\}, is a finite suite of elements such that a=x0<x1<⋯<xn=ba = x_0 < x_1 < \dots < x_n = b. The "pas" or mesh of this subdivision, denoted by δ(σ)\delta(\sigma), is defined as the maximum length of the sub-intervals: δ(σ)=max⁡0≤i≤n−1∣xi+1−xi∣\delta(\sigma) = \max_{0\leq i \leq n-1} |x_{i+1} - x_i|. A uniform subdivision of order n∈N∗n \in \mathbb{N}^* specifically divides the interval into nn sub-intervals of equal length, where δ(σn)=b−an\delta(\sigma_n) = \frac{b-a}{n} and each point is defined by the formula xi=a+ib−anx_i = a + i\frac{b-a}{n} for i=0,…,ni = 0, \dots, n. A subdivision σ\sigma is considered finer than another subdivision σ′\sigma' if σ⊂σ′\sigma \subset \sigma', which implies that δ(σ)≤δ(σ′)\delta(\sigma) \leq \delta(\sigma').

A function f:I→Rf: I \rightarrow \mathbb{R} is defined as a step function (fonction en escalier) if there exists a subdivision σ=(xi)0≤i≤n\sigma = (x_i)_{0\leq i \leq n} adapted to ff, such that ff remains constant on every open partial interval ]xi,xi+1[]x_i, x_{i+1}[. The choice of an adapted subdivision is not unique; if σ\sigma is adapted to ff, any finer subdivision σ′\sigma' is also adapted. For two step functions ff and gg with respective subdivisions σ\sigma and σ′\sigma', the union σ∪σ′\sigma \cup \sigma' serves as a subdivision adapted to both. The set of all step functions on [a,b][a,b], denoted E([a,b])\mathcal{E}([a,b]), forms a vector subspace of the set of all functions F([a,b])\mathcal{F}([a,b]), meaning that for any α,β∈R\alpha, \beta \in \mathbb{R} and f,g∈E(I)f, g \in \mathcal{E}(I), the combination αf+βg∈E(I)\alpha f + \beta g \in \mathcal{E}(I).

The Integral of Step Functions and Its Properties

For a step function f:[a,b]→Rf: [a,b] \rightarrow \mathbb{R} and an adapted subdivision σ=(xi)0≤i≤n\sigma = (x_i)_{0\leq i \leq n}, where f(x)=cif(x) = c_i on the interval ]xi,xi+1[]x_i, x_{i+1}[, the integral is defined by the expression I(σ,f)=∑i=0n−1ci(xi+1−xi)I(\sigma, f) = \sum_{i=0}^{n-1} c_i(x_{i+1} - x_i). Crucially, this value does not depend on the choice of the adapted subdivision σ\sigma. The integral of ff on [a,b][a,b] is denoted as ∫abf(x)dx\int_a^b f(x)dx. Geometrically, this value corresponds to the algebraic area between the x-axis, the graph of the function, and the lines x=ax = a and x=bx = b.

The integral of step functions follows several fundamental properties. Linearity dictates that ∫ab(αf(x)+βg(x))dx=α∫abf(x)dx+β∫abg(x)dx\int_a^b (\alpha f(x) + \beta g(x))dx = \alpha \int_a^b f(x)dx + \beta \int_a^b g(x)dx. Positivity ensures that if f≥0f \geq 0 on [a,b][a,b], then ∫abf(x)dx≥0\int_a^b f(x)dx \geq 0. Monotonicity implies that if f≥gf \geq g, then ∫abf(x)dx≥∫abg(x)dx\int_a^b f(x)dx \geq \int_a^b g(x)dx. The triangle inequality relates the absolute value: ∣∫abf(x)dx∣≤∫ab∣f(x)∣dx| \int_a^b f(x)dx | \leq \int_a^b |f(x)|dx. If ∣f(x)∣≤λ|f(x)| \leq \lambda for all x∈[a,b]x \in [a,b], then ∣∫abf(x)dx∣≤λ(b−a)| \int_a^b f(x)dx | \leq \lambda(b-a). Additionally, Chasles' relation states that for any c∈]a,b[c \in ]a,b[, ∫abf(x)dx=∫acf(x)dx+∫cbf(x)dx\int_a^b f(x)dx = \int_a^c f(x)dx + \int_c^b f(x)dx. Finally, if two step functions ff and gg are equal except for at a finite number of points, their integrals are identical.

Darboux Sums and the Riemann Integral Definition

To extend integration to broader classes of functions, we consider a bounded function f:[a,b]→Rf: [a,b] \rightarrow \mathbb{R}. For a uniform subdivision σn\sigma_n with pas b−an\frac{b-a}{n}, we define mni=inf⁡x∈]xi,xi+1[f(x)m_{ni} = \inf_{x \in ]x_i, x_{i+1}[} f(x) and Mni=sup⁡x∈]xi,xi+1[f(x)M_{ni} = \sup_{x \in ]x_i, x_{i+1}[} f(x). These generate two step functions, fn(x)=mnif_n(x) = m_{ni} and Fn(x)=MniF_n(x) = M_{ni} on ]xi,xi+1[]x_i, x_{i+1}[. The lower and upper Darboux sums are given by sn=∫abfn(x)dx=b−an∑i=0n−1mnis_n = \int_a^b f_n(x)dx = \frac{b-a}{n} \sum_{i=0}^{n-1} m_{ni} and Sn=∫abFn(x)dx=b−an∑i=0n−1MniS_n = \int_a^b F_n(x)dx = \frac{b-a}{n} \sum_{i=0}^{n-1} M_{ni}, respectively. Since mni≤Mnim_{ni} \leq M_{ni}, it is always true that sn≤Sns_n \leq S_n.

A bounded function ff is Riemann-integrable if lim⁡n→∞(Sn−sn)=0\lim_{n \to \infty} (S_n - s_n) = 0. The common limit of sns_n and SnS_n is the Riemann integral of ff on the segment [a,b][a,b], denoted ∫abf(t)dt\int_a^b f(t)dt. For example, the function f(x)=2xf(x) = 2x on [0,1][0,1] is Riemann-integrable with an integral of 1. Conversely, the Dirichlet function, defined as f(x)=1f(x) = 1 if x∈Qx \in \mathbb{Q} and 00 if x∉Qx \notin \mathbb{Q}, is not Riemann-integrable because on any interval, mni=0m_{ni} = 0 and Mni=1M_{ni} = 1, leading to sn=0s_n = 0 and Sn=1S_n = 1 for all n∈N∗n \in \mathbb{N}^*.

Families of Riemann-Integrable Functions and Sums of Riemann

Several specific classes of functions are guaranteed to be Riemann-integrable. These include monotone functions (whether increasing or decreasing) and continuous functions. For these functions, the integral can be calculated using Riemann sums: ∫abf(x)dx=lim⁡n→∞b−an∑i=0n−1f(a+ib−an)=lim⁡n→∞b−an∑i=1nf(a+ib−an)\int_a^b f(x)dx = \lim_{n \to \infty} \frac{b-a}{n} \sum_{i=0}^{n-1} f(a + i\frac{b-a}{n}) = \lim_{n \to \infty} \frac{b-a}{n} \sum_{i=1}^{n} f(a + i\frac{b-a}{n}). The Mean Value Theorem for integrals states that for a continuous function f∈C([a,b])f \in \mathcal{C}([a,b]), there exists a point c∈[a,b]c \in [a,b] such that 1b−a∫abf(x)dx=f(c)\frac{1}{b-a} \int_a^b f(x)dx = f(c).

Riemann-integrable functions share the same properties as step function integrals regarding linearity, positivity, growth, absolute value, and Chasles' relation. By convention, we define ∫abf(x)dx=−∫baf(x)dx\int_a^b f(x)dx = -\int_b^a f(x)dx and ∫aaf(x)dx=0\int_a^a f(x)dx = 0, which allows the Chasles relation ∫uvf(x)dx=∫uwf(x)dx+∫wvf(x)dx\int_u^v f(x)dx = \int_u^w f(x)dx + \int_w^v f(x)dx to hold for any u,v,w∈[a,b]u, v, w \in [a,b]. Like step functions, changing a Riemann-integrable function at a finite number of points does not change the value of its integral.

Primitives and the Fundamental Theorem of Analysis

A function F:I→RF: I \rightarrow \mathbb{R} is a primitive of ff if FF is derivable and F′=fF' = f on II. Primitives are unique up to an additive constant. The Fundamental Theorem of Analysis establishes that derivation and integration are reciprocal operations. Specifically, if ff is Riemann-integrable and has a primitive FF, then ∫abf(x)dx=F(b)−F(a)\int_a^b f(x)dx = F(b) - F(a). Furthermore, if ff is continuous on [a,b][a,b] and we define F(x)=∫x0xf(t)dtF(x) = \int_{x_0}^x f(t)dt, then FF is the unique primitive of ff that vanishes at x0x_0.

Indefinite integrals, denoted ∫f(x)dx\int f(x)dx, represent any primitive of ff. Integration techniques include integration by parts: ∫abf(x)g′(x)dx=[f(x)g(x)]ab−∫abf′(x)g(x)dx\int_a^b f(x)g'(x)dx = [f(x)g(x)]_a^b - \int_a^b f'(x)g(x)dx, provided ff and gg are of class C1\mathcal{C}^1. Another key method is the change of variable theorem: if ϕ:[α,β]→[a,b]\phi: [\alpha, \beta] \rightarrow [a,b] is of class C1\mathcal{C}^1, then ∫αβf(ϕ(x))ϕ′(x)dx=∫ϕ(α)ϕ(β)f(u)du\int_{\alpha}^{\beta} f(\phi(x))\phi'(x)dx = \int_{\phi(\alpha)}^{\phi(\beta)} f(u)du. If ϕ\phi is also a bijection, ∫abf(x)dx=∫ϕ−1(a)ϕ−1(b)f(ϕ(u))ϕ′(u)du\int_a^b f(x)dx = \int_{\phi^{-1}(a)}^{\phi^{-1}(b)} f(\phi(u))\phi'(u)du.

Advanced Integration of Complex Functions

Rational fractions P(x)Q(x)\frac{P(x)}{Q(x)} are integrated by decomposing them into a polynomial and simple elements of two types: γ(x−x0)k\frac{\gamma}{(x-x_0)^k} or αx+β(ax2+bx+c)k\frac{\alpha x + \beta}{(ax^2 + bx + c)^k} with b2−4ac<0b^2 - 4ac < 0. Specific results include ∫dxx−x0=ln⁡∣x−x0∣+C\int \frac{dx}{x-x_0} = \ln|x-x_0| + C and ∫dx(x−x0)k=(x−x0)−k+1−k+1+C\int \frac{dx}{(x-x_0)^k} = \frac{(x-x_0)^{-k+1}}{-k+1} + C for k≥2k \geq 2. For quadratic denominators, transformations into ln⁡\ln or arctan⁡\arctan forms are required.

Trigonometric functions ∫F(cos⁡(x),sin⁡(x))dx\int F(\cos(x), \sin(x))dx can be solved using the change of variable t=tan⁡(x2)t = \tan(\frac{x}{2}), where cos⁡(x)=1−t21+t2\cos(x) = \frac{1-t^2}{1+t^2}, sin⁡(x)=2t1+t2\sin(x) = \frac{2t}{1+t^2}, and dx=2dt1+t2dx = \frac{2dt}{1+t^2}. Alternatively, Bioche's rules are used: if ω(x)=f(x)dx\omega(x) = f(x)dx, then t=cos⁡(x)t = \cos(x) is used if ω(−x)=ω(x)\omega(-x) = \omega(x), t=sin⁡(x)t = \sin(x) if ω(π−x)=ω(x)\omega(\pi-x) = \omega(x), and t=tan⁡(x)t = \tan(x) if ω(π+x)=ω(x)\omega(\pi+x) = \omega(x). For rational functions involving exponentials, ∫F(ex)dx\int F(e^x)dx, the change of variable t=ext = e^x allows substitution via dx=dttdx = \frac{dt}{t}.

Nature and Properties of Generalized Integrals

A function ff is locally integrable on II (denoted Rloc(I)R_{loc}(I)) if it is Riemann-integrable on every segment of II. A generalized (or improper) integral ∫abf(t)dt\int_a^b f(t)dt is defined for intervals that are not necessarily closed or for functions that are not bounded. If F(x)=∫axf(t)dtF(x) = \int_a^x f(t)dt, the integral is convergent if lim⁡k→bF(x)\lim_{k \to b} F(x) is a finite limit LL. Otherwise, it is divergent. Examples include ∫0+∞dt1+t2=π2\int_0^{+\infty} \frac{dt}{1+t^2} = \frac{\pi}{2} (convergent) and ∫e+∞dtt=+∞\int_e^{+\infty} \frac{dt}{t} = +\infty (divergent).

If the integral has two "uncertain points," it is doublement généralisée. For f∈Rloc(]a,b[)f \in R_{loc}(]a,b[), the integral ∫abf(t)dt\int_a^b f(t)dt converges if there exists a c∈]a,b[c \in ]a,b[ such that both ∫acf(t)dt\int_a^c f(t)dt and ∫cbf(t)dt\int_c^b f(t)dt converge. Their sum provides the final value. An integral is "falsely generalized" if ff has a finite limit at the uncertain endpoint, in which case it is always convergent.

Convergence Criteria for Positive and Variable Functions

For positive functions (f≥0f \geq 0), ∫abf(t)dt\int_a^b f(t)dt converges if and only if the integral function F(x)F(x) is bounded. The comparison theorem states that if 0≤f≤g0 \leq f \leq g, the convergence of ∫g\int g implies the convergence of ∫f\int f, and the divergence of ∫f\int f implies the divergence of ∫g\int g. The equivalents theorem allows replacing ff with an equivalent function gg (f∼bgf \sim_b g) to determine nature. Standard scales include the Riemann Scale: ∫a+∞dxxα\int_a^{+\infty} \frac{dx}{x^{\alpha}} converges if and only if α>1\alpha > 1, whereas ∫0adxxα\int_0^a \frac{dx}{x^{\alpha}} and ∫badx(x−b)α\int_b^a \frac{dx}{(x-b)^{\alpha}} converge if and only if α<1\alpha < 1. The Bertrand scale ∫e+∞dxxαln⁡β(x)\int_e^{+\infty} \frac{dx}{x^{\alpha} \ln^{\beta}(x)} converges if α>1\alpha > 1 or if α=1\alpha = 1 and β>1\beta > 1.

For functions with variable signs, an integral is absolutely convergent if ∫ab∣f(x)∣dx\int_a^b |f(x)|dx converges. Absolute convergence implies convergence. If an integral converges but is not absolutely convergent, it is semi-convergente. The Cauchy criterion provides a theoretical basis for convergence without computing the limit. The Abel Rule for integrals suggests that if ff is decreasing towards zero and the integral of gg is bounded, then ∫f(t)g(t)dt\int f(t)g(t)dt converges.

Numerical Series and General Convergence Properties

A series ∑un\sum u_n is defined as the suite of partial sums Sn=∑k=0nukS_n = \sum_{k=0}^n u_k. The series converges if SnS_n has a finite limit as n→∞n \to \infty. The sum is denoted ∑n=0+∞un\sum_{n=0}^{+\infty} u_n, and the remainder is Rn=S−SnR_n = S - S_n. A geometric series ∑qn\sum q^n converges if and only if ∣q∣<1|q| < 1, with a sum of 11−q\frac{1}{1-q}. Telescopic series ∑(an+1−an)\sum (a_{n+1} - a_n) converge if and only if the sequence (an)(a_n) converges to a limit ll, yielding a sum of l−a0l - a_0.

A necessary condition for the convergence of any series is lim⁡n→∞un=0\lim_{n \to \infty} u_n = 0. If this limit is non-zero or nonexistent, the series is said to "diverge grossièrement." However, this is not a sufficient condition, as seen with the harmonic series ∑1n\sum \frac{1}{n}, which diverges despite the limit being zero. Series obey linearity: if ∑un\sum u_n and ∑vn\sum v_n converge, then ∑(αun+βvn)\sum (\alpha u_n + \beta v_n) converges to α∑un+β∑vn\alpha \sum u_n + \beta \sum v_n.

Convergence Criteria for Series with Positive and Arbitrary Terms

For series with positive terms (un≥0u_n \geq 0), convergence is equivalent to partial sums being bounded. The Comparison Criterion states that if un≤vnu_n \leq v_n, then convergence of ∑vn\sum v_n implies convergence of ∑un\sum u_n. The Equivalents Criterion holds that if un∼vnu_n \sim v_n, the series share the same nature. Positive series can also be compared to generalized integrals of decreasing functions: ∑f(n)\sum f(n) is of the same nature as ∫0+∞f(t)dt\int_0^{+\infty} f(t)dt. Riemann series ∑1nα\sum \frac{1}{n^{\alpha}} converge for α>1\alpha > 1, and Bertrand series ∑1nαln⁡β(n)\sum \frac{1}{n^{\alpha} \ln^{\beta}(n)} converge if α>1\alpha > 1 or (α=1\alpha = 1 and β>1\beta > 1).

Additional tests include D'Alembert's Ratio Test: if lim⁡un+1un=l\lim \frac{u_{n+1}}{u_n} = l, the series converges if l<1l < 1 and diverges if l>1l > 1. Cauchy's Root Test: if lim⁡∑(un)1/n=l\lim ∑(u_n)^{1/n} = l, the series converges for l<1l < 1 and diverges for l>1l > 1. For arbitrary terms, absolute convergence (∑∣un∣\sum |u_n|) implies convergence. The Leibniz Criterion for alternating series ∑(−1)nan\sum (-1)^n a_n states that if ana_n is a positive, decreasing sequence tending to zero, the series converges, and the remainder satisfies ∣Rn∣≤∣un+1∣|R_n| \leq |u_{n+1}|.

Equations Différentielles: First and Second Order

An ordinary differential equation (ODE) of order nn involves an unknown function y(x)y(x) and its derivatives up to y(n)y^{(n)}. First-order equations include those with separable variables, which take the form g(y)y′=f(x)g(y)y' = f(x). These are solved by calculating primitives: ∫g(y)dy=∫f(x)dx\int g(y)dy = \int f(x)dx, leading to G(y)=F(x)+CG(y) = F(x) + C. Linear first-order equations are defined as a(x)y′+b(x)y=c(x)a(x)y' + b(x)y = c(x). One first solves the homogeneous equation a(x)y′+b(x)y=0a(x)y' + b(x)y = 0, which has the solution y0=Ke−F(x)y_0 = K e^{-F(x)} where F(x)=∫b(x)a(x)dxF(x) = \int \frac{b(x)}{a(x)} dx. A particular solution ypy_p is then found via the Variation of Constants method, setting yp=K(x)e−F(x)y_p = K(x) e^{-F(x)}. The general solution is y=y0+ypy = y_0 + y_p.

Second-order linear ODEs with constant coefficients are of the form ay′′+by′+cy=f(x)ay'' + by' + cy = f(x). The associated characteristic equation is ar2+br+c=0ar^2 + br + c = 0. If the discriminant Δ>0\Delta > 0, the homogeneous solution is y0=λer1x+μer2xy_0 = \lambda e^{r_1 x} + \mu e^{r_2 x}; if Δ=0\Delta = 0, then y0=(λ+μx)er0xy_0 = (\lambda + \mu x) e^{r_0 x}; if Δ<0\Delta < 0, the solution is y0=eαx(λcos⁡(βx)+μsin⁡(βx))y_0 = e^{\alpha x} (\lambda \cos(\beta x) + \mu \sin(\beta x)) where the roots are r=α±iβr = \alpha \pm i\beta. Particular solutions are found based on the form of f(x)f(x). For f(x)=eαxPn(x)f(x) = e^{\alpha x} P_n(x), one looks for yp=eαxxmQn(x)y_p = e^{\alpha x} x^m Q_n(x), where mm depends on whether α\alpha is a root of the characteristic equation. For trigonometric RHS (eαxPn(x)cos⁡(βx)e^{\alpha x} P_n(x) \cos(\beta x)), the choice relates to whether α+iβ\alpha + i\beta is a root.