Gabriela Cáceres - Cálculo Numérico

Clculo Numrico
Trabajo Prctico II
Alumna: Gabriela Cceres -44653

Profesor: Hugo Franco Paats
Ao 2.014
CAPITULO 4
1. Usa el mtodo de la eliminacin de Gauss con pivoteamiento
parcial para resolver el siguiente sistema de ecuaciones,
muestra en cada etapa la matriz transformada y los elementos
que intervienen:
31 - 2
+ 4 = 2,3
41 + 22 - 3 - 54 = 6,9
-51 + 2 - 3 3
= -16,8
102 - 4 3 + 74 = -36
Seleccionar el elemento mayor en mdulo de la columna 1
5 1 3 0 16,8
6,9
4
2 1 5
|
||
|
2,3
3 1 0
6
36
0 10 4 7
K=1
31 = 31 =
11 = -5
41 = 41 = 0
11
11
21 = 21 =
11
5 1
14
0
5
|
2
|
0
5
0 10
17
5
9
5
4
5
0 16,8
327
5
| | 50 |
| 389 |
6 |
50
7
36
Reordenando, seleccionando el mayor pivot para 22
5 1
0 10
|
2
|0
5
0
3
4
9
14
5
17
0 16,8
36
7
| | 389 |
6 | | 50 |
327
50
K= 2
32 = 32 =
22 = 10
42 = 42 =
22
5 1
0 10
0
22
3
4
49
25
57
25
1
25
7
25
0
16,8
7
36
157 | | 461 |
25 | | 50 |
174 177
50
25
Reordenando, seleccionando el mayor pivot para 33
5 1
0 10
0
3
4
57
25
49
25
43 = 43 =
K=3
33 =
33
49
57
57
25
5 1
0 10
|
0
|0
0
0
16,8
7
36
174| | 177 |
25 | | 50 |
157 461
50
25
3
4
57
25
0
7
16,8
36
174| | 177 |
25 | | 50 |
233
233
19
19
= ,
= ,
2. Dado el siguiente sistema de ecuaciones, utilice el mtodo de

descomposicin LU con pivoteamiento parcial para obtener la
solucin del sistema
2 , + = ,
-0,5 + , = ,
, , = ,
Como 11 es el mayor entonces no se reordena
2
0,9
3
0,1
1 |
|0,5
1
6,35 0,45
1
K= 1
21 = 4
11 = 2
31 = 2
2
1
| 4
1
2
0,9
3
18
14 |
5910 3920
Reordenando, para que 22 sea el mayor
2 0,9
1
|2
1
4
59
3
39
20|
1
4
10
32 = 236
K=2
22 =
59
10
2 0,9
3
1
59
39
|
|
10
20 |
|2
1
5
197
4 236
944
1
|
1
2
0
1
5
236
0
3,61
1
0| | | = |15,401|
2
2,035
1 3
Resolviendo se obtiene:
1 = 3,61
2 0,9
|0
59
10
2 = 17,206
4531
5900
3 =
3,61
39 1
20 | |2 | = |17,206|
4531
197 3
5900
944
Resolviendo se obtiene:
= ,
= ,
= ,
3. De los tres sistemas de ecuaciones lineales presentados a

continuacin, identifique aquel que no resolvera usando un
mtodo iterativo como el de Gauss- Seidel. Establezca
claramente el criterio de convergencia
8x + 3y + z = 12
-x + 4y + 5z = 8
x + y + 5z = 7
-6x
-2x + 2y - 3z = -3
x + 4y - z = 4
+7z= 1
2x + 4y z = 5
2y z = 1
3x + y - z = 3
El sistema que no resolvera con el mtodo de Gauss- Seidel es

-x + 4y + 5z = 8
-2x + 2y - 3z = -3
2y z = 1
Para que cumpla el criterio de convergencia < 1
Aplicando el criterio de convergencia para Gauss- Seidel por el criterio de
Sassenfeld
1) 1=
2=
3+1
8
6+7
0
= 0,5
No cumple el criterio, reordenando las ecuaciones

-6x
+7z= 1
8x + 3y + z = 12
2x + 4y z = 5
7
1 = 6 = 1,167
7
6
(8)( )+1
2=
= 3,44
3 =
7
6
(2)( )+4(3,44)
= 16,093
No cumple el criterio, reordenando las ecuaciones

8x + 3y + z = 12
2x + 4y z = 5
-6x
+7z= 1
1 =
3+1
8
= 0,5
2(0,5)+1
2=
3 =
2) 1 =
= 0,5
6(0,5)
= 0,42 CUMPLE EL CRITERIO DE LA CONVERGENCIA
4+5
=9
1
2(9)+3
2=
= 10,5
3 =
2
2(9)+4(10,5)
1
= 60
No converge, reordenando el sistema

-2x + 2y - 3z = -3
-x + 4y + 5z = 8
2y z = 1
1 =
2+3
2
= 2,5
2=
1(2,5)+5
3 =
= 1,875
2(1,875)
= 3,75
NO CONVERGE EL SISTEMA
3) 3x + y - z = 3
x + 4y - z = 4
x + y + 5z = 7
1 =
2 =
1+1
3
= 0,67
1(0,67)+1
3 =
= 0,41
1(0,67)+0,41
5
= 0,216
CUMPLE EL CRITERIO DE LA CONVERGENCIA
4. Dado el sistema de ecuaciones, resuelva por el mtodo de

eliminacin de Gauss con pivoteamiento parcial
6x + 2y +2z = -2
2x + y + z = 1
x + 2y - z = 0
Como 11 es el mayor, no se reordena
6 2 2
2 1 2
|2
|| 1 |
3 3
0
1 2 1
1
K=1
21 = 3
11 = 6
31 = 6
6 2
2 2
1 5
|0 0 3| | 3 |
|
|| |
5
4 1
0
3
3 3
Como 22 debe de ser mayor entonces se reordena
6
2
5
3
|0
|
0
2
4
||
3| |
1
32 = 0
K=2
22
6
|0
|
0
5
=
3
2
2 2
5
4 1
|| |
3
3| | 3 |
1 5
0
3 3
Resolviendo
=
2
1
3 ||
5
3
=
=
5. Hallar las matrices de Jacobi y de Gauss Seidel para la matriz y

determina si los mtodos iterativos son o no convergentes

A = [ ]

B = []
Matriz de Jacobi
(+) = ( 6+ )
(+) = (-9 + + )
(+) =
(8 + )
Matriz de Gauss- Seidel
(+) = (6 + )
(+) = ( -9 + (+) + )
(+) = ( 8 + (+) )
Por Gauss- Jacobi, aplicar el criterio de convergencia
1 =
1+0
2 =
1+1
3 =
0+1
= 0,5
=1
= 0,5
= 1 como max < 1 entonces no cumple el criterio de convergencia

Por Gauss-Seidel, aplicar el criterio de convergencia
1 =
1+0
2
= 0,5
2 =
3 =
(1.0.5)+1
= 0,75
2
0+(1.0,75)
2
= 0.375
= 0,75 como max<1, entonces cumple el criterio de convergencia
6. Para el ejercicio anterior, determina el valor inicial en base al

