UNIVERSIDAD TCNICA DE MANAB

INSTITUTO DE CIENCIAS BSICAS

DEPARTAMENTO DE MATEMTICAS Y ESTADSTICA

POTAFOLIO DE MATODO NUMERICO

CUARTO SEMESTRE

PARALELOK

NOMBRE DEL ESTUDIANTE:

JOSE ANGEL ANDRADE TOALA

PROFESOR:
ING. YANDRI GUERRERO

PORTOVIEJO, JUNIO DE 2017

PERODO ACADMICO:
MAYO 2017-SEPTIEMBRE 2017

1
 NDICE

CONTENIDOS Pgs.

CARTA DE PRESENTACION.22

MISION Y VISION.3

SYLLABUS...5

MATRIZ DE SEGUIMIENTO AL SYLLABUS13

ACTA DE ACUERDO DE COMPROMISOS ETICOS Y

SOCIALIZACION DEL SYLLABO15

AUTORETRATO DEL ESTUDIANTE...16

EVIDENCIAS DE ACTIVIDADES DE APRENDIZAJE...17

EVALUACION DE PROCESO DE APRENDIZAJE..45

2
 CARTA DE PRESENTACIN

Este portafolio presenta mi trayectoria en el curso de:

METODO NUMERICO, este curso tiene como objetivos desarrollar las destrezas de
anlisis, razonamiento y comunicacin durante este medio ciclo pude conocer sobre varios
procedimientos matemticos como: porcentaje de error, grafica de funciones, races
mltiples, matrices, mtodo de newton, sistemas de ecuaciones, falsa posicin entre otros
mtodos realizas en el software Matlab y en Excel aplicndolas en el campo de la
ingeniera. Las tcnicas presentadas por el profesor me ayudaron a mejorar como futuro
profesional de la Ingeniera qumica.

Los tpicos ms dificultosos del curso fueron: el mtodo de falsa posicin, sistemas de
ecuaciones y mtodo de gauss.

MISION Y VISION DE LA UNIVERSIDAD TECNICA DE MANABI

Misin:

Formar acadmicos, cientficos y profesionales responsables, humanistas, ticos y

solidarios, comprometidos con los objetivos del desarrollo nacional, que contribuyan a la
solucin de los problemas del pas como universidad de docencia con investigacin,
capaces de generar y aplicar nuevos conocimientos, fomentando la promocin y difusin
de los saberes y las culturas, previstos en la Constitucin de la Repblica del Ecuador.

Visin:

Ser institucin universitaria, lder y referente de la educacin superior en el Ecuador,

promoviendo la creacin, desarrollo, transmisin y difusin de la ciencia, la tcnica y la
cultura, con reconocimiento social y proyeccin regional y mundial.

MISION Y VISION DEL INSTITUTO DE CIENCIAS BASICAS

Misin:

Ser una unidad acadmica que brinda una informacin de excelencia en el campo de las
ciencias matemticas, fsicas y qumicas a profesionales y acadmicos, a travs de la
transicin, creacin y difusin del conocimiento, que propendan y conduzca procesos de

3
cambios e innovacin en su desempeo profesional para contribuir al desarrollo y buen
vivir del pas.

Visin:

Ser una unidad acadmica reconocida por su excelencia, lder en la formacin de

profesionales y acadmicos con conocimientos de las ciencias matemticas, fsicas y
qumicas, en igualdad de oportunidades, con responsabilidad social y comprometida con
el desarrollo cientfico y tecnolgico de la sociedad brindando una educacin de calidad y
calidez, vinculada a la investigacin y la oferta de servicio.

Mision y Visin del departamento de Matemticas y Estadsticas

Misin:

Desarrollar herramientas que permitan transformar e interpretar la realidad del hombre

con el fin de solucionar problemas.

Al formar con el lenguaje de la matemtica se pretende crear un sistema de modelacin

del pensamiento estratgico, crtico y creativo a travs del cual se puede abordar,
analizar, argumentar racional y humanamente el mundo capacitndolo laboralmente y
profesionalmente

Visin:

Proyectar a los estudiantes como un matematizado de la realidad, capacitndolo para

resolver problemas de la matemtica, la ciencia y la cotidianidad.

Ser un rea que forme estudiantes de excelencia acadmica, creativos, dinmicos,

responsables con el espritu crtico del liderazgo y que al terminar posean un alto nivel
de competibidad en su campo laboral y profesional.

4
 Syllabus

5
6
7
8
9
10
11
12
 UNIVERSIDAD TCNICA DE MANAB
DEPARTAMENTO DE MATEMTICAS
Y ESTADSTICA MATRIZ
AUTOEVALUACIN ESTUDIANTE
ESTUDIANTE: NIVEL: FCHA:

PROFESOR: SEMESTRE: PUNTAJE:

