A contact lens manufacturing method comprising: a step of setting optical characteristics for a target subject; a step of setting a shape of an optical zone according to the set optical characteristics; a step of acquiring statistical information of a corneal surface shape obtained by a population of target eyes having the set optical characteristics; a definition step of defining a corneal surface shape corresponding to the population of the target eyes from the statistical information according to Expression 1 below which totals a mathematical expression that defines a rotationally symmetrical aspherical shape and a mathematical expression that defines a rotationally asymmetric aspherical shape, performing a fitting of the rotationally symmetrical aspherical shape by setting values of a radius of curvature and a conic constant according to Expression 2 below, and fitting a rotationally asymmetric element which is a residual difference during the fitting of the rotationally symmetrical aspherical shape using a Zernike polynomial shown in Expression 3 below which specifies a Zernike surface, the Zernike polynomial including at least Zernike terms where i = 1 to 20, 22 to 26, 30 to 33, 39 to 41, 49 to 50, and 60, wherein a lens rear surface shape corresponding to the corneal surface shape set at the definition step is used only for a rear surface of a peripheral area of the contact lens: Conic Formula (rotationally symmetrical aspherical surface formula) R _ average − 2 σ 1 ≤ R ≤ R _ average + 2 σ 1 C. C. _ average − 2 σ 2 ≤ C. C. ≤ C. C. _ average + 2 σ 2 Where, R: Radius of curvature of the contact lens, R_average: Average radius of curvature from the statistical information C.C.: Conic constant of the contact lens C.C._average: Average conic constant from the statistical information σ1: Standard deviation of radii from the statistical information σ2: Standard deviation of conic constants from the statistical information W ρ, θ = ∑ i = o ∞ N i P i Z i Where, W (ρ, θ): Zernike surface N i: i-th normalization coefficient P i: i-th standard Zernike coefficient Z i: i-th Zernike function. a step of setting optical characteristics for a target subject; a step of setting a shape of an optical zone according to the set optical characteristics; a step of acquiring statistical information of a corneal surface shape obtained by a population of target eyes having the set optical characteristics; a definition step of defining a corneal surface shape corresponding to the population of the target eyes from the statistical information according to Expression 1 below which totals a mathematical expression that defines a rotationally symmetrical aspherical shape and a mathematical expression that defines a rotationally asymmetric aspherical shape, performing a fitting of the rotationally symmetrical aspherical shape by setting values of a radius of curvature and a conic constant according to Expression 2 below, and fitting a rotationally asymmetric element which is a residual difference during the fitting of the rotationally symmetrical aspherical shape using a Zernike polynomial shown in Expression 3 below which specifies a Zernike surface, the Zernike polynomial including at least Zernike terms where i = 1 to 20, 22 to 26, 30 to 33, 39 to 41, 49 to 50, and 60, wherein a lens rear surface shape corresponding to the corneal surface shape set at the definition step is used only for a rear surface of a peripheral area of the contact lens: Conic Formula (rotationally symmetrical aspherical surface formula) R _ average − 2 σ 1 ≤ R ≤ R _ average + 2 σ 1 C. C. _ average − 2 σ 2 ≤ C. C. ≤ C. C. _ average + 2 σ 2 Where, R: Radius of curvature of the contact lens, R_average: Average radius of curvature from the statistical information C.C.: Conic constant of the contact lens C.C._average: Average conic constant from the statistical information σ1: Standard deviation of radii from the statistical information σ2: Standard deviation of conic constants from the statistical information W ρ, θ = ∑ i = o ∞ N i P i Z i Where, W (ρ, θ): Zernike surface N i: i-th normalization coefficient P i: i-th standard Zernike coefficient Z i: i-th Zernike function. R: Radius of curvature of the contact lens, R_average: Average radius of curvature from the statistical information C.C.: Conic constant of the contact lens C.C._average: Average conic constant from the statistical information σ1: Standard deviation of radii from the statistical information σ2: Standard deviation of conic constants from the statistical information W ρ, θ = ∑ i = o ∞ N i P i Z i Where, W (ρ, θ): Zernike surface N i: i-th normalization coefficient P i: i-th standard Zernike coefficient Z i: i-th Zernike function. W (ρ, θ): Zernike surface N i: i-th normalization coefficient P i: i-th standard Zernike coefficient Z i: i-th Zernike function.