2Eva_IT2010_T3 EDP elíptica, Placa no rectangular

2da Evaluación I Término 2010-2011. 31/Agosto/2010. ICM00158

Tema 3. La placa plana mostrada en la figura está construida con cierto metal, y se ha determinado que la temperatura en los bordes de la placa es la que se indica en la figura.

Ademas de tiene que el término no homogéneo asociado a la ecuación elíptica respectiva es f(x,y)=20

\frac{\delta^2 u}{\delta x^2} + \frac{\delta^2 u}{\delta y^2} = f

El problema consiste en determinar la temperatura en los puntos del interior de la placa en la malla que se muestra en la figura.

a. Determinar el algoritmo en diferencias finitas que resuelve el problema

b. Plantear el sistema de ecuaciones lineas que resuelve el problema

c. Utilice el método de Gauss para resolver el sistema de ecuaciones generado

2Eva_IT2012_T3_MN EDO Taylor 2 Contaminación de estanque

2da Evaluación I Término 2012-2013. 28/Agosto/2012. ICM02188 Métodos Numéricos

Tema 3. (30 puntos) Suponga un estanque de cierto tamaño con agua, la cual está siendo contaminada por una corriente que ingresa constantemente.

En la siguiente ecuación s representa la cantidad de contaminación en el tiempo t:

s'- \frac{26s}{200-t} - \frac{5}{2} = 0 0\leq t \lt 200

Con la condición inicial s(0) = 0, la cual significa que inicialmente el agua está limpia.

Determine la cantidad de contaminación s(t) para

t =  [0.1, 0.2, 0.3, 0.4]

usando la fórmula de Euler, es decir los dos primeros términos de la Serie de Taylor.

2Eva_IT2010_T2 Movimiento angular

2da Evaluación I Término 2010-2011. 31/Agosto/2010. ICM00158

Tema 2. La ecuación de un movimiento angular está dada por

y'' + 10 \sin (y) =0 0\leq t \leq 1 y(0)=0, y'(0)=0.1

Empleando el método de Runge-Kutta de 4to orden generalizado y un paso de 0.25, aproximar la solución de la ecuación en t=0.50


Referencia:  Chapra 28.4 p842 pdf 866

https://nitanperdida.com/2017/12/24/banos-y-el-columpio-del-fin-del-mundo/
BAÑOS DE AGUA SANTA Y EL COLUMPIO DEL FIN DEL MUNDO

2Eva_IT2012_T2_MN Altura del cable teleférico

2da Evaluación I Término 2012-2013. 28/Agosto/2012. ICM02188 Métodos Numéricos

Tema 2. (30 puntos) Utilice los datos del tema anterior para encontrar el valor aproximado de la altura del cable teleférico cuando x = 0.4. Use el polinomio de diferencias finitas de grado 3 y estime el error en la interpolación.

2Eva_IT2012_T1_MN Longitud de teleférico

2da Evaluación I Término 2012-2013. 28/Agosto/2012. ICM02188 Métodos Numéricos

Tema 1. (40 puntos)

La trayectoria de un teleférico está definida por una curva que tiene los puntos (x, f(x)) segun la tabla que se muestra a continuación:

 

x    = [0.00, 0.25, 0.50, 0.75, 1.00]
f(x) = [25.0, 22.0, 32.0, 51.0, 75.0] 

Para calcular la longitud de dicha curva se debe usar la integral:

L = \int_a^b \sqrt{1+[f'(x)]^2} \delta x

a. Aproxime el valor de la derivada f‘(x) en todos los puntos de la tabla con fórmulas de orden 2.

b. Aproxime el valor de la longitud del cable del teleférico entre 0 y 1 con la fórmula de Simpson

c. Aproxime el error de la longitud calculada.

2Eva_IT2012_T3 EDP elíptica

2da Evaluación I Término 2012-2013. 28/Agosto/2012. ICM00158

Tema 3. (20 puntos) Aproxime la solución de la ecuación diferencial parcial:

\frac{\delta ^2 u}{\delta x^2} + \frac{\delta ^2 u}{\delta y^2} = (x^2 + y^2) e^{xy} 0\lt x \lt 1, 0\lt y \lt 0.5

Con las condiciones de frontera:

u(0,y) = 1, u(1,y) = e^y , 0 \leq y \leq 0.5 u(x,0) = 1, u(x,0.5) = \sqrt{e^x} , 0 \leq x \leq 1

Usando un tamaño de paso hx = hy = 0.25

2Eva_IT2012_T2 Modelo de clima

2da Evaluación I Término 2012-2013. 28/Agosto/2012. ICM00158

Tema 2. (20 puntos) El meteorólogo Edward Lorenz propuso inicialmente el siguiente sistema para predecir el comportamiento del clima:

\begin{cases} x'(t) = \alpha (y(t) - x(t)) \\ y'(t) = \rho x(t) - y(t) - x(t) z(t) \\ z'(t) = -\beta z(t) + x(t) y(t) \end{cases}

Para su estudio eligió los parámetros α = 10, β = 8/3 , ρ=28 con las condiciones iniciales x(0) = 10, y(0) = 7, z(0) = 7

Use el método de Runge-Kutta clásico de cuarto orden con h = 0.25 para calcular la solución cuando t=1


Referencia: Chapra 28.2 p833 pdf857

 

https://en.wikipedia.org/wiki/Lorenz_system

2Eva_IT2012_T1 Longitud de teleférico

2da Evaluación I Término 2012-2013. 28/Agosto/2012. ICM00158

Tema 1. (20 puntos) La trayectoria de un teleférico está definida por una curva que tiene los puntos (x, f(x)) segun la tabla que se muestra a continuación:

x    = [ 0.00, 0.25, 0.50, 0.75, 1.00]
f(x) = [25.00,   22,   45,   62,   75  ] 

Para calcular la longitud de dicha curva se debe usar la integral:

\int_0^1 \sqrt{1+[f'(x)]^2} \delta x

a. Aproxime el valor de f'(x) pra cada uno de los valores de x de la tabla

b. Aproxime el valor de la longitud del cable usando el método de Simpson

2Eva_IIT2011_T3_MN Trazador cúbico

2da Evaluación II Término 2011-2012. 31/Enero/2012. ICM02188 Métodos Numéricos

Tema 3. (40 puntos) Dados los puntos (x,y): (2,3), (4,4), (5,6), (6,7), (8,5).

Use el trazador cúbico natural para determinar el valor de y cuando x=3.

Use el método iterativo de Gauss-Seidel para resolver el sistema de ecuaciones que se produce al aplicar la formulación del trazador cúbico.

Comience con un vector solución nulo e itere hasta obtener tres decimales exactos.


xi = [ 2, 4, 5, 6, 8]
yi = [ 3, 4, 6, 7, 5]

2Eva_IIT2011_T2_MN Profundidad en lago

2da Evaluación II Término 2011-2012. 31/Enero/2012. ICM02188 Métodos Numéricos

Tema 2. (30 puntos) Utilizando los 9 datos de las cuadrículas centrales de la tabla anterior, calcule aproximadamente la profundidad del lago en el punto de coordenadas x = 250, y = 125

Utilice en ambas direcciones el polinomio de Lagrange o el polinomio de diferencias finitas y estime el error en la interpolación.

 y\x  0 100 200 300 400
 0  0  0  4  6  0
 50  0  3  5  7  3
 100  1  5  6  9  5
 150  0  2  3  5  1
 200  0  0  1  2  0

profundidad= [[ 0, 0, 4, 6, 0],
              [ 0, 3, 5, 7, 3],
              [ 1, 5, 6, 9, 5],
              [ 0, 2, 3, 5, 1],
              [ 0, 0, 1, 2, 0]]

x = [ 0, 100, 200, 300, 400]
y = [ 0,  50, 100, 150, 200]