Universidad ECCI.

Aplicacin Mtodos numricos

APLICACIN DE MTODO NUMRICO

Caso I Aplicada al Crecimiento Poblacional

Barreto Natalia
(Cd.: 11114)
Universidad ECCI
Resumen
Este tiene como finalidad presentar la aplicacin de una ecuacin diferencial
de primer orden y aplicar el mtodo de Runge Kutta y en este sentido
comparar los resultados del problema.
Palabras clave: Modelo matemtico, mtodo, ecuacin, interrogante,
poblacin, tiempo, crecimiento.
Abstract
This is to present the application of a differential equation of the first order
and develop the method Runge Kuta to compare the results of Problem.
Keywords: Mathematical model , method , equation, question , people , time,
growth.

yeininatalia@gmail.com

Universidad ECCI.

Aplicacin Mtodos numricos

1. INTRODUCCIN
A diario los estudiantes
se
preguntan para que ver, modelos
matemticos, este es uno de los
tantos casos en que se puede
Implementar un mtodo numrico y
con base en el resultado se ayuda
por ejemplo: en la toma de
decisiones frente a un proyecto
urbanstico, control de riesgos por
sobrepoblacin, cobertura en salud,
escasez de empleo, abastecimiento
en servicios pblicos, entre otros.
Tambin
puede
aplicarse
al
crecimiento bacteriano.
Los mtodos numricos son una
herramienta
muy
valiosa
e
importante para la resolucin de
problemas prcticos de Ingeniera y
otras reas , los mtodos numricos
son tcnicas mediante las cuales es
posible formular problemas de
manera que puedan resolverse
utilizando operaciones aritmticas, o
tambin se puede interpretar como
el diseo y anlisis de algoritmos
necesarios para resolver problemas
de ciencia e ingeniera por ello el
objetivo de este trabajo es presentar
de manera sencilla el mtodo
numrico de Runge-Kutta y la
implementacin de los sistemas
numricos en el computador.

2.

OBJETIVO GENERAL

Plantear, aplicar y resolver una

ecuacin diferencial para
dar
respuesta a un interrogante directo
2

de la cotidianidad o vivencia que

surge en la vida real.

Objetivos especficos

Obtener una aproximacin de cuanto

aumenta
la
poblacin,
en
determinado lapso de tiempo.
Emplear el mtodo de Runge Kutta
para aproximar soluciones de
ecuaciones diferenciales.

3. DISEO EXPERIMENTAL

Ejercicio de aplicacin
El crecimiento de una ciudad es
proporcional al nmero de
habitantes, si la poblacin inicial es
de 4 habitantes y al cabo de 3 aos
es de 4,5.
a. Cunto tiempo tardar en
duplicarse
la
poblacin
inicial?
b. Qu poblacin habr en 10
aos?

Datos iniciales
Po = 4 H
t = 3 Aos
Pf en tres aos = 4,5H
2(poblacin inicial) = 8 H.
a. T=?
P=2 (400)
b. T=10 aos
P=?

yeininatalia@gmail.com

Universidad ECCI.

Aplicacin Mtodos numricos

Para resolver el interrogante b. tenemos

Solucin

Por lo tanto tenemos una

diferencial del siguiente tipo:

( )

ecuacin

( )
(

Calcular k

Habitantes

Aplicacin del Mtodo Runge Kutta

Para la aplicacin del mtodo de runge kutta

se tienen en cuenta los siguientes datos:

Ecuacin diferencial:

Y (4)=3 h=0,01 find Y (5)

Siendo Xo= La poblacin
Yo= Tiempo

k1 4*(e^0,039210*3)

Entonces tenemos la ecuacin de la forma:

( )
Para resolver el interrogante a. tenemos
( )

( )

( )

4,49931139

k2 Xo+h/2
Yo+h/2*K1

x
y

4,05
3,22496557

k3 Xo+h/2
Yo+h/2*K2

x
y

4,05
3,22979573

k4 Xo+h
Yo+h*K3

x
y

4,1
3,45967851

y1 Yo+h/6(K1+2K2+2K3+K4) 3,459673027

)
( )

( )

yeininatalia@gmail.com

Universidad ECCI.

Aplicacin Mtodos numricos

k2

k3

k4

KUTTA

0,1 4,4993114

4,5959146

4,5967851

4,6956708

3,460

4,1

3,4596730

0,1 4,6956697

4,7968909

4,7978429

4,9015459

3,939

4,2

3,9394511

0,1 4,9015447

5,0077882

5,0088314

5,1177789

4,440

4,3

4,4403271

0,1 5,1177774

5,2294937

5,2306392

5,3453077

4,963

4,4

4,9633830

0,1 5,3453060

5,4629985

5,4642592

5,5851819

5,510

4,5

5,5097997

0,1 5,5851799

5,7094130

5,7108038

5,8385800

6,081

4,6

6,0808696

0,1 5,8385776

