математическое ожидание и дисперсия

 #define MODMULT(a, b, c, m, s) q = s/a; s = b*(s-a*q)-c*q; if (s<0) s+=m;
using namespace std;
 
 
 
static long s1 = 9;
static long s2 = 5;
 long z;
/* MODMULT (a, b, c, m, s) рассчитывает s*b mod m при условии, что m = a*b+c и 0<=c<m
#define MODMULT(a, b, c, m, s) q = s/a; s = b*(s-a*q)-c*q; if (s<0) s+=m;
/* Возвращает действительное псевдослучайное значение в диапазоне (0,1) */
 
 
double CombLCG (void)
{
  long q;
 
  MODMULT(53668, 40014, 12211, 2147483563L, s1)
  MODMULT(52774, 40692, 3791, 2147483399L, s2)
  z = s1 - s2;
  if(z<1)
    z += 2147483562;
  return z*4.656613e-10;
    
}
void InitLCG(long InitS1, long InitS2)
{  s1 = InitS1; s2 = InitS2;
 
} 
 
    
    
 
 
int _tmain(int argc, _TCHAR* argv[]){setlocale(LC_CTYPE, "rus");
{   InitLCG(s1,s2);
    CombLCG();
       cout << "s1 = " << s1 << endl;
        cout << "s2 = " << s2 << endl;
        cout << "z = " << z << endl;
 
        
 
 cout << "Введите количество элементов в массиве  >1:";
    int N;cin >> N;
 
    long*  z=new long[N];
    long M=0;
    for (int i = 0; i < N; i++) {
       // cout << "Введите X[" << i << "]:";
        cin >> z[i];
        M += z[i];
    }
    cout << "Мат.ожидание = " << M / N;
    double D=0;
    for (int i = 0; i < N; i++) {
        D += (z[i] - M)*(z[i] - M);
    }
    cout << "\nДисперсия = " << D / (N - 1)<<endl;
    cout << "Среднеквадр. отклонение = " << sqrt(D / (N - 1))<<endl;
    delete[] z;
 
 
 
    getch();
    return 0;
}
}