CÁLCULO NUMÉRICO, MAT1123

M.J. Andrade, D. Sepúlveda, E.Cariaga

Primer Semestre 2020.

INFORME 1

1. En ciertas ocasiones, los ingenieros aeroespaciales deben calcular las

trayectorias de proyectiles, como cohetes. Un problema parecido tiene
que ver con la trayectoria de una pelota que se lanza. Dicha trayectoria
está definida por las coordenadas (x, y). La trayectoria se modela con
la ecuación
g
y = x tan θ0 − 2 2
x2 + y0 .
2v0 cos θ0
Se pide calcular el ángulo inicial θ0 apropiado si la velocidad inicial
v0 = 20[m/s] y la distancia x al catcher es de 40[m]. La pelota sale
de la mano del lanzador con una elevación y0 = 1.8[m], y el catcher
la recibe a 1[m]. Para g utilice un valor de 9.81[m/s2 ]. Exprese su
resultado final en grados sexagesimales.

2. En ingenierı́a oceanográfica, la ecuación de una ola estacionaria reflejada

en un puerto está dada por:
h 2πx 2πtv
− e−x = sin( ) cos( ),
h0 λ λ
con λ = 16, t = 12, y v = 48. Resuelva para el menor valor positivo de
x, si h = 0.4h0 .

3. La ecuación de Ergun

∆P Dp 3 µ
ρ 2
= 150 (1 − ) + 1.75,
Go L 1 −  Dp Go

modela el flujo de un lı́quido a través de un cauce comprimido, en

donde, ∆P es la caı́da de presión, ρ es la densidad del fluido, Go es la
velocidad música (el cociente del flujo de masa dividido entre el área
de la sección transversal), Dp es el diámetro de las partı́culas dentro
del cauce, µ es la viscocidad del fluido, L es la longitud del cauce, y 
es la fracción vacı́a del cauce. Se pide calcular la fracción vacı́a  del
cauce, si Dp Go µ−1 = 1000 y ∆P ρDp G−2 o L
−1
= 10.

1
INSTRUCCIONES

ˆ Fecha de entrega: se confirmará oportunamente, pero no será

antes del 5 de mayo.
ˆ Formato de archivo final: PDF.
ˆ Individual.
ˆ Tipeado en Word o Tex.
ˆ Gráficos hechos con software.
ˆ Gráficos en Anexo.
ˆ La respuesta final debe estar debidamente fundamentada.
ˆ Dependiendo del número de iteraciones requeridas puede informar
sólo valores representativos.
ˆ En Anexo los programas .m en Octave.

2
 CÁLCULO NUMÉRICO, MAT1123
M.J. Andrade, D. Sepúlveda, E.Cariaga
Primer Semestre 2020.

INFORME 2

PROBLEMA 1. Considere el sistema de ecuaciones lineales de dimensión

8 × 8 que surge luego de realizar un balance de fuerzas en una estructura
metálica:

√
2
−F1 + f1 + f2 = d1
√ 2 √
2 3
− f1 + f4 = 0
2 2
−f + f5 = d2
√ 2
3
− f4 − f5 = 0
√2
2
f1 − F 2 = d3
√ 2
2 1
− f1 − f3 − f4 = d4
2 2
f3 − 10000 = d5
1
f4 − F 3 = d6
2
con vector incógnita 8 × 1:

x = [F1 , F2 , F3 , f1 , f2 , f3 , f4 , f5 ]0 .

Se pide resolver utilizando los métodos iterativos de Jacobi, Gauss-Seidel y

Relajación (note que quizás alguno podrı́a ser divergente). En el caso del
método de Relajación responda con el mejor valor de ω calculado. Utilice
como medida de error y tolerancia ||A · xn − b||∞ < 10−10 . Note que di , i =
1, 2, 3, 4, 5, 6, son los seis primeros dı́gitos de su carné de identidad.

3
INSTRUCCIONES

ˆ Fecha de entrega: 31 de mayo.

ˆ Formato de archivo final: PDF.

ˆ Individual.

ˆ Tipeado en Word o Tex.

ˆ La respuesta final debe estar debidamente fundamentada.

ˆ Dependiendo del número de iteraciones requeridas puede informar sólo

valores representativos.

ˆ En Anexo los programas .m en Octave.

4
 CÁLCULO NUMÉRICO, MAT1123
