BANCO DE PREGUNTAS MÉTODOS NUMÉRICO

Pregunta 1:
1. La concentración x de una bacteria contaminante en un lago decrece según la expresión
−2t −0 ,5 t
x (t )=80 e +20 e siendo t el tiempo en horas. Determinar el tiempo que se necesita para
que el número de bacterias se reduzca 7.

Solución:
Por el método de Newton Raphson

t0 t1
1.100000 1.670399
1.670399 2.121097
2.121097 2.307707
2.307707 2.328861
2.328861 2.329088

El tiempo que se necesita para que se reduzca a 7 el número de bacterias es t= 2.32909 horas
con un margen de error de ε =0.005

Pregunta 2:

2.- Demostrar el método de secante

Demostración:

Sean P0=( x 0 , f ( x 0) ) , P1=( x1 , f ( x 1 ) ) y P 2=( x 2 , f ( x 2) ) puntos próximos a la solución

f ( x 0 ) −f ( x1 )
Pendiente entre los puntos P0 y P1 es:m= (1)
x 0−x 1

f ( x 1 )−f (x 2 )
Pendiente entre los puntos P1 y P2 es:m= (2)
x 1−x 2

Igualando pendiente (1) y (2) se tiene:

f ( x 0 )−f (x 1 ) f ( x 1 ) −f (x 2 )
=
x 0−x 1 x 1−x 2
Como la recta intersecta al eje x entonces la imagen en el punto x 2 es igual a cero es decir:
f ( x 2 ) =0

Por tanto:

f ( x 0 )−f (x 1 ) f ( x1)
=
x 0−x 1 x 1−x 2

Despejando x 2 se tiene:

x0 f ( x 1 ) −x1 f (x 0)
x 2=
f ( x 1 )−f ( x 0 )

Generalizando la fórmula se tiene:

x n−1 f ( x n )−x n f (x n−1 )

x n+1=
f ( x n ) −f (x n−1 )

Pregunta 3:

3.- Demostrar el método de Newton Raphson

Demostración:

Sea f (x) una función general:

La ecuación de recta tangente entre los puntos ( x 0 , f (x 0 )¿ y (x 1 , f (x 1)) dado por:

( y−f ( x 0 ) )=m(x−x 0 )

Como la pendiente es la derivada de la curva o ecuación no lineal en x 0 se tiene: f ' ( x0 )

'
( y−f ( x 0 ) )=f (x 0 )(x−x 0 )
Despejando la variable x se tiene:
f ( x0 )
x=x 0− '
f ( x0 )

Pregunta 4:

4.- Resolver por el método de punto fijo la siguiente ecuación:

−x
sen ( x ) e −1=0

Solución:

sin (x ) exp (x )-1

20

-20

-40

-60

-80

-100

-120

-140

-160

-180

-3 -2 -1 0 1 2 3 4 5 6
x

De la gráfica sea x 0=0 punto inicial o partida.

Despejando la variable x de la ecuación no lineal se tiene:

−1
x=si n (exp ⁡(x))

X0 X1
0.000000 1.570796
1.570796 1.570796
1.570796 -0.671664
-0.671664 -0.307546
-0.307546 0.687631
0.687631 1.159083
1.159083 0.226780
0.226780 -0.237983
-0.237983 0.332233
0.332233 0.774119
0.774119 0.506934
0.506934 0.066379
0.066379 0.213348
0.213348 0.523535
0.523535 0.525511
0.525511 0.273133
0.273133 0.228195
0.228195 0.393307
0.393307 0.470195
0.470195 0.365349
0.365349 0.283888
0.283888 0.343027
0.343027 0.415021
0.415021 0.390288
0.390288 0.330137
0.330137 0.335279
0.335279 0.379785
0.379785 0.387090
0.387090 0.355781
0.355781 0.343229
0.343229 0.363174
0.363174 0.377099
0.377099 0.365772
0.365772 0.352831
0.352831 0.358146
0.358146 0.368853
0.368853 0.367513
0.367513 0.359309
0.359309 0.358378
0.358378 0.364168
0.364168 0.366224
0.366224 0.362455

La solución es: 0.362455

Pregunta 5:

5.- Diseñar el código del método de punto fijo

Solución:

syms x;
fx=input('lea la ecuación no lineal f(x):','s');
gx=input('lea la ecuación no lineal g(x):','s');
x0=input('lea el punto inicial:');
e=input('lea el margen de error:');
tol=1;
fprintf('\n %7s %7s \n', 'X0','X1');
while tol>e
x=x0;
x1=eval(gx);
tol=abs(x1-x0);
fprintf('\n %7f %7f \n', x0,x1);
if abs(x0-x1)<e
break
else
x0=x1;
end
end
fprintf('\n la solución es:%7f\n',x1);
ezplot(fx)
grid on

Pregunta 6:

6.- Resolver el sistema AX=B un sistema de ecuaciones lineales aplicando el código del
método de L .U .
 x+ 3 y−z ¿4
x+ 2 y + 4 z ¿ 2
3 x+ 7 y−2 z ¿ 5
Solución:

Sea:

[ ] [] [ ]
1 3 −1 x 2
A= 1 2 4 X= y B= 0
3 7 −2 z −2

Aplicando el programa L.U. se tiene:

Lea A : [1 3 -1;1 2 4;3 7 -2]↲

Lea B:[2 0 -2]↲

x= -10.2222
y= 4.2222
z= 0.4444

Pregunta 7:

7.- Hallar el polinomio de interpolación y su gráfica de los siguientes datos a través del
método de Diferencias Divididas de Newton

X 2 4 5
Y -1 3 2
Solución:

El polinomio de interpolación es:

2
p ( x )=−x + 8 x −13
 8 x - x 2 - 13

-20

-40

-60

-80

-100

-6 -4 -2 0 2 4 6
x

Pregunta 8

Diseñar código para hallar el área de la región acotada por y−x 2=0; x 2−3=0 por trapecios

Solución:

Sea la gráfica de la región acotada:

∫ f ( x ) dx= h2 [ y 0 +2 y1 +2 y 2 +…+ y 0 ]
a

Hallando a , b

−√ 3 √3 por tanto a= −√3 ; b= √3

x 1= ; x 2=
√2 √2 √2 √2
syms x;
fx= x 2;
gx=3- x 2 ;

a=
−√ 3 ;
√2
b=
√ 3
√2
n=input('lea número de particiones n=');
h=(b-a)/n;
int=0;
for i=1:n-1
xn=i*h;
x=xn;
fxn=eval(fx);
gxn=eval(gx);
img=abs(fxn-gxn);
int=int+img;
I=(h/2)*(a+(2*img)+b);
end
fprintf('\n sol=%7f\n',I)

Pregunta 9:

dy
Resolver la siguiente ecuación diferencial =xy con valor inicial y ( 0 )=1 aproximar el
dx
valor de y (1), con h=0.5 a través del método de Runge-Kutta de cuarto orden.

Solución:
El valor inicial y ( 0 )=1 significa x 0=0 , y 0=1 , además h=0.5, f ( x , y )=xy

Primera Iteración:

1
y 1= y 0 + ( k 1 +2 k 2 +2 k 3 +k 4 ) ,se halla ( k 1 , k 2 , k 3 y k 4 )
6

k 1=hf ( x 0 , y 0 ) como f (x 0 , y 0 )=0

k 1=0.5 ( 0 )=0

( 1
2
1
)
k 2=hf x 0 + h , y 0 + k 1 Entonces:
2
1
2 ( 1
k 2=0.5 (0+ 0.5)(1+ 0)
2 )
k 2=0.5
1 1
=
4 8 ()
( 1
2
1
)
k 3=hf x 0 + h , y 0 + k 2 Entonces:
2

( )
1
∗1
1 2
k 3=0.5 (0+ 0.5)(1+ )
2 8
 ( )
1
∗17
4 17
k 3=0.5 =
16 128

k 4=hf ( x 0 +h , y 0 +k 3 ) Entonces:
k 4=0.5 ¿
k 4=0.5( 0.5∗145 145
128
= )
512
Portanto:
y 1=1+
1
6(0+
2∗1 2∗17 145 3481
8
+ + =
128 512 3072 )
Si n=1, x 1=x 0 +h=0+0.5=0.5

Segunda Iteración:

1
y 2= y 1 + ( k 1 +2 k 2 +2 k 3 +k 4 )hallando( k 1 , k 2 , k 3 y k 4 ) se tiene:
6

k 1=hf ( x 1 , y 1 )se tiene:

k 1=0.5 ( 0.5∗3481
3072 )=
3481
12288

( 1
2
1
)
k 2=hf x 1 + h , y 1+ k 1 se tiene:
2

(( )( ))
1
∗3481
1 3481 2
k 2=0.5 0.5+ 0.5 ∗ + =0.4780
2 3072 12288

( 1
2
1
)
k 3=hf x 1 + h , y 1+ k 2 se tiene:
2

((
1
k 3=0.5 0.5+ 0.5
2
3481 1
)(
+ ∗0.4780 =0.5146
3072 2 ))
k 4=hf ( x 1 +h , y 1+ k 3 )se tiene:

(
k 4=0.5∗ ( 0.5+0.5 )∗ ( 3481
3072 )
+0.5146 ) =0.8239
Portanto:
1
y 2= y 1 + ( k 1 +2 k 2 +2 k 3 +k 4 )
6

y 2= + (
3481 1 3481
3072 6 12288
+ 2∗0.4780+ 2∗0.5146 +0.8239 =1.6485 )
Si n=2, x 2=x 1+ h=0.5+ 0.5=1

Tercera Iteración:
 1
y 3= y2 + ( k 1 +2 k 2 +2 k 3 +k 4 )hallando( k 1 , k 2 , k 3 y k 4 ) se tiene:
6

k 1=hf ( x 2 , y 2 )se tiene:

k 1=0.5 ( 1∗1.6485 ) =0.8243

( 1
2
1
)
k 2=hf x 2 + h , y 2+ k 1 se tiene:
2
1
2 (( 1
)(
k 2=0.5 1+ 0.5 ∗ 1.6485+ ∗0.8243 =1.2879
2 ))
( 1
2
1
)
k 3=hf x 2 + h , y 2+ k 2 se tiene:
2

((
1
2 )( 1
k 3=0.5∗ 1+ 0.5 1.6485+ ∗1.2879 =1.4328
2 ))
k 4=hf ( x 2 +h , y 2 +k 3 )se tiene:
k 4=0.5∗( ( 1+0.5 )∗( 1.6485+1.4328 ) )=2.2310
Por tanto:
1
y 3= y2 + ( k 1 +2 k 2 +2 k 3 +k 4 )
6

1
y 3=1.6485+ ( 0.8243+2∗1.2879+2∗1.4328+ 2.2310 )=3.0646
6

Pregunta 10:

Demostrar el método de Euler

y i= y i−1 +hf ( xi , y i)i=1 ,2 , 3 , … n

Demostración:

El proceso para resolver mediante éste método consiste en encontrar ecuaciones de rectas
tangentes en cada intervalo, los intervalos son construidos a partir de puntos espaciados
x n−x 0
separados a distancias iguales por la siguiente fórmula: x i=x i−1+ h ; donde h=
n

f(x1,
x1)

f(x0, x0)

h h

x0 X1 X2
La primera ecuación de recta tangente que aproxima a la solución de la ecuación
dy
diferencial =f ( x , y )en un punto inicial y ( x 0 )= y 0 es:
dx
'
y ( x )− y ( x 0 )= y (x)(x −x 0)

' dy '
Como y ( x )= =f ( x , y ) en el punto y ( x0 ) = y 0 la pendiente es y ( x )=f (x 0 , y 0 ) por tanto
dx
'
reemplazando éste último valor en y ( x )− y ( x 0 )= y (x)(x −x 0)se tiene:

y ( x )− y ( x 0 )=f (x 0 , y 0 )( x−x 0 )

Primera ecuación de recta tangente en el punto ( x 0 , y 0 )

y ( x )= y ( x 0 )+ f ( x0 , y 0 )( x−x 0)

 La primera aproximación a la solución para el punto x=x 1=x 0 +h es:

y 1= y 0 + f (x 0 , y 0)(h)
'
Para el punto ( x 1 , y 1) la pendiente es: y ( x )=f ( x 1 , y 1) y la ecuación de recta tangente es:

y ( x )= y ( x1 ) + f (x 1 , y 1 )(x−x 1)

 La segunda aproximación a la solución para el punto x=x 2=x 1 +h es:

y 2= y 1 + f (x 1 , y 1 )(h)
'
Para el punto ( x 2 , y 2) la pendiente es: y ( x )=f ( x 2 , y 2) y la ecuación de recta tangente es:

y ( x )= y ( x 2) + f (x 2 , y 2)(x−x 2)

 La tercera aproximación a la solución para el punto x=x 3=x 2 +h es:

y 3= y2 + f ( x 2 , y 2 ) ( h )

Así sucesivamente se tiene:

y i= y i−1 +hf ( x i−1 , y i−1 )