かぜまち - 経済統計練習帳

package mon;
 
 
public class pro {
 
    double[] p=new double[11];
    int s;
    int[][] count=new int[10001][11];
     double[] y=new double[10001];
 
    public static void main(String[] args) {
        pro t=new pro();
    }
    pro(){
 
        for(s=1;s<11;s++){
    p[s]=0.1;        
        }
 
 
        int tr=1;
        for(tr=1;tr<10001;tr++){
        for(s=1;s<11;s++){
        count[tr][s]=0;        
        }
 
        for(s=1;s<101;s++){
        int n=randx();    
        count[tr][n]=count[tr][n]+1;
        }
 
        }
 
      for(tr=1;tr<10001;tr++){
        double z=0;
        for(s=1;s<11;s++){
        double z1=(0.01*count[tr][s]-p[s])*(0.01*count[tr][s]-p[s]);  
        z=z+z1/p[s];
        }
        y[tr]=z;
 
      }
 
       for(tr=1;tr<101;tr++){
       System.out.println(y[tr]);
       }
 
 
 
 
 
    }
    int randx(){
 
        double[] px=new double[11];
        int s;
        px[1]=p[1];
        for(s=2;s<11;s++){
        px[s]=px[s-1]+p[s];    
        }
 
        double r1=px[10]*Math.random();
 
        int n=1;
        for(s=1;s<10;s++){
        if(r1>px[s])n=s+1;
        }
 
        return n;
    }
 
 
}