Muchísimas gracias a vos también (no lo ví por tener otro asunto el mensaje... ), si, aparentemente la condición de estabilidad y convergencia se cumple... Muchas gracias a vos también, saludos..!!! (ya está aprobado... ahora lo voy a poder ver más tranquilo...) Saludoss! Gracias!
<br><br><div><span class="gmail_quote">El día 9/03/07, <b class="gmail_sendername">Agustín Cota Reguero</b> <<a href="mailto:clom2002@yahoo.es">clom2002@yahoo.es</a>> escribió:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hola, un pequeño comentario a lo que escribio Steven, si no recuerdo mal de alguna vez que tuve que enfrentarme a este tipo de problemas, el hacer lo que dijo Steven no garantiza que la solución que obtengas sea la solución del problema. Asi pues antes de programar algo deberias asegurarte de que para tu problema la solucion del algoritmo converge a la del problema real.
<br><br><b><i>Steven Davidson <<a href="mailto:steven@conclase.net" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">steven@conclase.net</a>></i></b> escribió:<blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; padding-left: 5px;">
 Hola Diego,<br><br>El pasado 2007-03-07 21:04:52, D1e6o! escribió:<br><br>D> Hola amigos.. Bueno, estoy en un gran problema, debo programar el método de<br>D> diferencias finitas para aproximar la ecuación de la onda
<br>D> (d^2T/dt^2=alfa*d^T/dx^2) y creo que mi algoritmo está mal, ya que no me da<br>D> la solución que debería... lo debo hacer en maple, pero si por ahí alguien<br>D> lo tiene al algoritmo o
 el programa en c me serviría muchísimo..! ... Perdon<br>D> por molestar con esto, pero es que debo entregarlo y estoy en los últimos<br>D> días sin nada de información....! Y en el amigo google no encuentro nada de
<br>D> nada.. Un saludo..!<br><br>Bueno, sería conveniente si nos mostraras el algoritmo que tienes y compararlo con el que escribiste. Sospecho que esto ya lo hiciste por tu cuenta, por lo que el otro detalle a analizar es el algoritmo o fórmula en sí.
<br><br>Al hablar del método de diferencias finitas, realmente estamos hablando de un conjunto de métodos que son: diferencia posterior (o avanzada, o adelantada), diferencia anterior (o retrógrada), y diferencia central o centrada. Para la ecuación diferencial parcial a considerar, sugiero la diferencia central. La "versión" finita y discreta es la siguiente ecuación finita parcial (efp):
<br><br>d^2T / dt^2 ~ ( T[j,m+1] - 2T[j,m] + T[j,m-1] ) / (Dt)^2<br>d^2T / dx^2 ~ ( T[j+1,m] - 2T[j,m] + T[j-1,m] ) /
 (Dx)^2<br><br>donde,<br><br>Dt := intervalo constante a "discretizar" para el tiempo, t<br>Dx := intervalo constante a "discretizar" para la longitud, x<br><br>x[] := {0,Dx,2Dx,3Dx,...,J*Dx}, donde tenemos J cantidad de intervalos
<br>t[] := {0,Dt,2Dt,3Dt,...,M*Dt}, donde tenemos M cantidad de intervalos<br><br>x[j] := j*Dx<br>t[m] := m*Dt<br><br>T[j,m] := T(x[j],t[m]), la solución aproximada para x[j] y t[m]<br><br>La solución viene a ser:<br><br>
T[j,m+1] = 2T[j,m] - T[j,m-1] + K ( T[j+1,m] - 2T[j,m] + T[j-1,m] ), donde<br><br>K = c^2 / (Dx/Dt)^2, donde c suele representar la velocidad de la onda<br><br>(Nota: En tu enunciado, has usado 'alfa' en lugar de c^2).
<br><br><br>Obviamente, necesitamos saber las condiciones iniciales, que típicamente son:<br><br>1. T(x,0) = f(x)<br>2. dT(x,0)/t = g(x)<br><br>Usando la efp, las CC.II. finitas son:<br><br>1. T[j,0] = f(x[j]) = f(j*Dx)<br>
2. ( T[j,1] - T[j,-1] ) / (2Dt) = g(x[j]) = g(j*Dx)<br><br><br>Espero que todo esto te
 ayude.<br><br>Steven<br>_______________________________________________<br>Lista de correo Cconclase <a href="mailto:Cconclase@listas.conclase.net" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">
Cconclase@listas.conclase.net</a><br><a href="http://listas.conclase.net/mailman/listinfo/cconclase_listas.conclase.net" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">http://listas.conclase.net/mailman/listinfo/cconclase_listas.conclase.net
</a><br>Bajas: <a href="http://listas.conclase.net/index.php?gid=2&mnu=FAQ" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">http://listas.conclase.net/index.php?gid=2&mnu=FAQ</a><br></blockquote>
<span class="ad"><br><p> 
                </p><hr size="1"><br><font face="Verdana" size="-2">LLama Gratis a cualquier PC del Mundo.<br>Llamadas a fijos y móviles desde 1 céntimo por minuto.<br><a href="http://us.rd.yahoo.com/mail/es/tagline/messenger/*http://es.voice.yahoo.com/" target="_blank" onclick="return top.js.OpenExtLink(window,event,this)">
http://es.voice.yahoo.com</a></font><p></p></span><br>_______________________________________________<br>Lista de correo Cconclase <a onclick="return top.js.OpenExtLink(window,event,this)" href="mailto:Cconclase@listas.conclase.net">
Cconclase@listas.conclase.net</a><br><a onclick="return top.js.OpenExtLink(window,event,this)" href="http://listas.conclase.net/mailman/listinfo/cconclase_listas.conclase.net" target="_blank">http://listas.conclase.net/mailman/listinfo/cconclase_listas.conclase.net
</a><br>Bajas: <a onclick="return top.js.OpenExtLink(window,event,this)" href="http://listas.conclase.net/index.php?gid=2&mnu=FAQ" target="_blank">http://listas.conclase.net/index.php?gid=2&mnu=FAQ</a><br></blockquote>
</div><br>