public override void Init_h()
{
h[1] = 0d;
h[2] = p * Geom.Get_dS(X);
h[3] = 0d;
for (int i = 0; i < mixture.powders.Count; i++)
{
h[i + 4] = q[1] * Geom.Square(X) * Math.Pow(p, mixture.powders[i].nu) / mixture.powders[i].Ik;
}
}
