雅克比矩阵程序

⑴ C++ 牛顿迭代解非线性方程组雅克比矩阵出现奇异

给你一个kernel做参考吧
void Jacobi(double* A , double* x , double lagma , unsigned int n)
{
unsigned int i , j , k=0 ;
unsigned int L=20 ;
double temp;
//double mu = 1e-12;//精度
//double Fanshu1; // x的一范数
double* AA=(double*)malloc((n+1)*(n+1));
for(i=1;i<=n;i++)
{
*(x+i)=1;//迭代初始值设置为1
*(AA+i*n+i) = *(A+i*n+i)-lagma; // AA=A-lagma*I
}
while(1)
{
for(i = 1 ; i<= n ; i++)
{
temp = 0;
//Fanshu1=0;
for(j = 1 ; j <= n ; j++)
if(j!=i)
temp -= *(AA+i*n+j) * *(x+j);
*(x+i) = temp/ *(AA+i*n+i);
}

k++;
if(k==L)
{
//printf("Jacobi法求解Ax=b时 迭代次数超过最大值");
break;
}
}
free(AA);
}