fd1d_heat_implicit_test, a FORTRAN90 code which calls fd1d_heat_implicit(), which solves the time-dependent 1D heat equation, using the finite difference method in space, and an implicit version of the method of lines to handle integration in time.
The computer code and data files described and made available on this web page are distributed under the MIT license
fd1d_heat_implicit, a FORTRAN90 code which uses the finite difference method (FDM) and implicit time stepping to solve the time dependent heat equation in 1D.
x = load ( 'x.txt' ); t = load ( 't.txt' ); u = load ( 'u.txt' ); [ xg, tg ] = meshgrid ( x, t ); mesh ( xg, tg, u );