高斯消元法Java、C、Fortran实现代码

Java实现：

package zz;

public class Gauss {

	// 高斯消元法
	public float[] guass(float[][] arr) {

		float[] result = solution(arr);
		for (int i = 0; i < arr.length; i++) {
			result[i] = arr[i][arr[0].length - 1];
		}
		
		System.out.println("高斯消元之后：");
		printArr(arr);

		for (float f : result) {
			System.out.print(f + " ");
		}
		return result;
	}

	/*
	 * 结果在行列式中 所以 列式比行数多1
	 */
	public float[] solution(float[][] arr) {
		int row = arr.length; // 行数
		int col = arr[0].length; // 列数
//		System.out.println("row = " + row + " col = " + col); 
		float[] result = new float[row]; // 有多少个方程就有多少个未知数， 所以结果的个数为 行数

		for (int i = 0; i < row; i++) {
			int maxRow = SelectIndex(arr, i);

			// 寻找最大行、交换
			if (maxRow != i) {
				swapRow(arr, i, maxRow);
			}

			// 归一化 当前行 每一个元素都除以当前 arr[i][i] 上的元素，归一化
			float divNum = arr[i][i];
			for (int k = i; k < col; k++) { // 当前行 的每一列元素都除以这个数
				arr[i][k] /= divNum;
			}

			// 用这个1， 消去下面所有行的数
			arr = elimination(arr, i);

		}

//		System.out.println("回代啦");
		// 回代 从最后一行开始
		backProgram(arr);

		return result;
	}

	/*
	 * 从最后一行回代，得到结果
	 */
	public float[][] backProgram(float[][] arr) {

		int row = arr.length; // 行
		int col = arr[0].length; // 列
		for (int i = row - 1; i >= 0; i--) {
			for (int j = i - 1; j >= 0; j--) {
//				System.out.println("arr[" + j + "][" + (col - 1) + "] = " + arr[j][col - 1]);
//				System.out.println("arr[" + j + "][" + i + "] = " + arr[j][i]);
//				System.out.println("arr[" + i + "][" + i + "] = " + arr[i][i]);
//				System.out.println("arr[" + j + "][" + i + "] * arr[" + i + "][" + (col - 1) + "] = "
//						+ arr[j][i] * arr[i][col - 1]);
//				printArr(arr);
				arr[j][col - 1] -= arr[j][i] * arr[i][col - 1];
				arr[j][i] = 0;

//				printArr(arr);
			}

		}

		return arr;
	}

	/*
	 * 用当前行消去下面的所有元素
	 */
	public float[][] elimination(float[][] arr, int