Unidad 3: Tarea 3 - Diferenciación e Integración Numérica y EDO

DANIA CAROLINA GONZÁLEZ

Código: 1.075.313.018

Tutor: Martín Gómez Orduz

Curso: 100401A_761

UNIVERSIDAD NACIONAL ABIERTA Y A DISTANCIA – UNAD

ESCUELA DE CIENCIAS BÁSICAS, TECNOLOGÍA E INGENIERÍA ECBTI
MÉTODOS NUMÉRICOS
MAYO 2020
NEIVA – HUILA
 Problema 5
La siguiente tabla representa datos físicos tomados igualmente espaciados:

i 0 1 2 3 4 5 6
x 0 0,5 1 1,5 2 2,5 3
f(X) 1 0,8 0,2 0,25 0,31 0,38 0,44
0,5 0,5 0,5 0,5 0,5 0,5

Evaluando los puntos dados en la anterior tabla en geogebra nos da que la función
en la siguiente:

f(x)=−0,27 x 6+
2,61 x5 −9,88 x 4 +17,69 x 3 - 14,43 x 2 +3,47x+1
f'(x)=-1,62x5+13,05x4-39,52x3+53,07x2-28,86x+3,47
 x0 1,5
h 0,39
x f(x)
0,16034
f(x0-h) 1,11 4
0,16796
f(x0) 1,5 9
0,01185
f(x0+h) 1,89 1
error rel
1292,35
Adelante -0,4003 2
0,01955 168,005
Atrás 2 8
-
Central 0,02896 0,71655
-
Exacta 0,02875

0.3

0.2

0.1

0
1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2 Puntos
hacia adelante
hacia atras
central
exacta

Análisis

Como se puede notar la mejor aproximación de la primera derivada se obtiene a través de la

Primera Diferencia atrás ya que en este se obtiene un 53% de error.
el problema es el tamaño de h ya que con un paso de 0,25 el error por cada uno de los métodos es
muy grande. si h se modifica a 0,35 el error disminuye en la diferencia finita central este queda
con 32% de error.

El mejor método es la diferenciación finita central porque esta es mucho más precisa que la otras
dos.

los otros dos métodos se utilizan cuando no se tiene información de los extremos.

INTEGRACION NUMERICA

= 0,6035

a 1
b 3
n 6
h 0,33333333

k a+k+*h f(a+hk)
0 1 0
1 1,33333333 0,21576155
2 1,66666667 0,30649537
3 2 0,34657359
4 2,33333333 0,36312765
5 2,66666667 0,36781097
6 3 0,3662041

Método Valor Error rel

Trapecio
Simple 0,3662041 39,3199509
Trapecio
compuesto 0,65532441 8,58730955
Simpsom 1/3
Simple 0,95037025 57,4764289
Simpsom 1/3
Compuesto 1,17137612 94,09712
Simpsom 3/8
Simple 0,31078211 48,503378
Simpsom 3/8
Compuesto 0,60236724 0,1876981
"Exacta" 0,6035
 1.4

1.2
Funcion evaluada en cada uno de
1 los paso
Trapecio SImple
0.8 Trapecio compuesto
Simpson 1/3 simple
0.6 simpson compuesto
simpson 3/8 simple
0.4
simpson 3/8 compuesta
exacta
0.2

0
0.5 1 1.5 2 2.5 3 3.5

trapecio Compuesta

trapecio simple

Simpson 1/3 Simple

Simpson 1/3 Compuesta

Simpson 3/8 Simple

Simpson 3/8 Compuesta

 EDO METODO DE EULER Y RK GRADO 2 Y 4

Ejercicio 5

y ´=− y (t)
solucion exacta

y=e−t+C 1

t0 0
yt0 0
h 0,1
a 0,1
b 0,5

Formula Euler

Solución

RK Orden 4 RK orden 2
n ti Yexacta Yeuler K1 k2 k3 k4 YRK orden 4 k2 YRK orden 2
0 0 1 1 - - - - 1 - 1
-
0,90483 - 0,1005012 0,996656
1 0,1 7 1 0 -0,05 0,050125 5 2 0,075 1,005
-
0,81873 - - 0,2024708 0,972380 -
2 0,2 1 0,99 0,0997 -0,9974 0,156979 3 9 0,026247 0,999928
- -
0,74081 - 0,9748 - 0,3054994 0,940142 -
3 0,3 8 0,9702 0,1945 1 0,459505 1 7 0,127903 0,984919
- -
0,9410 0,9450 - 0,3879058 0,908288 -
4 0,4 0,67032 9 -0,282 8 0,296219 4 6 0,226822 0,960396
- -
0,60653 0,9034 - 0,9164 0,4756118 0,871876 -
5 0,5 1 5 0,3633 6 -0,42935 2 3 0,320058 0,926948
- -
0,54881 0,8582 - 0,8838 - 0,5533560 0,832259 -
6 0,6 2 8 0,4359 6 0,503838 5 6 0,404822 0,885429
 - -