ESCLA: EXCELENTE (3) MUY BUENO (1) DEFICIENTE MARQUE: X

BUENO (2) (0)
N ASPECTO A 3 2 1 0 OBSERVACIONES
AVALUAR

DISCIPLINA
1 PUNTUALIDAD
EN LAS CLASES
2 MANTENER EL
SALON
ORDENADO
3 ETICA EN SUS
EVALUACIONES

RELACIN CON LOS COMPAEROS

4 PARTICIPACIN
ACTIVAMENTE
EN EL TRABAJO
COLABORATIVA

ATENCIN EN CLASE
5 PRESTAR
ATENCIN A L
CLASE
6 APLICAR LAS
INSTRUCCIONES
DEL DOCENTE

INTERS EN CLASE
7 UTILIZAR LOS
TICs PARA EL
APRENDIZAJE
8 TERMINO LAS
ACTIVIDADES
ASIGNADAS
PARA REALIZAR
EN EL AULA
9 ME ESFUERZO
EN LA
REALIZACIN Y
ENTREGA
PUNTUAL DE
LAS TAREAS
10 ORGANIZO Y
ESCRIBO L
INFORMACIN
DE LA CLASE
COMO
ESTRATEGIA.
11 REALIZO
PREGUNTAS
SOBRE LO QUE
NO ENTIENDO
DE LA CLSE

PRTICIPACIN EN CLASE
12 PARTICIPO DE
FORM ACTIVA Y
REGULARMENT
E EN LOS
TALLERES
13 REALIZO
PREGUNTAS
CLARAS Y
OPORTUNA

13
 14 EVITO
HABLAR
SOBRE TEMAS
QUE NO SON
15 PARTE DE LA
EXPRESO
MI
CLASE
S OPINIONES
DE
MANERA
RAZONABLE
Y
COHERENTE

14
15
 AUTORRETRATO

Mi nombre es Jos ngel Andrade Toala soy estudiante de la asignatura de Mtodos

Numricos , actualmente curso el tercer semestre en la Escuela de Ingeniera Civil, Facultad
de Ciencias Matemticas Fsicas y Qumicas de la Universidad Tcnica de Manab. Soy una
persona responsable, honesta, organizada.

Mis metas son convertirme en un profesional de la Ingeniera Civil en la Universidad Tcnica

de Manab y ejercer mi carrera de forma eficiente y correcta para as ser un ejemplo y para
eso tengo que adquirir conocimientos de todo mbito como lo estoy haciendo con la materia
de Mtodos Numricos.

16
 EVIDENCIA DE ACTIVIDADES DE APRENDIZAJE

CLASE N1: MTODOS NUMRICOS.

UNIDAD I: TEORIA DE ERRORES.

17
EJERCICIO EN CLASE:

X XA ERROR EA ER
2 5 -3 3 150
2 5 -3 3 150
6,00 5 1 1 16,6666667
5,00 56 -51 51 1020
4,00 6 -2 2 50

18
UNIDAD I: TEORA DE ERRORES.

CLASE N| 2: TRABAJO CON NMEROS APROXIMADOS.

ASUNTO. CIFRAS SIGNIFICATIVAS DE UN NMERO.

CIFRAS EXACTAS DE UN NMERO APROXIMADO.

REGLAS DE REDONDEO.

19
CLASE 3: UNIDAD N 2. OBTENCIN DE RACES DE FUNCIONES

Se trata de encontrar las races reales de una ecuacin con una incgnita del tipo
f(x)=0. Estos problemas surgen con mucha frecuencia en ingeniera. Encontrar
los nmeros reales que en cierto intervalo que sean solucin de una ecuacin
polinmica o trascendente.

Mtodo Grfico:
El mtodo grfico permite encontrar los intervalos donde se encuentran las
races de las ecuaciones estudiadas. Aunque el mtodo no permite obtener el
resultado de la raz, si localizar las mismas y precisar los intervalos ms ptimos
que deben introducirse para obtener la raz aplicando otros mtodos. Por ello
siempre comenzaremos el trabajo aplicando el mtodo grfico.

Ejercicio 1.
2
Sea la funcin f(x)= lnx + 4x x - 2, encuentre un intervalo donde se ubiquen
sus posibles ceros.

>> f='log(x)+4*x-x^2-2';

>> ezplot(f,[0.1,5]);grid on
log(x)+4 x-x 2-2

-1

-2

-3

-4

-5

-6
0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
x

>> format long

>> x=fzero('log(x)+4*x-x^2-2',0.5);

x = 0.711653712812500

20
EJERCICIO EN CLASE

CLASE 4 MTODO DE BISECCIN.

Sea f(x) una funcin real, definida en un intervalo [a,b] cerrado, si f(a) y f(b)
tienen signos opuestos, entonces f(a)*f(b)<0, entonces existe al menos una
raz en algn punto x perteneciente al intervalo.

