ku - 経済統計練習帳

package n;
 
 
public class pro {
 
  double[][] p=new double[11][11];
   double[] n=new double[11];
 
    double[] lb=new double[11];
    double[] lbs=new double[11];
     double[] lc=new double[11];
    double[] w=new double[11];
     double[] r=new double[11];
 
  double alpha,beta;  
  int s,sx,number;
 
    public static void main(String[] args) {
       pro test=new pro();
    }
 
    pro(){
 
 
        alpha=0.7;
        beta=0.6;  
        number=10;
 
        for(s=1;s<number+1;s++){
        for(sx=1;sx<number+1;sx++){
        p[s][sx]=0.01;
        }
        }
 
 
 for(s=1;s<number+1;s++){        
for(sx=1;sx<number+1;sx++){
n[sx]=p[s][sx]+n[sx];
}
}        
 
 
 
 
  for(s=1;s<11;s++){
   lb[s]=0.5;
   }
 
  int mode;
 
  mode=0;
 
  while(mode<50){
 
    for(s=1;s<number+1;s++){
    w[s]=beta*Math.pow(n[s],beta-1)*Math.pow(lb[s],1-beta);
    r[s]=(1-beta)*Math.pow(n[s],beta)*Math.pow(lb[s],-beta);
    }
 
     for(s=1;s<number+1;s++){
     lc[s]=0;
     for(sx=1;sx<number+1;sx++){
     lc[s]=lc[s]+p[s][sx]*(1-alpha)*w[sx]/r[s];
     }}
 
 
     for(s=1;s<number+1;s++){
     lbs[s]=seekl(s);
     }
 
 
for(s=1;s<11;s++){
 lb[s]=lbs[s];
 }
 
mode=mode+1;
 
  }     
 
 
 for(s=1;s<number+1;s++){
 w[s]=beta*Math.pow(n[s],beta-1)*Math.pow(lb[s],1-beta);
 r[s]=(1-beta)*Math.pow(n[s],beta)*Math.pow(lb[s],-beta);
 }
 
 
 
 
    }
double seekl(int s){
 
        double l1,l2,l3;
        double x1,x2;
        int mode;
 
        l1=0.1;
        l2=0.9;
        l3=0;
        mode=0;
 
        while(mode<50){
        l3=(l1+l2)/2;
        x1=lx(s,l3)+l3;
        if(x1>1)l2=l3;
        if(x1<1)l1=l3;
 
 
        mode=mode+1;    
        }
        return l3;
    }
 
    double lx(int s ,double lp){
 
        double ws,rs;
 
    ws=beta*Math.pow(n[s],beta-1)*Math.pow(lp,1-beta);
    rs=(1-beta)*Math.pow(n[s],beta)*Math.pow(lp,-beta);
 
    int sx;
    double lx;    
    lx=0;
     for(sx=1;sx<11;sx++){
    lx=lx+p[s][sx]*(1-alpha)*w[sx]/r[s];
     }
 
     return lx;
 
 
 
    }
 
 
 
}