52 C = 0.1 + -((
theFragZ-50.)/20.)*0.02;
56 C = 0.08 + -((
theFragZ-70.)/20.)*0.02;
90 result = 3.0/(P*(P-1.0)*(P-2.0)*(P-3.0)) *
91 (H*(H-1.0)*(H-2.0)*(H-3.0)*r*r*(r-1.0)*(r-1.0) +
92 2.0*H*(H-1.0)*(H-2.0)*(Pa*r*(1.0-
r)*(1.0-r)+Na*r*r*(1.0-
r)) +
93 H*(H-1.0)*(Pa*(Pa-1.0)*(1.0-
r)*(1.0-r)+4.0*Na*Pa*r*(1.0-
r)+Na*(Na-1.0)*r*
r) +
94 2*H*(Pa*Na*(Na-1.0)*r+Pa*(Pa-1.0)*Na*(1.0-
r)) +
95 Pa*(Pa-1.0)*Na*(Na-1.0));
97 result /= 6.0*r*r*(1. -
r) *(1. - r);