In the series of lectures, we will discuss probability laws of random points, curves, and surfaces. Starting from a brief review of the notion of martingales, one-dimensional Brownian motion (BM), and the D-dimensional Bessel processes, BES D , D ≥ 1, first we study Dyson's Brownian motion model with parameter β > 0, DYS β , which is a one-parameter family of repulsively interacting N Brownian motions on R, N ∈ N := {1, 2, . . . }, and is regarded as multivariate extensions of BES D with the relation β = D − 1. In particular, the determinantal structures are proved for DYS 2 , which is realized as the eigenvalue process of Hermitian-matrixvalued BM studied in random matrix theory. Next, using the reproducing kernels of Hilbert function spaces, the Gaussian analytic functions (GAFs) are defined on a unit disk D and an annulus A q := {z ∈ C : q < |z| < 1}, q ∈ (0, 1), which provide models of random surfaces. As zeros of the GAFs, determinantal point processes and permanental-determinantal point processes are obtained, which have symmetry and invariance associated with conformal transformations. Then, the Schramm-Loewner evolution with parameter κ > 0, SLE κ , is introduced, which is driven by a BM on R and generates a family of conformally invariant probability laws of random curves on the upper half complex plane H. We regard SLE κ as a complexification of BES D with the relation κ = 4/(D − 1). The last topic of lectures is the construction of the multiple SLE κ , which is driven by the N -particle DYS β on R and generates N interacting random curves in H. There, we define the Gaussian free field (GFF) and its generalization called the imaginary surface with parameter χ, which are considered as the distribution-valued random fields on H. Under the relation χ = 2/ √ κ − κ/ √ 2, we characterize the SLE/GFF coupling studied by Dubédat, Sheffield, and Miller by its temporal stationarity, and extend it to multiple cases. We prove that the multiple SLE/GFF coupling is established, if and only if the driving N -particle process on R is identified with DYS β with the relation β = 8/κ. This relation between parameters is very simple, but highly nontrivial, since if we simply combine the relations β = D − 1 and κ = 4/(D − 1) mentioned above, we will have a different relation β = 4/κ. Under the present multiple SLE/GFF coupling with β = 8/κ, we can prove that the multiple SLE driven by DYS β exhibits the phase transitions at κ = 4 and κ = 8 in the similar way to the original SLE κ with a single SLE curve.