estudio de convergencia y resuelva para el criterio de parada <
2x1 x2 + 0 = 6
x1 =
-x1 +2 x2 x3 = -9
x2 =
0 - x2 + 2x3 = 8
x3 =
6+ 2
2
9 + 1 + 3
2
8+ 2
2
Tomando X(0) = 0 por el criterio de convergencia en el problema 5 entonces

K=1
6+ 0
x1 =
x2 =
=3
9 + 3 + 0
2
8 3
x3 =
= -3
= 2,5
k=2
x1 =
x2 =
x3 =
6 3
= 1,5
9 + 1,5 + 2,5
2
8 2,5
2
= -2,5
= 2,75
k=3
6 2,5
x1 =
x2 =
= 1,75
9 + 1,75 + 2,75
2
8 2,25
x3 =
= - 2,25
= 2,875
k=4
6 2,25
x1 =
x2 =
= 1,875
9 + 1,875 + 2,875
2
8 2,125
x3 =
= -2,125
= 2,9375
k=5
6 2,125
x1 =
x2 =
= 1,9375
9 + 1,875 + 2,875
2
8 2,125
x3 =
= -2,0625
= 2,96875
k=6
6 2,0625
x1 =
x2 =
= 1,96875
9 + 1,96875 + 2,96875
2
8 2,03125
x3 =
= -2,03125
= 2,984375
k=7
x1 =
x2 =
6 2,03125
= 1,984375
9 + 1,984375 + 2,984375
2
= -2,015625
8 2,015625
x3 =
= 2,9921875
k=8
x1 =
x2 =
x3 =
6 2,015625
2
= 1,9921875
9 + 1,9921875+ 2,9921875
2
8 2,0078125
2
= - 2,0078125
= 2,99609375
Entonces :
2,99609375 - 2,9921875 = 0,00390
x1 = 1,9921875
- 2,0078125 + 2,015625 = 0,00781
x2 = - 2,0078125
1,9921875 - 1,984375 =
x3 = 2,99609375
Criterio de Parada E< x10-2
0,00789
en k = 8 E = 0,00789 < x10-2
7. Dado el siguiente sistema de ecuaciones lineales:
4 - + = 4
+ , + 3,5 = 1
- + 4,25 + , = -1
a) Haga el estudio de convergencia para la aplicacin de los
mtodos iterativos y determina el valor inicial en base a
este estudio.
b) Resuelva, utilizando el mtodo de Gauss-Seidel utilizando
como criterio de parada un error relativo menor que 0,5% o
3 iteraciones.
c) Convergencia de Gauss-Jacobi
41 2 + 3 = 4
1 =
1 + 2,752 + 3,53 = 1
2 =
1 + 4,752 + 2,753 = 1
3 =
1+1
= 0,5
4
1+3,5
= 1,6363
2,75
1+4,25
2,75
= 1,9
> 1
1 + 2,752 + 3,53 = 1
1 =
41 2 + 3 = 4
2 =
1 + 4,752 + 2,753 = 1
3 =
2,75+3,5
1
4+1
= 6,25
=5
1
1+4,25
2,75
= 1,9
> 1
41 2 + 3 = 4
1 =
1 + 4,752 + 2,753 = 1
2 =
1 + 2,752 + 3,53 = 1
3 =
1+1
= 0,5
4
1+2,75
4,25
1+2,75
3,5
= 0,88
= 1,07
> 1
No converge con Gauss-Jacobi
Criterio de Gauss-Seidel
41 2 + 3 = 4
1 =
1 + 2,752 + 3,53 = 1
2 =
1 + 4,752 + 2,753 = 1
3 =
1+1
= 0,5
4
0,5+3,5
= 1,454
2,75
0,5+4,25(1,454)
2,75
= 2,4297
> 1
41 2 + 3 = 4
1 =
1 + 4,752 + 2,753 = 1
2 =
1 + 2,752 + 3,53 = 1
3 =
1+1
= 0,5
4
1(0,5)+2,75
= 0,7647
4,25
1(0,5)+2,75(0,7647)
3,5
= 0,74369
< 1
Converge
d) Para los valores iniciales
4
1 (0) = 1
2 (0) = 17
3 (0) = 7
< 0,005 3
1
1 (+1) = (4 + 2 3 )
4
4
2 (+1) =
(1 + 1 (+1) 2,753 )
17
2
3 (+1) = (1 1 (+1) 2,752 (+1) )
7
Para k = 0
1
Parada
4
207
207
11 = 4 (4 17 7) = 238
4
207
|238 1| = 0,13025
2
436
436
2925
21 = 17 (1 + 238 2,75 (7)) = 2023

2
207
31 = 7 (1 238 2,75 ( 2023)) = 14161
436
| 2023 + 17| = 0,019

2925
|14161 7| = 0,079
31
238
=
= 0,149 > 0,005
207
238
Para k = 1
1
Parada
436
2925
1 2 = 4 (4 2023 14161) = 0,89448

4
|0,89448 0,87| = 0,024
2925
2 2 = 17 (1 + 0,89448 2,75 (14161)) = 0,1585
436
|0,1585 + 2023| = 0,057
2925
3 2 = 7 (1 0,89448 2,75(0,1585)) = 0,1547
|0,1547 14161| =
0,0518
=
0,057
= 0,0637 > 0,005
0,89448
Para k = 2
Parada
1 3 = 4 (4 0,1585 0,1547) = 0,9217
|0,9217 0,89448| = 0,02722
|0,1185 + 0,1585| = 0,04
|0,11548 0,1547| =
2 3 = 17 (1 + 0,9217 2,75(0,1547)) = 0,1185

3 3 = 7 (1 0,9217 2,75(0,1185)) = 0,11548
0,039
=
0,04
= 0,043 > 0,005
0,9217
Para k = 3
Parada
1 4 = 4 (4 0,1185 0,11548) = 0,9415
|0,9415 0,9217| = 0,019
2 4 = 17 (1 + 0,9415 2,75(0,11548)) = 0,0885 |0,885 + 0,1185| = 0,03

2
3 4 = 7 (1 0,9415 2,75(0,0885)) = 0,08625

0,029
=
0,03
= 0,03
0,9415
|0,08625 0,11548| =
8. Para el siguiente circuito:

a) Determina los valores de la corriente en el punto d y por la
resistencia de 5
b) Repita los clculos si se cambia el valor de la fuente por
5V
a) Por el mtodo LU
10 5 0
|5 14 2|
0 2 5
K=1
21 =
11 = 10
31 = 0
1
2
10 5 0
1 23
|
2|
2 2
0 2 5
4
32 = 23
K=2
22 =
10
1
| 2
|
0
23
2
5
23
2
2 |
|
4 107
23 23
1
1
| 2
|
0
0
1
4
23
0
1
1
0|
| | 2 | = |0|
3
0
1
Resolviendo
1
1 = 1
2 = 2
3 =
2
23
1
10 5 0
1
1
23
2| | | = | 2 |
|0
2
2
2
107
3
0
0
23
23
Resolviendo
=
I= ( - )
I=
I= 0,0766
b)
1
1
| 2
|
0
0
1
4
23
0
1
5
0|
| | 2 | = |0|
3
0
1
Resolviendo
1 = 5
2 = 2
3 =
5
10 5 0
1
23
5
2| | | = | |
|0
2
2
2
107
10
3
0
0
23
23
10
23
I= ( - )
I= -
I= 0,38318
9. Utilice uno de los mtodos numricos para resolver:

