ytju - プログラム倉庫

package simple;
 
 
public class pro {
 
    double a,beta;
double k[]=new double[101]; 
double h;    
double cx[]=new double[101];
int n;
 
 
    public static void main(String[] args) {
        pro test=new pro();
    }
 
 pro(){
 
 fasttime();
 
 fast();
 
 for(n=10;n<91;n++){
 System.out.println(cx[n]);
 }
 
}
void fasttime(){
 
double ks;
 
a=0.33; 
beta=0.95;
 
 
 
ks=Math.pow((1 / beta - 1) / a , 1 / (a - 1)); 
h=2*ks/100; 
for (n=1;n<101;n++){ 
k[n]=n*h;
cx[n]=Math.pow(k[n],a);   
} 
 
 
 
 
}
 
 
 
 
 
 
void fast(){
 
double cxs[]=new double[101]; 
    int t1;
 
double c1,l1,k1,uc,wedge,r1,w1;
double n1,s1,ep;
t1=0; 
while(t1<100){ 
 
for (n=10;n<91;n++){ 
c1=con(n);
k1=k[n]+Math.pow(k[n],a)-c1; 
r1=a*Math.pow(k1,a-1); 
uc=(beta*(1+r1))/c1;   
cxs[n]=1/uc; 
} 
 
 
ep=0; 
for (n=10;n<91;n++){ 
ep=ep+Math.pow(cx[n]-cxs[n],2); 
}  
 
for (n=10;n<91;n++){ 
cx[n]=cxs[n];  
} 
 
if (ep<0.0001)t1=1000;   
t1=t1+1; 
} 
 
 
 
}
 
double con(double nx){
double c1;
int n1,n2,s1,s2;
n1=(int)nx;
n2=n1+1;
if(n1>90)n1=90;
if(n2>90)n2=90;
c1=cx[n1]+(nx-n1)*(cx[n2]-cx[n1]);    
return c1;    
}
 
 
 
}