El intervalo se subdividir en intervalos ms pequeos, siempre que se cumpla

la condicin del cambio de signo en los extremos. El proceso se repite
logrando cada vez, una mejor aproximacin a la raz.

Pasos de solucin.

Paso 1: Elija valores iniciales inferior a y superior b que encierren la raz de

manera que la funcin cambie de signo en el intervalo.
+
Paso 2: Una aproximacin X1 de la raz se hallar mediante 1 = .
2

Paso 3: Realizar las siguientes evaluaciones para determinar en qu sub

intervalo se encuentra la raz.

Si f(a)*f(xi)< 0 entonces la raz se encuentra en el subntralo de la

izquierda, haga b=X1.

21
 y vuelva a realizar el paso 2.
Si f(a)*f(x1)>0, entonces la raz se encuentra en el subntralo de la
derecha y se hace a=x1
Y se vuelve al paso 2.
Si f(a)*f(b)=0 la raz es x1 y termina el clculo.

Criterio de paro y estimacin de errores: Se puede obtener el error relativo

mximo de la siguiente forma.

= 100. Cuando este error sea menor o igual que el

error relativo mximo fijado para encontrar la raz debemos parar.

Si se utiliza el criterio del error absoluto entonces ABS (X nueva X

anterior)<EA.

Entre las ventajas del mtodo, estn las siguientes, que lo hacen muy atractivo
p a r a el trabajo con computadoras , siempre que el tiempo a emplear no
sea un factor demasiado importante:
1. Las condiciones q u e se requieren para la convergencia (ver las
hiptesis) son mnimas.

2. La rapidez de la convergencia es independiente de la ecuacin a

resolver.

3. El trabajo manual previo a la solucin en la computadora es mnimo y

en la prctica consiste solamente en separar la raz deseada.

4. La acotacin del error es muy simple y segura.

5. La cantidad de memoria que requiere es insignificante.

EJERCICIO EN CLASE:
1. Encuentre la raz positiva de f(x)=x^4-8x^3-35x^2+450x-1001. Exprese
el resultado con 4 cifras significativa despus de la coma decimal.
Utilice el mtodo de falsa posicin

22
f='x^4-8*x^3-35*x^2+450*x-1001';

>> ezplot(f);grid on

>> ezplot(f, [2,6]);grid on

>> ezplot(f, [4,6]);grid on

>> x1=fzero(f, 5.603)

x1 =

5.6098

a b x f(a) f(b) f(x) CONDICION TOLERANCIA

5 6 5,125 -1 7 -1,054443359 0,0000005
5,125 6 5,23955018 -1,054443359 7 -1,118671487 SEGUIR
5,23955018 6 5,344332537 -1,118671487 7 -1,088857828 SEGUIR
5,344332537 6 5,432593286 -1,088857828 7 -0,929214254 SEGUIR
5,432593286 6 5,499086937 -0,929214254 7 -0,691596706 SEGUIR
5,499086937 6 5,544126977 -0,691596706 7 -0,459487624 SEGUIR
5,544126977 6 5,572207728 -0,459487624 7 -0,281664847 SEGUIR
5,572207728 6 5,588755323 -0,281664847 7 -0,164026742 SEGUIR

23
CLASE 5: RACES DE FUNCIONES TRASCENDENTES.

Se toma dos puntos tales que

f(xl)*f(xu)<0, se traza la lnea secante que une los dos
Criterio de parada. Se toma como xn-xn- puntos, esta lnea corta el eje x en el punto x1, la
1< Error. ecuacin de la secante ser.
Y=f(xl)+[(f(xu )-f(xl))/(xu-xl)] (x-xl).
Creacin de una plantilla de trabajo. Haciendo y=0, se despeja x obtenindose
x1=xl-[f(xl)(xu-xl)]/[f(xu)-f(xl)]. Se obtienen otras
Introducir a y b. Cercanos a la raz aproximaciones , como
buscada x2=x1-[f(x1)(xu-x1)]/[f(x1)-f(x1)].
Debe notar que en cada caso se tomen los puntos
() ()
= . cuyo producto
() () f(xn)*f(xn-1)<0, para asegurar que la raz se encuentre
en el intervalo utilizado.
si(f(a)*f(x) < 0 ;dejar a ; poner a = x).
si(f(b)*f(x) < 0 ;dejar b ; poner b = x ).

Criterio de parada. Si (abs(xnueva-xanterior) < EAM ; detener; seguir).

24
EJERCICIO EN CLASE

a b x f(a) f(b) f(x) condicion toleracia

