ひでよし - 経済統計練習帳

package pocaso;
 
 
public class pro {
 
    int number;
    int[] x=new int[101];
    int[] y=new int[101];
    double[] avg=new double[101];
 
    double[] mux=new double[11];
     double[] sigx=new double[11];
 
     double[] muy=new double[11];
     double[] sigy=new double[11];
 
     double[] p=new double[11];
 
     double[] px=new double[11];
 
     double[][] rate=new double[101][11];
 
    public static void main(String[] args) {
        pro test=new pro();
    }
    pro(){
 
        int s;
        for(s=1;s<11;s++){
        mux[s]=25+5*s;
        muy[s]=25+5*s;
        sigx[s]=10;
        sigy[s]=10;
     p[s]=0.1;
       }
 
 
 
        int n;
        for(n=1;n<101;n++){
        double x1=40;
        double y1=52;
        for(s=1;s<11;s++){
 
    double p1=f(x1,mux[s],sigx[s]);    
    double p2=f(y1,muy[s],sigy[s]);    
    px[s]=p1*p2*p[s];    
 
        }
 
        double p1=0;
        for(s=1;s<11;s++){
        p1=p1+px[s];    
        }
 
 
        for(s=1;s<11;s++){
        rate[n][s]=px[s]/p1;    
        }
 
        }
 
 
 
 
 
    }
 
     double f(double x1,double mux,double sigx){
double pi,f1,f2;
 
        pi=Math.PI;
 
        f2=Math.sqrt(2*pi*sigx*sigx);
        f1=-(x1-mux)*(x1-mux)/(2*sigx*sigx);
 
        return Math.exp(f1)/f2;
}
 
 
 
 
}