７ - 経済統計練習帳

package land2;
 
public class pro {
 
  double[] u=new double[11];
   int[] p=new int[11];
  int number,s;
   double theta;
 
 
    public static void main(String[] args) {
        pro test=new pro();
    }
 
   pro(){ 
 
 
 
        number=5;
 
        p[1]=10;
        p[2]=10;
        p[4]=10;
        p[5]=10;
        p[3]=10;
 
   int tr;
 
        for(tr=1;tr<10;tr++){
 
        ren sub=new ren();
        sub.p=p;
        sub.number=number;
        sub.makedata();
        u=sub.u;
 
 
        double maxu,minu;
        int s1,s2;
 
        s1=0;
        s2=0;
 
        maxu=0;
        for(s=1;s<number+1;s++){
          if(u[s]>maxu)s1=s;
        if(u[s]>maxu)maxu=u[s];
        }
 
 
        minu=9999;
        for(s=1;s<number+1;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<number+1;s++){
        System.out.println(p[s]);
        }
 
 
 
   }
}
 
 
 
package land2;
 
 
public class ren {
 
    double alpha,beta;
   double n; 
   double y,lb,lc;
   double rb,w;
 
   double[] c=new double[11];
   int[] p=new int[11];
    double[] r=new double[11];
    double[] nx=new double[11];
     double[] u=new double[11];
   int number,s;
   double theta;
 
    void makedata(){
 
        alpha=0.5;
        beta=0.7;
        theta=0.9;
 
 
    n=0;
   for(s=1;s<number+1;s++){    
   n=n+Math.pow(theta,s)*p[s];
   }
 
        lb=1;
        lc=1;
 
        y=Math.pow(n,beta)*Math.pow(lb,1-beta);
        rb=(1-beta)*Math.pow(n,beta)*Math.pow(lb,-beta);
        w=beta*Math.pow(n,beta-1)*Math.pow(lb,1-beta);
 
 
        for(s=1;s<number+1;s++){
        c[s]=alpha*w*Math.pow(theta,s);
        }
 
        for(s=1;s<number+1;s++){
        r[s]=p[s]*(1-alpha)*w*Math.pow(theta,s)/(lc);
        }
 
        for(s=1;s<number+1;s++){
        u[s]=Math.pow(c[s],alpha)*Math.pow((double)1/p[s],1-alpha);
        }
 
 
 
 
 
 
    }
}