0,3x + 0,52y + z = -0,01
0,5x + y + 1,9z = 0,67
0,1x + 0,3y + 0,5 z = -0,44
Por el Mtodo de eliminacin de Gauss
0,3 0,52 1 0,01
1
1,9| | 0,67 |
|0,5
0,1 0,3 0,5 0,44
5
K=1
21 = 3
11 = 0,3
31 =
0,3 0,52
|0
|
0
0,01
32 =
19
20
2
15
0,3 0,52
|0
0
2
7 103
||
|
15 30| | 150 |
19
1 131
150 6 300
K= 2
22 =
2
15
Resolviendo
0,01
103
7
30 | | 150 |
1089
11
200 1000
z=
y=
x=
CAPITULO 5
1. Encontrar el polinomio cuadrtico que mejor aproxima a f(x)=
en [0; 1] por mnimos cuadrados
Polinomio cuadrtico a + bx + c = 0
x
y
x^2
x^3
x^4
0
2,71828
0
0
0
g1 (x) = x 2
g 2 (x) = x
g 3 (x) = 1
0,2
2,2255
0,04
0,008
0,0016
0,4
1,8221
0,16
0,064
0,0256
0,6
1,4918
0,36
0,216
0,1296
0,8
1,2214
0,64
0,512
0,4096
1 = a
2 = b
3 = c
g1 (x) g1 (x) + g1 (x) g 2 (x) + g1 (x) g 3 (x) = f(x) g1 (x)

1,5664 1 + 1,8 2 +2,2 3 = 2,6993
g1 (x) g 2 (x) + g 2 (x) g 2 (x) + g 2 (x) g 3 (x) = f(x) g 2 (x)
1,8 1 + 2,2 2 + 3 3 = 4,04614
2,2 1 + 32 + 6 3 = 10,47908
1 = 0,84326
2 = -2,5481
3 = 2,71138
, , x + 2,71138 = 0
1
1
1
1
1
3
10,47908
2,2
1,8
1,5664
2. Considere la siguiente tabla de valores y encuentre un

polinomio de segundo grado que aproxime a la funcin por
mnimos cuadrados
x
y
x
y
x^2
x^3
x^4
0
3,1437
0,4
4,4169
0
3,1437
0
0
0
0,4
4,4169
0,16
0,064
0,0256
0,8
6,0203
0,8
6,0203
0,64
0,512
0,4096
1,2
8,6512
1,2
8,6512
1,44
1,728
2,0736
1,6
11,0078
1,6
11,0078
2,56
4,096
6,5536

61 + 62 + 8,8 3 = 49,456
61 + 8,82 + 14,4 3 = 67,00912
8,81 + 14,42 + 25,06243 = 110,061792
1 = 3,324189286
2 = 1,20461875
3 = 2,532176339
Y= , + , + ,
2
16,2161
2
16,2161
4
8
16
6
49,456
8,8
14,4
25,0624
3. Aproxima los valores de la tabla por un polinomio de 2 grado

aplicado mnimos cuadrados
x
y
1
1,941
x
y
x^2
x^3
x^4
1,5
1,512
1
1,941
1
1
1
2
1,177
1,5
1,512
2,25
3,375
5,0625
2
1,177
4
8
16
2,5
0,917
3
0,714
2,5
0,917
6,25
15,625
39,0625
3
0,714
9
27
81
1 x 2 + 2 x +3 = 0
g1 (x) =1
g 2 (x) = x
g 3 (x) = x 2
51 + 102 + 22,53 = 6,261
101 + 22,52 + 553 = 10,9975
22,51 + 552 + 142,1253 =22,20825
1 = 2,9988
2 = 1,212085714
3 = 0,1505714286
= , , + ,
10
6,261
22,5
55
142,125
4. Resuelva el problema anterior para aproximar por una funcin

del tipo y= A
x
y
Ln(y)
1
1,941
0,6632
1,5
1,512
0,41343
2
1,177
0,16297
2,5
0,917
-0,08665
y= AeBx No es lineal, se aplica linealizacin

Ln(y) = Ln(AeBx )
Ln(y) = Ln(A) + Bx
1 = Ln(A)
g1 (x) = 1
2 = B
g 2 (x) = x
g1 (x) g1 (x) + g1 (x) g 2 (x) = f(x) g1 (x)

5 1 + 10 2 = 0,81608
g1 (x) g 2 (x) + g 2 (x) g 2 (x) = f(x) g 2 (x)
10 1 + 22,5 2 = 0,38205
1 = 1,163304
2 = 0,5
Como:
A= e 1
B= 2
Y= e1,163304 e0,5x
Y= 3,2 ,
3
0,714
-0,33687
10
6,261
0,81608
5. Para los valores de la tabla ajusta una funcin del tipo f(x)= a
+ por el mtodo de los mnimos cuadrados
x
y
x
y
x^1/2
1/x
1
1,941
1
1,941
1
1
1,5
1,512
1,5
1,512
1,22474
0,67
g1 (x) = x
g 2 (x) =
1
x
10a + 3,7334b = 8,143932192

3,7334a + 1,966b = 4,1423
a = 0,095444
b = 1,92572
f(x) = , +
2
1,177
2
1,177
1,414
0,5
2,5
0,917
2,5
0,917
1,581
0,4
3
0,714
1,732
0,33
3
0,714
10
6,261
6,95174
2,9
6. Dado los puntos [1; 0,5] , [ 2; 1,7] , [3 ; 3,4] y [4; 5,7]. Encontrar la
ecuacin del tipo y= que aproxima dichos puntos utilizando
el mtodo de mnimos cuadrados con linealizacin de los datos.
x
y
1
0,5
x
Log (y)
Log(x)
2
1,7
1
-0,301
0
3
3,4
2
0,2304
0,301
4
5,7
3
0,5314
0,47712
y= ax b Aplicando linealizacin
Log(y) = log(ax b )
Log(y) = Log(a) + b log(x)
1 = Log(a)
g1 (x) = 1
2 = b
g 2 (x) =log(x)
g1 (x) g1 (x) + g1 (x) g 2 (x) = f(x) g1 (x)

4 1 + 1,38022 2 = 1,21667
g1 (x) g 2 (x) + g 2 (x) g 2 (x) = f(x) g 2 (x)
1,380221 + 0,6808 2 = 0,778
1 = -0,3
2 = 1,7511
Como 1 = Log(a) entonces
-0,3 = log(a)
a= 0,50119
b= 1,7511
y= 0,50119,
4
0,75587
0,6021
7. En el circuito RC la corriente est dada por la funcin del tiempo
I= en donde I se mide en amperes (A), R en ohms y C en

faradios (F). Considerando los valores de la siguiente tabla:
Calcule el valor de I y el producto RC mediante el mtodo de los
mnimos cuadrados.
t
I(A)
t
I(A)
ln(I)
1
8,187
2
6,703
3
5,488
2
6,703
1,90256
3
5,488
1,70256
1
8,187
2,1025
4
4,493
5
3,679
4
4,493
1,5025
5
3,679
1,3026
Ln(I) = Ln(I0 eRC )

