Cosmological N-Body Initial Condition Maker

Given a power spectrum, generate the 3D coordinates for a lattice of perturbed particles and their velocities.

For computational cosmology and N-Body simulations that want a MATLAB initial condition generator. Given the number of softwares, no particular parameter file is generated other than two saved MAT files for position and velocity in cgs unit system.

Cosmological parameters used from
"Planck 2018 results. VI. Cosmological parameters"
https://arxiv.org/abs/1807.06209

Practical Zel'Dovich Approximation from
Brian W. O'Shea (April 15, 2003)
University of California San Diego
https://enzo.readthedocs.io/en/enzo-2.1.0/_downloads/makeics.ps

Default Analytic Power Spectrum from
"The Statistics of Peaks of Gaussian Random Fields" (Bardeen et al 1986)
https://ui.adsabs.harvard.edu/abs/1986ApJ...304...15B/abstract