WYPiWYG hyperelasticity is a data-driven, model-free computational procedure for finite element analysis of soft materials. The procedure does not assume the shape of the stored energy function and does not employ material parameters, predicting accurately any smooth prescribed behavior from a complete set of experimental tests. However, fuzzy experimental data may yield useless highly oscillatory, unstable stored energy functions, and classical curvature smoothing gives frequently unsatisfactory results. Aside, the possibility of having experimental data from different specimens for the same test was not considered in previous procedures. In this work we present a novel technique based on spline regression and smoothing penalization using stability conditions. In general, this procedure reduces noisy experimental data or data from multiple specimens for ulterior determination of the stored energy. The procedure only needs the solution of a linear system of equations. Instead of classical curvature-based smoothing, we employ a novel stability-based smoothing, determining for each branch of the uniaxial stress-strain curve the most restrictive stability condition during uniaxial and equibiaxial tests. The resulting stored energy functions are smooth and stable. The procedure has little sensitivity to the number of spline segments or to the choice of the penalization parameter, which are computed automatically.