空間４ - 経済統計練習帳

package sj;
 
public class pro {
 
  int[] p=new int[101];
     double[] u=new double[101];
  int s;
     double maxu,minu;
 
 
    public static void main(String[] args) {
        pro test=new pro();
    }
 
   pro(){ 
 
       for(s=1;s<101;s++){
       p[s]=100;
       }
 
       int t;
 
       for(t=1;t<10000;t++){
 
       city sub=new city();
       sub.p=p;
       sub.makedata();
       u=sub.u;
 
       int s1,s2;
       s1=0;
       maxu=0;
       for(s=1;s<101;s++){
       if(u[s]>maxu)s1=s;
       if(u[s]>maxu)maxu=u[s];
       }
 
       s2=0;
       minu=9999;
       for(s=1;s<101;s++){
       if(u[s]<minu)s2=s;
       if(u[s]<minu)minu=u[s];    
       }
 
       p[s1]=p[s1]+1;
       p[s2]=p[s2]-1;
 
 
       }
 
       for(s=1;s<101;s++){
       System.out.println(p[s]);
       }
 
   }
}
 
 
 
 
package sj;
 
 
public class city{
 
    double nx,y;
  double alpha,beta;
  double rs,w;  
  double[] c=new double[101];
  double[] a=new double[101];
   double[] r=new double[101];
   double[] q=new double[101];
   double z;
  int[] p=new int[101];
     double[] u=new double[101];
  int s;
 
 
 
 
   void makedata(){ 
 
 
 
       for(s=1;s<101;s++){
       q[s]=1-0.005*s;
       nx=nx+p[s]*q[s];
       }
 
 
       alpha=0.9;
       beta=0.8;
 
 
 
       int t;
 
       z=0;
 
     for(t=1;t<10;t++){  
 
 
 
 
       y=Math.pow(nx,alpha)*Math.pow(1,1-alpha);
       w=alpha*Math.pow(nx,alpha-1)*Math.pow(1,1-alpha);
       rs=(1-alpha)*Math.pow(nx,alpha)*Math.pow(1,-alpha);
 
     for(s=1;s<101;s++){  
     a[s]=q[s]*w+z;
     c[s]=beta*a[s];
     r[s]=(1-beta)*a[s]*p[s];
     }
 
   z=rs;
    for(s=1;s<101;s++){
    z=z+r[s];
    }
    z=z/10000; 
 
 
   }
 
 
 
 int s;
 double p1;    
 
   for(s=1;s<101;s++){ 
   p1=(double)1/p[s];
   u[s]=Math.pow(c[s],beta)*Math.pow(p1,1-beta);
   }
 
 
 
 
 
   }
}