0,5 1 0,795686101 0,229425539 -0,158529015 0,081227532 SEGUIR 0,0000005
0,795686101 1 0,864905957 0,081227532 -0,158529015 0,012971946 SEGUIR
0,864905957 1 0,875124166 0,012971946 -0,158529015 0,001780781 SEGUIR
0,875124166 1 0,876511333 0,001780781 -0,158529015 0,000239268 SEGUIR
0,876511333 1 0,876697434 0,000239268 -0,158529015 3,20553E-05 SEGUIR
0,876697434 1 0,876722361 3,20553E-05 -0,158529015 4,29285E-06 SEGUIR
0,876722361 1 0,876725699 4,29285E-06 -0,158529015 5,7487E-07 DETENER

CLASE 6: NEWTON RHAPSON

f='x^4-8*x^3-35*x^2+450*x-1001';

>> ezplot(f);grid on

>> ezplot(f, [2,6]);grid on

>> ezplot(f, [4,6]);grid on

>> x1=fzero(f, 5.603)

x1 =

5.6098

25
a b X f(a) f(b) f(x) CONDICION TOLERANCIA
5 6 5,125 -1 7 -1,054443359 0,0000005
5,125 6 5,23955018 -1,054443359 7 -1,118671487 SEGUIR
5,23955018 6 5,344332537 -1,118671487 7 -1,088857828 SEGUIR
5,344332537 6 5,432593286 -1,088857828 7 -0,929214254 SEGUIR

26
CLASE 7: RACES MLTIPLES.

TRABAJO EN CLASE:

F(X)=4X^2-7.2X+3.24

27
>> f='4*X^2-7.2*X+3.24';

>> ezplot(f);grid on

>> ezplot(f, [-2,2]);grid on

>> d1=diff(f)

d1 =

8*X - 7.2

>> d2=diff(d1)

d2 =

8.0

28
A f(Xo) f'(Xo) f''(Xo) X CONDICION Tolerancia
0,8268 0,02143296 -0,5856 8 0,9 DETENER 0,000000005
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER

CLASE 8: SOLUCIN DE SISTEMAS DE ECUACIONES NO

LINEALES.

La solucin de este sistema consta de un conjunto de valores xi que

simultneamente hacen que todas las ecuaciones sean iguales a cero.
En la unidad tres, presentaremos los mtodos, para el caso en que las ecuaciones
simultneas son lineales, es decir, que se puedan expresar en la forma general
f(x) = a1x1 + a2x2 + + anxn b = 0 donde la b y las a son constantes. A las
ecuaciones algebraicas y trascendentes que no se pueden expresar de esta forma
se les llama ecuaciones no lineales. Por ejemplo
x2 + xy = 10
y + 3xy2 = 57
son dos ecuaciones simultneas no lineales con dos incgnitas, x y y, las cuales
se expresan en la forma de la ecuacin.
u(x, y) = x2 + xy 10 = 0
v(x, y) = y + 3xy2 57 = 0

As, la solucin seran los valores de x y de y que hacen a las funciones u(x, y)
y v(x, y) iguales a cero. La mayora de los mtodos para determinar tales
soluciones son extensiones de los mtodos abiertos para resolver ecuaciones
simples. En esta seccin presentaremos el de Newton - Raphson.
Recuerde que el mtodo de Newton-Raphson se utiliz empleando la derivada
(al evaluar, es la pendiente de la recta tangente) de una funcin, para calcular
su interseccin con el eje de la variable independiente; esto es, la raz (figura
6.5). Dicho clculo se bas en la expansin en serie de Taylor de primer orden
,
f(xi + 1) = f(xi) + (xi+1 xi) (xi)

29
donde xi es el valor inicial de la raz y xi+1 es el valor en el cual la recta tangente
intersecta el eje x. En esta interseccin, f(xi + 1) es, por definicin, igual a cero.
TRABAJO EN CLASE:

CLASE 9: RACES DE POLINOMIOS.

En este captulo estudiaremos los mtodos para encontrar las races de

ecuaciones polinmicas. De la forma general fn(x) = a0 + a1x + a2x2 +... + anxn
donde n es el grado del polinomio y las a son los coeficientes del polinomio.
Aunque los coeficientes pueden ser nmeros reales o complejos, este estudio se
limitar a los casos en que son reales. Entonces las races del polinomio pueden
ser reales y/o complejas.

Las races de los polinomios cumplen estas reglas:

1. En una ecuacin de grado n, hay n races reales o complejas. Se debe notar
que esas races no necesariamente son distintas.
2. Si n es impar, hay al menos una raz real.
30
3. Si existen races complejas, stas se encuentran por pares conjugados (es
decir, R +
i y R i), donde i = 1 .

MTODO DE MLLER
Recuerde que el mtodo de la secante obtiene una aproximacin de la raz
dirigiendo
una lnea recta hasta el eje x con dos valores de la funcin (figura 7.3a). El
mtodo de
Mller es similar; pero se construye una parbola con tres puntos (figura
7.3b).
El mtodo consiste en obtener los coeficientes de la parbola que pasa por los
tres
puntos. Dichos coeficientes se sustituyen en la frmula cuadrtica para
obtener el valor
donde la parbola interseca al eje x; es decir, la raz estimada. La
aproximacin se facilita
al escribir la ecuacin de la parbola en una forma conveniente.

