/*CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C                                                                             C
C  COMPUTE MERTENS FUNCTION (using Deleglise and Rivat's algorithm)           C
C  10/16/15 (DKC)							      C
C                                                                             C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC*/
#include <math.h>
double nuloc(unsigned int x1, unsigned int x2, unsigned int u, char *mob, int *M);
double nuhic(unsigned int x1, unsigned int x2, unsigned int u, char *mob, int *M);
int nuriv(unsigned int x1, unsigned int x2, unsigned int u, char *mobb, int *M) {
double f1,f2,f3;
f1=nuloc(x1,x2,u,mobb,M);
f2=nuhic(x1,x2,u,mobb,M);
f3=(double)M[u-1]-f1-f2;
return((int)f3);
}