RAICES DE ECUACIONES
MARIA JOSE LOAIZA ANDRADE
MARIA ALEXANDRA PERDOMO TOVAR
UNIVERSIDAD SURCOLOMBIANA – INGENIERIA CIVIL
1
Asignatura: Métodos Numéricos
RAICES DE ECUACIONES
Presenta:
Maria Jose Loaiza - Código: 20182171150
Maria Alexandra Perdomo - Código: 20182171109
Docente
Ing. YAMIL ARMANDO CERQUERA ROJAS MSc.
Neiva, agosto 05 de 2020
2
Tabla de contenido
PRIMER PROBLEMA.
1. PLANTEAMIENTO DEL PROBLEMA.................................................................................................................4
2. ANÁLISIS DE LA SITUACIÓN PLANTEADA.......................................................................................................4
3. MÉTODO DE NEWTON RAPHSON PARA SOLUCIONAR EL PROBLEMA...........................................................7
4. SOLUCIÓN DEL PROBLEMA:...........................................................................................................................8
5. ANALISIS DE RESULTADO...............................................................................................................................8
6. CODIGO USADO.............................................................................................................................................9
7. CONCLUSIONES...........................................................................................................................................10
3
1. PLANTEAMIENTO DEL PROBLEMA
Se plantea el siguiente problema:
Un modelo de crecimiento poblacional está dado por:
𝑓(𝑡) = 5𝑡 + 2𝑒 ˆ0.1𝑡
en donde n es el número de habitantes, 𝑡 es tiempo en años.
a) Calcule el número de habitantes que habrá en el año 25
b) Encuentre el tiempo para el cual la población es 200
2. ANÁLISIS DE LA SITUACIÓN PLANTEADA
La medición del crecimiento poblacional es muy importante porque con ella se conocen cifras
las cuales contribuyen en la planeación de lo que se necesita o se va a necesitar (por ejemplo,
cuántas escuelas, hospitales, transportes o fuentes de trabajo requiere un lugar), por eso se
realizan censos periódicos en cada país. Además, esta medición ayuda a predecir el
crecimiento de la población esto es de suma importancia ya que, con estos análisis
demográficos y organización urbanística, se tiene una idea de las causas que generan
el crecimiento o decrecimiento de la población, tomando en cuenta factores como lo son:
Tasas de nacimiento.
Figura 1. Crecimiento poblacional
Fuente: [Link]
despues-de-tres-siglos-1276640994/”
4
Como se puede apreciar en el planteamiento del problema, es decir, lo que se pregunta en el
problema, es el número de habitantes que habrá en el año 25 (n) y tiempo para que la
población sea igual a 200. Esto se puede interpretar de la siguiente manera:
Primer momento:
El tiempo esta dado desde cero al igual que el número de habitantes.
Segundo momento:
El crecimiento de la población cumple una función ascendente, al paso de los años el numero
de habitantes crece, se podría decir que su crecimiento es exponencial.
Tercer momento:
Cuando haya pasado 25 años se desea saber el número de habitantes (n) con los que cuenta,
además de eso se desea evaluar el tiempo (t) necesario para que la población sea igual a 200
personas.
Haciendo uso del programa Scilab se puede obtener una idea del aspecto de dicha gráfica:
f ( t )=5∗t+ 2 e0.1∗t
El código utilizado es el siguiente:
t=[Link];
y=5*t+(2*exp(0.1*t));
plot(t,y);xgrid
5
Figura 2. Forma de la ecuación
B
A
A: 25 años
B: 200 habitantes
Para hallar t cuando n(t)= 200, entonces la función f(t) será igual a:
f ( x )=5∗t +2(e 0.1∗t )−200
De la misma manera Haciendo uso del programa Scilab se puede obtener una idea del aspecto
de la anterior gráfica:
6
Figura 3. Imagen con el punto como raíz
En la Figura 3 se puede apreciar fácilmente el intervalo en el que probablemente se encuentre
la raíz [13 a 14], dicho intercepto sería el valor que toma t cuando la función n(t)=200
3. MÉTODO DE NEWTON RAPHSON PARA SOLUCIONAR EL PROBLEMA
Para el modelo de Newton la ecuación es:
f ( x0 )
x 1=x 0−
f ´ ( x0 )
f ( x +h )−f ( x)
f ' ( x )=lim
h→ 0 h
PLANTEAMIENTO DE LA SOLUCIÓN:
El problema se solucionará por el método de Newton Raphson teniendo como ecuación lo
siguiente:
A. Para hallar el numero de habitantes a los 25 años se usa esta ecuación:
f ( x )=5∗25+2( e0.1∗25)−t
B. Para hallar “t” cuando los habitantes sean igual a 200 se usa esta ecuación:
0.1∗t
f ( x )=5∗t +2(e )−200
7
4. SOLUCIÓN DEL PROBLEMA:
A. Para iniciar con la aplicación para encontrar el numero de habitantes que habrá a los 25 años es
necesario dar un valor inicial tal como lo exige, este valor será 25.
Con este valor se produce la siguiente salida.
Digite un valor para x0 :25
Metodos numericos
Metodo de Newton Raphson
x0 f(x0) x1 f(x1)
25.000000 124.364988 149.364988 -0.000000
La raiz es : 149.364988235398
B. Para iniciar con la aplicación del modelo es necesario dar un valor inicial tal como lo exige, este valor
será 15, de esta manera se encontrara el tiempo que es necesario para que los habitantes sean 200.
Con este valor se produce la siguiente salida.
Digite un valor para x0 :15
Metodos numericos
Metodo de Newton Raphson
x0 f(x0) x1 f(x1)
15.000000 -116.036622 34.679439 37.538664
34.679439 37.538664 31.390655 3.117849
31.390655 3.117849 31.066435 0.024004
31.066435 0.024004 31.063900 0.000001
31.063900 0.000001 31.063900 0.000000
La raiz es : 31.063899925758
5. ANALISIS DE RESULTADO
Si bien es cierto con el modelo se ha encontrado una raíz “Aparente”, lo real es
que al tiempo de 25 años habrán 149.364988235398 habitantes, además se
encontró que en t=31.063899925758 la población corresponderá a un total de
200 personas.
8
Raíz aparente
B
A
200
6. CODIGO USADO
A. CODIGO PARA HALLAR LA RAIZ APARENTE A 25 AÑOS:
function y=ff(n)
y=(5*25+2*exp(0.1*25))-n;
endfunction
h=1e-6; 1*10^(-6)
clc
x0 = input('Digite un valor para x0 :')
printf("Metodos numericos\nMetodo de Newton Raphson\n");
printf(" x0 f(x0) x1 f(x1)\n")
while abs(ff(x0))>h
x1=x0-ff(x0)*h/(ff(x0+h)-ff(x0))
printf("%9.6lf %9.6lf %9.6lf %9.6lf \n",x0,ff(x0),x1,ff(x1));
x0=x1;
end
printf("La raiz es : %15.12lf",x1);
9
B. CODIGO PARA HALLAR LA RAIZ APARENTE CUANDO LOS HABITANTES SEAN 200:
function y=ff(t)
y=(5*t+2*exp(0.1*t))-200;
endfunction
h=1e-6; 1*10^(-6)
clc
x0 = input('Digite un valor para x0 :')
printf("Metodos numericos\nMetodo de Newton Raphson\n");
printf(" x0 f(x0) x1 f(x1)\n")
while abs(ff(x0))>h
x1=x0-ff(x0)*h/(ff(x0+h)-ff(x0))
printf("%9.6lf %9.6lf %9.6lf %9.6lf \n",x0,ff(x0),x1,ff(x1));
x0=x1;
end
printf("La raiz es : %15.12lf",x1);
7. CONCLUSIONES
Mediante el método de Newton Rhason y con ayuda de la aplicación scilab se
pudo hallar de una forma muy fácil la raíz aparente, en la cual los datos con un
determinado tiempo tienen una determinada cantidad de habitantes, además se
concluye que con el numero de iteraciones utilizados en el proceso descrito
consiguió clasificar la población en grupos estables.
El análisis numérico que se utiliza en el método de Newton Rhason con el cual fue
resuelto el problema en compañía de un algoritmo que facilito la manera de
encontrar la respectiva solución de forma efectiva y precisa la raíz, el método fue
igualar la función a cero y aproximar de manera sucesiva, para así partir de un
valor inicial.
10