M.J. Andrade, D. Sepúlveda, E.Cariaga
Primer Semestre 2020.

INFORME 3. AJUSTE DE CURVAS.

PROBLEMA. Considere los datos:

(xi−1 , yi−1 ) = (1 − i, ln(di + 2)), i = 1, ...6

en donde: di representa el dı́gito i-ésimo de su RUN. Por ejemplo, si su RUN

es: 16.932.578-K entonces: d1 = 1, d2 = 6, d3 = 9, ... y ası́ sucesivamente.
Para estos datos se pide:

1. Construir un único polinomio de interpolación, de grado mı́nimo, y que

pase por todos los puntos, haciendo uso de notación matricial (matriz de
Vandermonde) para su cálculo. Grafique simultáneamente el polinomio
calculado con los datos.

2. Construir un único polinomio de interpolación, de grado mı́nimo, y que

pase por todos los puntos, utilizando un comando idóneo de Matlab u
Octave. Compare este polinomio con el calculado en el punto anterior:
¿Son distintos?, ¿son iguales?, ¿Por qué?.

3. Construir un interpolante spline lineal. Presente su respuesta en forma

gráfica. Utilize esta función para interpolar todos los puntos medios de
los nodos dados. Utilice una tabla para informar sus resultados.

4. Construir un interpolante spline cúbico. Presente su respuesta en forma

gráfica. Utilize esta función para interpolar todos los puntos medios de
los nodos dados. Utilice una tabla para informar sus resultados.

5. Utilice los interpolantes construidos en los puntos 1., 3., y 4. anteriores

para interpolar la respuesta Y para X = −4.2; −1.1; −0.5; −3.8. Utilice
una tabla para informar sus resultados.

6. Construya un polinomio de regresión de grado 4 para los datos dados.

Presente sus resultados en forma algebraica (expresión funcional) y
gráfica (polinomio y datos juntamente).

5
INSTRUCCIONES

ˆ Fecha de entrega: domingo 14 Junio.

ˆ Formato de archivo final: PDF.

ˆ Individual.

ˆ Tipeado en Word o Tex.

ˆ En Anexo los programas .m en Octave.

ˆ Sus desarrollos y conclusiones deben estar explicados con claridad, y

debidamente fundamentados.

6
 CÁLCULO NUMÉRICO, MAT1123
M.J. Andrade, D. Sepúlveda, E.Cariaga
Primer Semestre 2020.

INFORME 4. SOL NUMÉRICA DE EDO

PROBLEMA. Sean d1 , d2 , d3 , d4 , d5 , d6 los 6 primeros dı́gitos de su

RUN. Por ejemplo, si su RUN es 24.365.895 − k, entonces d1 = 2, d2 = 4,
d3 = 3, d4 = 6, d5 = 5, d6 = 8. Considere la siguiente ecuación diferencial
ordinaria o PVI
dy
= A · eB·x − C · y
dx
y(0) = D

donde A, B, C y D son constantes definidas como: A := d1 + d2 + d3 , B :=

d4 + d5
, C := d1 + d2 + d3 + d4 + d5 y D := d3 + d4 + d5 + d6 . Se pide:
d1 + d2 + d3
1. Obtener la solución analı́tica del PVI dado. Exprese la solución en
términos literales (A, B, C y D) y con sus respectivos valores numéricos.

2. Resolver numéricamente con los métodos de Euler, Taylor de orden 2,

y Runge-Kutta clásico de orden 4, con tamanho de paso h = 1/10 y
h = 1/5.

3. Construir una tabla resumen que informe el valor exacto y numérico

(para los tres métodos numéricos en el item anterior), en los puntos
x = 0, x = 1/10, x = 1/5, y x = 3/10. Esto para h = 1/10 y h = 1/5.

4. Para h = 1/5 construya tres gráficos. Uno para cada método numérico
mencionado en los puntos anteriores. Para cada caso incluya la curva
interpolante lineal de la solución numérica, y la curva de la solución
exacta del PVI dado.

5. Para h = 1/5 y h = 1/10 calcule los errores asociados a cada uno de los
tres métodos numéricos según la norma discreta ||x̂ − xE ||∞ , en donde
x̂ es la solución numérica, y xE la solución exacta.

7
INSTRUCCIONES

ˆ Fecha de entrega: VIERNES 3 DE JULIO.

