X := mu *(p[2] + p[55] + p[56] + p[57] + p[58] + p[59] + p[252] + p[323] +
     p[350] + p[389] + p[443] + p[453] + p[484] + p[508] + p[534] + p[542]
     + p[572] + p[606] + p[615] + p[627] + p[638] + p[650] + p[661] +
     p[673] + p[687] + p[696] + p[708] + p[720] + p[731] + p[742] + p[748]
     + p[764]);