最近在做一个关于施密特正交的问题:
源程序编写在附件中,

有时程序会运行成功,但是有时会有如下错误,
"Floatint point error:Domain"
请问一下大侠,这是什么原因,应该怎么样修改程序,使运行不再出这种错误!

希望大家帮忙啊

虽然我不明白施密特正交，但我可以建议逐个调试一下，多插入printf()去看看。

不晓得这个原因啊,期待达人解答.
贴一下我的程序,写的很粗糙,扩展性做的不好,也没有考虑溢出.见谅哈~
你修改下,加上clrscr()和getch(),还有注释就OK了.

#include<stdio.h>
#include<stdlib.h>

// 向量空间维数
#define VR 3
// 向量长度
#define VN 3

// 求内积
float neiji(float x1[], float x2[], int n)
{
float neiji = 0.0;
for(int i=0; i<n; i++)
{
neiji += (x1[i] * x2[i]);
}
return neiji;
}

// 向量相减
void xiangjian(float x1[], float x2[], int n)
{
for(int i=0; i<n; i++)
{
x1[i] -= x2[i];
}

}

// 数乘向量
void xiangcheng(float x1[], float x2[], float m, int n)
{
for(int i=0; i<n; i++)
{
x1[i] = m * x2[i];
}
}

// 施密特正交化
void zhengjiaohua(float (*b)[VN], float (*a)[VN])
{
int i,j;
float temp;
float tempNeiji[VR];
float tempXiangcheng[VN];

for(i=0; i<VR; i++)
{
for(j=0; j<VN; j++)
{
*(b[i]+j) = *(a[i]+j);
}
}

for(i=1; i<VR; i++)
{
tempNeiji[i-1] = neiji(b[i-1], b[i-1], VN);
for(j=0; j<i; j++)
{
temp = neiji(b[j], a[i], VN);
xiangcheng(tempXiangcheng, b[j], temp/tempNeiji[j], VN);
xiangjian(b[i], tempXiangcheng, VN);
}
}
}

int main()
{
int i,j;
float a[VR][VN];
float b[VR][VN];

for(i=0; i<VR; i++)
{
printf("please input a[%d] which has %d numbers: \n", i, VN);
for(j=0; j<VN; j++)
{
scanf("%f", &a[i][j]);
}
}

printf("\nthe numbers of inputed arr: \n");

for(i=0; i<VR; i++)
{
for(j=0; j<VN; j++)
{
printf("%.2f ", a[i][j]);
}
printf("\n");
}

zhengjiaohua(b, a);

printf("\nthe numbers of formated arr: \n");

for(i=0; i<VR; i++)
{
printf("\n");
for(j=0; j<VN; j++)
{
printf("%.2f ",b[i][j]);
}
}
}