ˆ Formato de archivo final: PDF.

ˆ Individual.

ˆ Tipeado en Word o Tex.

ˆ En Anexo los programas .m en Octave.

ˆ Sus desarrollos y conclusiones deben estar explicados con claridad, y

debidamente fundamentados.

ˆ Debe escribir explı́citamente las fórmulas para cada uno de los tres
métodos numéricos utilizados, esto según el PVI dado.

8
 PAUTAS DE CORRECCIÓN
Cálculo Numérico, MAT1123, Primer Semestre 2020

ˆ Informe 1

ˆ Informe 2

ˆ Informe 3

ˆ Informe 4

ˆ Informe 5

9
PAUTA INFORME 1
1. En ciertas ocasiones, los ingenieros aeroespaciales deben calcular las
trayectorias de proyectiles, como cohetes. Un problema parecido tiene
que ver con la trayectoria de una pelota que se lanza. Dicha trayectoria
está definida por las coordenadas (x, y). La trayectoria se modela con
la ecuación
g
y = x tan θ0 − 2 x2 + y0 .
2v0 cos2 θ0
Se pide calcular el ángulo inicial θ0 apropiado si la velocidad inicial
v0 = 20[m/s] y la distancia x al catcher es de 40[m]. La pelota sale
de la mano del lanzador con una elevación y0 = 1.8[m], y el catcher
la recibe a 1[m]. Para g utilice un valor de 9.81[m/s2 ]. Exprese su
resultado final en grados sexagesimales.

Solución: como función auxiliar se recurre a

g
f (θ0 ) = x tan θ0 − 2 x2 + y0 − y.
2v0 cos2 θ0
Ahora se trata de aproximar θ0 para el cual f (θ0 ) = 0. Para aplicar el
método de Newton-Raphson calculamos f 0 :

f 0 (θ0 ) = x sec2 θ0 − x2 (g/v02 ) sec2 θ0 tan θ0 .

Con esto las iteraciones del método de Newton-Raphson, para n =

0, 1, 2, ... se ejecutan como:

f (θ0,n )
θ0,n+1 = θ0,n − .
f 0 (θ0,n )
El código Octave para la función auxiliar f está dado por:

function u=p1inf1a(z)
% z denota al \’angulo $\theta_{0}$.
v0=20;
x=40;
y0=1.8;
y=1;
g=9.81;
%
u=x*tan(z)-g*((2*v0*v0*(cos(z))^2)^(-1))*(x^2)+(y0-y);
endfunction

10
Con la finalidad de identificar un candidato para θ0,0 , esto es, para el
valor inicial para las iteraciones se grafica la función auxiliar según la
instrucción

>> fplot(’p1inf1a’,[0 1.2])

en la lı́nea de comandos de Octave obteniéndose

Este plot se hizo sobre el intervalo [0; 1.2] pues θ0 se considera en

radianes, y tan(θ0 ) se indetermina en π2 ≈ 1.5708. Además, tratándose
de un ángulo de lanzamiento se asume que 0 < θ0 < π2 . Por otro
lado, se observa de la gráfica que existen dos puntos de corte de la
curva con el eje horixontal, y ambos con sentido fı́sico. Uno cerca de
0.6[rad](≈ 36.38[grados]) y el otro cerca de 0.9[rad](≈ 52.47[grados].)
Por lo tanto, se aplicará el método de NR para aproximar ambos ceros.
Ahora calculamos las iteraciones según el programa p1inf1 el cero cerca
de 0.6[rad]:

function u=p1inf1(z0,m)
% z denota al \’angulo $\theta_{0}$.
z=z0;
%
v0=20;
x=40;

11
y0=1.8;
y=1;
g=9.81;
%
for it=1:m
w=x*tan(z)-g*((2*v0*v0*(cos(z))^2)^(-1))*(x^2)+(y0-y);
dw=x*(sec(z))^2-(g/(v0*v0))*(sec(z))^2*tan(z)*(x^2);
z=z-w/dw;
endfor
u=z;
%
endfunction

cuya aplicación en la lı́nea de comandos permite obtener los siguientes

resultados en modo format long

