RAICES DE ECUACIONES

KAREN MOLINA MAZABEL

JUAN MANUEL CASTRO CHAVARRO
UNIVERSIDAD SURCOLOMBIANA – INGENIERIA CIVIL Y AGROINDUSTRIAL

1
 Asignatura: Métodos Numéricos

RAICES DE ECUACIONES

Presenta:
Karen Molina Mazabel; Código: 20181169736
Juan Manuel Castro Chavarro; Código: 20191179167

Docente
Ing. YAMIL ARMANDO CERQUERA ROJAS MSc.

Neiva, junio 30 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 Y STEFFENSEN PARA SOLUCIONAR EL PROBLEMA..................................7
4. SOLUCIÓN DEL PROBLEMA:...........................................................................................................................7
5. ANALISIS DE RESULTADO...............................................................................................................................8
6. CODIGO USADO.............................................................................................................................................8
7. CONCLUSIONES.............................................................................................................................................9

3
 1. PLANTEAMIENTO DEL PROBLEMA

Se plantea el siguiente problema:

Una partícula se mueve en el plano XY (la escala está en metros) con una trayectoria descrita
por la función:
u ( t )= ( x ( t ) , y ( t ) )
Teniendo que:
( x ( t ) , y ( t ) )=(2 et t +∜ t∗2 t3 )
Donde:
𝑡 ∈ [0,1]
a) Grafique la trayectoria u(t)
Como u(t) = (x(t), y(t)); por lo tanto:

u(t )=(2e t t+ ∜ t∗2t 3)

b) Encuentre la posición de la partícula cuando t=3.
c) En que instante la partícula se encuentra a una distancia de 4 metros del origen.
La partícula ha recorrido 4 metros, pero no es concreto el instante en que recorrió esta
distancia, por lo tanto, el problema pide que se determine el momento en el que la partícula
realizó su recorrido es decir u(t) = 4.

4=(2 e t t +∜ t∗2 t 3 )
Infiriendo la ecuación en términos del tiempo tendríamos:
f (t)=2 e t t+ ∜ t∗2t 3−4

2. ANÁLISIS DE LA SITUACIÓN PLANTEADA

El movimiento de una partícula se representa por medio de una ecuación y ésta puede ser
graficada para una mejor comprensión de su comportamiento. En este caso en particular la
gráfica se muestra en un plano XY, en el cual el eje X representa el tiempo y el eje Y
representa la posición de la partícula en metros.
A continuación, se muestran un ejemplo del movimiento de una partícula representados en
planos XY:

4
 Figura 1. Ejemplo de movimiento de una partícula en un plano XY

Fuente: “http://centros5.pntic.mec.es/ies.de.melilla/Movmto_rectil.htm ”

En este caso en el instante 0 la partícula se encuentra en el punto cero, y en una primera

instancia su posición sube hasta 3 para luego bajar hasta una posición de -3 y así
sucesivamente.
En el caso del ejercicio planteado, en un instante de tiempo t=0 la partícula se encuentra en la
posición Y=0, y a medida que aumenta el tiempo la posición incrementa teniendo un
comportamiento exponencial.
Haciendo uso del programa Scilab se logra obtener una representación gráfica con la ecuación
ya especificada:

u(t )=(2e t t+ ∜ t∗2t 3)

A. El código utilizado es el siguiente:

t=0:0.1:3;
x=2*t.*exp(t)+t.^ (1/4) *2. *t.^ (3)
plot (t, x); xgrid;

5
 Figura 2. Forma de la ecuación

B. Para determinar el valor de la posición en la que se encontrará la partícula cuando haya

transcurrido 3 horas se utilizó la función F(t):
f (t)=2 e t t+ ∜ t∗2t 3−4

Figura 3: imagen con punto indicando la posición que obtiene la partícula en un tiempo de 3 hr

En la represantacion grafica se indica con un punto la posicion que posee la particula al llevar
un recorrido de 3 horas, la distancia que logra obtener en este lapso de tiempo es de
aproximadamente 191.59 metros.
C. Cuando la particula posee una distancia de 4 metros tarda un cierto tiempo el cual se
6
determina gracias a su ecuacion definida anteriormente:

Figura 3. Imagen con punto señalando la posición de la partícula cuando recorre 4 m.

En la Figura 3 se puede apreciar fácilmente el tiempo que tarda la partícula al recorrer 4

metros, en el eje de las ordenas se especifica el tiempo y en el eje de las abscisas la distancia,
el punto que se denota en la gráfica posee una coordenada de (0.75, 4), esto indica que la
partícula en 4 metros gastó 0.75 horas.

3. MÉTODO DE NEWTON RAPHSON Y MÉTODO DE STEFFENSEN PARA SOLUCIONAR EL

PROBLEMA

Para el modelo de Newton la ecuación es:

f ( x0 )
x 1=x 0−
f ´ ( x0 )

Para el modelo de Steffensen la ecuación es:

f ( x 0 ) x f ( x0 )
x 1=x 0−
f ( x 0 + f ( x 0 ) ) −f ( x0 )

PLANTEAMIENTO DE LA SOLUCIÓN:

Con el objetivo de comparar los resultados obtenidos, el problema se solucionará por los
métodos de Newton Raphson y Steffensen con la siguiente fórmula:
u(t )=(2e t t+ ∜ t∗2t 3)

4. SOLUCIÓN DEL PROBLEMA:

7
Para el punto b: El problema desea que se halle la posición de una partícula en un tiempo
determinado, en este caso de 3 horas. Para esto, se usa la siguiente fórmula:
u(t )=(2e t t+ ∜ t∗2t 3) , al reemplazar t para encontrar la posición de la partícula:
u ( 3 )=( 2 e3 3+ √4 3∗2 t 3) =191,58121 metros
Para el punto c: A diferencia del punto b, este plantea una distancia la cual posee como
incógnita el tiempo que tarda la partícula en recórrela, la distancia propuesta es de 4 metros,
por lo tanto:
Se debe dar un valor inicial para x0, que va ser de 1; porque es la raíz que uno supone y con la
que el método va a empezar a trabajar.
Con este valor se produce la siguiente salida.

Método Newton Raphson:

Digite un valor para x0 : 1
Métodos numéricos
Método de Newton Raphson
x0 f(x0) x1 f(x1)
1.000000 3.436564 0.802193 0.555555
0.802193 0.555555 0.755889 0.024747
0.755889 0.024747 0.753627 0.000057
0.753627 0.000057 0.753622 0.000000
el instante de tiempo es: 0.753621759754 horas

Método Steffensen:

Métodos numéricos
Método de Steffensen
Digite un valor para x0: 1

x0 f(x0) x1 f(x1)
1.000000 3.436564 0.988139 3.232661
0.988139 3.232661 0.974794 3.008380
0.974794 3.008380 0.959691 2.760969
0.959691 2.760969 0.942506 2.487583
8
0.942506 2.487583 0.922873 2.185622
0.922873 2.185622 0.900421 1.853480
0.900421 1.853480 0.874870 1.492085
0.874870 1.492085 0.846287 1.108035
0.846287 1.108035 0.815696 0.719737
0.815696 0.719737 0.786260 0.367375
0.786260 0.367375 0.764238 0.116865
0.764238 0.116865 0.754895 0.013880
0.754895 0.013880 0.753641 0.000210
0.753641 0.000210 0.753622 0.000000

el instante de tiempo cuando la partícula se encuentra a 4 metros del origen será:

0.753621764137 horas

Figura 4: imagen adaptación del ejercicio

La gráfica anterior se adapta al ejercicio planteado ya que los métodos que se utilizaron para
dar solución encuentran una raíz, eso se logra cuando el dominio es igual a cero, por lo tanto
la raíz se encuentra cuando el recorrido pasa por 0 en el dominio, es decir en el eje de las
ordenadas (eje x), para ello la ecuación u(t) con u = 4, se replantea en función del tiempo
como se muestra a continuación:
u(t )=(2e t t+ ∜ t∗2t 3) siendo u(t) = 4

4=(2 e t t +∜ t∗2 t 3 ) por lo tanto:

f (t)=2 e t t+ ∜ t∗2t 3−4

9
 5. ANALISIS DE RESULTADO

Para el punto b: En este punto se evidencia que para un valor del dominio de x=3 (tiempo en
este caso) le corresponde un valor del rango equivalente a 191,58121.

Para el punto c: Se comparan los resultados obtenido en los dos métodos:

NEWTON: 0.753621759754
STEFFENSEN: 0.753621764137

Como es evidente los resultados son muy similares, la diferencia es de unos cuantos
decimales; esto indica que ambos métodos son viables a la hora de desarrollar problemas de
este tipo, su margen de error es mínimo.

6. CODIGO USADO

Para el punto a:
t=0:0.1:3;
y=2*t.*exp(t)+t.^(1/4)*2.*t.^3;
plot(t,y);xgrid;
Para el punto c:
Método Newton Raphson:
function u=ff(t)
u=2*t.*exp(t)+t.^(1/4)*2.*t.^3 - 4
endfunction
h=1e-5;
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("el instante de tiempo es:%15.12lf",x1);
printf(" horas")

10
 Método Steffensen:
function y=fx(t)
y=2*t.*exp(t)+t.^(1/4)*2.*t.^3-4;
endfunction
error=1e-6;
printf("Método numéricos\n Método de steffensen\n");
x0=input("Digite un valor para x0 : ")
printf( " x0 f(x0) x1 f(x1)\n")
while abs(fx(x0))>error
x1=x0-fx(x0)^2/(fx(x0+fx(x0))-fx(x0))
printf("%9.6lf %9.6lf %9.6lf %9.6lf\n",x0,fx(x0),x1,fx(x1));
x0=x1;
end
x2=x1*60,

printf("el instante de tiempo cuando la particula se encuentra a 4 metros del

origen será: %15.12lf ",x1);
printf(" horas")

7. CONCLUSIONES

1. El resultado obtenido en los dos modelos es muy similar, sin embargo, hay una diferencia
muy pequeña entre ellos, pero a final de cuentas la utilización de cualquiera de los dos
métodos es viable para dar solución a este tipo de problemas.
2. Estos métodos son viables cuando se tienen estados específicos de operación, como es el
caso de este; en el cual se desea encontrar el tiempo que ha transcurrido mientras una
partícula se mueve a una posición en específico.
3. Los métodos anteriormente utilizados suelen ser similares a la hora de programar, pero se
debe tener en cuenta que mientras un método trabaja con la derivada de una función en
el caso del método de Newton y para el método de steffesen se utiliza el límite de la
función para dar su respectiva solución.

11