Java编写的姐希尔伯特矩阵

里面顺便带着楚列斯基分解的程序，两个程序是独立的

import java.util.*;

public class CholeskyYiwei {
   private double b[];
   private int ijkn=3m;
   public static void main（String ss[]）{
	   CholeskyYiwei aa=new CholeskyYiwei（）;
	   aa.input（）;
	   aa.Calculate（）;
	   aa.outPut（）;
   }
   void input（）{
	   System.out.println（“输入数组维数：“）;
	   Scanner in=new Scanner（System.in）;
	   n=in.nextInt（）;
	   b=new double[n*（n+1）/2];
	   for（i=1;i<=n;i++）{
		   for（j=1;j<=i;j++）{
			   System.out.println（“a[“+i+“][“+j+“]“）;
			   b[i*（i-1）/2-1+j]=in.nextDouble（）;
		   }
	   }
   }
   void Calculate（）{
	   for（k=1;k<=n;k++）{
		   for（i=1;i<=k-1;i++）{
			  b[k*（k+1）/2-1]= b[k*（k+1）/2-1]-b[i*（i+1）/2-1]*b[k*（k-1）/2-1+i]*b[k*（k-1）/2-1+i];
		   }
		   for（m=k+1;m<=n;m++）{
			   for（j=1;j<=k-1;j++）{
				   b[m*（m-1）/2-1+k]= b[m*（m-1）/2-1+k]- b[m*（m-1）/2-1+j]* b[k*（k-1）/2-1+j]*b[k*（k+1）/2-1];
			   }
			   b[m*（m-1）/2-1+k]= b[m*（m-1）/2-1+k]/b[k*（k+1）/2-1];
		   }
	   }
   }
   void outPut（）{
	   for（i=1;i<=n;i++）{
		   for（j=1;j<=i;j++）{
			   System.out.print（b[i*（i-1）/2-1+j]+“\t“）;	 
		   }
		   System.out.println（）;
		  }
   }
   
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       3114  2010-04-16 06:31  Hilbert.java

     文件       1143  2010-04-14 14:28  CholeskyYiwei.java

----------- ---------  ---------- -----  ----

                 4257                    2