计算稀疏矩阵最大特征值及对应特征向量的算法

package MATRIXTR;

public class LeadingEigenvalue {


public static void main(String[] args) {
// int matrix[][]={{0,1,0},{1,1,0},{1,0,0}};
int matrix[][]={{2,1},{1,2}};
// int matrix[][] = { { 0, 1, 0, 1, 1, 0, 0 }, { 1, 0, 1, 1, 0, 0, 0 },
// { 0, 1, 0, 0, 0, 0, 0 }, { 1, 1, 0, 0, 0, 1, 0 },
// { 1, 0, 0, 0, 0, 1, 0 }, { 0, 0, 0, 1, 1, 0, 1 },
// { 0, 0, 0, 0, 0, 1, 0 } };
// float[][]matrix={{1 ,  2 , 1  , 4 , 6},
//        {1/2 , 1 ,1/2 , 2,  3},
//        {1/2 , 2,  1,   1,  1},
//        {1/4, 1/2 ,1   ,1 , 1},
//        {1/6, 1/3, 1,   1,  1}} ;
//double [][]matrix={{1, 1/3,1/3,1/5,1/9},{3, 1, 1, 1/2,1/3},{3, 1, 1, 1/2,1/3},{5, 2, 2, 1, 1/2},{9, 3, 3, 2, 1}};
double LeadingVector[]=new double[matrix.length];
for(int i=0;i<matrix.length;i++)
{
LeadingVector[i]=1;
}

double vector[]=new double[matrix.length];
for(int x=0;x<100000;x++)
{
for(int i=0;i<matrix.length;i++)
{
double sum=0;