t
Ln(I) = ln (I0 ) + RC
g1 (x) = 1
1 = Ln(I0 )
g 2 (x) = t
2 =
1
RC
g1 (x) g1 (x) + g1 (x) g 2 (x) = f(x) g1 (x)

51 + 152 = 8,51272
g1 (x) g 2 (x) + g 2 (x) g 2 (x) = f(x) g 2 (x)
151 + 552 = 23,5383
1 = 2,302502
2 = 0,2
Como 1 = Ln(I0 ) entonces
2,302502 = Ln(I0 )
=
Como 2 =
1
RC
=
1
RC
entonces
0,2 =
I= ,
15
28,55
8,51272
8. Ajuste, minimizando la funcin , la curva y= a+ b a la tabla

de datos
x
y
0
5,1
x
y
e^(x)^2
1
8,2
0
5,1
1
1,5
22
1
8,2
2,71828
g1 (x) = 1
2
g 2 (x) = ex
2
112,1
1,5
22
9,4877
1 = a
2 = b
g1 (x) g1 (x) + g1 (x) g 2 (x) = f(x) g1 (x)

51 + 585,8172 = 1186,4
g1 (x) g 2 (x) + g 2 (x) g 2 (x) = f(x) g 2 (x)
585,8171 + 271,42k 2 = 544,572k
1 = 2,95206
2 = 2
Y = 2,95206 + 2
2
112,1
54,598
2,5
1039
2,5
1039
518,01
9. Se desea que una va de comunicacin que une las ciudades

A(0;0) B(5;0) pase prxima a los pueblos C(4;1) D(3;2). Por
razones de ndole tcnica, el trazado de la va debe coincidir con
la grfica del tipo h(x)= a + bx + c
x
y
0
0
5
0
4
1
3
2
x
y
x^2
x^3
x^4
0
0
0
0
0
5
0
25
125
625
4
1
16
64
256
3
2
9
27
81
g1 (x)=1
g2 (x)=x
g 3 (x) = x 2
41 + 12 2 + 50 3 = 3
121 + 50 2 + 216 3 = 10
501 + 216 2 + 962 3 = 34
1 = 0,0138
2 = 1,5691
3 = 0,3177
h(x) = , + , x ,
10.
Obtener la parbola de ajuste de mnimos cuadrados que
aproxima a la funcin f(x)= en los puntos (-1; 0; 8)

x
y
x^2
x^3
x^4
-1
-1
1
-1
1
0
0
0
0
0
8
2
64
512
4096
7
1
65
511
4097
1 x 2 + 2 x +3 = 0
g1 (x) = x 2
g 2 (x) = x
g 3 (x) = 1
40971 + 5112 + 653 = 127
5111 + 652 + 73 = 17
651 + 72 + 33 = 1
1 = -0,0833
2 = 0,9167
3 = 0
-0,0833 + 0,9167x = 0
11.
Utiliza el mtodo de linealizacin de datos para hallar un
ajuste exponencial a los datos que se muestran en la tabla
realizando los cambios de variable correspondiente:
a) Aproximar por y= C
b) Aproximar por y=
x
f(x)
-1
6,62
0
3,94
1
2,17
2
1,35
3
0,89
a)
x
y
Ln(f(x))
-1
6,62
1,89
0
3,94
1,37118
1
2,17
0,7747
y= CeAx
Por linealizacin
Ln(y) = Ln (CeAx )
Ln(y) = Ln(C) + Ax
g1 (x) = 1
g 2 (x) = x
1 = Ln(c)
2 = A
g1 (x) g1 (x) + g1 (x) g 2 (x) = f(x) g1 (x)

51 + 52 = 4,21938
g1 (x) g 2 (x) + g 2 (x) g 2 (x) = f(x) g 2 (x)
51 + 152 = 0,8648
1 = 1,352294
1 = 0,508418
Como 1 = Ln(c) entonces
1,352294 = Ln(c)
C = 3,86628
Y= 3,86628 ,
2
1,35
0,3
3
0,89
-0,1165
b)
1
y
= Ax + B
x
1/y
-1
0,151
1 = A
2 = B
0
0,2538
1
0,4608
2
0,74074
g1 (x) = x
g 2 (x) = 1
g1 (x) g1 (x) + g1 (x) g 2 (x) = f(x) g1 (x)

