平方根法求解方程组编程

用平方根法求解方程组 #include#include#includeusing namespace std;void choleskty(int n,double a10,double b10)int k,p,i;for(k=1;k=n;k+)for(p=1;pk;p+)akk-=akp*akp;akk=sqrt(akk);for(i=k+1;i=n;i+)for(p=1;pk;p+)aik-=aip*akp;aik/=akk;for(i=1;i=n;i+)for(k=1;k=1;i-)for(k=i+1;k=n;k+)bi-=aki*bk;bi/=aii;int main()int n,i,j;double a1010,b10;cout请输入方程组未知数的个数：n;cout请输入方程组的下三角元素：endl;for(i=1;i=n;i+)for(j=1;jaij;cout请输入方程组的最后一列：endl;for(i=1;ibi;for(i=1;i=n;i+)for(j=i+1;j=n;j+)aij=aji;cout此方程组的增广矩阵为：endl;for(i=1;i=n;i+)for(j=1;j=n;j+)coutsetw(10)aij;coutsetw(10)bi;coutendl;choleskty(n,a,b);cout此方程组的解是：endl;for(i=1;i=n;i+)coutsetw(10)bi;coutendl;return 0;