Constructing realistic Szekeres models from initial and final data

The Szekeres family of inhomogeneous solutions, which are defined by six arbitrary metric functions, offers a wide range of possibilities for modelling cosmic structure. Here we present a model construction procedure for the quasispherical case using given data at initial and final times. Of the six arbitrary metric functions, the three which are common to both Szekeres and Lemaître-Tolman models are determined by the model construction procedure of Krasinski and Hellaby. For the remaining three functions, which are unique to Szekeres models, we derive exact analytic expressions in terms of more physically intuitive quantities — density profiles and dipole orientation angles. Using MATLAB, we implement the model construction procedure and simulate the time evolution.