ｊｈｏ - 経済統計練習帳

package docomo;
 
 
public class city {
 
   double alpha,beta,theta;
   int[] p=new int[101];
   int[] q=new int[101];
   double[] n=new double[11];
    double[] w=new double[11];
      double[] r=new double[11];
     double[] lp=new double[101];
     double[] cp=new double[101];
  double[] u=new double[101];
     int s;
 
 
 
   void makedata(){ 
 
 
 
    alpha=0.5;
    beta=0.7;
    theta=0.9;
 
 
    int sx,s1,s2;
 
    for(s=1;s<101;s++){
    sx=p[s]-q[s];
    if(sx<0)sx=-sx;    
    s1=q[s];
    n[s1]=n[s1]+Math.pow(theta,sx);    
    }
 
 
    for(s=4;s<7;s++){
    w[s]=beta*Math.pow(n[s],beta-1)*Math.pow(1,1-beta);
    r[s]=(1-beta)*Math.pow(n[s],beta)*Math.pow(1,-beta);
    }
 
 
    double lx;
    int t;
 
 
    for(s=1;s<4;s++){
    lx=0;
    for(sx=1;sx<101;sx++){
    s1=q[sx];
    t=p[sx]-q[sx];
    if(t<0)t=-t;
    if(s==p[sx])lx=lx+Math.pow(theta,t)*(1-alpha)*w[s1];
    }
    r[s]=lx;
    }
 
 
    for(s=7;s<11;s++){
    lx=0;
    for(sx=1;sx<101;sx++){
    s1=q[sx];
    t=p[sx]-q[sx];
    if(t<0)t=-t;
    if(s==p[sx])lx=lx+Math.pow(theta,t)*(1-alpha)*w[s1];
    }
    r[s]=lx;
    }
 
 
 
 
     for(s=1;s<101;s++){
    s1=p[s];
    s2=q[s];     
    t=s1-s2;
    if(t<0)t=-t;
    cp[s]=alpha*Math.pow(theta,t)*w[s2];     
   lp[s]=(1-alpha)*Math.pow(theta,t)*w[s2]/r[s1];   
     }
 
 
      for(s=1;s<101;s++){
      u[s]=Math.pow(cp[s],alpha)*Math.pow(lp[s],1-alpha);
      }
 
 
 
 
 
 
 
 
 
 
 
 
 
 
   }
   int trial_p(){
 
       int mode,n;
       mode=0;
       n=0;        
       while(mode<50){
       n=trial();
       if(n<4)mode=100;    
       if(n>6)mode=100;    
       }
 
   return n;
   }
 
 
   int trial_q(){
 
       int mode,n;
       mode=0;
       n=0;        
       while(mode<50){
       n=trial();
       if(n==4)mode=100;    
       if(n==5)mode=100;    
       if(n==6)mode=100;  
       }
 
   return n;
   }
 
 
 
int trial(){
   int n;
   n=(int)(10*Math.random())+1;
    return n;
 
}
}