def hseq_profile(w, z, temperature, mstar, mu=2.279):

    from hyperion.util.constants import G, m_h, k
    from ..util.integrate import integrate, integrate_subset

    # Compute the integrand
    integrand = z / temperature / (w ** 2 + z ** 2) ** 1.5

    # Compute the integral for all cells
    # TODO - inefficient to compute integral from scratch - optimize
    i = np.array([integrate_subset(z, integrand, 0., zmax) for zmax in z])
    i[z < 0] = -i[z < 0]

    # Compute the factor for the integrand
    factor = G * mstar * mu * m_h / k

    # Compute the profile
    density = np.exp(-i * factor) / temperature

    # Normalize the density profile
    density = density / integrate(z, density)

    return density


# The mean molecular weight of H2 + He is given by:
#
# mu = 4 * (X + 1) / (X + 2)
#
# where X is the mass fraction of Helium to Hydrogren. Assuming
#
# X = 0.32476319350473615
#
# gives:
#
# mu = 2.279
Unknown, Nodule or Tumor. This information will be saved in the GeometryTopologyData that
dataset/ETHPy150Open hyperion-rt/hyperion/hyperion/model/helpers.py/hseq_profile--crossover with--dataset/ETHPy150Open acil-bwh/SlicerCIP/Scripted/CIP_LesionModel/CIP_LesionModel.py/CIP_LesionModelWidget.lesionType