TRABAJO EN CLASE:

y=x^2+ln(x-1)

e^y=5*sen(x)

>> u='x^2+log(x-1)-y';

31
>> v='5+sin(x)-exp(y)';
>> ezplot(u);grid on;hold on;ezplot(v);
>> [x,y]=solve(u,v);
>> x=eval(x)
x=
-1.2154 + 1.0575i
>> y=eval(y)
y=
1.2571 + 0.1256i
>> syms x
>> d1=diff(u,x)
d1 =
2*x + 1/(x - 1)
5+sin(x)-exp(y) = 0
>> d3=diff(v,x) 6

d3 =
4
cos(x)
X: 1.539
X: -0.9053 Y: 1.792
>> syms y 2 Y: 1.438

>> d2=diff(u,y)
0
y

d2 =
-1 -2

>> d4=diff(v,y)
-4
d4 =
-exp(y)
-6
-6 -4 -2 0 2 4 6
x

32
 DETERMINA
Xi Yi d1 d2 d3 d4 U V X Y
NTE

- -
4,94679 0,04078 0,00971 1,54732 1,79174
1,53 1,79 -1 5,98945 0,08397 -29,58779325
245 501 548 807 009
247 827

- - -
1,54732 1,79174 4,92171 0,02346 1,54736 1,79171
-1 5,99988 0,00022 0,00015 -29,50624515
807 009 39 61 799 373
374 283 911

- -
1,54736 1,79171 4,92166 0,02342 -1,066E- 1,54736 1,79171
-1 5,99972 2,8815E -29,50518615
799 373 05 619 09 799 373
557 -09

CLASE 10: POLY, ROOTS, POLIVAL

33
CLASE 11: METODO DE GAUUS

El mtodo de Gauss consiste en realizar sobre el sistema dado de ecuaciones

una sucesin de transformaciones que, sin afectar la solucin, conduzcan a un
sistema especialmente simple, el cual es resuelto de inmediato. Estas
transformaciones, llamadas comnmente elementales, son:
1. Intercambiar dos ecuaciones del sistema.
2. Multiplicar ambos miembros de una ecuacin del sistema por un nmero
diferente de cero.
3. Sumar una ecuacin, miembro a miembro, a otra ecuacin del sistema
multiplicada por un nmero real cualquiera.
Es bastante evidente que la realizacin de cualquier nmero de
transformaciones elementales sobre un sistema conduce a uno equivalente, es
decir, a un sistema con las mismas incgnitas y las mismas soluciones.
El mtodo de Gauss consta de dos etapas, llamadas proceso directo y proceso
inverso.

Proceso directo
Consiste en la realizacin de transformaciones elementales que conduzcan a un
sistema triangular superior:

Proceso inverso
De la ltima ecuacin se puede obtener la incgnita xn; conocida esta, se
sustituye en la penltima ecuacin y se obtiene xn-1, etctera, hasta llegar a la
primera ecuacin, de la cual se halla x1
En el siguiente ejemplo se ilustra el mtodo para un sistema de 4x4.

34
35
TRABAJO EN CLASE.

Inversa de una matriz

CLASE 12: METODO DE GAUUS

36
REALIZACION DEL METODO DE GAUUS EN EXCEL

factor c1 c2 c3 r
2 1 -1 7
1 2 1 1 8
0,5 1 -3 1 -2

factor c1 c2 c3 r
2 1 -1 7
0 -3,5 1,5 -5,5
0 0 0 2 1

factor c1 c2 c3 r
2 1 -1 7 x1 2,857142857
0 -3,5 1,5 -5,5 x2 1,785714286
0 0 2 1 x3 0,5

37
CLASE 13: METODO DE JACOBI

Sea AX = b un sistema de n ecuaciones lineales. Se supone, como hasta ahora,

que este sistema es cuadrado y de solucin nica. En forma desarrollada el

sistema es:

38
39
CLASE 14: .METODO SEIDEL

El mtodo de Seidel, tambin llamado de Gauss-Seidel, constituye una

variacin del mtodo de Jacobi que logra, sin complicar el algoritmo, mejorar
la rapidez de convergencia en la mayora de los casos.
En el mtodo de Seidel se halla x1(k) de la misma forma que en Jacobi; pero,
una vez calculado, se le utiliza de inmediato para el clculo de las restantes
incgnitas; una vez hallado x2(k) se usa este en el clculo de las dems en lugar
de x2k-1) de la aproximacin anterior, y as con las dems incgnitas. En
trminos ms precisos:

40
41
42
CLASE15: UNIDAD IV: INTERPOLACION Y RECTIFICACION DE
CURVAS