>> u=p1inf1(0.6,5);
w = -6.375347113308372e-01
z = 6.330135022315971e-01
w = -3.330086995741577e-02
z = 6.349448666617670e-01
w = -1.263941388252032e-04
z = 6.349522533534134e-01
w = -1.860151810362254e-09
z = 6.349522534621271e-01
w = 6.439293542825908e-15
z = 6.349522534621268e-01
>>

En los resultados numéricos anteriores w denota a f (θ0,n ) y z a θ0,n .

Para mayor completitud y certeza numérica se informa la aplicación
del comando Octave fzero

>> fzero(’p1inf1a’,0.6)
ans = 6.349522534621277e-01
>>

el cual considera una combinación inteligente de varios algoritmos, y

una programación experta. Tal como puede apreciarse los resultados

12
coinciden exactamente, salvo los dos últimos dı́gitos decimales.

Conclusión: el ángulo inicial θ0 tiene un valor aproximado de

θ0 ≈ 0.6349522[rad] ≈ 36.380081[grados].

En términos prácticos basta decir que: θ0 ≈ 36.4[grados].

Respecto del otro punto de corte cercano a 0.9[rad] los resultados son
los siguientes:

>> u=p1inf1(0.9,5);
w = 4.297935554592656e-01
z = 9.175763214708446e-01
w = -5.224619527246399e-02
z = 9.158646688165071e-01
w = -5.360108545786613e-04
z = 9.158467413106386e-01
w = -5.832043536990739e-08
z = 9.158467393596194e-01
w = -2.553512956637860e-14
z = 9.158467393596185e-01
>> fzero(’p1inf1a’,0.9)
ans = 9.158467393596199e-01

Por lo tanto, para el segundo punto de corte cerca de 0.9[rad]:

θ0 ≈ 0.9158467[rad] ≈ 52.47415059[grados].

En términos prácticos basta decir que: θ0 ≈ 52.5[grados].

13
2. En ingenierı́a oceanográfica, la ecuación de una ola estacionaria reflejada
en un puerto está dada por:
h 2πx 2πtv
− e−x = sin( ) cos( ),
h0 λ λ
con λ = 16, t = 12, y v = 48. Resuelva para el menor valor positivo de
x, si h = 0.4h0 .

Solución: como función auxiliar se recurre a

2πx 2πtv h
f (x) = sin( ) cos( ) + e−x −
λ λ h0
Ahora se trata de aproximar x para el cual f (x) = 0. Para aplicar el
método de Newton-Raphson calculamos f 0 :

2π 2πx 2πtv
f 0 (x) = ( ) cos( ) cos( ) − e−x .
λ λ λ
Con esto las iteraciones del método de Newton-Raphson, para n =
0, 1, 2, ... se ejecutan como:

f (xn )
xn+1 = xn − .
f 0 (xn )
El código Octave para la función auxiliar f está dado por:

function u=p2inf1a(x)
%
L=16;
t=12;
v=48;
r=0.4;
%
u=sin(2*pi*x/L)*cos(2*pi*t*v/L)+exp(-x)-r;
%
endfunction

Note que r = h/h0 = 0.4.

Con la finalidad de identificar un candidato para x0 , esto es, para el
valor inicial para las iteraciones se grafica la función auxiliar según la
instrucción

14
>> fplot(’p2inf1a’,[0 20])

en la lı́nea de comandos de Octave obteniéndose

Notamos que existen al menos dos ceros positivos para la función f,

pero consideramos sólo el más cercano al origen por las condiciones
dadas en el enunciado, esto es, el valor cercano a x0 = 7.
Ahora calculamos las iteraciones según el programa p2inf1 a partir de
x0 = 7 :

function u=p2inf1(x0,m)
x=x0;
%
L=16;
t=12;
v=48;
r=0.4;
%
for it=1:m
y=sin(2*pi*x/L)*cos(2*pi*t*v/L)+exp(-x)-r;
dy=(2*pi/L)*cos(2*pi*x/L)*cos(2*pi*t*v/L)-exp(-x);
x=x-y/dy;
endfor
u=x;
%
endfunction

15
cuya aplicación en la lı́nea de comandos permite obtener los siguientes
resultados en modo format long

>> u=p2inf1(7,5);
y = -1.640468566935560e-02
x = 6.954897305203191
y = -5.993768752793338e-05
x = 6.954731292193685
y = -8.347711966116833e-10
x = 6.954731289881504
y = 3.330669073875470e-16
x = 6.954731289881505
y = -1.110223024625157e-16
x = 6.954731289881505
>>

