Laboratorio 6. Resultados
PROBLEMA 2a, Método de Jacobi
[ortiz@castor tp5]$ ./jacobi
Ingrese la tolerancia: 1.0e-11
Ingrese el número máximo de iteraciones: 100
Matrix A, vector b e iteración inicial x0:
4.0000e+00
-1.0000e+00 0.0000e+00 -1.0000e+00 0.0000e+00
0.0000e+00
1.0000e+00 0.0000e+00
-1.0000e+00
4.0000e+00 -1.0000e+00 0.0000e+00 -1.0000e+00
0.0000e+00 1.0000e+00
0.0000e+00
0.0000e+00
-1.0000e+00 4.0000e+00 0.0000e+00 0.0000e+00
-1.0000e+00
1.0000e+00 0.0000e+00
-1.0000e+00
0.0000e+00 0.0000e+00 4.0000e+00 -1.0000e+00
0.0000e+00
0.0000e+00 0.0000e+00
0.0000e+00
-1.0000e+00 0.0000e+00 -1.0000e+00 4.0000e+00
-1.0000e+00 0.0000e+00
0.0000e+00
0.0000e+00
0.0000e+00 -1.0000e+00 0.0000e+00 -1.0000e+00
4.0000e+00
0.0000e+00 0.0000e+00
iteración norm
1 7.50e-01
2 4.38e-01
3 2.66e-01
4 1.60e-01
5 9.67e-02
..................
..................
46 9.88e-11
47 5.96e-11
48 3.60e-11
49 2.17e-11
50 1.31e-11
51 7.91e-12
La iteración final es:
4.1614906832e-01
5.0931677018e-01
4.1614906832e-01
1.5527950310e-01
2.0496894410e-01
1.5527950310e-01
Problema 2a, Método de Gauss-Seidel
[ortiz@castor tp5]$ ./seidel
Ingrese la tolerancia: 1.0e-11
Ingrese el número máximo de iteraciones: 100
Matrix A, vector b e iteración inicial x0:
4.0000e+00
-1.0000e+00 0.0000e+00 -1.0000e+00 0.0000e+00
0.0000e+00
1.0000e+00 0.0000e+00
-1.0000e+00
4.0000e+00 -1.0000e+00 0.0000e+00 -1.0000e+00
0.0000e+00 1.0000e+00
0.0000e+00
0.0000e+00
-1.0000e+00 4.0000e+00 0.0000e+00 0.0000e+00
-1.0000e+00
1.0000e+00 0.0000e+00
-1.0000e+00
0.0000e+00 0.0000e+00 4.0000e+00 -1.0000e+00
0.0000e+00
0.0000e+00 0.0000e+00
0.0000e+00
-1.0000e+00 0.0000e+00 -1.0000e+00 4.0000e+00
-1.0000e+00 0.0000e+00
0.0000e+00
0.0000e+00
0.0000e+00 -1.0000e+00 0.0000e+00 -1.0000e+00
4.0000e+00
0.0000e+00 0.0000e+00
iteración norm
1 1.15e+00
2 4.31e-01
3 1.72e-01
4 6.49e-02
5 2.38e-02
..................
..................
22 8.35e-10
23 3.04e-10
24 1.11e-10
25 4.03e-11
26 1.47e-11
27 5.35e-12
La iteración final es:
4.1614906832e-01
5.0931677019e-01
4.1614906832e-01
1.5527950311e-01
2.0496894410e-01
1.5527950311e-01
PROBLEMA 2b, Método de Jacobi
[ortiz@castor tp5]$ ./jacobi
Ingrese la tolerancia: 1.0e-11
Ingrese el número máximo de iteraciones: 100
Matrix A, vector b e iteración inicial x0:
4.0000e+00
-1.0000e+00 0.0000e+00 -1.0000e+00 0.0000e+00
0.0000e+00
1.0000e+00 0.0000e+00
-1.0000e+00
4.0000e+00 -1.0000e+00 0.0000e+00 -1.0000e+00
0.0000e+00 1.0000e+00
0.0000e+00
0.0000e+00
-1.0000e+00 4.0000e+00 0.0000e+00 0.0000e+00
-1.0000e+00
1.0000e+00 0.0000e+00
-1.0000e+00
0.0000e+00 0.0000e+00 4.0000e+00 -1.0000e+00
0.0000e+00
1.0000e+00 0.0000e+00
0.0000e+00
-1.0000e+00 0.0000e+00 -1.0000e+00 4.0000e+00
-1.0000e+00 1.0000e+00
0.0000e+00
0.0000e+00
0.0000e+00 -1.0000e+00 0.0000e+00 -1.0000e+00
4.0000e+00
1.0000e+00 0.0000e+00
iteración norm
1 1.50e+00
2 8.75e-01
3 5.31e-01
4 3.20e-01
5 1.93e-01
...................
...................
47 1.19e-10
48 7.20e-11
49 4.34e-11
50 2.62e-11
51 1.58e-11
52 9.55e-12
La iteración final es:
5.7142857143e-01
7.1428571428e-01
5.7142857143e-01
5.7142857143e-01
7.1428571428e-01
5.7142857143e-01
Problema 2b, Método de Gauss-Seidel
[ortiz@castor tp5]$ ./seidel
Ingrese la tolerancia: 1.0e-11
Ingrese el número máximo de iteraciones: 100
Matrix A, vector b e iteración inicial x0:
4.0000e+00
-1.0000e+00 0.0000e+00 -1.0000e+00 0.0000e+00
0.0000e+00
1.0000e+00 0.0000e+00
-1.0000e+00
4.0000e+00 -1.0000e+00 0.0000e+00 -1.0000e+00
0.0000e+00 1.0000e+00
0.0000e+00
0.0000e+00
-1.0000e+00 4.0000e+00 0.0000e+00 0.0000e+00
-1.0000e+00
1.0000e+00 0.0000e+00
-1.0000e+00
0.0000e+00 0.0000e+00 4.0000e+00 -1.0000e+00
0.0000e+00
1.0000e+00 0.0000e+00
0.0000e+00
-1.0000e+00 0.0000e+00 -1.0000e+00 4.0000e+00
-1.0000e+00 1.0000e+00
0.0000e+00
0.0000e+00
0.0000e+00 -1.0000e+00 0.0000e+00 -1.0000e+00
4.0000e+00
1.0000e+00 0.0000e+00
iteración norm
1 2.04e+00
2 9.74e-01
3 4.27e-01
4 1.71e-01
5 6.31e-02
...................
...................
24 2.94e-10
25 1.07e-10
26 3.90e-11
27 1.42e-11
28 5.18e-12
La iteración final es:
5.7142857143e-01
7.1428571429e-01
5.7142857143e-01
5.7142857143e-01
7.1428571429e-01
5.7142857143e-01
PROBLEMA 3a Método de Jacobi:
[ortiz@castor tp5]$ ./jacobi
Ingrese la tolerancia: 1.0e-11
Ingrese el número máximo de iteraciones: 100
Matrix A, vector b e iteración inicial x0:
3.0000e+00 1.0000e+00 1.0000e+00 5.0000e+00 0.0000e+00
2.0000e+00 6.0000e+00 1.0000e+00 9.0000e+00 0.0000e+00
1.0000e+00 1.0000e+00 4.0000e+00 6.0000e+00 0.0000e+00
iteración norm
1 4.67e+00
2 2.60e+00
3 1.45e+00
4 8.08e-01
5 4.50e-01
..........................
40 5.93e-10
41 3.31e-10
42 1.84e-10
43 1.03e-10
44 5.73e-11
45 3.19e-11
46 1.78e-11
47 9.93e-12
La iteración final es:
1.0000000000e+00
1.0000000000e+00
1.0000000000e+00
PROBLEMA 3a, Métdo de Gauss-Seidel:
[ortiz@castor tp5]$ ./seidel
Ingrese la tolerancia: 1.0e-11
Ingrese el número máximo de iteraciones: 100
Matrix A, vector b e iteración inicial x0:
3.0000e+00 1.0000e+00 1.0000e+00 5.0000e+00 0.0000e+00
2.0000e+00 6.0000e+00 1.0000e+00 9.0000e+00 0.0000e+00
1.0000e+00 1.0000e+00 4.0000e+00 6.0000e+00 0.0000e+00
iteración norm
1 3.46e+00
2 7.90e-01
3 8.16e-02
4 7.64e-03
5 4.58e-04
6 6.50e-05
7 1.44e-05
8 2.09e-06
9 2.37e-07
10 2.36e-08
11 1.82e-09
12 1.09e-10
13 3.12e-11
14 5.28e-12
La iteración final es:
1.0000000000e+00
1.0000000000e-00
1.0000000000e-00
PROBLEMA 3b, Método de Jacobi: DIVERGE!!!!!!!
[ortiz@castor tp5]$ ./jacobi
Ingrese la tolerancia: 1.0e-6
Ingrese el número máximo de iteraciones: 100
Matrix A, vector b e iteración inicial x0:
5.0000e+00 7.0000e+00 6.0000e+00 5.0000e+00 2.3000e+01 0.0000e+00
7.0000e+00 1.0000e+01 8.0000e+00 7.0000e+00 3.2000e+01 0.0000e+00
6.0000e+00 8.0000e+00 1.0000e+01 9.0000e+00 3.3000e+01 0.0000e+00
5.0000e+00 7.0000e+00 9.0000e+00 1.0000e+01 3.1000e+01 0.0000e+00
iteración norm
1 1.42e+01
2 3.52e+01
3 8.71e+01
4 2.16e+02
5 5.34e+02
....................
....................
95 1.45e+38
96 3.59e+38
97 8.90e+38
98 2.20e+39
99 5.45e+39
100 1.35e+40
101 3.34e+40
La iteración final es:
7.7891263324e+39
5.4203503315e+39
5.4941191187e+39
5.1028189791e+39
PROBLEMA 3b, Método de Gauss Seidel: CONVERGE LENTAMENTE....
[ortiz@castor tp5]$ ./seidel
Ingrese la tolerancia: 1.0e-4
Ingrese el número máximo de iteraciones: 2000
Matrix A, vector b e iteración inicial x0:
5.0000e+00 7.0000e+00 6.0000e+00 5.0000e+00 2.3000e+01 0.0000e+00
7.0000e+00 1.0000e+01 8.0000e+00 7.0000e+00 3.2000e+01 0.0000e+00
6.0000e+00 8.0000e+00 1.0000e+01 9.0000e+00 3.3000e+01 0.0000e+00
5.0000e+00 7.0000e+00 9.0000e+00 1.0000e+01 3.1000e+01 0.0000e+00
iteración norm
1 5.49e+00
2 1.46e+00
3 8.64e-01
4 5.08e-01
5 2.96e-01
.........................
.........................
1433 1.02e-04
1434 1.02e-04
1435 1.01e-04
1436 1.01e-04
1437 1.01e-04
1438 1.01e-04
1439 1.00e-04
1440 9.99e-05
La iteración final es:
1.0160190995e+00
9.9037190708e-01
9.9593271014e-01
1.0023906762e+00