5,9699864

5,9715246

6,1068296

6,678

4,7

6,6780101

0,1 6,1068267

6,2461285

6,2478345

6,3914327

7,303

4,8

7,3027798

0,1 6,3914293

6,5394372

6,5413350

6,6940949

7,957

4,9

7,9568976

0,1 6,6940908

6,8517308

6,8538487

7,0167617

8,642

8,6422645

0,1 7,0167567

7,1850881

7,1874597

7,3616623

9,361

5,1

9,3609897

0,1 7,3616563

7,5418963

7,5445618

7,7313637

10,115

5,2

10,1154220

1,1 7,7313564

7,9249111

7,9279189

8,1288372

10,908

k1

RUNGE KUTTA

CRECIMIENTO DE LA
POBLACION
15,000
10,000
5,000
0,000

Ecuacin analtica

4 4,14,24,34,44,54,64,74,84,9 5 5,15,2
TIEMPO

5. CONCLUSIONES DEL CASO

1. En 10 aos crecer menos la

poblacin que en 17, 65 Aos,
teniendo sentido.

2. La poblacin crecer ms que la

poblacin inicial.

4. ANALISIS DE RESULTADOS
EXPERIMENTAL

3. A medida que el tiempo avanza la

poblacin crece.

CRECIMIENTO DE LA
POBLACION
10

8
6

0
3

10

POBLACION
17

yeininatalia@gmail.com

Universidad ECCI.

Aplicacin Mtodos numricos

Caso I Aplicada a Ley de Enfriamiento

Resumen
El experimento tiene como finalidad evidenciar un proceso de variacin de temperatura, para un
objeto de uso cotidiano como lo es el computador de mesa. Nos centraremos en la funcin del
disipador el cual es un aparato electrnico ubicado dentro del PC.

Lo que queremos reflejar son los resultados de la funcin como modelo matemtico de Ley de
enfriamiento de Newton y que esta a su vez sea un comparativo con nuestros datos
experimentales; para calcular la temperatura mxima a la que permite llegar la placa base antes
de que se encienda el sistema de ventilacin. Y con base a esto podamos realizar un proceso de
verificacin sobre la calidad del equipo o una debida programacin de un mantenimiento
preventivo.

Palabras Claves: Proceso, Ecuacin, Ecuacin Diferencial, Mantenimiento preventivo, Disipadores,

Ley de enfriamiento.

Abstract
The experiment aims to demonstrate a process of temperature variation to an object of daily use
such as the computer table. We will focus on the role of the sink which is an electronic device
located inside the PC.

What we want are the results reflect the function as a mathematical model of "Newton's Law of
Cooling" and that this in turn is a comparison with our experimental data; to calculate the
maximum temperature that can reach the motherboard before the ventilation system is turned
on. And based on this we can make a verification process on the quality of the equipment or
programming proper preventive maintenance.

Keywords: Process, Equation, Differential Equation, Preventive Maintenance, sinks, cooling Act.

yeininatalia@gmail.com

Universidad ECCI.

Aplicacin Mtodos numricos

1. INTRODUCION

Este artculo tiene como finalidad ver como un proceso fsico como lo es la variacin de
temperatura se puede reflejar en un mtodo numrico y de ella realizar una base para la
implementacin de un mantenimiento preventivo, este proceso de mejora lo vamos a reflejar
hacia el rendimiento de los disipadores de un PC, utilizando una funcin derivada de la ecuacin
planteado por el gran fsico Newton conocida como la ley de enfriamiento.
2. MARCO TEORICO
Uno de los principales riesgos que se tienen
al adquirir un PC o un porttil es el riego que
se obtiene al daarse el disipador o
ventilador, ya que al adquirirse una
temperatura mayor a la que resiste nuestra
tarjeta madre, placa base o mini board, se
llega a un apagado total, y solo nos llega a
mostrar un mensaje de sobrecalentamiento
de dicha tarjeta, lo que deseamos expresar
con ejemplos tanto prcticos como
matemticos la temperatura mxima que
pueda llegar a adquirir dicha tarjeta para
que no ocurra este tipo de dao al PC,
nuestro modelo matemtico lo vamos a
trabajar con una ecuacin diferencial de
temperatura ms conocida como Ley de
enfriamiento,
plasmada por el Fsico
Newton la cual es:

En donde encontramos
una variacin de la temperatura respecto al
tiempo con una K como constante que
difiere respecto al material que estemos
usando, en nuestro proyecto esta K ira
positiva ya que es un sistema de
calentamiento y no de enfriamiento,
multiplicada por la diferencia entre nuestra
temperatura inicial con nuestra temperatura
del medio, luego realizamos una separacin
de variables y se procede a integrar instante
t=0, la temperatura del cuerpo es T0.

Despus de realizar la integracin,

