ぽうぇｒ - プログラム倉庫

var ax=new Array();
var ap=new Array();
var bx=new Array();
var b=new Array();
var z;
var lam=new Array();
var v=new Array();
y2=new Array();
x2=new Array();
y1=new Array();
x1=new Array();
 
function matx(){
 
var avx,avy,sigx,sigy;
 
avx=0;
for(s=1;s<artistnumber+1;s++){
avx=avx+x[s];
}
avx=avx/artistnumber;
 
avy=0;
 
for(s=1;s<artistnumber+1;s++){
avy=avy+y[s];
}
avy=avy/artistnumber;
 
sigx=0;
for(s=1;s<artistnumber+1;s++){
sigx=sigx+(x[s]-avx)*(x[s]-avx);
}
 
sigx=Math.sqrt(sigx/artistnumber);
 
sigy=0;
for(s=1;s<artistnumber+1;s++){
sigy=sigy+(y[s]-avy)*(y[s]-avy);
}
sigy=Math.sqrt(sigy/artistnumber);
 
 
for(s=1;s<artistnumber+1;s++){
x1[s]=(x[s]-avx)/sigx;
y1[s]=(y[s]-avy)/sigy;
}
 
 
for(s=1;s<3;s++){
ap[s]=new Array();
}
 
z=0;
for(s=1;s<artistnumber+1;s++){
z=z+x1[s]*x1[s];
}
ap[1][1]=z;
 
z=0;
for(s=1;s<artistnumber+1;s++){
z=z+y1[s]*y1[s];
}
ap[2][2]=z;
 
z=0;
for(s=1;s<artistnumber+1;s++){
z=z+x1[s]*y1[s];
}
ap[1][2]=z;
ap[2][1]=z;
 
}
 
function power(){
 
var s;
 
 
for(s=1;s<3;s++){
v[s]=new Array();
}
 
for(s=1;s<3;s++){
ax[s]=new Array();
}
 
ax=ap;
 
seekbs();
 
lam[1]=z;
 
for(s=1;s<3;s++){
v[1][s]=b[s];
}
 
 
for(s=1;s<3;s++){
for(sx=1;sx<3;sx++){
ax[s][sx]=ax[s][sx]-z*b[s]*b[sx];
}
}
 
seekbs();
 
lam[2]=z;
 
for(s=1;s<4;s++){
v[2][s]=b[s];
}
 
 
for(s=1;s<number+1;s++){
x2[s]=x1[s]*v[1][1]+y1[s]*v[1][2];
y2[s]=x1[s]*v[2][1]+y1[s]*v[2][2];
}
 
 
 
}
 
function seekbs(){
 
var mode;
 
b[1]=1;
b[2]=0;
 
mode=0;
 
while(mode<10){
 
seekb();
 
var e;
      e=0;
for(s=1;s<3;s++){
    e=e+ (bx[s]-b[s])*(bx[s]-b[s]);
     }
if(e<0.01)mode=100;
 
      for(s=1;s<3;s++){
      b[s]=bx[s];
      }
 
      mode=mode+1;
 
}
 
}
 
 
function seekb(){
 
var s,sx;
var b1;
 
    for(s=1;s<3;s++){
    b1=0;
    for(sx=1;sx<3;sx++){
    b1=b1+ax[s][sx]*b[sx];
    }
    bx[s]=b1;
    }
 
   z=0;
     for(s=1;s<3;s++){
     z=z+bx[s]*bx[s];
     }
 
z=Math.sqrt(z);
 
 
 for(s=1;s<3;s++){
     bx[s]=bx[s]/z;
     }
 
}