求高手 给看一下我的平方根算法解正定矩阵代码 请给改一下

#include<iostream.h>
#include<math.h>
#define N 10
void main()
{
double a[N][N],b[N];
int i,j,k,h,n;
cout<<"please input n :"<<endl;
cin>>n;
cout<<"please input a["<<n<<"]["<<n<<"]"<<endl;
for(i=1;i<=n;i++)
{
for(j=1;j<=n;j++)
{
cin>>a[i][j];
}
}
cout<<"please input b["<<n<<"]"<<endl;
for(i=1;i<=n;i++)
{
cin>>b[i];
}
for(j=1;j<=n;j++)
{
double temp=0;
for(k=1;k<j;k++)
{
temp=temp+a[j][k]*a[j][k];
}
a[j][j]=sqrt(a[j][j]-temp);
for(i=j+1;i<=n;i++)
{
double remp=0;
for(h=1;h<i-1;h++)
{
remp=a[j][h]*a[i][h]+remp;
}
a[i][j]=(a[i][j]-remp)/a[j][j];
}
}
for(i=1;i<=n;i++)
{
for(j=1;j<=i;j++)
{
a[j][i]=a[i][j];
}
}
for(i=1;i<=n;i++)
{
double dp=0;
for(j=1;j<i;j++)
{
dp=dp+a[i][j]*b[j];
}
b[i]=(b[i]-dp)/a[i][i];
}
for(i=n;i>=1;i--)
{
double gemp;
for(j=i+1;j>=n;j++)
{
gemp=gemp+a[i][j]*b[i];
}
b[i]=(b[i]-gemp)/a[i][i];
}
for(i=1;i<=n;i++)
{
cout<<"x["<<i<<"]:"<<b[i]<<endl;
}
}

你好哦楼主~

很高兴看到你的问题。
但是又很遗憾到现在还没有人回答你的问题。也可能你现在已经在别的地方找到了答案，那就得恭喜你啦。
可能是你问的问题有些专业了，没人会。或者别人没有遇到或者接触过你的问题，所以帮不了你。建议你去问题的相关论坛去求助，那里的人通常比较多，也会比较热心，能快点帮你解决问题。

希望我的回答能够帮到你！
祝你好运。。