INTERACION LINEAL

Sea una funcin g : [a, b] R tal que en a y b se conoce g(a) y g(b), entonces
se puede interpolar (valor intermedio) g(x) por una funcin lineal f(x) = mx +
n.

Veamos:
F=m*x+n
()() () ()
g(a) = f(a) = ma + n = =

g(b) = f(b) = mb + n
8
que equivale a
Sean los datos empricos obtenidos de una investigacin:

X=[0 1 2 3 4 5]; y=[15 10 9 6 2 0];

Cada tramo es unido por una lnea recta interpolar es encontrar el valor de f(x)
para una x que este comprendida en el dominio.

Como realizar esta interpolacin lineal en Matlab.

15

Sean los datos empricos obtenidos de una

investigacin:
Realice la interpolacin lineal en matlab. 10

Como realizar esta interpolacin lineal en

Matlab.
5

0
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5

43
TRABAJO EN CLASE

44
 TAREAS

TAREA#5-6 (RAICES MULTIPLES Y SOLUCIONES DE

ECUACIONES)

RAICES MULTIPLES

F(X)=4X^2-7.2X+3.24

>> f='4*X^2-7.2*X+3.24';

>> ezplot(f);grid on

>> ezplot(f, [-2,2]);grid on

>> d1=diff(f)

d1 =

8*X - 7.2

>> d2=diff(d1)

d2 =

8.0

45
A f(Xo) f'(Xo) f''(Xo) X CONDICION Tolerancia
0,8268 0,02143296 -0,5856 8 0,9 DETENER 0,000000005
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER
0,9 0 -1,0658E-14 8 0,9 DETENER

46
 A f(Xo) f'(Xo) f''(Xo) X CONDICION Tolerancia
1,002 0,041616 0,816 8 0,9 DETENER 0,000000005
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER
0,9 0 2,1316E-14 8 0,9 DETENER

SOLUCIONES DE ECUACIONES

y=x^2+ln(x-1)

e^y=5*sen(x)

>> u='x^2+log(x-1)-y';

>> v='5+sin(x)-exp(y)';

>> ezplot(u);grid on;hold on;ezplot(v);

>> [x,y]=solve(u,v);

>> x=eval(x)

x=

-1.2154 + 1.0575i

>> y=eval(y)

y=

1.2571 + 0.1256i

>> syms x

47
 5+sin(x)-exp(y) = 0
>> d1=diff(u,x)
6
d1 =

2*x + 1/(x - 1)
4
>> d3=diff(v,x)
X: 1.539
X: -0.9053 Y: 1.792
d3 = Y: 1.438
2
cos(x)

>> syms y
0
y
>> d2=diff(u,y)

d2 =
-2
-1

>> d4=diff(v,y)
-4
d4 =

-exp(y)
-6
-6 -4 -2 0 2 4 6
x

Xi Yi d1 d2 d3 d4 U V DETERMINANTE X Y
- -
1,53 1,79 4,94679245 -1 0,04078501 0,00971548 -29,58779325 1,54732807 1,79174009
5,98945247 0,08397827
- - -
1,54732807 1,79174009 4,9217139 -1 0,0234661 -29,50624515 1,54736799 1,79171373
5,99988374 0,00022283 0,00015911
- -2,8815E-
1,54736799 1,79171373 4,9216605 -1 0,02342619 -1,066E-09 -29,50518615 1,54736799 1,79171373
5,99972557 09
-
1,54736799 1,79171373 4,9216605 -1 0,02342619 0 0 -29,50518613 1,54736799 1,79171373
5,99972557
-
1,54736799 1,79171373 4,9216605 -1 0,02342619 0 0 -29,50518613 1,54736799 1,79171373
5,99972557
-
1,54736799 1,79171373 4,9216605 -1 0,02342619 0 0 -29,50518613 1,54736799 1,79171373
5,99972557
-
1,54736799 1,79171373 4,9216605 -1 0,02342619 0 0 -29,50518613 1,54736799 1,79171373
5,99972557
-
1,54736799 1,79171373 4,9216605 -1 0,02342619 0 0 -29,50518613 1,54736799 1,79171373
5,99972557
-
1,54736799 1,79171373 4,9216605 -1 0,02342619 0 0 -29,50518613 1,54736799 1,79171373
5,99972557
-
1,54736799 1,79171373 4,9216605 -1 0,02342619 0 0 -29,50518613 1,54736799 1,79171373
5,99972557

48
TAREA # 3-4(FALSA POSICION Y NWETON RHASON)
MTODO DE NEWTON RHASON

1. Determine las races de f(x)=8sen(x)e^(-x)-1. Con seis cifras exactas, despus de la coma decimal,
tolerancias=0.0000005; en el intervalo entre (-6.5;-2)

>> f='8*sin(x)*exp(-x)-1';

