28 May 2024 02:32:35 PM

demo1():
  C version
  Solve heat equation with exact solution
  u(x,t) = e^(-pi^2 t/2) * cos ( pi/2 x )

  Final time T = 1
  Number of spatial nodes n = 10
  Number of time steps m = 50

  Explicit finite difference scheme:
  r = 0.5
  Geomview script written to 'demo1_explicit.gv'
  Maple script written to 'demo1_explicit.mpl'
  Matlab script written to 'demo1_explicit.m'
  Absolute error = 0.0061635

  Implicit finite difference scheme:
  r = 0.5
  Geomview script written to 'demo1_implicit.gv'
  Maple script written to 'demo1_implicit.mpl'
  Matlab script written to 'demo1_implicit.m'
  Absolute error = 0.0118496

  Crank Nicolson linear equation solver:
  r = 0.5
  Geomview script written to 'demo1_crank_nicolson.gv'
  Maple script written to 'demo1_crank_nicolson.mpl'
  Matlab script written to 'demo1_crank_nicolson.m'
  Absolute error = 0.00295428

  Seidman sweep scheme:
  r = 0.25
  Geomview script written to 'demo1_seidman_sweep.gv'
  Maple script written to 'demo1_seidman_sweep.mpl'
  Matlab script written to 'demo1_seidman_sweep.m'
  Absolute error = 0.0048706

demo1():
  Normal end of execution.
28 May 2024 02:32:35 PM