예제 #1
0
파일: Practica6.c 프로젝트: Nigsia/MN1
int resol(int n, double **A, double *b, int *p){
	int res = 0;
	//Primer permutam el vector b.
	permutarB(p, b, n);

	//Despres resolem el sistema Lz = Pb
	res = triinf(n, A, b, TOL);
	//Despres resolem el sistema Ux = z
	res += trisup(n, A, b, TOL);

	//Si la funcio ens retorna un valor > 1 vol dir que haura fallat en alguna part dels algorismes anteriors.
	//Haurem de parar els algorismes.
	return res;

}
예제 #2
0
파일: 2ex.c 프로젝트: TresLos/Universitat
int main (void)
{
	int n;
	double tol = 1e-14;
	double **A, *b;

	for (n = 1; n <= 15; n++)
	{
		A = GMH (n);
		b = GV1 (n, A);

		gauss (n, A, b, tol);
		trisup (n, A, b, tol);

		printf ("%3d:\t%.30le\n", n, Norm (n, b)-1);

		FM (A, n, n);
		FV (b, n);
	}

	return 0;
}