Random-sample-analysis.zip_random

#include "math.h" #include "stdio.h" void rhis(x,n,x0,h,m,l,dt,g,q) int n,m,l,g[],q[]; double x[],x0,h,dt[]; { int i,j,k,z; double s; char a[50]; dt[0]=0.0; for (i=0; i<=n-1; i++) dt[0]=dt[0]+x[i]/n; dt[1]=0.0; for (i=0; i<=n-1; i++) dt[1]=dt[1]+(x[i]-dt[0])*(x[i]-dt[0]); dt[1]=dt[1]/n; dt[2]=sqrt(dt[1]); for (i=0; i<=m-1; i++) { q[i]=0; s=x0+(i+0.5)*h-dt[0]; s=exp(-s*s/(2.0*dt[1])); g[i]=n*s*h/(dt[2]*2.5066); } s=x0+m*h; for (i=0; i<=n-1; i++) if ((x[i]-x0)>=0.0) if ((s-x[i])>=0.0) { j=(x[i]-x0)/h; q[j]=q[j]+1; } if (l==0) return; printf("\n"); printf("n=%d\n",n); printf("\n"); printf("x0=%e h=%e m=%d\n",x0,h,m); printf("\n"); printf("xa=%e s=%e t=%e\n",dt[0],dt[1],dt[2]); printf("\n"); k=1; z=0; for (i=0; i<=m-1; i++) if (q[i]>z) z=q[i]; while (z>50) { z=z/2; k=2*k;} for (i=0; i<=m-1; i++) { s=x0+(i+0.5)*h; for (j=0; j<=49; j++) a[j]=' '; j=q[i]/k; for (z=0; z<=j-1; z++) a[z]='X'; j=g[i]/k; if ((j>0)&&(j<=50)) a[j]='*'; printf("%e %7d ",s,q[i]); for (j=0; j<=49; j++) printf("%c",a[j]); printf("\n"); } printf("\n"); printf("k=%d\n",k); printf("\n"); return; }