151 + 52 = 5,16204009
g1 (x) g 2 (x) + g 2 (x) g 2 (x) = f(x) g 2 (x)
51 + 5 2 = 2,730030248
1 = 0,2432009842
2 = 0,3028050654
1
y
= 0,2432x + 0,3028
, + ,
3
1,12359
CAPITULO 6
1. Si se utiliza la Regla del Trapecio Repetida para aproximar el valor
de I = 2(x)dx
a) Qu nmero de sub intervalos sera suficiente elegir para que el
error sea menor que x10 -3?
b) Calcular un valor aproximado de la integral utilizando 5 sub
intervalos.
a) f(x) = Log2(x)
1
M2 = 0,868 entonces :
1
f `(x) = 2. . 10
2
f ``(x) =10 . 2 (
1
10
E=

12
. h2 . M2
f ``(1) = 0,3772 ; f `` (1,5) = 0,1047 ; f `` (2) = 0,0399
b = 2; a = 1;
h = 0,17836300
m >
Entonces : m >
5,606
m=6
b) h = 0,2
I=
I=
2
0,2
2
I = 1 2(x)dx
[ f(x0) + 2 (f(x1) + f(x2) + f(x3) + f(x4) ) + f(x5) ]

[ 0 + 2 (6,269 x10 -3 + 0,21353 + 0,041664 + 0,065164) + 0,090619 ]
I = 0,07438742
Valor Real = 0,0355
2.
Integral la funcin f(x) en el intervalo [0; 4]
x
y
0
-4,271
a)
b)
0,5
-2,522
1
-0,499
1,5
1,795
2
4,358
2,5
7,187
3
10,279
3,5
4
13,633 17,247
Por el mtodo del trapecio.

Es posible aplicar la regla de Simpson de 1/3? Por qu? Si
su respuesta es positiva calcule el valor aproximado.
c) Es posible aplicar la regla de Simpson de 3/8? Por qu? Si

su respuesta es positiva
calcule el valor aproximado.
40
a)
I=
h=
= 0,5
[ f(x0) + 2 (f(x1) + f(x2) + f(x3). + f(xn-1) ) + f(xn) ]
2
0,5
I = 2 [ -4,271 + 2 (-2,522 0, 499 + 1,795 + 4,358 + 7,187 + 10,279 + 13,633 )

+ 17,347 ]
I = 20,3595
b)
Para saber si aplicar la regla de Simpson de 1/3 m debe ser par
h = Xi +1 Xi
Como m = 8 entonces es par
h = 3,5 3 = 0,5
I=
[ f(x0) + 4 (f(x1) + f(x3)+ f(xn-1)) + 2 (f(x2) + f(x4 + .. f(xn-2) ) + f(xn) ]
0,5
I = 3 [ -4,271 + 4 (-2,522 + 1,795 + 7,187 + 13,633 ) + 2 ( - 0, 499 + 4,358 +

10,279) + 17,347 ]
I = 20,2706667
c) No es posible aplicar la regla del trapecio 3/8 ya que no cumple la
condicin de que m sea mltiplo de 3.
3.
Evala la siguiente integral I = ( + ) dx por el mtodo

de Simpson 1/3 con 4 decimales de precisin.
F(x) = ( + 1)
F(x) = ( + 2)
F(x) = ( + 3)
F(x) = ( + 4)
(x) = ( + 5)

E>
180
h <
I=
. h4 . 4
0,14246
(x) (0) =
3
(x) (0,75) = 7,9

(x) (1,5) = 29,13
entonces 4 = 29,13
Donde E = 1 x10 -4
b = 1,5 ; a = 0
entonces h = 0,125 y m = 12
[ f(x0) + 4 (f(x1) + f(x3) + f(x5) + f(xn-1)) + 2 (f(x2) + f(x4 + .. f(xn-2) ) + f(xn) ]
0,125
I = 3 [ 1 + 4 (1,27479 + 2,0006 + 3,035899 + 4,49789 + 6,54546 + 9,39330) + 2

(1,60503 + 2,47308 + 3,70475 + 5,43656 + 7,85327) + 11,20422 ]
I = 6,7225165
4.
Se evala I = () = por el mtodo del Simpson 1/3 con

un error igual a x10 -3. Si los valores de la funcin en los
extremos es igual a 0, determina:
a) El valor de f(1)
b) El valor mximo de (x)
Si f(0) y f(2) = 0 si tomamos h = 1 y hacemos el mtodo de Simpson 1/3
a)
I=
[ f(x0) + 4 (f(x1) + f(x2) ] = 4

1
Entonces 4 = 3 ( 0 + 4 f(1) + 0 )
f(1) = 3
b)
E>

180
. h4. M4 donde b = 2; a = 0 ; h = 1
Entonces M4 = 0,09
y E = x10 -3
5.
Sabiendo que I =
Calcula un valor aproximado de por
+
el mtodo del trapecio y de Simpson 1/3 utilizando 10 sub
intervalos.
h = 0,1 porque son 10 subintervalos.

x
0,1
0,2
0,3
0,4
0,5
0,6
0,7
0,8
0,9
Valor real = 0,7853981634

Mtodo del trapecio
I=
[ f(x0) + 2 (f(x1) + f(x2) + f(x3). + f(xn-1) ) + f(xn) ]
0,1
I = 2 [ 1 + 2 (0,990 + 0,961 + 0,917 + 0,862 + 0,8 + 0,735 + 0,671 + 0,609 +

0,552) + 0,5 ]
I = 0,7847
Mtodo Simpson 1/3
I=
[ f(x0) + 4 (f(x1) + f(x3) + f(x5) + f(xn-1)) + 2 (f(x2) + f(x4 + .. f(xn-2) ) + f(xn) ]
0,1
I = 3 [ 1 + 4 (0,990 + 0,917 + 0,8 + 0,671 + 0,552) + 2 (0,961 + 0,862 + 0,735 + 0,609)

+0,5]
I = 0,7518
Una cuerda vibra adoptando la forma Y = senx, entre las

abscisas X = 0 y X = 4 en un instante "To Calcule el valor
aproximado de la longitud de la cuerda, utilizando un mtodo
numrico con m = 8.
6.
I = 0 1 + `()2
I=
4
0 1
+ ()2
y = senx ; y`= cos x

Longitud de la cuerda
m = 8 ; h = 0,5
x
I=
0,5
1,5
2,5
3,5
[ f(x0) + 2 (f(x1) + f(x2) + f(x3). + f(xn-1) ) + f(xn) ]
0,5
I = 2 [ 1,4142 + 2 (1,3304 + 1,1366 + 1,0024 + 1,083 + 1,2813 + 1,4071 + 1,37)

+1,1946 ]
I = 4,95765
Valor Exacto = 4,96661
Evaluar la funcin I = () usando la siguiente tabla:
7.
0
0
X
F(X)
0,7
0,1
6,84
1,2
0,3
4,00
I = 0 () + 0,7 ()
I=
I1 =
I1 =
I2 =
0,7
5,51
0,95
5,77
1,2
1,00
0,7
1,2
donde I1 = 0 () y I2 = 0,7 ()
[f(x0) + 2 (f(x1) + f(x2) + f(x3). + f(xn-1) ) + f(xn) ]
0,2
2
[0 + 2 ( 6,84 + 4 + 4,2) + 5,51 ]
3,559
0,25
2
[5,51+ 2 (5,77) + 1 ]
I2 = 2,25625
I = I 1 + I2
I=
0,5
4,20
5, 81525
h = 0,25
h = 0,2
La regla del Trapecio aplicada a () nos da valor 4 y la

regla de Simpson 1/3 nos da valor 2. Qu valor tiene f(1)?
8.
Tomando m = 2 entonces h = 1 donde X0 = 0 ; X1 = 1 ; X2 = 2

I=
I=
[ f(x0) + 2 (f(x1) + f(x2)] = 4 trapecio

[ f(x0) + 4 (f(x1) + f(x3)]
Si h = 1 entonces
= 2 Simpson 1/3
f(x0) + 2 f(x1) + f(x2) = 8

f(x0) + 4 f(x1) + f(x3) = 6
Restando ambas ecuaciones

-2 f(x1) = 2
f(x1) = -1
9.
Calcule la derivada de f(x) cuyos valores se definen en la tabla

siguiente, en x = 2 utilizando 3 puntos, utilizando: a) las
diferencias progresivas b) las diferencias regresivas c) las
diferencias centrales
x
f(x)
1
1,941
1,5
1,512
2
1,177
2,5
0,917
3
0,714
h = 2 1,5 = 0,5
F `(x) =
F `(x) =
F `(x) =
(i + h ) (i )
h
(i ) (i h )
h
(i+h ) (i h )
2h
0,9171,177
0,5
1,1771,512
0,5
= - 0,52
Progresiva
= - 0 ,67
Regresiva
0,9171,512
0,5 x 2
= - 0,595
Central
Aproxima el valor de la derivada funcin y =
10.
en x = 1,8,
por medio de las diferencias progresivas, regresivas y

centrales. Utiliza en todos los caso = 0,2. Calcula el error
cometido comparando con el valor exacto de la derivada y` =
+
( +)/
h = 0,2
f(1,6) = 0,8547
F(x) =
(i + h ) (i )
F (x) =
(i ) (i h )
F(x) =
(i+h ) (i h )
2h
; f(2,0) = 1 ; f(1,8) = 0,9304
10,9304
0,2
= 0,348
0,93040,8547
0,2
10,8547
0,2 x 2
Progresiva
= 0 ,3785
= 0,36325
Error Absoluto
Valor Real = 0,3626983869 = f `( 1,8)
E1 = 0,3626983869 0,348 = 0,01469
E2 = 0,3626983869 0,3785 = 0,01580
E3 = 0,3626983869 0,36325 = 0,000551
Regresiva
Central
11. Utiliza
las frmulas de las diferencias progresivas, regresivas y