En los resultados numéricos anteriores y denota a f (xn ) y x a xn .

Para mayor completitud y certeza numérica se informa la aplicación
del comando Octave fzero

>> fzero(’p2inf1a’,7)
ans = 6.954731289881509

el cual considera una combinación inteligente de varios algoritmos, y

una programación experta. Tal como puede apreciarse los resultados
coinciden exactamente, salvo los dos últimos dı́gitos decimales.

Conclusión:

x ≈ 6.95473128.

En términos prácticos basta decir que: x ≈ 6.95.

16
3. La ecuación de Ergun

∆P Dp 3 µ
ρ = 150 (1 − ) + 1.75,
G2o L 1 −  Dp Go

modela el flujo de un lı́quido a través de un cauce comprimido, en

donde, ∆P es la caı́da de presión, ρ es la densidad del fluido, Go es la
velocidad música (el cociente del flujo de masa dividido entre el área
de la sección transversal), Dp es el diámetro de las partı́culas dentro
del cauce, µ es la viscocidad del fluido, L es la longitud del cauce, y 
es la fracción vacı́a del cauce. Se pide calcular la fracción vacı́a  del
cauce, si Dp Go µ−1 = 1000 y ∆P ρDp G−2 o L
−1
= 10.

Solución: como función auxiliar se recurre a

∆P Dp 3 µ
f () = ρ − 150 (1 − ) − 1.75,
G2o L 1 −  Dp Go

Ahora se trata de aproximar  para el cual f () = 0. Para aplicar el

método de Newton-Raphson calculamos f 0 :

f 0 () = A(32 (1 − ) + 3 (−1)) − 150B(−1),

Dp µ
con A = ρ ∆P
G2 L
yB= Dp Go
.
o

Con esto las iteraciones del método de Newton-Raphson, para n =

0, 1, 2, ... se ejecutan como:

f (n )
n+1 = n − .
f 0 (n )
El código Octave para la función auxiliar f está dado por:

function y=p3inf1a(x)
%
A=10;
B=1/1000;
%
y=A*x^3*((1-x)^(-1))-150*B*(1-x)-1.75;
%
endfunction

17
Con la finalidad de identificar un candidato para x0 , esto es, para el
valor inicial para las iteraciones se grafica la función auxiliar según la
instrucción

>> fplot(’p3inf1a’,[0.01 0.8])

en la lı́nea de comandos de Octave obteniéndose

Ahora calculamos las iteraciones según el programa p3inf1 a partir de

0 = 0.4 :

function y=p3inf1(x0,m)
x=x0;
%
A=10;
B=1/1000;
%
for it=1:m
y=A*x^3*((1-x)^(-1))-150*B*(1-x)-1.75;
dy=A*(3*x^2*((1-x)^(-1))+x^3*((1-x)^(-2)))+150*B;
x=x-y/dy;
endfor
y=x;
%
endfunction

18
cuya aplicación en la lı́nea de comandos permite obtener los siguientes
resultados en modo format long

>> u=p3inf1(0.4,5);
y = -7.733333333333331e-01
x = 4.778959149412423e-01
y = 2.621492073051530e-01
x = 4.627224799836014e-01
y = 1.341775633821629e-02
x = 4.618589090568725e-01
y = 4.054922745955203e-05
x = 4.618562834416719e-01
y = 3.734170750391286e-10
x = 4.618562834174922e-01
>>

En los resultados numéricos anteriores y denota a f (n ) y x a n .

Para mayor completitud y certeza numérica se informa la aplicación
del comando Octave fzero

>> fzero(’p3inf1a’,0.4)
ans = 4.618562834174916e-01

el cual considera una combinación inteligente de varios algoritmos, y

una programación experta. Tal como puede apreciarse los resultados
coinciden exactamente, salvo los dos últimos dı́gitos decimales.

Conclusión:

 ≈ 0.46185628.

En términos prácticos basta decir que:  ≈ 0.46.

19
 CÁLCULO NUMÉRICO, MAT1123
M.J. Andrade, D. Sepúlveda, E.Cariaga
Primer Semestre 2020.

INFORME 5

Considere el sistema de ecuaciones no lineales con incógnita (x, y, z, u) :

4x − y + z = xu
−x + 3y − 2z = yu
x − 2y + 3z = zu
x2 + y 2 + z 2 = 1

