An optimal trans-atmospheric vehicle and trajectory design are presented to simultaneously maximize the change in inclination angle and minimize total ΔV for an atmospheric skip entry maneuver. As an alternative to contemporary optimization techniques such as pseudospectral methods, a Design of Experiments approach featuring orthogonal arrays of experiments is implemented to solve the multi-objective maneuver optimization problem of atmospheric reentry. Conducting both Pareto front and main effects analysis, the optimal solution is determined by evaluating the relative performance of six design variables, to include mass, planform area, aerodynamic coefficients, perigee altitude, and bank angle. Depending on the chosen re-circularization altitude, the optimal design performing an out-of-plane skip entry maneuver can achieve an inclination change of 19.91 deg with 50-85% less ΔV than a simple plane change. Nomenclature a = semi-major axis, m C D , C L = drag and lift coefficients, respectively D = aerodynamic drag force, N g = gravitational acceleration, m∕s 2 h = altitude, m i = inclination angle, deg L = aerodynamic lift force, N m = vehicle mass, kg _ Q = heat flux, kW∕m 2 r = geocentric radial distance, m S = planform area, m 2 V = velocity, m∕s β = inverse scale height, 1∕m γ = flight-path angle, deg θ = longitude, deg μ = gravitational parameter, m 3 ∕s 2 ρ = atmospheric density, kg∕m 3 σ = bank angle, deg ϕ = latitude, deg ψ = heading angle, deg ω· = planetary rotation rate, rad∕s Subscripts f = final condition i = initial condition I = measured with respect to an inertial frame R = measured with respect to a rotating frame SL = sea-level condition = conditions for the Earth