jacpra.c

雅克比迭代法

#include "stdio.h"
#include "jacobi.h"

int main(){
	double x[6][6];
	double b[6];
	double x1[6];

	x[1][1]=8;
	x[1][2]=6;
	x[1][3]=0;
	x[1][4]=0;
	x[1][5]=1;
	x[2][1]=1;
	x[2][2]=9;
	x[2][3]=2;
	x[2][4]=0;
	x[2][5]=-1;
	x[3][1]=0;
	x[3][2]=1;
	x[3][3]=15;
	x[3][4]=4;
	x[3][5]=2;
	x[4][1]=-1;
	x[4][2]=0;
	x[4][3]=1;
	x[4][4]=10;
	x[4][5]=6;
	x[5][1]=1;
	x[5][2]=3;
	x[5][3]=2;
	x[5][4]=1;
	x[5][5]=-15;

	b[1]=1;
	b[2]=0;
	b[3]=2;
	b[4]=-1;
	b[5]=1;
	
	x1[1]=0;
	x1[2]=0;
	x1[3]=0;
	x1[4]=0;
	x1[5]=0;

	jacobi(x,5,b,x,0.0000000001);

	return 0;
}