>> ezplot(f);grid on;

>> ezplot(f,[-4,2]);grid on;

>> ezplot(f,[-3.5,2.5]);grid on;

>> diff(f)

ans =

(8*cos(x))/exp(x) - (8*sin(x))/exp(x)

>> x1=fzero(f,-3.1)

x1 =

-3.146965474469272

>> x2=fzero(f,0.1)

x2 =

0.145014812533098

>> x3=fzero(f,2.0)

x3 =

1.989221390478615

49
 8 sin(x) exp(-x)-1

40

20
X: -3.151 X: 0.1046 X: 2.013
Y: 0.7135 Y: -0.2476 Y: -0.03401
0

-20

-40

-60

-3 -2 -1 0 1 2
x

A f(x) f'(x) x1 condicin tolerancia

-3,1 -8,38404415 -170,045983 -3,14930457 0,0000005
-3,14930457 0,43871118 -187,991897 -3,1469709 SEGUIR
-3,1469709 0,00101519 -187,122197 -3,146965474 SEGUIR
-3,14696547 5,4782E-09 -187,120178 -3,146965474 DETENER
-3,14696547 -9,4369E-15 -187,120178 -3,146965474 DETENER

A f(x) f'(x) x2 condicin tolerancia

0,1 -0,27733591 6,479871911 0,142799598 0,0000005
0,142799598 -0,01298697 5,87781672 0,145009086 SEGUIR
0,145009086 -3,3485E-05 5,847519523 0,145014812 SEGUIR
0,145014812 -2,2453E-10 5,847441102 0,145014813 DETENER
0,145014813 0 5,847441101 0,145014813 DETENER

A f(x) f'(x) x1 condicin tolerancia

2 -0,0155198 -1,435035 1,989185071 0,0000005
1,989185071 5,2471E-05 -1,44471715 1,98922139 SEGUIR
1,98922139 5,8655E-10 -1,44468485 1,98922139 DETENER
1,98922139 0 -1,44468485 1,98922139 DETENER

50
 2. El polinomio f(x)=0.0074x^4-0.284x^3+3.355x^2-12.183x+5 tiene una raz real entre 15 y 20. Aplique el
mtodo de newton-raphson a dicha funcin. Encuentre el resultado con error absoluto mximo menor o igual
a 0.0000005.

syms x

>> f='0.074*x^4-0.284*x^3+3.355*x^2-12.183*x+5';

>> diff (f)

ans =

0.296*x^3 - 0.852*x^2 + 6.71*x - 12.183

>> format long

>> ezplot(f);grid on

x1=fzero(f, 0.5949)

x1 =

0.468828633652398

a f(x) f'(x) x condicin

0,46 0,07842591 -9,24787174 0,46848043
0,46848043 -6,2916E-06 -9,19605361 0,46847974 seguir
0,46847974 1,8741E-08 -9,19605779 0,46847975 detener
0,46847975 -5,5822E-11 -9,19605778 0,46847975 detener
0,46847975 1,6609E-13 -9,19605778 0,46847975 detener
0,46847975 0 -9,19605778 0,46847975 detener
0,46847975 0 -9,19605778 0,46847975 detener
0,46847975 0 -9,19605778 0,46847975 detener
0,46847975 0 -9,19605778 0,46847975 detener
0,46847975 0 -9,19605778 0,46847975 detener
0,46847975 0 -9,19605778 0,46847975 detener
0,46847975 0 -9,19605778 0,46847975 detener
0,46847975 0 -9,19605778 0,46847975 detener
0,46847975 0 -9,19605778 0,46847975 detener
0,46847975 0 -9,19605778 0,46847975 detener
0,46847975 0 -9,19605778 0,46847975 detener

51
 0,46847975 0 -9,19605778 0,46847975 detener
0,46847975 0 -9,19605778 0,46847975 detener
0,46847975 0 -9,19605778 0,46847975 detener
0,46847975 0 -9,19605778 0,46847975 detener
0,46847975 0 -9,19605778 0,46847975 detener

3. Dada f(x)=-2x^5-1.5x^4+10x+2 use el mtodo de falsa posicin para determinar el mximo de esta funcin.
Haga una aproximacin hasta la tercera cifra exacta despus del punto decimal.

>> f='-2*x^5-1.5*x^4+10*x+2';

>> ezplot(f);grid on

>> ezplot(f, [-2,2]);grid on

>> x1=fzero(f, 1.39)

x1 =

1.388215842381694

