Practico 2 - Ej7

Practico 2 - Ej7

de Esteban Normey Rieta -
Número de respuestas: 3

Buenas,

Estuve intentando hacer la matriz tridiagonal y resolver el problema como se pide en la letra del ejercicio. Sin embargo, no me quedan soluciones congruentes con lo que debería ( y=sin(x) ), y por más que intenté ver donde está mi error, no consigo verlo.
Me podrían dar una mano e indicarme que estoy haciendo mal?
Este es el código:

 

En respuesta a Esteban Normey Rieta

Re: Practico 2 - Ej7

de Gustavo Rama -
Buenas Esteban,
hay un par de problemas con el código. Primero, la letra pide el intervalo [0, 5]
y no estás teniendo en cuenta el valor inicial y final de y(x). Cuando hacés
la primer ecuación lineal y la última te tienen que aparecer. Y ahí se lo tenés que
restar al vector independiente.

Otra cosa, la próxima manda el código de manera que se pueda copiar y pegar
para poder verificar más rápido.

Saludos,
Gustavo.
En respuesta a Gustavo Rama

Re: Practico 2 - Ej7

de Esteban Normey Rieta -
Corregí lo del intervalo. Sin embargo, sigo sin entender muy bien que tengo que hacer.
Además, en la letra no dice que se debe crear una matriz A de [N-1 * N-1] y resolver utilizando esta matriz? Tipo, al ser de N-1 asumo que no contamos el valor inicial y final en esta matriz.
Otra pregunta que tengo es: los coeficientes que puse en al matriz y en el vector b son correctos (sin contar al valor inicial y final)?
Disculpa la molestia!
En respuesta a Esteban Normey Rieta

Re: Practico 2 - Ej7

de Gustavo Rama -
Buenas,
es correcto, la matriz A no cambia, lo que deberías corregir es el vector b.
Cuando planteas la ecuación y_0+(g_1h^2-2)y_1+y_2 = f_1h^2 tenés que pasar
el y_0 como constante a la derecha, o sea (g_1h^2-2)y_1+y_2 = f_1h^2-y_0.

Lo mismo pasa cuando tomás i=N-1, te queda y_{N-2} + (g_{N-1}h^2-2)y_{N-1} = f_{N-1}h^2 - y_{N}.

En el ejemplo del ejercicio tenés que y_0 = y(0) = 0, y_N = y(5) = \sin(5).

Con eso te debería quedar listo.

Saludos,
Gustavo.