0,49658 0,8067 - 0,8484 - 0,6223798 0,789946 -
7 0,7 5 8 0,4994 9 0,568545 8 8 0,478728 0,836868
- -
0,44932 0,7503 0,8106 - 0,6817863 0,745475 -
8 0,8 9 1 -0,553 8 0,622861 2 3 0,539838 0,782447
-
0,6902 - 0,7708 - 0,699399
9 0,9 0,40657 8 0,5964 2 0,666414 -0,730905 9 -0,58674 0,723452
0,36787 0,6281 - - -
10 1 9 6 0,6295 -0,7293 0,699072 -0,7693071 0,652281 0,618595 0,66123
- -
0,33287 0,5653 - 0,6865 - 0,7968122 0,604671 -
11 1,1 1 4 0,6523 3 0,720938 1 7 0,635151 0,597144
- -
0,30119 0,5031 - 0,6429 0,8134868 -
12 1,2 4 5 0,6651 2 -0,73234 6 0,557107 0,636732 0,532524
- -
0,27253 0,4427 - 0,5988 - 0,8196349 0,510092 -
13 1,3 2 7 0,6685 9 0,733814 4 5 0,624196 0,468626
- -
0,24659 0,3852 - 0,5548 - 0,8157803 0,464095 -
14 1,4 7 1 0,6631 5 0,726077 7 3 0,598868 0,406598
-
0,3312 - 0,8026430 0,419535 -
15 1,5 0,22313 8 0,6497 -0,5112 -0,71 4 7 0,562459 0,347443
- -
0,20189 0,2815 - 0,4683 - 0,7811089 0,376780 -
16 1,6 7 9 0,6293 1 0,686574 8 8 0,516961 0,292002
- -
0,18268 0,2365 - 0,4265 - 0,7521962 0,336140 -
17 1,7 4 4 0,6028 2 0,656876 8 2 0,464545 0,240937
- -
0,16529 0,1963 - 0,3861 - 0,7170182 -
18 1,8 9 3 0,5714 4 0,622033 2 0,297863 0,407452 0,194726
- -
0,14956 0,1609 - 0,3474 - 0,6767451 0,262137 -
19 1,9 9 9 0,5362 6 0,583186 4 3 0,347887 0,153662
0,13533 - - 0,229090 -
20 2 5 0,1304 0,4981 -0,3107 0,541461 -0,6325668 9 0,287931 0,117864
- -
0,12245 0,1043 - 0,2760 - 0,5856565 0,198793 -
21 2,1 6 2 0,4582 5 0,497932 1 8 0,229456 0,087294
- -
0,11080 0,0824 - 0,2436 - 0,5371386 0,171262 -
22 2,2 3 1 0,4175 7 0,453601 3 5 0,174069 0,061774
23 2,3 0,10025 0,0642 - - - - 0,146464 - 0,04101
 0,2136 0,4880602
9 8 0,3768 5 0,409376 4 8 0,123072
- -
0,09071 - 0,1860 - 0,4393680 0,124325 -
24 2,4 8 0,0495 0,3369 5 0,366053 5 8 0,077433 0,024619
- -
0,08208 0,0376 - 0,1608 - 0,3918909 0,104734 -
25 2,5 5 2 0,2984 8 0,324306 6 9 0,037794 0,012153
- -
0,07427 0,0282 - 0,1381 - 0,3463285 0,087552 -
26 2,6 4 1 0,2618 2 0,284684 4 1 0,004478 0,003127
- -
0,06720 0,0208 - 0,1177 0,3032452 0,072615
27 2,7 6 8 0,2276 1 -0,24761 6 3 0,022481 -0,00296
- -
0,0152 - 0,0995 - 0,2630702 0,059747
28 2,8 0,06081 4 0,1961 7 0,213384 8 2 0,043296 -0,00661
- -
0,05502 0,0109 - 0,0835 0,2261020
29 2,9 3 7 0,1673 9 -0,18219 3 0,048761 0,05838 -0,0083

1.2

0.8

0.6

0.4

0.2

0
0 0.2 0.4 0.6 0.8 1 1.2 1.4

Exacta Euler
RK orden 4 RK orden 2

Formulas RK orden 4
Formulas RK orden 2