j.h

来自「本实例通过ＶＣ＋＋实现了JACOBI 和SOR 算法」· C头文件 代码 · 共 53 行

#include<iostream.h>

void Jacobi(double **a,double *b,int n)//本函数以课本习题六.3要求为例，即初值和限制次数条件确定。
{
	double *crux=new double[n];
	double *nextx=new double[n];
	double *chax=new double[n];

	int i;
	
	for(i=1;i<=n;i++)
	{
		crux[i-1]=0;
		chax[i-1]=1;
	}


	double max=1;
    while(max>0.0001)
	{
		for(i=1;i<=n;i++)
		{
			double temp1=0;
			double temp2=0;

			for(int j=1;j<=i-1;j++)
				temp1=temp1+a[i-1][j-1]*crux[j-1];
		
			for(int k=i+1;k<=n;k++)
				temp2=temp2+a[i-1][k-1]*crux[k-1];

			nextx[i-1]=(b[i-1]-temp1-temp2)/a[i-1][i-1];
			
		}
		
		for(i=1;i<=n;i++)
			chax[i-1]=nextx[i-1]-crux[i-1];


		max=chax[0];
		for(i=2;i<=n;i++)
			if(max<chax[i-1])
				max=chax[i-1];

		//给crux赋值
		for(i=1;i<=n;i++)
			crux[i-1]=nextx[i-1];
	}

	for(i=1;i<=n;i++)
		cout<<nextx[i-1]<<" ";
	cout<<endl;
}