1. Demuestre que si (x, y, z, u) es una solución de este sistema, entonces

(−x, −y, −z, u) también lo es.

2. Demuestre que los siguientes 6 puntos son soluciones analı́ticas del

sistema dado:
x y√ z√ u
0 −1/√ 2 −1/√ 2 1
p0 1/ √2 1/√2 1
−p 2/3 −1/√ 6 1/ √6 3
2/3
√ 1/√6 −1/√6 3
−1/√ 3 1/ √3 −1/√ 3 6
1/ 3 −1/ 3 1/ 3 6
Obs.: puede resolverlo manualmente, y/o a través de un software idóneo.
La evaluación directa del punto en el sistema sólo puede utilizarla como
verificación, pero no como demostración.

3. Seleccione 1 de los 6 puntos. Luego aplique el método de Newton para

calcular una aproximación Pn+1 del mismo. Su aproximación debe
satisfacer la tolerancia ||Pn+1 − Pn || < 10−8 . Ayuda: recuerde que P0
debe ser un punto cercano a la solución analı́tica seleccionada.

———————–

20
INSTRUCCIONES

ˆ Plazo máximo de entrega: Miércoles 5 de Agosto.

ˆ Este Informe 5 será ponderado como un Examen según Reglamento de su

carrera o facultad.

ˆ Individual.

ˆ Archivo final en PDF.

ˆ Puede venir escrito en .doc o .tex.

ˆ Porcentaje mı́nimo para una calificación aprobatoria: 70%.

ˆ Todo resultado que no está debidamente respaldado por su respectiva captura

de pantalla de Octave, o de otro, no será considerado en la calificación.

ˆ CRITERIOS DE CALIDAD
- El informe debe contar con una Portada con la información básica usual.
- Cada problema debe contar con su Enunciado, Desarrollo, y Conclusión.
- La subsección Desarrollo debe contener todos los detalles, argumentos y
fundamentos.
- Debe indicar explı́citamente, en caso de que ası́ sea, el enunciado del
teorema que está utilizando. Este enunciado debe estar en Anexos.
- En caso de utilizar notación o resultados matemáticos no tratados durante
el curso debe informarlos explı́citamente, junto con su fuente bibliográfica.
- En cada Desarrollo se espera una secuencia lógica y ordenada de pasos.
- El informe finaliza con la Bibliografı́a, seguida de los Anexos.
- En Anexos debe incluir: códigos Octave .m utilizados, copia o muestra
de los resultados obtenidos en la lı́nea de comandos de Octave, captura de
pantalla de los resultados de cualquier software utilizado.

21
 PAUTA INFORME 5. Considere el sistema de ecuaciones no lineales
con incógnita (x, y, z, u) :

4x − y + z = xu
−x + 3y − 2z = yu
x − 2y + 3z = zu
x2 + y 2 + z 2 = 1

1. Demuestre que si (x, y, z, u) es una solución de este sistema, entonces

(−x, −y, −z, u) también lo es.

Solución: sea (x, y, z, u) una solución del sistema dado. Multiplicando

por (−1) las tres primeras ecuaciones, y reescribiendo la cuarta:

4(−x) − (−y) + (−z) = (−x)u

−(−x) + 3(−y) − 2(−z) = (−y)u
(−x) − 2(−y) + 3(−z) = (−z)u
(−x)2 + (−y)2 + (−z)2 = 1.

Con lo anterior se muestra que (−x, −y, −z, u) también es solución del
sistema original. Por lo tanto, queda demostrado que si (x, y, z, u) es
una solución del sistema dado, entonces (−x, −y, −z, u) también lo es.

2. Demuestre que los siguientes 6 puntos son soluciones analı́ticas del

sistema dado:
x y√ z√ u
0 −1/√ 2 −1/√ 2 1
p0 1/ √2 1/√2 1
−p 2/3 −1/√ 6 1/ √6 3
2/3
√ 1/√6 −1/√6 3
−1/√ 3 1/ √3 −1/√ 3 6
1/ 3 −1/ 3 1/ 3 6
Obs.: puede resolverlo manualmente, y/o a través de un software idóneo.
La evaluación directa del punto en el sistema sólo puede utilizarla como
verificación, pero no como demostración.
Solución: digitando la instrucción en la barra de Wolfgramalpha (ver
captura de pantalla) se obtienen las 5 primeras soluciones analı́ticas de

