高斯消元模板啦啦啦..

来源:互联网 发布:三菱plc编程入门pdf 编辑:程序博客网 时间:2024/05/24 07:22

放一下高斯消元模板啦啦啦..

    for ( i = 1; i <= n; i ++ ){        for ( j = i; j <= n; j ++ ){            if ( _abs ( fc[j][i] ) > eps ){                for ( k = i; k <= n+1; k ++ ){                    swap ( fc[j][k], fc[i][k] );                }                break;            }        }        for ( j = 1; j <= n; j ++ ){            if ( i != j && _abs (fc[j][i]) > eps ){                double tmp = fc[j][i]/fc[i][i];                for ( k = i; k <= n+1; k ++ ) fc[j][k] -= tmp*fc[i][k];            }        }    }    for ( i = 1; i <= n; i ++ ){        for ( j = 1; j <= n; j ++ ){            if ( _abs (fc[j][i]) > eps ) ans[i] = fc[j][n+1]/fc[j][i];        }    }
0 0
原创粉丝点击