Hola Steven..! Muchas Gracias de nuevo..! je veo que estás en todos los temás...y bien..!!!!! Es ese el algoritmo que tengo. Estaba debiendo un trabajo de ecuaciones diferenciales en derivadas parciales y se me complicó en ese punto... creí que estaba mal mejor dicho, me trajo muchas dudas. El algoritmo es ese que tenés vos, solo falta ahí la primera aproximación por tratarse de condiciones en derivadas... Pero aparentemente los datos obtenidos por taylor están bien... Capaz que tenga mal las condiciones de estabilidad... (c^2*t^2/h^2<=1 para que sea estable... pero eso es lo que dicen los libros donde h es el paso en x y t en el tiempo) En este momento no te lo paso porque no puedo abrir Maple pero el archivo lo tengo (no se pueden adjuntar acá, y no es texto plano el formato mws. voy a tener que formatear... o al menos reparar... je ni la mayor de las seguridades sirve con windows a veces... No hay firewall, antispyware ni antivirus que valga ya... cuestion de mala suerte :-( ) y muchas gracias por todo!!! información que me sirvio realmente de mucho sino seguiría en la duda del algoritmo. El trabajo ya fué visto y aprobado, lo único que me dijeron que era un problema de constantes y no como debería ser la gráfica... pero aprobado está, estaba entre otros algoritmos para el mismo trabajo y los demás estabn bien... por el momento me voy a poner a pensar puede estar mal porque quiero saberlo y numéricamente me da igual... pero no tengo seguridad si lo está por no saber el resultado que debería dar. Apenas solucione el tema de la pc lo posteo acá... Igualmente ya esta entregado, supongo que la profesora me va a saber decir cuando le vuelva a preguntar... Así que no es tanto problema AHORA mas que un capricho por saber...

<br><br><br> Un saludo!!!!! O mas de uno....!!!!! GRRAAACCIIAASS!!<br><br><div><span class="gmail_quote">El día 9/03/07, <b class="gmail_sendername">Steven Davidson</b> <<a href="mailto:steven@conclase.net">steven@conclase.net

</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 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

D> días sin nada de información....! Y en el amigo google no encuentro nada de D> nada.. Un saludo..! 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">

Cconclase@listas.conclase.net</a><br><a href="http://listas.conclase.net/mailman/listinfo/cconclase_listas.conclase.net">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">

http://listas.conclase.net/index.php?gid=2&mnu=FAQ</a><br></blockquote></div><br>