def alpha_from_gamma(gamma):
    """
    Spectral index from adiabatic index, Sedov-Taylor expansion.

    Parameters
    ----------
    gamma : adiabatic index
    """
    alpha = ((5./4)*gamma - 1.)/2.
    return alpha