Question

如何用fortran编程求解四元线性方程组

方程如下：
http://hi.baidu.com/%BD%E8%B9%FDcsu/album/%E9%BB%98%E8%AE%A4%E7%9B%B8%E5%86%8C
从文件中读入m1、m2和θ的值，求a1、a2、N1 和N2的值。其中g取9.8，输入θ时以角度为单位。
要求：
（1）分别用两种方法（例如高斯消去法、矩阵求逆法、三角分解法、追赶法等），定义求解线性方程组Ax=b的子程序，要求该子程序能求解任意线性方程组。
（2）在主程序中分别调用上面定义的两个子程序，并对求解结果进行对比分析。

Answer 1

在MKL或是IMSL里有现成的算法，您直接调用就行了。求逆是用.i.（在IMSL中）。