久久久久久久999_99精品久久精品一区二区爱城_成人欧美一区二区三区在线播放_国产精品日本一区二区不卡视频_国产午夜视频_欧美精品在线观看免费

標題: n階(<20)行列式計算器 (有bug) [打印本頁]

作者: xiaos 時間: 2015-4-10 17:28
標題: n階(<20)行列式計算器 (有bug)

最近線代的計算量越來越大了網上的計算器只能算低階的行列式 .所以我自己就寫了一個最多可以計算20階的行列式的代碼計算結果總是會出問題,我發在這里希望有誰能給我建議 .基于C語言主要使用遞歸算法完成.

復制代碼

#include "stdafx.h"
float p[][20] = { 0 };
void T(float(*a)[20], float(*p2)[20], int y, int x, int n)//把a中的數據去掉第 y 行第x列后給p2
{
int k = 0, l = 0;
for (int i = 0; i<n; i++)
{
if (i != y)
{
for (int j = 0; j<n; j++)
{
if (j != x)
{
p2[k][l] = a[i][j];
if (++l == n - 1)
{
l = 0;
}
}
}
k++;
}
}
}
int cc(int i, int j)
{
if ((i + j) % 2 == 0)
{
return 1;
}
else
return -1;
}
float Deter(float(*a)[20], int n)//計算行列式的主要函數
{
float s = 0;
if (n == 2)
{
s = a[0][0] * a[1][1] - a[1][0] * a[0][1];
}
else
{
for (int i = 0; i<n; i++)
{
T(a, p, i, 0, n);
s += a[i][0] * Deter(p, n - 1)*cc(i + 1, 1);
}
}
return s;
}

復制代碼

歡迎光臨 (http://m.zg4o1577.cn/bbs/)