22
 la tabla. La número 6 se obtiene aplicando la proposición del punto 1.
También se acepta como solución válida una deducción algebraica, la
cual también es posible en este caso.

UPGRADE TO PRO APPS TOUR Sign in

Balance chemical
reactions like a pro.
Unlock Step-by-Step

solve system of 4 equations

Extended Keyboard Upload Examples Random

Assuming "solve system of 4 equations" refers to a computation | Use as a general topic instead

Computational Inputs:

» equation 1: 4x-y+z-xu=0

» equation 2: -x+3y-2z-yu=0

» equation 3: x-2y+3z-zu=0

» equation 4: x^2+y^2+z^2-1=0

Compute

Input interpretation:

Results: Approximate forms More solutions

Download Page POWERED BY THE WOLFRAM LANGUAGE

Related Queries:

integral (-x + 3 y - u y - 2 z) du wiki popularity Labrador Retriever vs Cocker Spani…

3. Seleccione 1 de los 6 puntos. Luego aplique el método de Newton para

calcular una aproximación Pn+1 del mismo. Su aproximación debe
satisfacer la tolerancia ||Pn+1 − Pn || < 10−8 . Ayuda: recuerde que P0
debe ser un punto cercano a la solución analı́tica seleccionada.
Solución: sólo a modo ilustrativo se selecciona el primer punto de la
tabla, esto es,
1 1
(0, − √ , − √ , 1)
2 2
como solución analı́tica a aproximar. Note que
1
− √ = −7.071067811865475e − 01.
2
Ejecutando el programa Octave exnr.m con semilla arbitrario
P0 = (x0 , y0 , z0 , u0 ) = (1/2, −1/2, −1/2, 1/2),
se obtiene para m = 5

23
>> [xm ym zm um f1m f2m f3m f4m]=exnr(x0,y0,z0,u0,5)
err_inf = 2.786539472499072e-07
xm = 1.316045500713324e-14
ym = -7.071067811865547e-01
zm = -7.071067811865676e-01
um = 1.000000000000044
f1m = -1.941407085718596e-07
f2m = 1.091167005640159e-07
f3m = 2.849596157927792e-07
f4m = 3.172858851208105e-07

Es decir, ||P5 − P4 || = 2.786539472499072e − 07. Y para m = 6 :

>> [xm ym zm um f1m f2m f3m f4m]=exnr(x0,y0,z0,u0,6)

err_inf = 4.352074256530614e-14
xm = 1.713926170109991e-17
ym = -7.071067811865476e-01
zm = -7.071067811865476e-01
um = 1.000000000000000
f1m = 2.658552927444679e-14
f2m = 4.329869796038111e-14
f3m = 1.820765760385257e-14
f4m = 3.841371665203042e-14

Es decir, ||P6 − P5 || = 4.352074256530614e − 14 < e − 8.

Es decir, esta última ejecución con m=6 satisface por primera vez la
tolerancia exigida en el enunciado.

24
ANEXO:

function [xm ym zm um f1m f2m f3m f4m]=exnr(x0,y0,z0,u0,m)

%
P(1)=x0;
P(2)=y0;
P(3)=z0;
P(4)=u0;
%
for it=1:m
x=P(1);
y=P(2);
z=P(3);
u=P(4);
%
f1=4*x-y+z-x*u;
f2=-x+3*y-2*z-y*u;
f3=x-2*y+3*z-z*u;
f4=x^2+y^2+z^2-1;
%
f1x=4-u;
f1y=-1;
f1z=1;
f1u=-x;
%
f2x=-1;
f2y=3-u;
f2z=-2;
f2u=-y;
%
f3x=1;
f3y=-2;
f3z=3-u;
f3u=-z;
%
f4x=2*x;
f4y=2*y;
f4z=2*z;
f4u=0;
%
J=[f1x f1y f1z f1u;f2x f2y f2z f2u;f3x f3y f3z f3u;f4x f4y f4z f4u];

25
b=-[f1;f2;f3;f4];
deltaP=J\b;
Pa=[P(1);P(2);P(3);P(4)];
P=Pa+deltaP;
endfor
%
err_inf=norm(P-Pa,’inf’)
%
xm=P(1);
ym=P(2);
zm=P(3);
um=P(4);
%
f1m=f1;
f2m=f2;
f3m=f3;
f4m=f4;
%
endfunction