despejamos la temperatura quedndonos
como resultado , en la que tendremos una
temperatura m, la cual ser la temperatura
del medio en el que nos encontremos,
sumando entre la diferencia entre un
tiempo inicial con un tiempo final el cual
ser reflejado en el momento que el
disipador se encienda para refrigerar
nuestra tarjeta y esta se multiplica por el
exponencial, el cual est elevado a la
constante del material que estamos
midiendo el cual sera la placa

base por el tiempo en cual se realiza el

cambio de la temperatura, todas estas
operaciones nos dan como resultado:
la temperatura
mxima a la cual se puede someter dicha
tarjeta antes de que llegue a tener un
sobrecalentamiento , en nuestra parte
experimental se va a realizar una tabulacin
por medio de una termo cupla y un
multmetro en donde tendremos como
variable dependiente el tiempo que
transcurre antes de que llegue a activarse el
sistema de enfriamiento al igual tendremos
yeininatalia@gmail.com

Universidad ECCI.

Aplicacin Mtodos numricos

nuestra variable dependiente que sera la

temperatura a la cual sube el sistema al
transcurrir el tiempo .

3. METODOLOGIA

- Cronometro del celular.

6. ANALISIS DE RESULTADOS

La funcin que tomamos como base para la

identificacin del proceso fue la de
derivacin de la Ley de enfriamiento o
calentamiento de Newton.

4. HIPOTESIS O RESULTADOS
ESPERADOS:

Para realizar la comparacin del modelo de

mtodos numricos utilizamos la ecuacin
diferencial de "ley de enfriamiento de
Newton" vrs datos experimentales. La
ecuacin de enfriamiento se adapt con las
condiciones del proceso experimental que
fueron:

Temperatura ambiental: 16 C

Temperatura inicial al tiempo cero: 18C

La ecuacin utilizada en el proyecto fue:

Que los datos experimentales vs los

datos reales reflejan un nivel de
similitud en concordancia con los dos
procedimientos.

Determinacin de factores que afectan

de forma directa e indirecta la toma de
los datos, e influyeron en los resultados
del procedimiento.

Llegar a una temperatura exacta a la

cual el procesador no sea afectado.

Llegar a una temperatura final de un

procesador

5. MATERIALES
1 - Multimetro UT33C.

= K (T Tm)
Siendo T (Temperatura), Tm (Temperatura
Ambiental) y K una constante de ritmo.
Reemplazamos los valores de las
condiciones iniciales de las ecuaciones,
realizamos separacin de variables e
integramos. Despejamos T (Temperatura) y
para hallar nuestra C (constante) igualamos
a las condiciones iniciales y creamos la
FUNCIN.

Ln (T-16)= Kt + C

T=C

+ 16

18 = C

+ 16 C=2

1 - Termocupla.
Para hallar valor de C (constante)
reemplazamos en uno de los intervalos
medios del ejercicio y se despeja.

2 - Camara de Video.
7

yeininatalia@gmail.com

Universidad ECCI.

47 = 2

K=

Aplicacin Mtodos numricos

7. FORMULACION DE TABLA
COMPARATIVA
+ 16

= 0.1958

Posteriormente se tabula la FUNCION.

T=C

t
(Minutos)

C
(Experimental)

C
(Teorico)

C
(Taylor)

18

18

18,4

23

19

20,7

28

20

24,1

31

22

30,7

38

26

43,5

10

43

30

69,5

+ 16

Para el mtodo de Taylor se utiliza la funcin

anterior y con las condiciones iniciales del
experimento.
METODO DE TAYLOR

Funcin: T = C

+ 16

Derivada:
75
65

H=2

55
45
35

Yn + Y(Xn)*H

25
15
0

(2*

+ 16) + 2 (

C (Experimental)

C (Teorico)

C (Taylor)

Lineal (C (Taylor))

18,3916

yeininatalia@gmail.com

10

Universidad ECCI.

Aplicacin Mtodos numricos

8. ERRORES ABSOLUTOS

Los errores de la serie de Taylor se

encuentran por encima del valor real.

Experimental

Teoricos

0,4

0,4

2,3

1,7

10.
REFERENCIAS
BIBLIOGRAFICAS

3,9

4,1

0,3

8,7

5,5

17,5

26,5

39,5

Ecuaciones
diferenciales
y
problemas con valores en la frontera
Nagle,saff,Snider- Cuarta edicin.
http://doctoradosociales.com.ar/wp
-content/uploads/Lineamientosgenerales-para-escribir-un-paperDOC-27.pdf

9. CONCLUSIONES DEL CASO

La ecuacin diferencial junto con el

modelo de mtodos numricos nos
mostr valores crecientes en un rango
alto debido a que no se tuvo en cuenta
la participacin del disipador en el
sistema que a diferencia del modelo
experimental llega a una temperatura
constante despus de un punto.

La temperatura crece de manera

acelerada y exponencialmente.

La serie de Taylor con respecto al valor

terico tiene un error relativo ms
grande que el experimental.

yeininatalia@gmail.com