Robust trajectory optimization for Mars ascent vehicles was achieved using polynomial chaos expansion and Picard iteration.