うぇｍｘ - 経済統計練習帳

package p4;
 
 
 
public class pro {
 
   int[][] x=new int[7][7];
 
 int power[]=new int[7];
   int p[]=new int[7];
 
 
    public static void main(String[] args) {
        pro t=new pro();
    }
 
    pro(){
 
 
 
 
        x[1][1]=0;
        x[1][2]=12;
        x[1][3]=12;
        x[1][4]=14;
        x[1][5]=13;
        x[1][6]=17;
 
 
        x[2][1]=13;
        x[2][2]=0;
        x[2][3]=16;
        x[2][4]=10;
        x[2][5]=14;
        x[2][6]=15;
 
 
 
        x[3][1]=13;
        x[3][2]=9;
        x[3][3]=0;
        x[3][4]=8;
        x[3][5]=16;
        x[3][6]=14;
 
 
        x[4][1]=11;
        x[4][2]=15;
        x[4][3]=15;
        x[4][4]=0;
        x[4][5]=9;
        x[4][6]=10;
 
 
        x[5][1]=11;
        x[5][2]=10;
        x[5][3]=9;
        x[5][4]=15;
        x[5][5]=0;
        x[5][6]=10;
 
         x[6][1]=8;
        x[6][2]=10;
        x[6][3]=11;
        x[6][4]=15;
        x[6][5]=15;
        x[6][6]=0;
 
 
        power[1]=50;
        power[2]=50;
        power[3]=50;
        power[4]=50;
        power[5]=50;
        power[6]=50;
 
        int s;
 
 
        int trial;
 
       for(trial=1;trial<3;trial++){ 
        int k,tr;
        for(tr=1;tr<7;tr++){
 
        for(s=1;s<7;s++){
        p[s]=power[s];   
        }    
 
        double minx=999999;
        int ks;
        ks=0;
        for(k=1;k<101;k++){
        p[tr]=k;
        double g1=g();
        if(g1<minx)ks=k;
        if(g1<minx)minx=g1;
        }
        power[tr]=ks;
 
        }
 
 
 
 
 
 
 
       }
 
 
 
        for(s=1;s<7;s++){
       System.out.println(s+","+power[s]);  
       }
 
 
 
    }
    double  g(){
 
        int s,sx;
        double g1=0; 
 
        for(s=1;s<7;s++){
        for(sx=s+1;sx<7;sx++){
        g1=g1+f(s,sx);
        }
        }
 
    return g1;
    }
 
 
double f(int s1,int s2){
 
 
    double z1=p[s1]+p[s2];
    double z2=100*p[s1];
    double z3=z2/z1;
 
    double x1=x[s1][s2]+x[s2][s1];
    double x2=100*x[s1][s2];
    double x3=x2/x1;
 
    double f1=(z3-x3)*(z3-x3);
 
    return f1;
 
 
 
}
 
 
 
 
}