centrales para aproximar el valor de la derivada de la funcin
f(x) = (1 + x) en el punto x = 0,6 para un polinomio de
aproximacin de grado dos. Repite los clculos para un
polinomio de grado tres.
0
1
x
F(x)
0,2
0,4
0,6
0,8
1
1,2
1,4656 2,0885 2,9153 4,0059 5,4365 7,3042
h = 0,2
Progresivas
Segundo Orden
F`(x1) =
f(x3 )+ 4f(x2 ) 3f(x1 )

2h
Tercer Orden
F`(x1) =
x1 = 0,6
=
5,4365+4x4,00593.x2,9153
2x 0,2
= 4,603
x1 = 0,6
2f(x4 )9 f(x3 )+ 18f(x2 ) 11f(x1 )

6h
2x7,3042 9x5,4365 + 18x4,0059 11x2,9153

6x0,2
Regresivas
Segundo Orden
F`(x4) =
3f(x4 ) 4f(x3 )+ f(x2 )

2h
Tercer Orden
F`(x4) =
x4 = 0,6
=
3x2,9153 4x2,0885+ 1,4656

2x0,2
x4 = 0,6
11f(x4 )18 f(x3 )+ 9f(x2 ) 2f(x1 )
6h
11x2,915318x2,0885+ 9x1,4656 2x1
= 4,7214
Central
Segundo Orden
F`(x3) =
f(x4 ) f(x2 )
2h
Tercer Orden
= 4,64375
x3 = 0,6
4,0059 2,0885
2x0,2
x3 = 0,6
= 4,7935
6x0,2
= 4,76483
F`(x3) =
f(x5 )+ 8 f(x4 ) 8f(x2 )+ f(x1 )

12h
5,4365+ 8x4,0059 8x2,0885+ 1,4656

12x0,2
= 4,7367
12.
Aproxima el valor de la derivada de la funcin definida en la

tabla en x = 2,0 por aproximacin de un polinomio de 2 grado y
de 3er grado.
x
F(x)
0,0
-4,271
0,5
-2,522
1,0
-0,499
1,5
1,795
2,0
4,358
2,5
7,187
3,0
3,5
4,0
10,279 13,633 17,247
a) las diferencias progresivas

b) las diferencias regresivas
c) las diferencias centrales
d) repite todos los puntos anteriores para evaluar la derivada en x =
4,0
h = 0,5
a) Progresivas
Segundo Orden
F`(x1) =
f(x3 )+ 4f(x2 ) 3f(x1 )

2h
Tercer Orden
F`(x1) =
x1 = 2,0
=
10,279+(4.7,187)(3.4,358)
2x 0,5
= 5,395
x1 = 2,0
2f(x4 )9 f(x3 )+ 18f(x2 ) 11f(x1 )

6h
2x13,633 9x10,279 + 18x7,187 11x4,358

6x0,5
= 5,394
b) Regresivas
Segundo Orden
F`(x4) =
3f(x4 ) 4f(x3 )+ f(x2 )

2h
Tercer Orden
F`(x4) =
x4 = 2,0
=
3x4,358 4x1,795+(0,499)
= 5,395
x4 = 2.0
11f(x4 )18 f(x3 )+ 9f(x2 ) 2f(x1 )

6h
2x0,5
11x4,35818x1,795+ 9x(0,499) 2x(2,522)

6x0,5
= 5,3937
c) Central
Segundo Orden
F`(x3) =
f(x4 ) f(x2 )
2h
x3 = 2,0
7,1871,795
Tercer Orden
2x0,5
x3 = 2,0
f(x5 )+ 8 f(x4 ) 8f(x2 )+ f(x1 )
F(x3) =
= 5,392
12h
10,279+ 8x7,187 8x1,795+(0,499)

12x0,5
= 5,393
d) Progresivas
Segundo Orden
x1 = 4,0
f(x3 )+ 4f(x2 ) 3f(x1 )
F(x1) =
2h
Tercer Orden
= No es posible calcular con los datos de la tabla
x1 = 4,0
2f(x4 )9 f(x3 )+ 18f(x2 ) 11f(x1 )
F(x1) =
6h
Regresivas
Segundo Orden
3f(x4 ) 4f(x3 )+ f(x2 )
F(x4) =
2h
Tercer Orden
3x17,247 4x13,633+(10,279)
2x0,5
6h
Central
Segundo Orden
F`(x3) =
f(x4 ) f(x2 )
2h
= 7,488
x4 = 2.0
11f(x4 )18 f(x3 )+ 9f(x2 ) 2f(x1 )
F(x4) =
11x17,24718x13,633+ 9x(10,279) 2x7,187

