///Doolittle分解
function Doolittle(myMatrix:Array, n:Number) {
 //將myMatrix分解成一個(gè)單位下三角矩陣和一個(gè)上三角矩陣的乘積
 // var temp:Number;
 for (var i = 1; i<=n; i++) {
  //計(jì)算U的第i行
  for (var j = i; j<=n; j++) {
   var sum:Number = 0;
   for (var r = 1; r<=i-1; r++) {
    sum = sum+myMatrix[i][r]*myMatrix[r][j];
   }
   myMatrix[i][j] = myMatrix[i][j]-sum;
  }
  //計(jì)算L的第i列
  for (var j = i+1; j<=n; j++) {
   var sum:Number = 0;
   for (var r = 1; r<=i-1; r++) {
    sum = sum+myMatrix[i][r]*myMatrix[r][i];
   }
   sum = myMatrix[j][i]-sum;
   myMatrix[j][i] = sum/myMatrix[i][i];
  }
 }
}
//Doolittle分解user52/comeonboy/upload/824171064.rar