In 1990, based on numerical and formal asymptotic analysis, Ori and Piran predicted the existence of selfsimilar spacetimes, called relativistic Larson-Penston solutions, that can be suitably flattened to obtain examples of spacetimes that dynamically form naked singularities from smooth initial data, and solve the radially symmetric Einstein-Euler system. Despite its importance, a rigorous proof of the existence of such spacetimes has remained elusive, in part due to the complications associated with the analysis across the so-called sonic hypersurface. We provide a rigorous mathematical proof. Our strategy is based on a delicate study of nonlinear invariances associated with the underlying non-autonomous dynamical system to which the problem reduces after a selfsimilar reduction. Key technical ingredients are a monotonicity lemma tailored to the problem, an ad hoc shooting method developed to construct a solution connecting the sonic hypersurface to the so-called Friedmann solution, and a nonlinear argument to construct the maximal analytic extension of the solution. Finally, we reformulate the problem in double-null gauge to flatten the selfsimilar profile and thus obtain an asymptotically flat spacetime with an isolated naked singularity.