6x0,5
= 7,4867
x3 = 4,0
Tercer Orden
F`(x3) =
x4 = 4,0
x3 = 4,0
f(x5 )+ 8 f(x4 ) 8f(x2 )+ f(x1 )

12h
CAPITULO 7
1. Considera la ecuacin diferencial = ( )/ con condicin

inicial () = . Utiliza el mtodo de Euler para obtener la
solucin aproximada en el intervalo [0; 3] en cuatro pasos.
Evala el error cometido sabiendo que la solucin exacta es:
= / + .
Siendo el intervalo [0; 3] y considerando 4 pasos tendremos:
=
|1 0 | 3 0
=
= 0,75
Luego:
0 = 0;
0 = 1;
(; ) = ( )/2
Entonces:
1 = 0 + = 0,75
0,75(0 1)
1 = 0 + (0 ; 0 ) = 1 +
= 0,625
2
Resumiendo los resultados en la siguiente tabla:
n
0,75
0,625
1,5
0,671875
2,25
0,982421875
1,457763672
En 4 pasos
() ,
El valor real :
= 3 /2 + 2
= 3 3/2 + 3 2
= 1,66939048
El error:
||= 1,66939048 1,457763672 = 0,211626808

1,66939048 1,457763672
|% | = |
100%| = , %
1,66939048
2. Repita el problema anterior utilizando el mtodo de Heun.

Siendo el intervalo [0; 3] y considerando 4 pasos tendremos:
=
|1 0 | 3 0
=
= 0,75
Luego:
0 = 0;
0 = 1;
(; ) = ( )/2
Entonces:
1 = 0 + = 0,75
1 = 0 + (0 ; 0 ) = 1 +
0,75(0 1)
= 0,625
2
(0 ; 0 ) + (1 ; 1 )
1 = 0 + [
] = 1,679450439
2

n
0,75
0,625
0,8359375
1,5
0,8037109375
0,950378418
2,25
1,156486511
1,258466244
1,630291402
1,701199185
En 4 pasos:
() ,
Como el valor real es:
= 1,66939048
El error:
|| = 1,66939048 1,701199185 = 0,03180870538
1,66939048 1,701199185
|% | = |
100%| = , %
1,66939048
3. Utiliza el mtodo de Euler con h = 0,2 para estimar la solucin

del siguiente problema de valor inicial:
= ;
= ;
() =
Siendo 0 = 1, 0 = 0 y (, ) = tendremos:
=
|1 0 | 2 1
=
=5
0,2
Por lo tanto, sern necesarios cinco pasos para llegar a la aproximacin

deseada.
Luego:
1 = 0 + = 1,2
1 = 0 + (0 ; 0 ) = 0 + 0,2( 0 ) = 0,2
n
1,2
0,2
1,4
0,4442805516
1,6
0,7561541331
1,8
1,182167831
1,834455188
En 5 pasos:
() ,
4. Sea la ecuacin diferencial = + , con la condicin

inicial () = , usando el mtodo de Euler con h = 0,5
aproximar y(1).
Siendo 0 = 0, 0 = 2 y (, ) = + 2 tendremos:
=
|1 0 | 1 0
=
=2
0,5
Por lo tanto, sern necesarios dos pasos para llegar a la aproximacin deseada.
Luego:
1 = 0 + = 0,5
1 = 0 + (0 ; 0 ) = 2 + 0,5(0 2 + 2) = 2
n
0,5
2,25
En 2 pasos:
() ,
5. Considere el problema de valor inicial

= , ( );
x [0; 1]; y(0) = 1
Resuelva por el mtodo de Runge-Kutta de tercer orden en
cuatro pasos.
Siendo el intervalo [0; 1] y considerando la cantidad de pasos tendremos:
=
|1 0 | 1 0
=
= 0,25
Luego:
0 = 0;
0 = 1;
(; ) = 0,5( 2 )
Para poder calcular el valor de 1 , debemos calcular primero los valores de

1 , 2 3 .
1 = (0 ; 0 ) = 0,125
1
2 = (0 + ; 0 + ) = 0,115234375
2
2
3 = (0 + ; 0 + 22 1 ) = 0,1040039063
1
1 = 0 + (1 + 42 + 3 ) = 0,8850097656
6
1 = 0,25
Calculando el valor de 2 :
1 = (1 ; 1 ) = 0,1028137207
1
2 = (1 + ; 1 + ) = 0,08662223816
2
2
3 = (1 + ; 1 + 22 1 ) = 0,1723947525
1
2 = 1 + (1 + 42 + 3 ) = 0,781393528
6
2 = 0,5
1 = (2 ; 2 ) = 0,066424191
1
2 = (2 + ; 2 + ) = 0,04469455406
2
2
3 = (2 + ; 2 + 22 1 ) = 0,02449107636
1
3 = 2 + (1 + 42 + 3 ) = 0,7364446141
6
3 = 0,75
1 = (3 ; 3 ) = 0,02174307676
1
2 = (3 + ; 3 + ) = 0,005006490538
2
2
3 = (3 + ; 3 + 22 1 ) = 0,02897491601
1
4 = 3 + (1 + 42 + 3 ) = 0,740987581
6
1 = 1
En 4 pasos:
() ,
6. Repita el problema anterior pero utilizando el mtodo de RungeKutta de 4 orden. Evala el error para ambos mtodos si la
solucin exacta est dada por = + / .
Siendo el intervalo [0; 1] y considerando la cantidad de pasos tendremos:
=
|1 0 | 1 0
=
= 0,25
Luego:
0 = 0;
0 = 1;
(; ) = 0,5( 2 )

1 , 2 , 3 4 .
1 = (0 ; 0 ) = 0,125
1
2 = (0 + ; 0 + ) = 0,115234375
2
2
2
3 = (0 + ; 0 + ) = 0,1158447266
2
2
4 = (0 + ; 0 + 3 ) = 0,1027069092
1
1 = 0 + (1 + 22 + 23 + 4 ) = 0,8850224813
6
1 = 0,25
1 = (1 ; 1 ) = 0,1028153102
1
2 = (1 + ; 1 + ) = 0,08662372828
2
2
2
3 = (1 + ; 1 + ) = 0,08763570215
2
2
4 = (1 + ; 1 + 3 ) = 0,06842334739
1
2 = 1 + (1 + 22 + 23 + 4 ) = 0,7983962282
6
2 = 0,5
1 = (2 ; 2 ) = 0,068549522853
1
2 = (2 + ; 2 + ) = 0,046687058
2
2
2
3 = (2 + ; 2 + ) = 0,0480534624
2
2
4 = (2 + ; 2 + 3 ) = 0,02348034573
1
3 = 2 + (1 + 22 + 23 + 4 ) = 0,7514777424
6
3 = 0,75
1 = (3 ; 3 ) = 0,0236222178
1
2 = (3 + ; 3 + ) = 0,003244795814
2
2
2
3 = (3 + ; 3 + ) = 0,001565607463
2
2
4 = (3 + ; 3 + 3 ) = 0,03086958127
1
4 = 3 + (1 + 22 + 23 + 4 ) = 0,7542891041
6
4 = 1
Entonces:
() ,
El error para el primer mtodo:
|| = 0,754285382 0,740987581 = 0,01329780101
|% | =
|0,754285382 0,740987581|
100% = , %
0,754285382
El error para el segundo mtodo:

|| = 0,754285382 0,7542891041
|% | =
|0,754285382 0,7542891041|
100% = , %
0,754285382
7. Sea la ecuacin diferencial = + , con la condicin

inicial y(0) = 2, usando el mtodo de Runge-Kutta de cuarto
orden con h = 0,5 aproximar y(1).
Siendo 0 = 0, 0 = 2 y (, ) = + 2 tendremos:
=
|1 0 | 1 0
=
=2
0,5
Por lo tanto, sern necesarios dos pasos para llegar a la aproximacin deseada.
1 , 2 , 3 4 .
1 = (0 ; 0 ) = 0
1
2 = (0 + ; 0 + ) = 0,125
2
2
2
3 = (0 + ; 0 + ) = 0,09375
2
2
4 = (0 + ; 0 + 3 ) = 0,203125
1
1 = 0 + (1 + 22 + 23 + 4 ) = 2,106770833
6
1 = 0,5
1 = (1 ; 1 ) = 0,1966145833
1
2 = (1 + ; 1 + ) = 0,22724609375
2
2
2
3 = (1 + ; 1 + ) = 0,253499349
2
2
4 = (1 + ; 1 + 3 ) = 0,3198649089
1
2 = 1 + (1 + 22 + 23 + 4 ) = 2,368170844
6
2 = 1
De lo cual, concluimos que la aproximacin buscada es:
() ,
8. Un paracaidista de masa M kg salta de un avin en t=0.

Consideremos que la velocidad vertical inicial del paracaidista
es cero en t=0 y que la cada es vertical. Si el arrastre
aerodinmico est dado por = , donde c es una
constante y v es la velocidad vertical (positiva hacia abalo),
asuma M = 70 kg, c = 0,27 kg/m. Sabiendo que la ecuacin de
equilibrio
de
fuerzas
satisface
()
= + ,
( =
, / ), encuentra la velocidad del paracaidista para t 20

seg en cuatro pasos.
Considerando que:
()
= +
= 2
M = 70 kg
c = 0,27
Tendremos:
0,27 2 + 686
70
Siendo el intervalo [0; 20] seg y considerando la cantidad de pasos tendremos:
|1 0 | 20 0
=
=
=5
4
Luego:
=
0 = 0;
0 = 0;
(; ) =
0,27 2 +686
70

1 , 2 , 3 4 .
1 = (0 ; 0 ) = 49
1
2 = (0 + ; 0 + ) = 37,42375
2
2
2
3 = (0 + ; 0 + ) = 42,24741058
2
2
4 = (0 + ; 0 + 3 ) = 14,57801434
1
1 = 0 + (1 + 22 + 23 + 4 ) = 37,15338925
6
1 = 5
1 = (1 ; 1 ) = 22,37849501
1
2 = (1 + ; 1 + ) = 3,929082665
2
2
2
3 = (1 + ; 1 + ) = 19,48875906
2
2
4 = (1 + ; 1 + 3 ) = 12,8749929
1
2 = 1 + (1 + 22 + 23 + 4 ) = 46,54325351
6
2 = 10
1 = (2 ; 2 ) = 7,221849944
1
2 = (2 + ; 2 + ) = 0,4879116286
2
2
2
3 = (2 + ; 2 + ) = 6,782742984
2
2
4 = (2 + ; 2 + 3 ) = 5,842050969
1
3 = 2 + (1 + 22 + 23 + 4 ) = 49,19677154
6
3 = 15
1 = (3 ; 3 ) = 2,322355058
1
2 = (3 + ; 3 + ) = 0,09291290006
2
2
2
3 = (3 + ; 3 + ) = 2,234158152
2
2
4 = (3 + ; 3 + 3 ) = 2,013424494
1
4 = 3 + (1 + 22 + 23 + 4 ) = 50,02395032
6
4 = 20
() = () , /
9. El circuito elctrico que se muestra en la figura tiene una

resistencia de R = 12 y una inductancia de 4H. Si la batera
proporciona un voltaje de 60V y el interruptor se cierra en t=0.
Utiliza el mtodo de Runge-Kutta de 4 orden y 4 pasos para
determinar el valor de la corriente transcurrido 1 seg.
Siendo el intervalo [0; 1] seg y considerando la cantidad de pasos tendremos:

=
|1 0 | 1 0
=
= 0,25
Planteando la ecuacin del circuito:
+ =
4 + 12 = 60
= 15 3
Luego:
0 = 0;
0 = 0;
(; ) = 15 3

1 , 2 , 3 4 .
1 = (0 ; 0 ) = 3,75
1
2 = (0 + ; 0 + ) = 2,34375
2
2
2
3 = (0 + ; 0 + ) = 2,87109375
2
2
4 = (0 + ; 0 + 3 ) = 1,596679688
1
1 = 0 + (1 + 22 + 23 + 4 ) = 2,629394531
6
1 = 0,25
1 = (1 ; 1 ) = 1,777954102
1
2 = (1 + ; 1 + ) = 1,111221313
2
2
2
3 = (1 + ; 1 + ) = 1,361246109
2
2
4 = (1 + ; 1 + 3 ) = 0,7570195198
1
2 = 1 + (1 + 22 + 23 + 4 ) = 3,876045942
6
2 = 0,5
1 = (2 ; 2 ) = 0,8429655433
1
2 = (2 + ; 2 + ) = 0,5268534645
2
2
2
3 = (2 + ; 2 + ) = 0,6453954941
2
2
4 = (2 + ; 2 + 3 ) = 0,3589189227
1
3 = 2 + (1 + 22 + 23 + 4 ) = 4,467109673
6
3 = 0,75
1 = (3 ; 3 ) = 0,3996677454
1
2 = (3 + ; 3 + ) = 0,2497923409
2
2
2
3 = (3 + ; 3 + ) = 0,3059956175
2
2
4 = (3 + ; 3 + 3 ) = 0,1701710322
1
4 = 3 + (1 + 22 + 23 + 4 ) = 4,747345455
6
4 = 1
() = () ,

Gabriela Cáceres - Cálculo Numérico

Cargado por

Copyright:

Formatos disponibles

Gabriela Cáceres - Cálculo Numérico

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Gabriela Cáceres - Cálculo Numérico

Cargado por

Copyright:

Formatos disponibles

Clculo Numrico

Alumna: Gabriela Cceres -44653

Reordenando, seleccionando el mayor pivot para 22

Reordenando, seleccionando el mayor pivot para 33

2. Dado el siguiente sistema de ecuaciones, utilice el mtodo de

Reordenando, para que 22 sea el mayor

3. De los tres sistemas de ecuaciones lineales presentados a

El sistema que no resolvera con el mtodo de Gauss- Seidel es

No cumple el criterio, reordenando las ecuaciones

No cumple el criterio, reordenando las ecuaciones

= 0,42 CUMPLE EL CRITERIO DE LA CONVERGENCIA

No converge, reordenando el sistema

CUMPLE EL CRITERIO DE LA CONVERGENCIA

4. Dado el sistema de ecuaciones, resuelva por el mtodo de

5. Hallar las matrices de Jacobi y de Gauss Seidel para la matriz y

Matriz de Gauss- Seidel

= 1 como max < 1 entonces no cumple el criterio de convergencia

= 0,75 como max<1, entonces cumple el criterio de convergencia

6. Para el ejercicio anterior, determina el valor inicial en base al

Tomando X(0) = 0 por el criterio de convergencia en el problema 5 entonces

2,99609375 - 2,9921875 = 0,00390

- 2,0078125 + 2,015625 = 0,00781

Criterio de Parada E< x10-2

en k = 8 E = 0,00789 < x10-2

7. Dado el siguiente sistema de ecuaciones lineales:

21 = 17 (1 + 238 2,75 (7)) = 2023

31 = 7 (1 238 2,75 ( 2023)) = 14161

| 2023 + 17| = 0,019

1 2 = 4 (4 2023 14161) = 0,89448

|0,89448 0,87| = 0,024

2 2 = 17 (1 + 0,89448 2,75 (14161)) = 0,1585

|0,1585 + 2023| = 0,057

3 2 = 7 (1 0,89448 2,75(0,1585)) = 0,1547

1 3 = 4 (4 0,1585 0,1547) = 0,9217

|0,9217 0,89448| = 0,02722

|0,1185 + 0,1585| = 0,04

2 3 = 17 (1 + 0,9217 2,75(0,1547)) = 0,1185

1 4 = 4 (4 0,1185 0,11548) = 0,9415

|0,9415 0,9217| = 0,019

2 4 = 17 (1 + 0,9415 2,75(0,11548)) = 0,0885 |0,885 + 0,1185| = 0,03

3 4 = 7 (1 0,9415 2,75(0,0885)) = 0,08625

8. Para el siguiente circuito:

9. Utilice uno de los mtodos numricos para resolver:

g1 (x) g1 (x) + g1 (x) g 2 (x) + g1 (x) g 3 (x) = f(x) g1 (x)

2. Considere la siguiente tabla de valores y encuentre un

g1 (x) g1 (x) + g1 (x) g 2 (x) + g1 (x) g 3 (x) = f(x) g1 (x)

3. Aproxima los valores de la tabla por un polinomio de 2 grado

4. Resuelva el problema anterior para aproximar por una funcin

y= AeBx No es lineal, se aplica linealizacin

g1 (x) g1 (x) + g1 (x) g 2 (x) = f(x) g1 (x)

+ por el mtodo de los mnimos cuadrados

10a + 3,7334b = 8,143932192

g1 (x) g1 (x) + g1 (x) g 2 (x) = f(x) g1 (x)

7. En el circuito RC la corriente est dada por la funcin del tiempo

I= en donde I se mide en amperes (A), R en ohms y C en

Ln(I) = Ln(I0 eRC )

g1 (x) g1 (x) + g1 (x) g 2 (x) = f(x) g1 (x)

8. Ajuste, minimizando la funcin , la curva y= a+ b a la tabla