Stohastic / Stohastic / Binary_destribution / Новый текстовый документ

#include <stdio.h>
#include <stdlib.h>
#define N 100
#define s 2
#define n 11

int rec(int e, int k)
{
    if(k>e)
        return 0;
    if(e==k)
        return 1;
    if(k==1)
        return e;
    return rec(e-1, k-1)+rec(e-1, k);
}


int X=200;
int myrand()
{
int m=0x7FF;
int a=214013,c=2531011;
    X = (a*X+c)%m;
return X;
}

int main()
{
    void srand(unsigned seed);
  int rand (void);

  int y=0,H[n];
 for(int j=0;j<n;j++)
  H[j]=0;

float x;
float Pk[n]={1./36.,3./36.,6./36.,10./36.,15./36.,21./36.,26./36.,30./36.,33./36.,35./36.,1};
double Pkk[n],p=0.5;
FILE *fp = fopen("test.txt","w");
//FILE *fp1= fopen("gyst.dat","w");

for(int i=0;i<N;i++)
{
    x=rand()/(RAND_MAX+1.);
    for(int j=1;j<n+1;j++)
        if((x>Pk[j-1])&&(x<Pk[j]))
        {
        H[j]=H[j]++;
      //printf("%f %d %d\n",x,y=j+2,H[j]);

        }
}
 for(int j=0;j<n;j++)
fprintf(fp,"%d\n", H[j]);

printf("%f ooooo\n",rec(36,4)*pow(p,4)*pow(1-p,32));


 fclose(fp);
// fclose(fp1);
    return 0;
}