最適線形い - プログラム倉庫

th=new Array();
 
function tax(){
 
var s;
var h;
 
h=0.01;
 
for(s=1;s<101;s++){
th[s]=s*h;
}
 
var tr,tl;
var sx;
var we,maxw,trs,tls;
 
maxw=-9999;
trs=0;
tls=0;
 
sx=10;
 
maxw=-999;
for(sx=10;sx<50;sx++){
 
tl=0.01*sx;
tr=g(tl);
we=wel(tl,tr);
if(we>maxw)tls=tl;
if(we>maxw)trs=tr;
if(we>maxw)maxw=we;
}
 
$("#memo").html(tls);
 
 
 
}
 
function g(tlx){
 
var mode;
var e,tr1,tr2,tr3,b1,b2;
 
mode=0;
 
tr1=0.01;
tr2=0.4;
 
while(mode<100){
 
b1=bud(tlx,tr1);
b2=bud(tlx,tr2);
tr3=tr2-b2*(tr2-tr1)/(b2-b1);
tr1=tr2;
tr2=tr3;
e=b2*b2;
if(e<0.001)mode=1000;
mode=mode+1;
 
}
 
return tr2;
}
 
 
function wel(tl,tr){
 
var s;
var b1,w1,c1,l1,w2;
 
w2=0;
for(s=1;s<101;s++){
w1=(1-tl)*th[s];
l1=lab(w1,tr);
c1=w1*l1+tr;
w2=w2+u(c1,1-l1);
}
 
return w2;
}
 
function u(c1,x1){
var u1;
u1=Math.log(c1)+Math.log(x1);
return u1;
}
 
function bud(tl,tr){
 
var s;
var b1,w1,c1,l1,y1;
b1=0;
 
for(s=1;s<101;s++){
w1=(1-tl)*th[s];
l1=lab(w1,tr);
y1=th[s]*l1;
c1=w1*l1+tr;
b1=b1+y1-c1;
}
 
return b1;
}
 
function lab(w1,tr){
 
var l1;
l1=(w1-tr)/(2*w1);
if(l1<0)l1=0;
return l1;
}