52
a b x f(a) f(b) f(x) CONDICION TOLERANCIA
1 2 1,11409396 8,5 -66 7,397327873 0,0000005
1,11409396 2 1,203379737 7,397327873 -66 5,841089843 SEGUIR
1,203379737 2 1,268149502 5,841089843 -66 4,242344667 SEGUIR
1,268149502 2 1,31235022 4,242344667 -66 2,888856243 SEGUIR
1,31235022 2 1,341186834 2,888856243 -66 1,879308466 SEGUIR
1,341186834 2 1,359426754 1,879308466 -66 1,185803185 SEGUIR
1,359426754 2 1,370732622 1,185803185 -66 0,733702798 SEGUIR
1,370732622 2 1,377651093 0,733702798 -66 0,448444594 SEGUIR
1,377651093 2 1,381851177 0,448444594 -66 0,272035441 SEGUIR
1,381851177 2 1,384388572 0,272035441 -66 0,164266724 SEGUIR
1,384388572 2 1,385916957 0,164266724 -66 0,098916251 SEGUIR
1,385916957 2 1,386835925 0,098916251 -66 0,059464588 SEGUIR
1,386835925 2 1,387387876 0,059464588 -66 0,035711788 SEGUIR
1,387387876 2 1,387719174 0,035711788 -66 0,021433932 SEGUIR
1,387719174 2 1,387917951 0,021433932 -66 0,0128598 SEGUIR
1,387917951 2 1,388037189 0,0128598 -66 0,007713862 SEGUIR
1,388037189 2 1,388108705 0,007713862 -66 0,004626502 SEGUIR
1,388108705 2 1,388151595 0,004626502 -66 0,002774594 SEGUIR
1,388151595 2 1,388177315 0,002774594 -66 0,001663894 SEGUIR
1,388177315 2 1,388192739 0,001663894 -66 0,000997791 SEGUIR
1,388192739 2 1,388201989 0,000997791 -66 0,000598338 SEGUIR
1,388201989 2 1,388207535 0,000598338 -66 0,000358797 SEGUIR
1,388207535 2 1,388210861 0,000358797 -66 0,000215153 SEGUIR
1,388210861 2 1,388212855 0,000215153 -66 0,000129017 SEGUIR
1,388212855 2 1,388214051 0,000129017 -66 7,73648E-05 SEGUIR
1,388214051 2 1,388214768 7,73648E-05 -66 4,63917E-05 SEGUIR
1,388214768 2 1,388215198 4,63917E-05 -66 2,78187E-05 DETENER
1,388215198 2 1,388215456 2,78187E-05 -66 1,66814E-05 DETENER
1,388215456 2 1,388215611 1,66814E-05 -66 1,0003E-05 DETENER

53
TAREA#7(MATRICES)

>> c=[3 -9 5 2 2 -15;9 -3 -8 -2 4 64;-5 4 4 2 6 -80;4 7 7 5 8 -112;4 4 5 1 5 -59]

c=

3 -9 5 2 2 -15

9 -3 -8 -2 4 64

-5 4 4 2 6 -80

4 7 7 5 8 -112

4 4 5 1 5 -59

>> rref(c)

ans =

1.0000 0 0 0 0 1.9281

0 1.0000 0 0 0 -2.9117

0 0 1.0000 0 0 -5.6708

0 0 0 1.0000 0 -4.9695

0 0 0 0 1.0000 -4.3484

factor c1 c2 c3 c4 c5 r
3 -9 5 2 2 -15
3 9 -3 -8 -2 4 64
-1,6666667 -5 4 4 2 6 -80
1,33333333 4 7 7 5 8 -112
1,33333333 4 4 5 1 5 -59

54
factor c1 c2 c3 c4 c5 r
3 -9 5 2 2 -15
0 24 -23 -8 -2 109
-0,4583333 0 -11 12,3333333 5,33333333 9,33333333 -105
0,79166667 0 19 0,33333333 2,33333333 5,33333333 -92
0,66666667 0 16 -1,66666667 -1,66666667 2,33333333 -39

factor c1 c2 c3 c4 c5 r
3 -9 5 2 2 -15
0 24 -23 -8 -2 109
0 0 1,79166667 1,66666667 8,41666667 -55,0416667
10,3488372 0 0 18,5416667 8,66666667 6,91666667 -178,291667
7,62790698 0 0 13,6666667 3,66666667 3,66666667 -111,666667

factor c1 c2 c3 c4 c5 r
3 -9 5 2 2 -15
0 24 -23 -8 -2 109
0 0 1,79166667 1,66666667 8,41666667 -55,0416667
0 0 0 -8,58139535 -80,1860465 391,325581
1,05420054 0 0 0 -9,04651163 -60,5348837 308,186047

factor c1 c2 c3 c4 c5 r
3 -9 5 2 2 -15
0 24 -23 -8 -2 109
0 0 1,79166667 1,66666667 8,41666667 -55,0416667
0 0 0 -8,58139535 -80,1860465 391,325581
0 0 0 0 23,99729 -104,349593
X1 1,92810088
X2 -2,91172596
X3 -5,67080745
X4 -4,96950875
X5 -4,34839074

55
56
57