Оценка псевдослучайной последовательности

#include <iostream>
#include <cmath>
using namespace std;
const int n=150, r=12;
double mass(double x, double y, double z, int a)
{
    double f=(abs(sin(x+y)+sin(z+a)))/2; 
    double s=(f*10)/1;
    double masss=(int)s%10/2.0;
    for (int i=1;i<10;i++) masss+=s/(((int)pow(10.0,10.0-i)%(i*10)))*pow(2.0,-i-1);
    z=y; y=x; x=masss;
    return masss;
}
void inter(double e[r], double&z)
{
    for (int k=0;k<r;k++) if ((z>=(double)k/r)&&(z<(double)(k+1)/r)) e[k]++;
}
void main()
{
    double x0,x1,x2,z,x[n],me,de;
    double e[r];
    x0=0.841;
    x1=0.126;
    x2=0.385;
    me=0;de=0;
    for (int i=0;i<=1;i++)
    {
        for (int j=0;j<n;j++)
        {
            x[j]=mass(x0,x1,x2,i);
            inter (e,x[j]);
        }
        for (int k=0;k<r;k++) me+=(double)(2*k+1)/(2*r)*e[k]/n;
        for (int k=0;k<r;k++) de+=(double)((2*k+1)/(2*r)-me)*((2*k+1)/(2*r)-me)*e[k]/n;
        cout<<"for a="<<i<<endl<<"e[k]: ";
        for (int k=0;k<r;k++) cout<<e[k]<<" ";
        cout<<endl<<"me="<<me<<endl<<"de="<<de<<endl;
    }
    system ("pause");
}