Instituto Tecnologico de Morelia

Lista 1: Introduccion
Analisis Numerico a los Metodos Numericos Agosto-Diciembre 2017
Indicaciones

Para resolver los siguientes ejercicios toma en cuenta las siguientes indicaciones:
Se puede realizar en equipos de a lo mas tres personas. Pero, se sugiere, que la distri-
bucion de los ejercicios sea de modo equitativo teniendo en cuenta que algunos de ellos
llevaran bastante tiempo en su resolucion.
En los ejercicios que se realicen de forma manual debes escribir todo el procedimien-
to; mientras que los ejercicios en los que utilices Matlab debes escribir los comandos
utilizados.
Cada ejercicio debera ir en hojas distintas y, si el ejercicio no se encuentra resuelto
completamente deberas escribir que la solucion es parcial (de otro modo no se tomara
en cuenta).
En total son 238 puntos distribuidos en 23 ejercicios (cada ejercicio muestra su valor),
de los cuales deberas entregar un maximo de 200 puntos que equivaldra al 60 % de la
tarea.
En todos los ejercicios en los que se haya utilizado Matlab debera anexar comentarios
sobre el tiempo que tardo en realizar dicha tarea.
En los ejercicios de numeros felices, narcicistas y perfectos debera anexar un grafico
que muestre la cantidad de numeros que cumplen dicha condicion.
En el ejercicio sobre las ternas de Cardano deberas ubicar todas las ternas en N3 .
Los ejercicios deberas entregarlos a mas tardar el da 25 de septiembre a las 2 : 30pm.
Despues, de la fecha se recibira con una amonestacion de 5 % por cada hora, o fraccion,
de retraso.
Es comprensible que algunos ejercicios se resuelven de forma similar pero si se descubre
que alguno realizo plagio se penalizara con un 10 % (plagiado y plagiador) por ejercicio.
Cualquier duda que se tenga podra resolverse en clase o podra escribir al correo
[email protected] especificando el numero de ejercicio en el asunto.
Cualquier error que se tenga en la redaccion favor de comunicarlo de inmediato en el
mismo correo.

Ultima modificacion: 7 de septiembre de 2017

 Instituto Tecnologico de Morelia
Lista 1: Introduccion
Analisis Numerico a los Metodos Numericos Agosto-Diciembre 2017

Nombre: Matrcula:
Nombre: Matrcula:
Nombre: Matrcula:

1. (Valor 10 puntos) Sea n un numero natural 1 . Sea dk . . . d2 d1 d0 su expan-

sion decimal; es decir,

n = dk 10k + + d2 102 + d1 10 + d0

con d0 , d1 , d2 , . . . , dk {0, 1, . . . , 9} y dk 6= 0.
Definamos la operacion : N N de la siguiente forma:

(n) = d2k + + d22 + d21 + d20 .

Por ejemplo, (367) = 32 + 62 + 72 = 94.

Diremos que un numero natural n es feliz si existe un numero natural m tal
que m (n) = 1. En caso contrario, diremos que el numero n no es feliz.
Por ejemplo, el numero natural 31 es feliz pues, 2 (31) = ((31)) =
(10) = 1.
a) Determina, de forma manual, si los numeros 4 y 91 son felices.
b) Elabora un script en Matlab para determinar si un numero natural
n es feliz.
c) Determina cuantos numeros felices existen que sean menores que
1000000.
d ) Determina cuantos numeros felices existen que sean menores que
1000000 y que terminen en 1.
e) Calcula la suma de todos los numeros felices menores que 1000000
que sean primos.
2. (Valor 10 puntos) Sea n un numero natural. Sea dk . . . d2 d1 d0 su expansion
decimal; es decir,

n = dk 10k + + d2 102 + d1 10 + d0

con d0 , d1 , d2 , . . . , dk {0, 1, . . . , 9} y dk 6= 0.
Definamos la operacion : N N de la siguiente forma:

(n) = dkk + + dk2 + dk1 + dk0 .

Por ejemplo, (153) = 13 + 53 + 33 = 153.

Diremos que un numero natural n es narcicista si (n) = n (153 es nar-
cicista).
1 En esta tarea consideraremos al conjunto de los numeros naturales N como {1, 2, 3, 4, . . . }.

Ultima modificacion: 7 de septiembre de 2017

 Instituto Tecnologico de Morelia
Lista 1: Introduccion
Analisis Numerico a los Metodos Numericos Agosto-Diciembre 2017

a) Elabora un script en Matlab para determinar si un numero natural

n es narcicista.
b) Completa la siguiente tabla.
Numero Cantidad de numeros Suma de los numeros
de cifras narcicistas narcicistas de tales cifras
1
2
3
4
5
6
7
8
9
10
11
12