26
 INFORME 1
Ecuaciones No Lineales
Métodos Numéricos, MAT-1123.
Segundo Semestre 20201

Los 4 problemas deben ser resueltos con el método de Newton-

Raphson (NR). En cada caso debe indicar explı́citamente: función
auxiliar f tal que f (x) = 0, fórmula de NR con la notación del
ejercicio, una tabla con n, xn , y f (xn ) con formato format long.
Debe garantizar, al menos, 5 decimales significativos. Debe graficar
f de tal modo que x0 sea identificable.

1. La ecuación de Bernoulli para un fluido que fluye en una canal abierto

con una pequeña estrechez está dada por

Q2 Q2
+ h0 = + h + H,
2gb2 h20 2gb2 h2

en donde Q = 1.2[m3 /s] razón volumétrica de flujo, g = 9.81[m/s2 ]

aceleración gravitacional, b = 1.8[m] ancho del canal, h0 = 0.6[m] nivel
de agua hacia arriba, H = 0.075[m] largo de la estrechez, y h nivel del
agua abajo de la estrechez. Calcular h.

2. La rapidez v de un proyectil en un vuelo vertical cerca de la superficie

de la Tierra puede ser aproximada por
 
M0
v = u ln − gt,
M0 − mt

en donde, u = 2510[m/s], M0 = 2.8 × 106 [kg], m = 13.3 × 103 [kg/s],

g = 9.81[m/s2 ], y t[s] es el tiempo. Calcular el tiempo t para el cual el
proyectil alcanza la velocidad de la luz v = 335[m/s].

3. Considere el ciclo termondinámico de un motor. La eficiencia de este

motor para un gas monoatómico se puede calcular como

ln[T2 /T1 ] − (1 − T1 /T2 )

η= ,
ln[T2 /T1 ] + (1 − T1 /T2 )/(γ − 1)

en donde T es la temperatura absoluta y γ = 5/3. Calcular la razón

T2 /T1 para la cual la eficiencia es de un 30%, es decir, η = 0.3.
1
Prof. Marı́a José Andrade, Emilio Cariaga.

27
 4. La energı́a libre de Gibb de un mol de hidrógeno a temperatura T está
dada por
G = −RT ln[(T /T0 )5/2 ],
en donde, R = 8.31441[J/K], y T0 = 4.44418[K]. Calcular la temperatura
a la cual G = −105 [J].
INSTRUCCIONES
ˆ Plazo máximo de entrega: 18 Octubre 2020. NO se reciben informes fuera
de plazo.
ˆ Formato de archivo final: PDF.
ˆ Puede venir firmado por 1, 2 o 3 estudiantes.
ˆ Tipeado en Word o Tex.
ˆ Gráficos hechos con software.
ˆ En Anexo los programas .m en Octave.
ˆ CRITERIOS DE CALIDAD
- El informe debe contar con una Portada con la información básica usual.
- Cada problema debe contar con Enunciado, Desarrollo, y Conclusión.
- La subsección Desarrollo debe contener todos los detalles, argumentos y
fundamentos.
- Debe indicar explı́citamente, en caso de que ası́ sea, el enunciado del
teorema que está utilizando. Este enunciado debe estar en Anexos.
- En caso de utilizar notación o resultados matemáticos no tratados durante
el curso debe informarlos explı́citamente, junto con su fuente bibliogrı́fica.
Todo esto en el Anexo.
- En cada Desarrollo se espera una secuencia lógica y ordenada de pasos.
- El informe finaliza con la Bibliografı́a, seguida de los Anexos.
- En Anexos debe incluir: códigos Octave .m utilizados, copia o muestra
de los resultados obtenidos en la lı́nea de comandos de Octave, captura de
pantalla de los resultados de cualquier software utilizado.
ˆ EVALUACIÓN
- Nivel de dificultad: 70%.
- Informes que contengan secciones copiadas textualmente serán calificadas
con la nota mı́nima 1.0
- Vencido el plazo de entrega se hará una selección aleatoria de estudiantes
los cuales serán interrogados oralmente. Si demuestran un desconocimiento
significativo de los desarrollos informados su informe será calificado con nota
máxima 3.9.
———————–

28