ｒｓｊｊｊ - 経済統計練習帳

package simple;
 
import org.apache.commons.math3.random.*;
import org.apache.commons.math3.random.GaussianRandomGenerator;
 
public class pro {
 
    int[] point=new int[500];
   int avg;
          double sig;
    double[] delta=new double[10001];
    int[] count=new int[101];      
 
    public static void main(String[] args) {
        pro te=new pro();
    }
 
 
    pro(){
 
        int s;
 
        for(s=1;s<101;s++){
            point[s]=25+(int)(50*Math.random());
        }
 
        int a1=0;
        for(s=1;s<101;s++){
        a1=a1+point[s];   
        }
        avg=a1/100;
 
        double s1=0;
         for(s=1;s<101;s++){
        s1=s1+(point[s]-avg)*(point[s]-avg);   
        }
        sig=s1/100;
        sig=Math.sqrt(sig);
        System.out.println(avg+","+sig);
 
        RandomGenerator rg = new JDKRandomGenerator();
 
        rg.setSeed(17399225432l); 
 
   GaussianRandomGenerator r = new GaussianRandomGenerator(rg);
 
 for(s=1;s<10001;s++){
   double x=r.nextNormalizedDouble();
   delta[s]=avg+sig*x;
    }
 
 for(s=0;s<101;s++){
 count[s]=0;    
 }
 
 for(s=1;s<10001;s++){
 int n=(int)delta[s];    
 if(n<0)n=0;
 if(n>100)n=100;
 count[n]=count[n]+1;
 }
 
 
 for(s=0;s<101;s++){
 System.out.println(s+","+count[s]);    
 }
 
 
 
 
 
 
 }
 
 
}