3. (Valor 15 puntos) Sean a, b, c numeros naturales. Decimos que la terna

(a, b, c) es de Cardano si

q q
3 3
a + b c + a b c = 1.

Por ejemplo, (2, 1, 5) es una terna de Cardano. Determina todas las ternas
de Cardano de modo que a + b + c = 1000.
4. (Valor 8 puntos) Sea {fn } la sucesion de Fibonacci; esto es,

f1 = 1
f2 = 1
fn+2 = fn+1 + fn para cada n 1

fn+1
Definamos la sucesion an = fn con n 1.

a) Demuestra que la sucesion converge.

b) Determina el valor de convergencia con 10 cifras significativas.
5. (Valor 9 puntos) Elabora un script en Matlab que convierta cualquier
numero racional positivo a sistema binario. (Sugerencia: Toma en cuenta
la periodicidad).
6. (Valor 20 puntos) Usando el script anterior, responde lo siguiente:
a) Cuantos dgitos son necesarios para representar n cifras significati-
vas del numero ?

Ultima modificacion: 7 de septiembre de 2017

 Instituto Tecnologico de Morelia
Lista 1: Introduccion
Analisis Numerico a los Metodos Numericos Agosto-Diciembre 2017

n Cantidad de dgitos
1
2
3
4
5
6
7
8
9
10
b) (Valor 5 puntos) Convierte a sistema binario cada uno de los siguien-
tes numeros.

1) 451 4) 13
2) 78912 5) 18.25
3) 13
16 6) 101.0101

7. (Valor 3 puntos) Tenemos que 11.0010010001dos . Halle el error de

dicha aproximacion si consideramos a = 3.14159265358979 . . .
8. (Valor 5 puntos) Use la tabla con mantisa de 4 cifras y n {3, 2, 1, 0, 1, 2, 3}
para calcular

3 1 1 3 1 1



a) 17 + 9 + 7 b) 17 + 9 + 7

Calcula el error relativo verdadero en cada uno de los incisos anteriores y

comenta a que se deben las diferencias en caso de que las hubiera.
9. (Valor 2 puntos) En cada uno de los siguientes casos, halle el error absoluto
y el error relativo; ademas, determina el numero de cifras significativas de
la aproximacion.
a) x = 2.71828182 y x = 2.7182
b) y = 98350 y y = 98000
c) z = 0.000068 y z = 0.00006
d ) z = 0.000068 y z = 0.00007
10. (Valor 6 puntos) Determina el tipo de error que se presenta si
1/4 1/4
x4 x6
Z Z  
x2 2
e dx 1+x + + dx = p
0 0 2! 3!

y compare su resultado con el valor de p = 0.2553074606.

Ultima modificacion: 7 de septiembre de 2017

 Instituto Tecnologico de Morelia
Lista 1: Introduccion
Analisis Numerico a los Metodos Numericos Agosto-Diciembre 2017

11. (Valor 5 puntos) Determina el resultado adecuado de la suma p + q y el

producto pq si consideramos los datos
a) p = 1.414 y q = 0.09125, que vienen dados con una precision de
cuatro cifras significativas.
b) p = 31.415 y q = 0.027182, que vienen dados con una precision de
cinco cifras significativas.
12. (Valor 10 puntos) Encuentre, en cada uno de los siguientes casos, una
formula equivalente a la dada que evite la perdida de cifras significativas.
Justifica.
a) ln(x + 1) ln(x) para x grande.

b) x2 + 1 x para x grande.
c) cos2 x sen2 x para x 4 .
q
1+cos x
d) 2 para x .

13. (Valor 7 puntos) Usando la aritmetica de punto flotante con tres cifras y
corte, calcule las siguientes sumas (sumando en el orden indicado).
6
1
P
a) 3k
k=1
6
1
P
b) 37k
k=1

14. (Valor 10 puntos) Sean p(x) = x3 3x2 + 3x 1, q(x) = ((x 3)x + 3)x 1
y r(x) = (x 1)3 .
a) Usando aritmetica de punto flotante con cuatro cifras significativas
y redondeo, calcule p(2.72), q(2.72) y r(2.72). En el calculo de p(x),
suponga que (2.72)3 = 20.12 y (2.72)2 = 7.398.
b) Usando aritmetica de punto flotante con cuatro cifras significativas y
redondeo, calcule p(0.975), q(0.975) y r(0.975). En el calculo de p(x),
suponga que (0.975)3 = 0.9268 y (0.975)2 = 0.9506.
Ademas, calcula el error relativo verdadero en cada uno de los incisos.
15. (Valor 5 puntos) La formula mejorada para la resolucion de la ecuacion de
segundo grado. Supongamos que a 6= 0 y que b2 4ac > 0 y consideremos la
ecuacion cuadratica ax2 + bx + c = 0. Sus races pueden hallarse mediante
la conocida formula

