Difusion_Unidimensional_
DownloadTélécharger
Actions
Vote :
ScreenshotAperçu
Informations
Catégorie :Category: nCreator TI-Nspire
Auteur Author: caxior
Type : Classeur 3.0.1
Page(s) : 1
Taille Size: 2.15 Ko KB
Mis en ligne Uploaded: 29/09/2024 - 21:02:01
Mis à jour Updated: 29/09/2024 - 21:07:56
Uploadeur Uploader: caxior (Profil)
Téléchargements Downloads: 1
Visibilité Visibility: Archive publique
Shortlink : http://ti-pla.net/a4228675
Type : Classeur 3.0.1
Page(s) : 1
Taille Size: 2.15 Ko KB
Mis en ligne Uploaded: 29/09/2024 - 21:02:01
Mis à jour Updated: 29/09/2024 - 21:07:56
Uploadeur Uploader: caxior (Profil)
Téléchargements Downloads: 1
Visibilité Visibility: Archive publique
Shortlink : http://ti-pla.net/a4228675
Description
Fichier Nspire généré sur TI-Planet.org.
Compatible OS 3.0 et ultérieurs.
<<
Define LibPub Difusion1D()= Prgm // Definir las constantes y parámetros iniciales Local k, L, Nx, Nt, dx, dt, S, A, B Local U, Q, x, t Local i, j, m // Parámetros del problema k := 0.01 // Coeficiente de difusión L := 10 // Longitud del dominio Nx := 20 // Número de puntos en el espacio Nt := 100 // Número de pasos de tiempo dx := L / (Nx - 1) // Tamaño del paso en x dt := 0.001 // Tamaño del paso en t S := k * dt / (dx^2) // Parámetro S para estabilidad // Condiciones de frontera y función fuente Func A(t) Return 0 // Condición en x = 0 EndFunc Func B(t) Return 0 // Condición en x = L EndFunc Func f(x) Return sin(À * x / L) // Condición inicial EndFunc Func Q(x, t) Return 0 // Función fuente (puedes modificarla si es necesario) EndFunc // Inicializar el arreglo U (Matriz de tamaño Nx x Nt) U := NewMat(Nx, Nt) // Definir la condición inicial For j, 1, Nx Do U[j, 1] := f((j - 1) * dx) // f(x) EndFor // Implementar el método explícito For m, 1, Nt - 1 Do // Aplicar las condiciones de frontera U[1, m + 1] := A(m * dt) U[Nx, m + 1] := B(m * dt) // Iterar en el espacio For j, 2, Nx - 1 Do U[j, m + 1] := U[j, m] + S * (U[j - 1, m] - 2 * U[j, m] + U[j + 1, m]) + dt * Q((j - 1) * dx, (m + 1) * dt) EndFor EndFor // Mostrar la solución para un tiempo específico (último paso de tiempo) DispMat(U) EndPrgm Made with nCreator - tiplanet.org
>>
Compatible OS 3.0 et ultérieurs.
<<
Define LibPub Difusion1D()= Prgm // Definir las constantes y parámetros iniciales Local k, L, Nx, Nt, dx, dt, S, A, B Local U, Q, x, t Local i, j, m // Parámetros del problema k := 0.01 // Coeficiente de difusión L := 10 // Longitud del dominio Nx := 20 // Número de puntos en el espacio Nt := 100 // Número de pasos de tiempo dx := L / (Nx - 1) // Tamaño del paso en x dt := 0.001 // Tamaño del paso en t S := k * dt / (dx^2) // Parámetro S para estabilidad // Condiciones de frontera y función fuente Func A(t) Return 0 // Condición en x = 0 EndFunc Func B(t) Return 0 // Condición en x = L EndFunc Func f(x) Return sin(À * x / L) // Condición inicial EndFunc Func Q(x, t) Return 0 // Función fuente (puedes modificarla si es necesario) EndFunc // Inicializar el arreglo U (Matriz de tamaño Nx x Nt) U := NewMat(Nx, Nt) // Definir la condición inicial For j, 1, Nx Do U[j, 1] := f((j - 1) * dx) // f(x) EndFor // Implementar el método explícito For m, 1, Nt - 1 Do // Aplicar las condiciones de frontera U[1, m + 1] := A(m * dt) U[Nx, m + 1] := B(m * dt) // Iterar en el espacio For j, 2, Nx - 1 Do U[j, m + 1] := U[j, m] + S * (U[j - 1, m] - 2 * U[j, m] + U[j + 1, m]) + dt * Q((j - 1) * dx, (m + 1) * dt) EndFor EndFor // Mostrar la solución para un tiempo específico (último paso de tiempo) DispMat(U) EndPrgm Made with nCreator - tiplanet.org
>>
