lp100 - プログラム倉庫

class pro{
public static void main(String args[]){
int s;
double th[]=new double[11];
double c[]=new double[11];
double y[]=new double[11];
double w[]=new double[11];
double u[]=new double[11];
double lam,la1,la2;
double mu[]=new double[11];
double b1,c1,y1,u1;
int t1;
double z,z1,z2,z3,p1,p2,p3;
for (s=1;s<11;s++){
th[s]=1+0.1*s;
}
b1=0;
la1=0.1;
la2=3;
t1=0;
lam=0;
while(t1<100){
lam=(la1+la2)/2;;  
b1=0;
for (s=1;s<11;s++){
c1=1/lam;
y1=th[s]-1/lam;
b1=b1+y1-c1;
}
if (b1>0)la2=lam;
if (b1<0)la1=lam;
if (b1*b1<0.0001)t1=1000;
t1=t1+1;
}
System.out.println(lam);
System.out.println(b1);

for (s=1;s<11;s++){
c1=1/lam;
y1=th[s]-1/lam;
u[s]=Math.log(c1)+Math.log(1-y1/th[s]);
}

for (s=1;s<10;s++){
c1=1/lam;
y1=th[s-1]-1/lam;
w[s]=Math.log(c1)+Math.log(1-y1/th[s]);
}

s=5;
mu[s]=0.01;
c1=(1-mu[s])/lam;
System.out.println(c1);

z1=0;
z2=th[s]-0.01;

p1=1/(z1-th[s])+lam+mu[s]/(th[s+1]-z1);
System.out.println(p1);

p2=1/(z2-th[s])+lam+mu[s]/(th[s+1]-z2);
System.out.println(p2);

t1=0;
z3=0;
while (t1<100){
z3=(z1+z2)/2;
p3=1/(z3-th[s])+lam+mu[s]/(th[s+1]-z3);
if (p3>0)z1=z3;
if (p3<0)z2=z3;
if (p3*p3<0.0001)t1=1000;
t1=t1+1;
}
y1=z3;
u1=Math.log(c1)+Math.log(1-y1/th[s+1]);
System.out.println(u1);





}
}