ｐｐｐｐ１００ - プログラム倉庫

class pro{
public static void main(String args[]){
int s,t1;
double th[]=new double[11];
double c[]=new double[11];
double y[]=new double[11];
double mu[]=new double[11];
double lam,la1,la2;
double b1,u1;
for (s=1;s<11;s++){
mu[s]=0.1;
}

lam=1;
la1=0.1;
la2=3;
t1=0;
while(t1<100){
lam=(la1+la2)/2;;
for (s=1;s<11;s++){
th[s]=1+0.1*s;
c[s]=0.5*th[s];
y[s]=0.5*th[s];
}
b1=0;
for (s=1;s<11;s++){
c[s]=1/lam;
y[s]=th[s]-1/lam;
b1=b1+y[s]-c[s];
}
System.out.println(b1);
if (b1>0)la2=lam;
if (b1<0)la1=lam;
if (b1*b1<0.0001)t1=1000;
t1=t1+1;
}
for (s=2;s<11;s++){
u1=Math.log(c[s])+Math.log(1-y[s]/th[s])+lam*(y[s]-c[s]);
u1=u1+mu[s]*(Math.log(c[s-1])+Math.log(1-y[s-1]/th[s])-Math.log(c[s])-Math.log(1-y[s]/th[s]));
System.out.println(u1);
}


}
}