b + b2 4ac b b2 4ac
x1 = y x2 = (1)
2a 2a
Pruebe que estas races pueden calcularse mediante las formulas equiva-
lentes
2c 2c
x1 = y x2 = (2)
b + b2 4ac b b2 4ac

Ultima modificacion: 7 de septiembre de 2017

 Instituto Tecnologico de Morelia
Lista 1: Introduccion
Analisis Numerico a los Metodos Numericos Agosto-Diciembre 2017

Nota: Cuando |b| b2 4ac, hay que proceder con cuidado para evitar
la perdida de precision por cancelacion. Si b > 0, entonces x1 debera ser
calculado con la formula (2) y x2 debera ser calculado con la formula (1);
mientras que, si b < 0, entonces x1 debera ser calculado con la formula
(1) y x2 debera ser calculado con la formula (2).
Explica el porque de la nota anterior.
16. (Valor 10 puntos) Elabora un script en Matlab que calcule las races de
una ecuacion de segundo grado de dos formas: usando la formula general
y usando la formula general mejorada.
17. (Valor 8 puntos) Usando el script enterior, calcula las races de las siguien-
tes ecuaciones de segundo grado.
a) x2 1000.001x + 1 = 0
b) x2 10000.0001x + 1 = 0
c) x2 100000.00001x + 1 = 0
d ) x2 1000000.000001x + 1 = 0
Comenta si hubo diferencia en los resultados obtenidos. En caso de que
haya habido diferencia calcula el error relativo porcentual.
18. (Valor 8 puntos) Convierte los siguientes numeros a sistema decimal.

a) (0.11110101)dos d ) 11.01011011100001dos
b) 410111cinco e) 0.11001dos
c) 11101.11011dos f ) 0.111011dos

19. (Valor 15 puntos) Una forma de calcular el valor de es el algoritmo

de Gauss-Legendre. Dicho algoritmo toma como datos iniciales a0 = 1,
b0 = 12 = 0.7071067812, t0 = 14 = 0.25 y p0 = 1. A partir de ellos
realizamos las siguientes operaciones:
an + bn
an+1 =
2
p
bn+1 = an bn
tn+1 = tn pn (an an+1 )2
pn+1 = 2pn
Entonces
(an + bn )2

4tn
a) Realiza un script en matlab que aproxime el valor de pidiendo el
numero de iteraciones necesarias n.

Ultima modificacion: 7 de septiembre de 2017

 Instituto Tecnologico de Morelia
Lista 1: Introduccion
Analisis Numerico a los Metodos Numericos Agosto-Diciembre 2017

b) Usa el script para estimar el valor de pi para n = 1, 2, 3, 4, 5. (Use

format long para tener mayor cantidad de cifras significativas.
c) Calcule el error relativo porcentual verdadero, en cada caso, conside-
rando = 3.141592653589793.

20. (Valor 5 puntos) Evalue e5 con el uso de dos metodos

x2 x3 x4
ex = 1 x + + ...
2 3! 4!
y
1 1
ex = = x2 3 x4
ex 1+x+ 2 + x3! + 4! + ...
y comparelo con el valor verdadero 6.737947 103 . Calcule el error por-
centual aproximado y el error porcentual verdadero.
21. (Valor 20 puntos) Un numero natural n es perfecto si la suma de sus
divisores propios (i.e. menores a el) da por resultado n. Por ejemplo, 6 es
un numero perfecto ya que sus divisores propios son 1, 2 y 3; y, 6 = 1+2+3.
a) Elabora un script en Matlab que calcule los primeros m numeros
perfectos; es decir, que pida un valor m e imprima dicha cantidad de
numeros perfectos.
b) Haciendo uso del script anterior encuentra los primeros 8 numeros
perfectos.
22. (Valor 12 puntos) El numero e, tambien conocido como numero de Euler
o constante de Napier, es uno de los mas importantes numeros reales que
aparece en diversas areas de las matematicas. Existen diversas formas de
poder calcular dicho numero y una de ellas es

X k3
e=
5(k!)
k=1

a) Realiza un script en Matlab que realice una aproximacion de e hasta

que el error porcentual verdadero sea menor a una tolerancia dada.
b) Hasta que valor de k es necesario realizar dicha suma de modo que
el error porcentual verdadero sea menor a 0.00001 %?

Considere a e = 2.71828182845945.
23. (Valor 30 puntos) Elabora un script en Matlab en el que se especifiquen
el numero de cifras significativas de un numero, y cuando se realice una
operacion aritmetica arroje el resultado con tal restriccion.

Ultima modificacion: 7 de septiembre de 2017