The purpose of this study is to apply a stochastic optimal control method to the guidance of an aircraft. The aircraft is given optimal control to make a final approach through a microburst, a hazardous weather phenomenon for lowflying aircraft that may cause accidents during landing. As the determination of the size and the strength of a microburst always involves uncertainty, its precise detection is difficult. To minimize the risk of accidents, we must account for system uncertainty, a difficult task using conventional optimization techniques, which cannot handle random variables. This paper aims to solve this problem by employing a stochastic optimization algorithm that incorporates the generalized polynomial chaos method into a conventional direct method. The generalized polynomial chaos method is a numerical technique for solving stochastic differential equations. To demonstrate the effectiveness of the proposed algorithm, we conduct a numerical simulation in which an aircraft flies through a microburst and attempts to land. The results of the simulation successfully verify its effectiveness. Nomenclature a = vector of static variables that are independent of timeacceleration, 9.81 m∕s i = index of random variables p j = index of set of collocation points and weight functions L = lift, N l = index of root of the Nth-degree Legendre orthogonal polynomial M = aircraft mass, kg m, n = indices indicating degree of orthogonal polynomials N = number of random variables N 0 = number of Legendre-Gauss collocation points P = maximum degree of N-variate orthogonal polynomial p = vector of random variables that are independent of time Q = total number of sets of collocation points and weight functions q i = number of sets of collocation points and weight functions of each random variable r p = microburst radius parameter, m S = wing surface, m 2 T = thrust, N t = time, s u = control variable u d = deterministic optimal control u m = microburst intensity parameter, m∕s u s = stochastic optimal control u = approximating function of control u = optimal control Var = variance v a = airspeed, m∕s w x = wind speed in x direction, m∕s w y = wind speed in y direction, m∕s w z = wind speed in z direction, m∕s x = horizontal position, m x c = microburst center parameter, m x = state variablẽ x = approximated state variable by generalized polynomial chaos methodx = approximatedx by Legendre-Gauss pseudospectral method z = altitude, m z m = microburst height parameter, where horizontal wind speed is strongest, m α = angle of attack, rad Γ = N-dimensional probability space γ = flight-path angle relative to air, rad γ g = flight-path angle relative to ground, rad ρ = probability density functions ρ air = density of air from standard atmosphere model, kg∕m 3 σ = standard deviation τ = normalized time τ T = time constant of thrust, s τ α = time constant of angle of attack, s Ψ = multidimensional orthogonal polynomials ψ = one-dimensional orthogonal polynomials Subscripts com = commanded value f = final nom = necessary value with no d...